Micro-colinearity between rice, Brachypodium, and Triticum monococcum at the wheat domestication locus Q

Genome-wide analysis of heavy metal ATPases (HMAs) in Poaceae species and their potential role against copper stress in Triticum aestivum

Plants require copper for normal growth and development and have evolved an efficient system for copper management based on transport proteins such as P_1B-ATPases, also known as heavy metal ATPases (HMAs). Here, we report HMAs in eleven different Poaceae species, including wheat. Furthermore, the possible role of wheat HMAs in copper stress was investigated. BlastP searches identified 27 HMAs in wheat, and phylogenetic analysis based on the Maximum Likelihood method demonstrated a separation into four distinct clades. Conserved motif analysis, domain identification, gene structure, and transmembrane helices number were also identified for wheat HMAs using computational tools. Wheat seedlings grown hydroponically were subjected to elevated copper and demonstrated toxicity symptoms with effects on fresh weight and changes in expression of selected HMAs TaHMA7, TaHMA8, and TaHMA9 were upregulated in response to elevated copper, suggesting a role in wheat copper homeostasis. Further investigations on these heavy metal pumps can provide insight into strategies for enhancing crop heavy metal tolerance in the face of heavy metal pollution.

The MYB transcription factor TaPHR3-A1 is involved in phosphate signaling and governs yield-related traits in bread wheat

Improved inorganic phosphate (Pi) use efficiency in crops will be important for sustainable agriculture. Exploring molecular mechanisms that regulate Pi uptake could provide useful information for breeding wheat with improved Pi use efficiency. Here, a TaPHR3-A1 (Gene ID: TraesCS7A02G415800) ortholog of rice OsPHR3 that functions in transcriptional regulation of Pi signaling was cloned from wheat chromosome 7A. Ectopic expression of TaPHR3-A1 in Arabidopsis and rice produced enhanced vegetative growth and more seeds. Overexpression in transgenic rice led to increased biomass, grain number, and primary panicle branching by 61.23, 42.12, and 36.34% compared with the wild type. Transgenic wheat lines with down-regulation of TaPHR3-A1 exhibited retarded growth and root hair development at the seedling stage, and showed yield-related effects at the adult stage when grown in both low- and sufficient Pi conditions, indicating that TaPHR3-A1 positively regulated tolerance to low Pi. Introgression lines further confirmed the effect of TaPHR3-A1 in improving grain number. The Chinese wheat mini core collection and a recombinant inbred line analysis demonstrated that the favorable allele TaPHR3-A1-A associated with higher grain number was positively selected in breeding. A TaPHR3-A1-derived cleaved amplified polymorphic sequence marker effectively identified haplotype TaPHR3-A1-A. Our results suggested that TaPHR3-A1 was a functional regulatory factor for Pi uptake and provided useful information for marker-assisted selection for high yield in wheat.

Genome-wide identification and expression profiling of cytosine-5 DNA methyltransferases during drought and heat stress in wheat (Triticum aestivum)

DNA methylation is a potential epigenetic mechanism that regulates genome stability, development, and stress mitigation in plants. It is mediated by cytosine-5 DNA methyltransferases (C5-MTases). We identified 52 wheat C5-MTases; and based on domain structure and phylogenetics, these 52 C5-MTases were classified into four sub-families including MET, CMT, DRM and DNMT2; and were distributed on 18 chromosomes. Cis-acting regulatory elements analysis identified abiotic stress-responsive, phytohormone-responsive, development-related and light-related elements in the promoters of TaC5-MTases. We also examined the transcript abundance of TaC5-MTases in different tissues, developmental stages and under abiotic stresses. Notably, most of the TaC5-MTases (TaCMT2, TaCMT3b, TaCMT3c, TaMET1, TaDRM10, TaDNMT2) showed differential regulation of their transcript abundance during drought and heat stress. Overall, the above results provide significant insights into the expression and the probable functions of TaC5-MTases and will also expedite future research programs to explore the mechanisms of epigenetic regulation in wheat.

Systematic analysis and comparison of the PHD-Finger gene family in Chinese pear (Pyrus bretschneideri) and its role in fruit development

PHD-finger proteins, which belongs to the type of zinc finger family, and that play an important role in the regulation of both transcription and the chromatin state in eukaryotes. Currently, PHD-finger proteins have been well studied in animals, while few studies have been carried out on their function in plants. In the present study, 129 non-redundant PHD-finger genes were identified from 5 Rosaceae species (pear, apple, strawberry, mei, and peach); among them, 31 genes were identified in pear. Subsequently, we carried out a bioinformatics analysis of the PHD-finger genes. Thirty-one PbPHD genes were divided into 7 subfamilies based on the phylogenetic analysis, which are consistent with the intron-exon and conserved motif analyses. In addition, we identified five segmental duplication events, implying that the segmental duplications might be a crucial role in the expansion of the PHD-finger gene family in pear. The microsynteny analysis of five Rosaceae species showed that there were independent duplication events in addition to the genome-wide duplication of the pear genome. Subsequently, ten expressed PHD-finger genes of pear fruit were identified using qRT-PCR, and one of these genes, PbPHD10, was identified as an important candidate gene for the regulation of lignin synthesis. Our research provides useful information for the further analysis of the function of PHD-finger gene family in pear.

Genome-wide analysis of wheat calcium ATPases and potential role of selected ACAs and ECAs in calcium stress

BACKGROUND

P₂- type calcium ATPases (ACAs-auto inhibited calcium ATPases and ECAs-endoplasmic reticulum calcium ATPases) belong to the P- type ATPase family of active membrane transporters and are significantly involved in maintaining accurate levels of Ca²⁺, Mn²⁺ and Zn²⁺ in the cytosol as well as playing a very important role in stress signaling, stomatal opening and closing and pollen tube growth. Here we report the identification and possible role of some of these ATPases from wheat.

RESULTS

In this study, ACA and ECA sequences of six species (belonging to Poaceae) were retrieved from different databases and a phylogenetic tree was constructed. A high degree of evolutionary relatedness was observed among P₂ sequences characterized in this study. Members of the respective groups from different plant species were observed to fall under the same clade. This pattern highlights the common ancestry of P_2- type calcium ATPases. Furthermore, qRT-PCR was used to analyse the expression of selected ACAs and ECAs from Triticum aestivum (wheat) under calcium toxicity and calcium deficiency. The data indicated that expression of ECAs is enhanced under calcium stress, suggesting possible roles of these ATPases in calcium homeostasis in wheat. Similarly, the expression of ACAs was significantly different in plants grown under calcium stress as compared to plants grown under control conditions. This gives clues to the role of ACAs in signal transduction during calcium stress in wheat.

CONCLUSION

Here we concluded that wheat genome consists of nine P_2B and three P_2A -type calcium ATPases. Moreover, gene loss events in wheat ancestors lead to the loss of a particular homoeolog of a gene in wheat. To elaborate the role of these wheat ATPases, qRT-PCR was performed. The results indicated that when plants are exposed to calcium stress, both P_2A and P_2B gene expression get enhanced. This further gives clues about the possible role of these ATPases in wheat in calcium management. These findings can be useful in future for genetic manipulations as well as in wheat genome annotation process.

The wheat Phs-A1 pre-harvest sprouting resistance locus delays the rate of seed dormancy loss and maps 0.3 cM distal to the PM19 genes in UK germplasm

The precocious germination of cereal grains before harvest, also known as pre-harvest sprouting, is an important source of yield and quality loss in cereal production. Pre-harvest sprouting is a complex grain defect and is becoming an increasing challenge due to changing climate patterns. Resistance to sprouting is multi-genic, although a significant proportion of the sprouting variation in modern wheat cultivars is controlled by a few major quantitative trait loci, including Phs-A1 in chromosome arm 4AL. Despite its importance, little is known about the physiological basis and the gene(s) underlying this important locus. In this study, we characterized Phs-A1 and show that it confers resistance to sprouting damage by affecting the rate of dormancy loss during dry seed after-ripening. We show Phs-A1 to be effective even when seeds develop at low temperature (13 °C). Comparative analysis of syntenic Phs-A1 intervals in wheat and Brachypodium uncovered ten orthologous genes, including the Plasma Membrane 19 genes (PM19-A1 and PM19-A2) previously proposed as the main candidates for this locus. However, high-resolution fine-mapping in two bi-parental UK mapping populations delimited Phs-A1 to an interval 0.3 cM distal to the PM19 genes. This study suggests the possibility that more than one causal gene underlies this major pre-harvest sprouting locus. The information and resources reported in this study will help test this hypothesis across a wider set of germplasm and will be of importance for breeding more sprouting resilient wheat varieties.

Toward a better understanding of the genomic region harboring Fusarium head blight resistance QTL Qfhs.ndsu-3AS in durum wheat

New molecular markers were developed and mapped to the FHB resistance QTL region in high resolution. Micro-collinearity of the QTL region with rice and Brachypodium was revealed for a better understanding of the genomic region. The wild emmer wheat (Triticum dicoccoides)-derived Fusarium head blight (FHB) resistance quantitative trait locus (QTL) Qfhs.ndsu-3AS previously mapped to the short arm of chromosome 3A (3AS) in a population of recombinant inbred chromosome lines (RICLs). This study aimed to attain a better understanding of the genomic region harboring Qfhs.ndsu-3AS and to improve the utility of the QTL in wheat breeding. Micro-collinearity of the QTL region with rice chromosome 1 and Brachypodium chromosome 2 was identified and used for marker development in saturation mapping. A total of 42 new EST-derived sequence tagged site (STS) and simple sequence repeat (SSR) markers were developed and mapped to the QTL and nearby regions on 3AS. Further comparative analysis revealed a complex collinearity of the 3AS genomic region with their collinear counterparts of rice and Brachypodium. Fine mapping of the QTL region resolved five co-segregating markers (Xwgc1186/Xwgc716/Xwgc1143/Xwgc501/Xwgc1204) into three distinct loci proximal to Xgwm2, a marker previously reported to be closely linked to the QTL. Four other markers (Xwgc1226, Xwgc510, Xwgc1296, and Xwgc1301) mapped farther proximal to the above markers in the QTL region with a higher resolution. Five homozygous recombinants with shortened T. dicoccoides chromosomal segments in the QTL region were recovered by molecular marker analysis and evaluated for FHB resistance. Qfhs.ndsu-3AS was positioned to a 5.2 cM interval flanked by the marker Xwgc501 and Xwgc510. The recombinants containing Qfhs.ndsu-3AS and new markers defining the QTL will facilitate utilization of this resistance source in wheat breeding.

Identification of candidate genes, regions and markers for pre-harvest sprouting resistance in wheat (Triticum aestivum L.)

BACKGROUND

Pre-harvest sprouting (PHS) of wheat grain leads to a reduction in grain yield and quality. The availability of markers for marker-assisted selection (MAS) of PHS resistance will serve to enhance breeding selection and advancement of lines for cultivar development. The aim of this study was to identify candidate regions and develop molecular markers for PHS resistance in wheat. This was achieved via high density mapping of single nucleotide polymorphism (SNP) markers from an Illumina 90 K Infinium Custom Beadchip in a doubled haploid (DH) population derived from a RL4452/'AC Domain' cross and subsequent detection of quantitative trait loci (QTL) for PHS related traits (falling number [FN], germination index [GI] and sprouting index [SI]). SNP marker sequences flanking QTL were used to locate colinear regions in Brachypodium and rice, and identify genic markers associated with PHS resistance that can be utilized for MAS in wheat.

RESULTS

A linkage map spanning 2569.4 cM was constructed with a total of 12,201 SNP, simple sequence repeat (SSR), diversity arrays technology (DArT) and expressed sequence tag (EST) markers. QTL analyses using Multiple Interval Mapping (MIM) identified four QTL for PHS resistance traits on chromosomes 3B, 4A, 7B and 7D. Sequences of SNPs flanking these QTL were subject to a BLASTN search on the International Wheat Genome Sequencing Consortium (IWGSC) database (http://wheat-urgi.versailles.inra.fr/Seq-Repository). Best survey sequence hits were subject to a BLASTN search on Gramene (www.gramene.org) against both Brachypodium and rice databases, and candidate genes and regions for PHS resistance were identified. A total of 18 SNP flanking sequences on chromosomes 3B, 4A, 7B and 7D were converted to KASP markers and validated with matching genotype calls of Infinium SNP data.

CONCLUSIONS

Our study identified candidate genes involved in abscissic acid (ABA) and gibberellin (GA) metabolism, and flowering time in four genomic regions of Brachypodium and rice respectively, in addition to 18 KASP markers for PHS resistance in wheat. These markers can be deployed in future genetic studies of PHS resistance and might also be useful in the evaluation of PHS in germplasm and breeding material.

Dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium distachyon: comparison with common wheat and Aegilops peregrina

BACKGROUND

Thorough understanding of seed starch biosynthesis and accumulation mechanisms is of great importance for agriculture and crop improvement strategies. We conducted the first comprehensive study of the dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium distachyon and compared the findings with those reported for common wheat (Chinese Spring, CS) and Aegilops peregrina.

RESULTS

Only B-granules were identified in Brachypodium Bd21, and the shape variation and development of starch granules were similar in the B-granules of CS and Bd21. Phylogenetic analysis showed that most of the Bd21 starch synthesis-related genes were more similar to those in wheat than in rice. Early expression of key genes in Bd21 starch biosynthesis mediate starch synthesis in the pericarp; intermediate-stage expression increases the number and size of starch granules. In contrast, these enzymes in CS and Ae. peregrina were mostly expressed at intermediate stages, driving production of new B-granules and increasing the granule size, respectively. Immunogold labeling showed that granule-bound starch synthase (GBSSI; related to amylose synthesis) was mainly present in starch granules: at lower levels in the B-granules of Bd21 than in CS. Furthermore, GBSSI was phosphorylated at threonine 183 and tyrosine 185 in the starch synthase catalytic domain in CS and Ae. peregrina, but neither site was phosphorylated in Bd21, suggesting GBSSI phosphorylation could improve amylose biosynthesis.

CONCLUSIONS

Bd21 contains only B-granules, and the expression of key genes in the three studied genera is consistent with the dynamic development of starch granules. GBSSI is present in greater amounts in the B-granules of CS than in Bd21; two phosphorylation sites (Thr183 and Tyr185) were found in Triticum and Aegilops; these sites were not phosphorylated in Bd21. GBSSI phosphorylation may reflect its importance in amylose synthesis.

Comparative high-resolution mapping of the wax inhibitors Iw1 and Iw2 in hexaploid wheat

The wax (glaucousness) on wheat leaves and stems is mainly controlled by two sets of genes: glaucousness loci (W1 and W2) and non-glaucousness loci (Iw1 and Iw2). The non-glaucousness (Iw) loci act as inhibitors of the glaucousness loci (W). High-resolution comparative genetic linkage maps of the wax inhibitors Iw1 originating from Triticum dicoccoides, and Iw2 from Aegilops tauschii were developed by comparative genomics analyses of Brachypodium, sorghum and rice genomic sequences corresponding to the syntenic regions of the Iw loci in wheat. Eleven Iw1 and eight Iw2 linked EST markers were developed and mapped to linkage maps on the distal regions of chromosomes 2BS and 2DS, respectively. The Iw1 locus mapped within a 0.96 cM interval flanked by the BE498358 and CA499581 EST markers that are collinear with 122 kb, 202 kb, and 466 kb genomic regions in the Brachypodium 5S chromosome, the sorghum 6S chromosome and the rice 4S chromosome, respectively. The Iw2 locus was located in a 4.1 to 5.4-cM interval in chromosome 2DS that is flanked by the CJ886319 and CJ519831 EST markers, and this region is collinear with a 2.3 cM region spanning the Iw1 locus on chromosome 2BS. Both Iw1 and Iw2 co-segregated with the BF474014 and CJ876545 EST markers, indicating they are most likely orthologs on 2BS and 2DS. These high-resolution maps can serve as a framework for chromosome landing, physical mapping and map-based cloning of the wax inhibitors in wheat.

The α-gliadin genes from Brachypodium distachyon L. provide evidence for a significant gap in the current genome assembly

Brachypodium distachyon, is a new model plant for most cereal crops while gliadin is a class of wheat storage proteins related with wheat quality attributes. In the published B. distachyon genome sequence databases, no gliadin gene is found. In the current study, a number of gliadin genes in B. distachyon were isolated, which is contradictory to the results of genome sequencing projects. In our study, the B. distachyon seeds were found to have no gliadin protein expression by gel electrophoresis, reversed-phase high-performance liquid chromatography and Western blotting analysis. However, Southern blotting revealed a presence of more than ten copies of α-gliadin coding genes in B. distachyon. By means of AS-PCR amplification, four novel full-ORF α-gliadin genes, and 26 pseudogenes with at least one stop codon as well as their promoter regions were cloned and sequenced from different Brachypodium accessions. Sequence analysis revealed a few of single-nucleotide polymorphisms among these genes. Most pseudogenes were resulted from a C to T change, leading to the generation of TAG or TAA in-frame stop codon. To compare both the full-ORFs and the pseudogenes among Triticum and Triticum-related species, their structural characteristics were analyzed. Based on the four T cell stimulatory toxic epitopes and two ployglutamine domains, Aegilops, Triticum, and Brachypodium species were found to be more closely related. The phylogenetic analysis further revealed that B. distachyon was more closely related to Aegilops tauschii, Aegilops umbellulata, and the A or D genome of Triticum aestivum. The α-gliadin genes were able to express successfully in E. coli using the functional T7 promoter. The relative and absolute quantification of the transcripts of α-gliadin genes in wheat was much higher than that in B. distachyon. The abundant pseudogenes may affect the transcriptional and/or posttranscriptional level of the α-gliadin in B. distachyon.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Arabidopsis Fused kinase and the Kinesin-12 subfamily constitute a signalling module required for phragmoplast expansion

The conserved Fused kinase plays vital but divergent roles in many organisms from Hedgehog signalling in Drosophila to polarization and chemotaxis in Dictyostelium. Previously we have shown that Arabidopsis Fused kinase termed TWO-IN-ONE (TIO) is essential for cytokinesis in both sporophytic and gametophytic cell types. Here using in vivo imaging of GFP-tagged microtubules in dividing microspores we show that TIO is required for expansion of the phragmoplast. We identify the phragmoplast-associated kinesins, PAKRP1/Kinesin-12A and PAKRP1L/Kinesin-12B, as TIO-interacting proteins and determine TIO-Kinesin-12 interaction domains and their requirement in male gametophytic cytokinesis. Our results support the role of TIO as a functional protein kinase that interacts with Kinesin-12 subfamily members mainly through the C-terminal ARM repeat domain, but with a contribution from the N-terminal kinase domain. The interaction of TIO with Kinesin proteins and the functional requirement of their interaction domains support the operation of a Fused kinase signalling module in phragmoplast expansion that depends upon conserved structural features in diverse Fused kinases.

The WRKY transcription factor family in Brachypodium distachyon

Background

A complete assembled genome sequence of wheat is not yet available. Therefore, model plant systems for wheat are very valuable. Brachypodium distachyon (Brachypodium) is such a system. The WRKY family of transcription factors is one of the most important families of plant transcriptional regulators with members regulating important agronomic traits. Studies of WRKY transcription factors in Brachypodium and wheat therefore promise to lead to new strategies for wheat improvement.

Results

We have identified and manually curated the WRKY transcription factor family from Brachypodium using a pipeline designed to identify all potential WRKY genes. 86 WRKY transcription factors were found, a total higher than all other current databases. We therefore propose that our numbering system (BdWRKY1-BdWRKY86) becomes the standard nomenclature. In the JGI v1.0 assembly of Brachypodium with the MIPS/JGI v1.0 annotation, nine of the transcription factors have no gene model and eleven gene models are probably incorrectly predicted. In total, twenty WRKY transcription factors (23.3%) do not appear to have accurate gene models. To facilitate use of our data, we have produced The Database of Brachypodium distachyon WRKY Transcription Factors. Each WRKY transcription factor has a gene page that includes predicted protein domains from MEME analyses. These conserved protein domains reflect possible input and output domains in signaling. The database also contains a BLAST search function where a large dataset of WRKY transcription factors, published genes, and an extensive set of wheat ESTs can be searched. We also produced a phylogram containing the WRKY transcription factor families from Brachypodium, rice, Arabidopsis, soybean, and Physcomitrella patens, together with published WRKY transcription factors from wheat. This phylogenetic tree provides evidence for orthologues, co-orthologues, and paralogues of Brachypodium WRKY transcription factors.

Conclusions

The description of the WRKY transcription factor family in Brachypodium that we report here provides a framework for functional genomics studies in an important model system. Our database is a resource for both Brachypodium and wheat studies and ultimately projects aimed at improving wheat through manipulation of WRKY transcription factors.

Fine mapping of the Bsr1 barley stripe mosaic virus resistance gene in the model grass Brachypodium distachyon

The ND18 strain of Barley stripe mosaic virus (BSMV) infects several lines of Brachypodium distachyon, a recently developed model system for genomics research in cereals. Among the inbred lines tested, Bd3-1 is highly resistant at 20 to 25°C, whereas Bd21 is susceptible and infection results in an intense mosaic phenotype accompanied by high levels of replicating virus. We generated an F_6∶7 recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21 and used the RILs, and an F₂ population of a second Bd21 × Bd3-1 cross to evaluate the inheritance of resistance. The results indicate that resistance segregates as expected for a single dominant gene, which we have designated Barley stripe mosaic virus resistance 1 (Bsr1). We constructed a genetic linkage map of the RIL population using SNP markers to map this gene to within 705 Kb of the distal end of the top of chromosome 3. Additional CAPS and Indel markers were used to fine map Bsr1 to a 23 Kb interval containing five putative genes. Our study demonstrates the power of using RILs to rapidly map the genetic determinants of BSMV resistance in Brachypodium. Moreover, the RILs and their associated genetic map, when combined with the complete genomic sequence of Brachypodium, provide new resources for genetic analyses of many other traits.

Application of Brachypodium to the genetic improvement of wheat roots

To meet the demands of a larger and more affluent global population, wheat yields must increase faster this century than last, with less irrigation, fertilizer, and land. Modelling and experiments consistently demonstrate a large potential for increasing wheat productivity by improving root systems; however, application of research to new varieties is slow because of the inherent difficulties associated with working underground. This review makes the case for the use of the model grass Brachypodium distachyon to simplify root research and accelerate the identification of genes underlying wheat root improvement. Brachypodium is a small temperate grass with many genomic, genetic, and experimental resources that make it a tractable model plant. Brachypodium and wheat have very similar root anatomies which are distinct from rice root anatomy that is specialized to help it overcome anaerobic conditions associated with submerged roots. As a dicotyledonous plant, Arabidopsis has an even more divergent root system that features a tap root system and cambia with secondary growth, both of which are lacking in the grasses. The major advantage of Brachypodium is its small stature that allows the adult grass root system to be readily phenotyped, unlike rice and maize. This will facilitate the identification of genes in adult roots that greatly influence yield by modulating water uptake during flowering and grain development. A summary of the advantages of Brachypodium for root studies is presented, including the adult root system architecture and root growth during grain development. Routes to translate discoveries from Brachypodium to wheat are also discussed.

Identification and comparative mapping of a powdery mildew resistance gene derived from wild emmer (Triticum turgidum var. dicoccoides) on chromosome 2BS

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is an important foliar disease of wheat worldwide. Wild emmer (Triticum turgidum var. dicoccoides) is a valuable genetic resource for improving disease resistance in common wheat. A powdery mildew resistance gene conferring resistance to B. graminis f. sp. tritici isolate E09 at the seedling and adult stages was identified in wild emmer accession IW170 introduced from Israel. An incomplete dominant gene, temporarily designated MlIW170, was responsible for the resistance. Through molecular marker and bulked segregant analyses of an F(2) population and F(3) families derived from a cross between susceptible durum wheat line 81086A and IW170, MlIW170 was located in the distal chromosome bin 2BS3-0.84-1.00 and flanked by SSR markers Xcfd238 and Xwmc243. MlIW170 co-segregated with Xcau516, an STS marker developed from RFLP marker Xwg516 that co-segregated with powdery mildew resistance gene Pm26 on 2BS. Four EST-STS markers, BE498358, BF201235, BQ160080, and BF146221, were integrated into the genetic linkage map of MlIW170. Three AFLP markers, XPaacMcac, XPagcMcta, XPaacMcag, and seven AFLP-derived SCAR markers, XcauG2, XcauG3, XcauG6, XcauG8, XcauG10, XcauG20, and XcauG25, were linked to MlIW170. XcauG3, a resistance gene analog (RGA)-like sequence, co-segregated with MlIW170. The non-glaucousness locus Iw1 was 18.77 cM distal to MlIW170. By comparative genomics of wheat-Brachypodium-rice genomic co-linearity, four EST-STS markers, CJ658408, CJ945509, BQ169830, CJ945085, and one STS marker XP2430, were developed and MlIW170 was mapped in an 2.69 cM interval that is co-linear with a 131 kb genomic region in Brachypodium and a 105 kb genomic region in rice. Four RGA-like sequences annotated in the orthologous Brachypodium genomic region could serve as chromosome landing target regions for map-based cloning of MlIW170.

Fine mapping and comparative genomics integration of two quantitative trait loci controlling resistance to powdery mildew in a Spanish barley landrace

The intervals containing two major quantitative trait loci (QTL) from a Spanish barley landrace conferring broad spectrum resistance to Blumeria graminis were subjected to marker saturation. First, all the available information on recently developed marker resources for barley was exploited. Then, a comparative genomic analysis of the QTL regions with other sequenced grass model species was performed. As a result of the first step, 32 new markers were added to the previous map and new flanking markers closer to both QTL were identified. Next, syntenic integration revealed that the barley target regions showed homology with regions on chromosome 6 of rice (Oryza sativa), chromosome 10 of Sorghum bicolor and chromosome 1 of Brachypodium distachyon. A nested insertion of ancestral syntenic blocks on Brachypodium chromosome 1 was confirmed. Based on sequence information of the most likely candidate orthologous genes, 23 new barley unigene-derived markers were developed and mapped within the barley target regions. The assessment of colinearity revealed an inversion on chromosome 7HL of barley compared to the other three grass species, and nearly perfect colinearity on chromosome 7HS. This two-step marker enrichment allowed for the refinement of the two QTL into much smaller intervals. Inspection of all predicted proteins for the barley unigenes identified within the QTL intervals did not reveal the presence of resistance gene candidates. This study demonstrates the usefulness of sequenced genomes for fine mapping and paves the way for the use of these two loci in barley breeding programs.

Duplication and partitioning in evolution and function of homoeologous Q loci governing domestication characters in polyploid wheat

The Q gene encodes an AP2-like transcription factor that played an important role in domestication of polyploid wheat. The chromosome 5A Q alleles (5AQ and 5Aq) have been well studied, but much less is known about the q alleles on wheat homoeologous chromosomes 5B (5Bq) and 5D (5Dq). We investigated the organization, evolution, and function of the Q/q homoeoalleles in hexaploid wheat (Triticum aestivum L.). Q/q gene sequences are highly conserved within and among the A, B, and D genomes of hexaploid wheat, the A and B genomes of tetraploid wheat, and the A, S, and D genomes of the diploid progenitors, but the intergenic regions of the Q/q locus are highly divergent among homoeologous genomes. Duplication of the q gene 5.8 Mya was likely followed by selective loss of one of the copies from the A genome progenitor and the other copy from the B, D, and S genomes. A recent V(329)-to-I mutation in the A lineage is correlated with the Q phenotype. The 5Bq homoeoalleles became a pseudogene after allotetraploidization. Expression analysis indicated that the homoeoalleles are coregulated in a complex manner. Combined phenotypic and expression analysis indicated that, whereas 5AQ plays a major role in conferring domestication-related traits, 5Dq contributes directly and 5Bq indirectly to suppression of the speltoid phenotype. The evolution of the Q/q loci in polyploid wheat resulted in the hyperfunctionalization of 5AQ, pseudogenization of 5Bq, and subfunctionalization of 5Dq, all contributing to the domestication traits.

Fine mapping and comparative genomics integration of two quantitative trait loci controlling resistance to powdery mildew in a Spanish barley landrace

The intervals containing two major quantitative trait loci (QTL) from a Spanish barley landrace conferring broad spectrum resistance to Blumeria graminis were subjected to marker saturation. First, all the available information on recently developed marker resources for barley was exploited. Then, a comparative genomic analysis of the QTL regions with other sequenced grass model species was performed. As a result of the first step, 32 new markers were added to the previous map and new flanking markers closer to both QTL were identified. Next, syntenic integration revealed that the barley target regions showed homology with regions on chromosome 6 of rice (Oryza sativa), chromosome 10 of Sorghum bicolor and chromosome 1 of Brachypodium distachyon. A nested insertion of ancestral syntenic blocks on Brachypodium chromosome 1 was confirmed. Based on sequence information of the most likely candidate orthologous genes, 23 new barley unigene-derived markers were developed and mapped within the barley target regions. The assessment of colinearity revealed an inversion on chromosome 7HL of barley compared to the other three grass species, and nearly perfect colinearity on chromosome 7HS. This two-step marker enrichment allowed for the refinement of the two QTL into much smaller intervals. Inspection of all predicted proteins for the barley unigenes identified within the QTL intervals did not reveal the presence of resistance gene candidates. This study demonstrates the usefulness of sequenced genomes for fine mapping and paves the way for the use of these two loci in barley breeding programs.

The domestication syndrome genes responsible for the major changes in plant form in the Triticeae crops

The process of crop domestication began 10,000 years ago in the transition of early humans from hunter/gatherers to pastoralists/farmers. Recent research has revealed the identity of some of the main genes responsible for domestication. Two of the major domestication events in barley were (i) the failure of the spike to disarticulate and (ii) the six-rowed spike. The former mutation increased grain yield by preventing grain loss after maturity, while the latter resulted in an up to 3-fold increase in yield potential. Here we provide an overview of the disarticulation systems and inflorescence characteristics, along with the genes underlying these traits, occurring in the Triticeae tribe.

Targeted mapping of Cdu1, a major locus regulating grain cadmium concentration in durum wheat (Triticum turgidum L. var durum)

Some durum wheat (Triticum turgidum L. var durum) cultivars have the genetic propensity to accumulate cadmium (Cd) in the grain. A major gene controlling grain Cd concentration designated as Cdu1 has been reported on 5B, but the genetic factor(s) conferring the low Cd phenotype are currently unknown. The objectives of this study were to saturate the chromosomal region harboring Cdu1 with newly developed PCR-based markers and to investigate the colinearity of this wheat chromosomal region with rice (Oryza sativa L.) and Brachypodium distachyon genomes. Genetic mapping of markers linked to Cdu1 in a population of recombinant inbred substitution lines revealed that the gene(s) associated with variation in Cd concentration resides in wheat bin 5BL9 between fraction breakpoints 0.76 and 0.79. Genetic mapping and quantitative trait locus (QTL) analysis of grain Cd concentration was performed in 155 doubled haploid lines from the cross W9262-260D3 (low Cd) by Kofa (high Cd) revealed two expressed sequence tag markers (ESMs) and one sequence tagged site (STS) marker that co-segregated with Cdu1 and explained >80% of the phenotypic variation in grain Cd concentration. A second, minor QTL for grain Cd concentration was also identified on 5B, 67 cM proximal to Cdu1. The Cdu1 interval spans 286 kbp of rice chromosome 3 and 282 kbp of Brachypodium chromosome 1. The markers and rice and Brachypodium colinearity described here represent tools that will assist in the positional cloning of Cdu1 and can be used to select for low Cd accumulation in durum wheat breeding programs targeting this trait. The isolation of Cdu1 will further our knowledge of Cd accumulation in cereals as well as metal accumulation in general.

A highly conserved gene island of three genes on chromosome 3B of hexaploid wheat: diverse gene function and genomic structure maintained in a tightly linked block

BACKGROUND

The complexity of the wheat genome has resulted from waves of retrotransposable element insertions. Gene deletions and disruptions generated by the fast replacement of repetitive elements in wheat have resulted in disruption of colinearity at a micro (sub-megabase) level among the cereals. In view of genomic changes that are possible within a given time span, conservation of genes between species tends to imply an important functional or regional constraint that does not permit a change in genomic structure. The ctg1034 contig completed in this paper was initially studied because it was assigned to the Sr2 resistance locus region, but detailed mapping studies subsequently assigned it to the long arm of 3B and revealed its unusual features.

RESULTS

BAC shotgun sequencing of the hexaploid wheat (Triticum aestivum cv. Chinese Spring) genome has been used to assemble a group of 15 wheat BACs from the chromosome 3B physical map FPC contig ctg1034 into a 783,553 bp genomic sequence. This ctg1034 sequence was annotated for biological features such as genes and transposable elements. A three-gene island was identified among >80% repetitive DNA sequence. Using bioinformatics analysis there were no observable similarity in their gene functions. The ctg1034 gene island also displayed complete conservation of gene order and orientation with syntenic gene islands found in publicly available genome sequences of Brachypodium distachyon, Oryza sativa, Sorghum bicolor and Zea mays, even though the intergenic space and introns were divergent.

CONCLUSION

We propose that ctg1034 is located within the heterochromatic C-band region of deletion bin 3BL7 based on the identification of heterochromatic tandem repeats and presence of significant matches to chromodomain-containing gypsy LTR retrotransposable elements. We also speculate that this location, among other highly repetitive sequences, may account for the relative stability in gene order and orientation within the gene island.Sequence data from this article have been deposited with the GenBank Data Libraries under accession no. GQ422824.

Fine mapping and syntenic integration of the semi-dwarfing gene sdw3 of barley

The barley mutant allele sdw3 confers a gibberellin-insensitive, semi-dwarf phenotype with potential for breeding of new semi-dwarfed barley cultivars. Towards map-based cloning, sdw3 was delimited by high-resolution genetic mapping to a 0.04 cM interval in a "cold spot" of recombination of the proximal region of the short arm of barley chromosome 2H. Extensive synteny between the barley Sdw3 locus (Hvu_sdw3) and the orthologous regions (Osa_sdw3, Sbi_sdw3, Bsy_sdw3) of three other grass species (Oryza sativa, Sorghum bicolor, Brachypodium sylvaticum) allowed for efficient synteny-based marker saturation in the target interval. Comparative sequence analysis revealed colinearity for 23 out of the 38, 35, and 29 genes identified in Brachypodium, rice, and Sorghum, respectively. Markers co-segregating with Hvu_sdw3 were generated from two of these genes. Initial attempts at chromosome walking in barley were performed with seven orthologous gene probes which were delimiting physical distances of 223, 123, and 127 kb in Brachypodium, rice, and Sorghum, respectively. Six non-overlapping small bacterial artificial chromosome (BAC) clone contigs (cumulative length of 670 kb) were obtained, which indicated a considerably larger physical size of Hvu_sdw3. Low-pass sequencing of selected BAC clones from these barley contigs exhibited a substantially lower gene frequency per physical distance and the presence of additional non-colinear genes. Four candidate genes for sdw3 were identified within barley BAC sequences that either co-segregated with the gene sdw3 or were located adjacent to these co-segregating genes. Identification of genic sequences in the sdw3 context provides tools for marker-assisted selection. Eventual identification of the actual gene will contribute new information for a basic understanding of the mechanisms underlying growth regulation in barley.

Characterization of the primary cell walls of seedlings of Brachypodium distachyon--a potential model plant for temperate grasses

The genome of Brachypodium distachyon, also known as purple false brome, was fully sequenced in 2008 largely in response to the demand for a model plant for temperate grasses. A comparative study of the primary cell walls of seedlings of B. distachyon, Hordeum vulgare and Triticum aestivum was carried out. The cell walls of the three species were characterized by similar relative levels of, and developmental changes in, hemicelluloses. The occurrence of (1,3;1,4)-beta-D-glucans was correlated with phases of growth involving cell elongation. Expression profiling of the genes involved in (1,3;1,4)-beta-D-glucan synthesis (cellulose synthase-like F family (CSLF), CSLH and a putative synthase gene CSLJ) did not show a transcriptional regulation that corresponded to the abundance of (1,3;1,4)-beta-D-glucans. CSLF6 transcripts were similarly highly expressed in all three grasses, and were much more abundant than any of the other transcripts. The CSLH transcript was relatively abundant in B. distachyon but almost undetectable in the other species. The deposition of arabinoxylans increased steadily during seedling growth in all three grasses, but they became less substituted and more cross-linked into the wall matrix during cell maturation. Moreover, arabinoxylans in B. distachyon differed from the two other grasses in having a lower degree of arabinose substitution, a higher percentage of ferulic acid in form of dimers and a larger proportion of ester-linked p-coumaric acid.

A BAC-based physical map of Brachypodium distachyon and its comparative analysis with rice and wheat

Background

Brachypodium distachyon (Brachypodium) has been recognized as a new model species for comparative and functional genomics of cereal and bioenergy crops because it possesses many biological attributes desirable in a model, such as a small genome size, short stature, self-pollinating habit, and short generation cycle. To maximize the utility of Brachypodium as a model for basic and applied research it is necessary to develop genomic resources for it. A BAC-based physical map is one of them. A physical map will facilitate analysis of genome structure, comparative genomics, and assembly of the entire genome sequence.

Results

A total of 67,151 Brachypodium BAC clones were fingerprinted with the SNaPshot HICF fingerprinting method and a genome-wide physical map of the Brachypodium genome was constructed. The map consisted of 671 contigs and 2,161 clones remained as singletons. The contigs and singletons spanned 414 Mb. A total of 13,970 gene-related sequences were detected in the BAC end sequences (BES). These gene tags aligned 345 contigs with 336 Mb of rice genome sequence, showing that Brachypodium and rice genomes are generally highly colinear. Divergent regions were mainly in the rice centromeric regions. A dot-plot of Brachypodium contigs against the rice genome sequences revealed remnants of the whole-genome duplication caused by paleotetraploidy, which were previously found in rice and sorghum. Brachypodium contigs were anchored to the wheat deletion bin maps with the BES gene-tags, opening the door to Brachypodium-Triticeae comparative genomics.

Conclusion

The construction of the Brachypodium physical map, and its comparison with the rice genome sequence demonstrated the utility of the SNaPshot-HICF method in the construction of BAC-based physical maps. The map represents an important genomic resource for the completion of Brachypodium genome sequence and grass comparative genomics. A draft of the physical map and its comparisons with rice and wheat are available at .

Control of flowering time and spike development in cereals: the earliness per se Eps-1 region in wheat, rice, and Brachypodium

The earliness per se gene Eps-A^m1 from diploid wheat Triticum monococcum affects heading time, spike development, and spikelet number. In this study, the Eps1 orthologous regions from rice, Aegilops tauschii, and Brachypodium distachyon were compared as part of current efforts to clone this gene. A single Brachypodium BAC clone spanned the Eps-A^m1 region, but a gap was detected in the A. tauschii physical map. Sequencing of the Brachypodium and A. tauschii BAC clones revealed three genes shared by the three species, which showed higher identity between wheat and Brachypodium than between them and rice. However, most of the structural changes were detected in the wheat lineage. These included an inversion encompassing the wg241-VatpC region and the presence of six unique genes. In contrast, only one unique gene (and one pseudogene) was found in Brachypodium and none in rice. Three genes were present in both Brachypodium and wheat but were absent in rice. Two of these genes, Mot1 and FtsH4, were completely linked to the earliness per se phenotype in the T. monococcum high-density genetic map and are candidates for Eps-A^m1. Both genes were expressed in apices and developing spikes, as expected for Eps-A^m1 candidates. The predicted MOT1 protein showed amino acid differences between the parental T. monococcum lines, but its effect is difficult to predict. Future steps to clone the Eps-A^m1 gene include the generation of mot1 and ftsh4 mutants and the completion of the T. monococcum physical map to test for the presence of additional candidate genes.

DArT markers: diversity analyses, genomes comparison, mapping and integration with SSR markers in Triticum monococcum

BACKGROUND

Triticum monococcum (2n = 2x = 14) is an ancient diploid wheat with many useful traits and is used as a model for wheat gene discovery. DArT (Diversity Arrays Technology) employs a hybridisation-based approach to type thousands of genomic loci in parallel. DArT markers were developed for T. monococcum to assess genetic diversity, compare relationships with hexaploid genomes, and construct a genetic linkage map integrating DArT and microsatellite markers.

RESULTS

A DArT array, consisting of 2304 hexaploid wheat, 1536 tetraploid wheat, 1536 T. monococcum as well as 1536 T. boeoticum representative genomic clones, was used to fingerprint 16 T. monococcum accessions of diverse geographical origins. In total, 846 polymorphic DArT markers were identified, of which 317 were of T. monococcum origin, 246 of hexaploid, 157 of tetraploid, and 126 of T. boeoticum genomes. The fingerprinting data indicated that the geographic origin of T. monococcum accessions was partially correlated with their genetic variation. DArT markers could also well distinguish the genetic differences amongst a panel of 23 hexaploid wheat and nine T. monococcum genomes. For the first time, 274 DArT markers were integrated with 82 simple sequence repeat (SSR) and two morphological trait loci in a genetic map spanning 1062.72 cM in T. monococcum. Six chromosomes were represented by single linkage groups, and chromosome 4Am was formed by three linkage groups. The DArT and SSR genetic loci tended to form independent clusters along the chromosomes. Segregation distortion was observed for one third of the DArT loci. The Ba (black awn) locus was refined to a 23.2 cM region between the DArT marker locus wPt-2584 and the microsatellite locus Xgwmd33 on 1Am; and the Hl (hairy leaf) locus to a 4.0 cM region between DArT loci 376589 and 469591 on 5Am.

CONCLUSION

DArT is a rapid and efficient approach to develop many new molecular markers for genetic studies in T. monococcum. The constructed genetic linkage map will facilitate localisation and map-based cloning of genes of interest, comparative mapping as well as genome organisation and evolution studies between this ancient diploid species and other crops.

DArT markers: diversity analyses, genomes comparison, mapping and integration with SSR markers in Triticum monococcum

Background

Triticum monococcum (2n = 2x = 14) is an ancient diploid wheat with many useful traits and is used as a model for wheat gene discovery. DArT (Diversity Arrays Technology) employs a hybridisation-based approach to type thousands of genomic loci in parallel. DArT markers were developed for T. monococcum to assess genetic diversity, compare relationships with hexaploid genomes, and construct a genetic linkage map integrating DArT and microsatellite markers.

Results

A DArT array, consisting of 2304 hexaploid wheat, 1536 tetraploid wheat, 1536 T. monococcum as well as 1536 T. boeoticum representative genomic clones, was used to fingerprint 16 T. monococcum accessions of diverse geographical origins. In total, 846 polymorphic DArT markers were identified, of which 317 were of T. monococcum origin, 246 of hexaploid, 157 of tetraploid, and 126 of T. boeoticum genomes. The fingerprinting data indicated that the geographic origin of T. monococcum accessions was partially correlated with their genetic variation. DArT markers could also well distinguish the genetic differences amongst a panel of 23 hexaploid wheat and nine T. monococcum genomes. For the first time, 274 DArT markers were integrated with 82 simple sequence repeat (SSR) and two morphological trait loci in a genetic map spanning 1062.72 cM in T. monococcum. Six chromosomes were represented by single linkage groups, and chromosome 4A^m was formed by three linkage groups. The DArT and SSR genetic loci tended to form independent clusters along the chromosomes. Segregation distortion was observed for one third of the DArT loci. The Ba (black awn) locus was refined to a 23.2 cM region between the DArT marker locus wPt-2584 and the microsatellite locus Xgwmd33 on 1A^m; and the Hl (hairy leaf) locus to a 4.0 cM region between DArT loci 376589 and 469591 on 5A^m.

Conclusion

DArT is a rapid and efficient approach to develop many new molecular markers for genetic studies in T. monococcum. The constructed genetic linkage map will facilitate localisation and map-based cloning of genes of interest, comparative mapping as well as genome organisation and evolution studies between this ancient diploid species and other crops.

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.

TriFLDB: a database of clustered full-length coding sequences from Triticeae with applications to comparative grass genomics

The Triticeae Full-Length CDS Database (TriFLDB) contains available information regarding full-length coding sequences (CDSs) of the Triticeae crops wheat (Triticum aestivum) and barley (Hordeum vulgare) and includes functional annotations and comparative genomics features. TriFLDB provides a search interface using keywords for gene function and related Gene Ontology terms and a similarity search for DNA and deduced translated amino acid sequences to access annotations of Triticeae full-length CDS (TriFLCDS) entries. Annotations consist of similarity search results against several sequence databases and domain structure predictions by InterProScan. The deduced amino acid sequences in TriFLDB are grouped with the proteome datasets for Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and sorghum (Sorghum bicolor) by hierarchical clustering in stepwise thresholds of sequence identity, providing hierarchical clustering results based on full-length protein sequences. The database also provides sequence similarity results based on comparative mapping of TriFLCDSs onto the rice and sorghum genome sequences, which together with current annotations can be used to predict gene structures for TriFLCDS entries. To provide the possible genetic locations of full-length CDSs, TriFLCDS entries are also assigned to the genetically mapped cDNA sequences of barley and diploid wheat, which are currently accommodated in the Triticeae Mapped EST Database. These relational data are searchable from the search interfaces of both databases. The current TriFLDB contains 15,871 full-length CDSs from barley and wheat and includes putative full-length cDNAs for barley and wheat, which are publicly accessible. This informative content provides an informatics gateway for Triticeae genomics and grass comparative genomics. TriFLDB is publicly available at http://TriFLDB.psc.riken.jp/.

Orthology between genomes of Brachypodium, wheat and rice

BACKGROUND

In the past, rice genome served as a good model for studies involving comparative genomics of grass species. More recently, however, Brachypodium distachyon genome has emerged as a better model system for genomes of temperate cereals including wheat. During the present study, Brachypodium EST contigs were utilized to resolve orthologous relationships among the genomes of Brachypodium, wheat and rice.

FINDINGS

Comparative sequence analysis of 3,818 Brachypodium EST (bEST) contigs and 3,792 physically mapped wheat EST (wEST) contigs revealed that as many as 449 bEST contigs were orthologous to 1,154 wEST loci that were bin-mapped on all the 21 wheat chromosomes. Similarly 743 bEST contigs were orthologous to specific rice genome sequences distributed on all the 12 rice chromosomes. As many as 183 bEST contigs were orthologous to both wheat and rice genome sequences, which harbored as many as 17 SSRs conserved across the three species. Primers developed for 12 of these 17 conserved SSRs were used for a wet-lab experiment, which resolved relatively high level of conservation among the genomes of Brachypodium, wheat and rice.

CONCLUSION

The present study confirmed that Brachypodium is a better model than rice for analysis of the genomes of temperate cereals like wheat and barley. The whole genome sequence of Brachypodium, which should become available in the near future, will further facilitate greatly the studies involving comparative genomics of cereals.

Structural characterization of Brachypodium genome and its syntenic relationship with rice and wheat

Brachypodium distachyon (Brachypodium) has been recently recognized as an emerging model system for both comparative and functional genomics in grass species. In this study, 55,221 repeat masked Brachypodium BAC end sequences (BES) were used for comparative analysis against the 12 rice pseudomolecules. The analysis revealed that approximately 26.4% of BES have significant matches with the rice genome and 82.4% of the matches were homologous to known genes. Further analysis of paired-end BES and approximately 1.0 Mb sequences from nine selected BACs proved to be useful in revealing conserved regions and regions that have undergone considerable genomic changes. Differential gene amplification, insertions/deletions and inversions appeared to be the common evolutionary events that caused variations of microcolinearity at different orthologous genomic regions. It was found that approximately 17% of genes in the two genomes are not colinear in the orthologous regions. Analysis of BAC sequences also revealed higher gene density (approximately 9 kb/gene) and lower repeat DNA content (approximately 13.1%) in Brachypodium when compared to the orthologous rice regions, consistent with the smaller size of the Brachypodium genome. The 119 annotated Brachypodium genes were BLASTN compared against the wheat EST database and deletion bin mapped wheat ESTs. About 77% of the genes retrieved significant matches in the EST database, while 9.2% matched to the bin mapped ESTs. In some cases, genes in single Brachypodium BACs matched to multiple ESTs that were mapped to the same deletion bins, suggesting that the Brachypodium genome will be useful for ordering wheat ESTs within the deletion bins and developing specific markers at targeted regions in the wheat genome.

Sixty million years in evolution of soft grain trait in grasses: emergence of the softness locus in the common ancestor of Pooideae and Ehrhartoideae, after their divergence from Panicoideae

Together maize, Sorghum, rice, and wheat grass (Poaceae) species are the most important cereal crops in the world and exhibit different "grain endosperm texture." This trait has been studied extensively in wheat because of its pivotal role in determining quality of products obtained from wheat grain. Grain softness protein-1 and Puroindolines A and B (grain storage proteins), encoded by Ha-like genes: Gsp-1, Pina, and Pinb, of the Hardness (Ha) locus, are the main determinants of the grain softness/hardness trait in wheat. The origin and evolution of grain endosperm texture in grasses was addressed by comparing genomic sequences of the Ha orthologous region of wheat, Brachypodium, rice, and Sorghum. Results show that the Ha-like genes are present in wheat and Brachypodium but are absent from Sorghum bicolor. A truncated remnant of an Ha-like gene is present in rice. Synteny analysis of the genomes of these grass species shows that only one of the paralogous Ha regions, created 70 My by whole-genome duplication, contained Ha-like genes. The comparative genome analysis and evolutionary comparison with genes encoding grain reserve proteins of grasses suggest that an ancestral Ha-like gene emerged, as a new member of the prolamin gene family, in a common ancestor of the Pooideae (Triticeae and Brachypoidieae tribes) and Ehrhartoideae (rice), between 60 and 50 My, after their divergence from Panicoideae (Sorghum). It was subsequently lost in Ehrhartoideae. Recurring duplications, deletions, and/or truncations occurred independently and appear to characterize Ha-like gene evolution in the grass species. The Ha-like genes gained a new function in Triticeae, such as wheat, underlying the soft grain phenotype. Loss of these genes in some wheat species leads, in turn, to hard endosperm seeds.

The 'inner circle' of the cereal genomes

Early marker-based macrocolinearity studies between the grass genomes led to arranging their chromosomes into concentric 'crop circles' of synteny blocks that initially consisted of 30 rice-independent linkage groups representing the ancestral cereal genome structure. Recently, increased marker density and genome sequencing of several cereal genomes allowed the characterization of intragenomic duplications and their integration with intergenomic colinearity data to identify paleo-duplications and propose a model for the evolution of the grass genomes from a common ancestor. On the basis of these data an 'inner circle' comprising five ancestral chromosomes was defined providing a new reference for the grass chromosomes and new insights into their ancestral relationships and origin, as well as an efficient tool to design cross-genome markers for genetic studies.