Assessment of adaptive evolution between wheat and rice as deduced from full-length common wheat cDNA sequence data and expression patterns

Understanding DNA Methylation Patterns in Wheat

The bread wheat genome is large (17 Gb), allohexaploid, and highly repetitive (80-90% of the genome), which makes genomic and epigenomic analyses expensive to conduct and a challenge to analyze. Here we provide an overview of recent bioinformatic and experimental methods that have been developed to understand DNA methylation patterns in the complex polyploid genome of wheat.

The non-random patterns of genetic variation induced by asymmetric somatic hybridization in wheat

BACKGROUND

Asymmetric somatic hybridization is an efficient crop breeding approach by introducing several exogenous chromatin fragments, which leads to genomic shock and therefore induces genome-wide genetic variation. However, the fundamental question concerning the genetic variation such as whether it occurs randomly and suffers from selection pressure remains unknown.

RESULTS

Here, we explored this issue by comparing expressed sequence tags of a common wheat cultivar and its asymmetric somatic hybrid line. Both nucleotide substitutions and indels (insertions and deletions) had lower frequencies in coding sequences than in un-translated regions. The frequencies of nucleotide substitutions and indels were both comparable between chromosomes with and without introgressed fragments. Nucleotide substitutions distributed unevenly and were preferential to indel-flanking sequences, and the frequency of nucleotide substitutions at 5'-flanking sequences of indels was obviously higher in chromosomes with introgressed fragments than in those without exogenous fragment. Nucleotide substitutions and indels both had various frequencies among seven groups of allelic chromosomes, and the frequencies of nucleotide substitutions were strongly negatively correlative to those of indels. Among three sets of genomes, the frequencies of nucleotide substitutions and indels were both heterogeneous, and the frequencies of nucleotide substitutions exhibited drastically positive correlation to those of indels.

CONCLUSIONS

Our work demonstrates that the genetic variation induced by asymmetric somatic hybridization is attributed to both whole genomic shock and local chromosomal shock, which is a predetermined and non-random genetic event being closely associated with selection pressure. Asymmetric somatic hybrids provide a worthwhile model to further investigate the nature of genomic shock induced genetic variation.

The Agropyron cristatum karyotype, chromosome structure and cross-genome homoeology as revealed by fluorescence in situ hybridization with tandem repeats and wheat single-gene probes

Fluorescence in situ hybridization with probes for 45 cDNAs and five tandem repeats revealed homoeologous relationships of Agropyron cristatum with wheat. The results will contribute to alien gene introgression in wheat improvement. Crested wheatgrass (Agropyron cristatum L. Gaertn.) is a wild relative of wheat and a promising source of novel genes for wheat improvement. To date, identification of A. cristatum chromosomes has not been possible, and its molecular karyotype has not been available. Furthermore, homoeologous relationship between the genomes of A. cristatum and wheat has not been determined. To develop chromosome-specific landmarks, A. cristatum genomic DNA was sequenced, and new tandem repeats were discovered. Their distribution on mitotic chromosomes was studied by fluorescence in situ hybridization (FISH), which revealed specific patterns for five repeats in addition to 5S and 45S ribosomal DNA and rye subtelomeric repeats pSc119.2 and pSc200. FISH with one tandem repeat together with 45S rDNA enabled identification of all A. cristatum chromosomes. To analyze the structure and cross-species homoeology of A. cristatum chromosomes with wheat, probes for 45 mapped wheat cDNAs covering all seven chromosome groups were localized by FISH. Thirty-four cDNAs hybridized to homoeologous chromosomes of A. cristatum, nine hybridized to homoeologous and non-homoeologous chromosomes, and two hybridized to unique positions on non-homoeologous chromosomes. FISH using single-gene probes revealed that the wheat-A. cristatum collinearity was distorted, and important structural rearrangements were observed for chromosomes 2P, 4P, 5P, 6P and 7P. Chromosomal inversions were found for pericentric region of 4P and whole chromosome arm 6PL. Furthermore, reciprocal translocations between 2PS and 4PL were detected. These results provide new insights into the genome evolution within Triticeae and will facilitate the use of crested wheatgrass in alien gene introgression into wheat.

Identification and expression profiling of Pht1 phosphate transporters in wheat in controlled environments and in the field

Phosphorus (P) is an important macronutrient with critical functions in plants. Phosphate (Pi) transporters, which mediate Pi acquisition and Pi translocation within the plant, are key factors in Pi deficiency responses. However, their relevance for adaptation to long-term Pi limitation under agronomic conditions, particularly in wheat, remains unknown. Here, we describe the identification of the complete Pi transporter gene family (Pht1) in wheat (Triticum aestivum). Gene expression profiles were compared for hydroponic and field-grown plant tissues of wheat at multiple development stages. Cis-element analysis of selected Pht1 promoter regions was performed. A broad range of expression patterns of individual TaPht1 genes was observed in relation to tissue specificity and the nutrient supply in the soil or in liquid culture, as well as an influence of the experimental system. The expression patterns indicate the involvement of specific transporters in Pi uptake, and in Pi transport and remobilisation within the plant, at different growth developmental stages. Specifically, the influence of Pi nutrition indicates a complex regulatory pattern of TaPht1 gene transcript abundances as a response to low Pi availability in different culture systems, correlating with the existence of different cis-acting promoter elements.

Comparative metabolite profiling of drought stress in roots and leaves of seven Triticeae species

BACKGROUND

Drought is a lifestyle disease. Plant metabolomics has been exercised for understanding the fine-tuning of the potential pathways to surmount the adverse effects of drought stress. A broad spectrum of morphological and metabolic responses from seven Triticeae species including wild types with different drought tolerance/susceptibility level was investigated under control and water scarcity conditions.

RESULTS

Significant morphological parameters measured were root length, surface area, average root diameter and overall root development. Principal Component Analysis, Partial Least-Squares-Discriminant Analysis and Hierarchical Cluster Analysis were applied to the metabolomic data obtained by Gas Chromatography-Mass Spectrometry technique in order to determine the important metabolites of the drought tolerance across seven different Triticeae species. The metabolites showing significant accumulation under the drought stress were considered as the key metabolites and correlated with potential biochemical pathways, enzymes or gene locations for a better understanding of the tolerance mechanisms. In all tested species, 45 significantly active metabolites with possible roles in drought stress were identified. Twenty-one metabolites out of forty-five including sugars, amino acids, organic acids and low molecular weight compounds increased in both leaf and root samples of TR39477, IG132864 and Bolal under the drought stress, contrasting to TTD-22, Tosunbey, Ligustica and Meyeri samples. Three metabolites including succinate, aspartate and trehalose were selected for further genome analysis due to their increased levels in TR39477, IG132864, and Bolal upon drought stress treatment as well as their significant role in energy producing biochemical pathways.

CONCLUSION

These results demonstrated that the genotypes with high drought tolerance skills, especially wild emmer wheat, have a great potential to be a genetic model system for experiments aiming to validate metabolomics-genomics networks.

UGT74S1 is the key player in controlling secoisolariciresinol diglucoside (SDG) formation in flax

BACKGROUND

Flax lignan, commonly known as secoisolariciresinol (SECO) diglucoside (SDG), has recently been reported with health-promoting activities, including its positive impact in metabolic diseases. However, not much was reported on the biosynthesis of SDG and its monoglucoside (SMG) until lately. Flax UGT74S1 was recently reported to sequentially glucosylate SECO into SMG and SDG in vitro. However, whether this gene is the only UGT achieving SECO glucosylation in flax was not known.

RESULTS

Flax genome-wide mining for UGTs was performed. Phylogenetic and gene duplication analyses, heterologous gene expression and enzyme assays were conducted to identify family members closely related to UGT74S1 and to establish their roles in SECO glucosylation. A total of 299 different UGTs were identified, of which 241 (81%) were duplicated. Flax UGTs diverged 2.4-153.6 MYA and 71% were found to be under purifying selection pressure. UGT74S1, a single copy gene located on chromosome 7, displayed no evidence of duplication and was deemed to be under positive selection pressure. The phylogenetic analysis identified four main clusters where cluster 4, which included UGT74S1, was the most diverse. The duplicated UGT74S4 and UGT74S3, located on chromosomes 8 and 14, respectively, were the most closely related to UGT74S1 and were differentially expressed in different tissues. Heterologous expression levels of UGT74S1, UGT74S4 and UGT74S3 proteins were similar but UGT74S4 and UGT74S3 glucosylation activity towards SECO was seven fold less than UGT74S1. In addition, they both failed to produce SDG, suggesting neofunctionalization following their divergence from UGT74S1.

CONCLUSIONS

We showed that UGT74S1 is closely related to two duplicated genes, UGT74S4 and UGT74S3 which, unlike UGT74S1, failed to glucosylate SMG into SDG. The study suggests that UGT74S1 may be the key player in controlling SECO glucosylation into SDG in flax although its closely related genes may also contribute to a minor extent in supplying the SMG precursor to UGT74S1.

Gene Overexpression Resources in Cereals for Functional Genomics and Discovery of Useful Genes

Identification and elucidation of functions of plant genes is valuable for both basic and applied research. In addition to natural variation in model plants, numerous loss-of-function resources have been produced by mutagenesis with chemicals, irradiation, or insertions of transposable elements or T-DNA. However, we may be unable to observe loss-of-function phenotypes for genes with functionally redundant homologs and for those essential for growth and development. To offset such disadvantages, gain-of-function transgenic resources have been exploited. Activation-tagged lines have been generated using obligatory overexpression of endogenous genes by random insertion of an enhancer. Recent progress in DNA sequencing technology and bioinformatics has enabled the preparation of genomewide collections of full-length cDNAs (fl-cDNAs) in some model species. Using the fl-cDNA clones, a novel gain-of-function strategy, Fl-cDNA OvereXpressor gene (FOX)-hunting system, has been developed. A mutant phenotype in a FOX line can be directly attributed to the overexpressed fl-cDNA. Investigating a large population of FOX lines could reveal important genes conferring favorable phenotypes for crop breeding. Alternatively, a unique loss-of-function approach Chimeric REpressor gene Silencing Technology (CRES-T) has been developed. In CRES-T, overexpression of a chimeric repressor, composed of the coding sequence of a transcription factor (TF) and short peptide designated as the repression domain, could interfere with the action of endogenous TF in plants. Although plant TFs usually consist of gene families, CRES-T is effective, in principle, even for the TFs with functional redundancy. In this review, we focus on the current status of the gene-overexpression strategies and resources for identifying and elucidating novel functions of cereal genes. We discuss the potential of these research tools for identifying useful genes and phenotypes for application in crop breeding.

Carotenoids in Staple Cereals: Metabolism, Regulation, and Genetic Manipulation

Carotenoids play a critical role in animal and human health. Animals and humans are unable to synthesize carotenoids de novo, and therefore rely upon diet as sources of these compounds. However, major staple cereals often contain only small amounts of carotenoids in their grains. Consequently, there is considerable interest in genetic manipulation of carotenoid content in cereal grain. In this review, we focus on carotenoid metabolism and regulation in non-green plant tissues, as well as genetic manipulation in staple cereals such as rice, maize, and wheat. Significant progress has been made in three aspects: (1) seven carotenogenes play vital roles in carotenoid regulation in non-green plant tissues, including 1-deoxyxylulose-5-phosphate synthase influencing isoprenoid precursor supply, phytoene synthase, β-cyclase, and ε-cyclase controlling biosynthesis, 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase and carotenoid cleavage dioxygenases responsible for degradation, and orange gene conditioning sequestration sink; (2) provitamin A-biofortified crops, such as rice and maize, were developed by either metabolic engineering or marker-assisted breeding; (3) quantitative trait loci for carotenoid content on chromosomes 3B, 7A, and 7B were consistently identified, eight carotenogenes including 23 loci were detected, and 10 gene-specific markers for carotenoid accumulation were developed and applied in wheat improvement. A comprehensive and deeper understanding of the regulatory mechanisms of carotenoid metabolism in crops will be beneficial in improving our precision in improving carotenoid contents. Genomic selection and gene editing are emerging as transformative technologies for provitamin A biofortification.

Metabolomic profiling and genomic analysis of wheat aneuploid lines to identify genes controlling biochemical pathways in mature grain

Metabolomics is becoming an increasingly important tool in plant genomics to decipher the function of genes controlling biochemical pathways responsible for trait variation. Although theoretical models can integrate genes and metabolites for trait variation, biological networks require validation using appropriate experimental genetic systems. In this study, we applied an untargeted metabolite analysis to mature grain of wheat homoeologous group 3 ditelosomic lines, selected compounds that showed significant variation between wheat lines Chinese Spring and at least one ditelosomic line, tracked the genes encoding enzymes of their biochemical pathway using the wheat genome survey sequence and determined the genetic components underlying metabolite variation. A total of 412 analytes were resolved in the wheat grain metabolome, and principal component analysis indicated significant differences in metabolite profiles between Chinese Spring and each ditelosomic lines. The grain metabolome identified 55 compounds positively matched against a mass spectral library where the majority showed significant differences between Chinese Spring and at least one ditelosomic line. Trehalose and branched-chain amino acids were selected for detailed investigation, and it was expected that if genes encoding enzymes directly related to their biochemical pathways were located on homoeologous group 3 chromosomes, then corresponding ditelosomic lines would have a significant reduction in metabolites compared with Chinese Spring. Although a proportion showed a reduction, some lines showed significant increases in metabolites, indicating that genes directly and indirectly involved in biosynthetic pathways likely regulate the metabolome. Therefore, this study demonstrated that wheat aneuploid lines are suitable experimental genetic system to validate metabolomics-genomics networks.

Delimitation of the Earliness per se D1 (Eps-D1) flowering gene to a subtelomeric chromosomal deletion in bread wheat (Triticum aestivum)

Earliness per se (Eps) genes account for the variation in flowering time when vernalization and photoperiod requirements are satisfied. Genomics and bioinformatics approaches were used to describe allelic variation for 40 Triticum aestivum genes predicted, by synteny with Brachypodium distachyon, to be in the 1DL Eps region. Re-sequencing 1DL genes revealed that varieties carrying early heading alleles at this locus, Spark and Cadenza, carry a subtelomeric deletion including several genes. The equivalent region in Rialto and Avalon is intact. A bimodal distribution in the segregating Spark X Rialto single seed descent (SSD) populations enabled the 1DL QTL to be defined as a discrete Mendelian factor, which we named Eps-D1. Near isogenic lines (NILs) and NIL derived key recombinants between markers flanking Eps-D1 suggest that the 1DL deletion contains the gene(s) underlying Eps-D1. The deletion spans the equivalent of the Triticum monoccocum Eps-A (m) 1 locus, and hence includes MODIFIER OF TRANSCRIPTION 1 (MOT1) and FTSH PROTEASE 4 (FTSH4), the candidates for Eps-A (m) 1. The deletion also contains T. aestivum EARLY FLOWERING 3-D1 (TaELF3-D1) a homologue of the Arabidopsis thaliana circadian clock gene EARLY FLOWERING 3. Eps-D1 is possibly a homologue of Eps-B1 on chromosome 1BL. NILs carrying the Eps-D1 deletion have significantly reduced total TaELF3 expression and altered TaGIGANTEA (TaGI) expression compared with wild type. Altered TaGI expression is consistent with an ELF3 mutant, hence we propose TaELF3-D1 as the more likely candidate for Eps-D1. This is the first direct fine mapping of Eps effect in bread wheat.

Asymmetric somatic hybridization induces point mutations and indels in wheat

Background

Allopolyploid genome needs wide structural variation to deal with genomic shock. The introgression line, generated via asymmetric somatic hybridization, is introgressed with a minimum of exogenous chromatin, which also leads to genomic shock to induce genetic variation. However, the extent of its genomic variation and its difference from allopolyploidies remains unknown.

Methods

Here, we explored this issue using the bread wheat cultivar SR3, a derivative of an asymmetric somatic hybrid between the cultivar JN177 and an accession of tall wheatgrass (Thinopyrum elongatum). The ESTs (expressed sequence taqs) were large-scale sequenced using the cDNA library constructed in each of SR3 and JN177. Point mutations and indels (insertions and deletions) of SR3 were calculated, and their difference from the genetic variation of bread wheat and its ancestors were compared, with aim to analyze the extent and pattern of sequence variation induced by somatic hybridization.

Results

Both point mutations and indels (insertions and deletions) were frequently induced by somatic hybridization in the coding sequences. While the genomic shock caused by allopolyploidization tends to favor deletion over insertion, there was no evidence for such a preference following asymmetric somatic hybridization. The GC content of sequence adjacent to indel sites was also distinct from what has been observed in allopolyploids.

Conclusions

This study demonstrates that asymmetric somatic hybridization induces high frequency of genetic variation in a manner partially different from allopolipoidization. Asymmetric somatic hybridization provides appropriate material to comprehensively explore the nature of the genetic variation induced by genomic shock.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1974-6) contains supplementary material, which is available to authorized users.

MOROKOSHI: transcriptome database in Sorghum bicolor

In transcriptome analysis, accurate annotation of each transcriptional unit and its expression profile is essential. A full-length cDNA (FL-cDNA) collection facilitates the refinement of transcriptional annotation, and accurate transcription start sites help to unravel transcriptional regulation. We constructed a normalized FL-cDNA library from eight growth stages of aerial tissues in Sorghum bicolor and isolated 37,607 clones. These clones were Sanger sequenced from the 5' and/or 3' ends and in total 38,981 high-quality expressed sequence tags (ESTs) were obtained. About one-third of the transcripts of known genes were captured as FL-cDNA clone resources. In addition to these, we also annotated 272 novel genes, 323 antisense transcripts and 1,672 candidate isoforms. These clones are available from the RIKEN Bioresource Center. After obtaining accurate annotation of transcriptional units, we performed expression profile analysis. We carried out spikelet-, seed- and stem-specific RNA sequencing (RNA-Seq) analysis and confirmed the expression of 70.6% of the newly identified genes. We also downloaded 23 sorghum RNA-Seq samples that are publicly available and these are shown on a genome browser together with our original FL-cDNA and RNA-Seq data. Using our original and publicly available data, we made an expression profile of each gene and identified the top 20 genes with the most similar expression. In addition, we visualized their relationships in gene co-expression networks. Users can access and compare various transcriptome data from S, bicolor at http://sorghum.riken.jp.

Characterization of the multigene family TaHKT 2;1 in bread wheat and the role of gene members in plant Na(+) and K(+) status

BACKGROUND

A member of the TaHKT2;1 multigene family was previously identified as a Na(+) transporter with a possible role in root Na(+) uptake. In the present study, the existing full-length cDNA of this member was used as a basis to query the International Wheat Genome Survey Sequence to identify all members of the TaHKT2;1 family. Individual TaHKT2;1 genes were subsequently studied for gene and predicted protein structures, promoter variability, tissue expression and their role in Na(+) and K(+) status of wheat.

RESULTS

Six TaHKT2;1 genes were characterized which included four functional genes (TaHKT2;1 7AL-1, TaHKT2;1 7BL-1, TaHKT2;1 7BL-2 and TaHKT2;1 7DL-1) and two pseudogenes (TaHKT2;1 7AL-2 and TaHKT2;1 7AL-3), on chromosomes 7A, 7B and 7D of hexaploid wheat. Variability in protein domains for cation specificity and in cis-regulatory elements for salt response in gene promoters, were identified amongst the functional TaHKT2;1 members. The functional genes were expressed under low and high NaCl conditions in roots and leaf sheaths, but were down regulated in leaf blades. Alternative splicing events were evident in TaHKT2;1 7AL-1. Aneuploid lines null for each functional gene were grown in high NaCl nutrient solution culture to identify potential role of each TaHKT2;1 member. Aneuploid lines null for TaHKT2;1 7AL-1, TaHKT2;1 7BL-1 and TaHKT2;1 7BL-2 showed no difference in Na(+) concentration between Chinese Spring except for higher Na(+) in sheaths. The same aneuploid lines had lower K(+) in roots, sheath and youngest fully expanded leaf but only under high (200 mM) NaCl in the external solution. There was no difference in Na(+) or K(+) concentration for any treatment between aneuploid line null for the TaHKT2;1 7DL-1 gene and Chinese Spring.

CONCLUSIONS

TaHKT2;1 is a complex family consisting of pseudogenes and functional members. TaHKT2;1 genes do not have an apparent role in controlling root Na(+) uptake in bread wheat seedlings under experimental conditions in this study, contrary to existing hypotheses. However, TaHKT2;1 genes or, indeed other genes in the same chromosome region on 7AL, are candidates that may control Na(+) transport from root to sheath and regulate K(+) levels in different plant tissues.

Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism

Plant metabolites are important to world food security in terms of maintaining sustainable yield and providing food with enriched phytonutrients. Here we report comprehensive profiling of 840 metabolites and a further metabolic genome-wide association study based on ∼6.4 million SNPs obtained from 529 diverse accessions of Oryza sativa. We identified hundreds of common variants influencing numerous secondary metabolites with large effects at high resolution. We observed substantial heterogeneity in the natural variation of metabolites and their underlying genetic architectures among different subspecies of rice. Data mining identified 36 candidate genes modulating levels of metabolites that are of potential physiological and nutritional importance. As a proof of concept, we functionally identified or annotated five candidate genes influencing metabolic traits. Our study provides insights into the genetic and biochemical bases of rice metabolome variation and can be used as a powerful complementary tool to classical phenotypic trait mapping for rice improvement.

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

Unlocking Triticeae genomics to sustainably feed the future

The tribe Triticeae includes the major crops wheat and barley. Within the last few years, the whole genomes of four Triticeae species-barley, wheat, Tausch's goatgrass (Aegilops tauschii) and wild einkorn wheat (Triticum urartu)-have been sequenced. The availability of these genomic resources for Triticeae plants and innovative analytical applications using next-generation sequencing technologies are helping to revitalize our approaches in genetic work and to accelerate improvement of the Triticeae crops. Comparative genomics and integration of genomic resources from Triticeae plants and the model grass Brachypodium distachyon are aiding the discovery of new genes and functional analyses of genes in Triticeae crops. Innovative approaches and tools such as analysis of next-generation populations, evolutionary genomics and systems approaches with mathematical modeling are new strategies that will help us discover alleles for adaptive traits to future agronomic environments. In this review, we provide an update on genomic tools for use with Triticeae plants and Brachypodium and describe emerging approaches toward crop improvements in Triticeae.

Large-scale collection and analysis of full-length cDNAs from Brachypodium distachyon and integration with Pooideae sequence resources

A comprehensive collection of full-length cDNAs is essential for correct structural gene annotation and functional analyses of genes. We constructed a mixed full-length cDNA library from 21 different tissues of Brachypodium distachyon Bd21, and obtained 78,163 high quality expressed sequence tags (ESTs) from both ends of ca. 40,000 clones (including 16,079 contigs). We updated gene structure annotations of Brachypodium genes based on full-length cDNA sequences in comparison with the latest publicly available annotations. About 10,000 non-redundant gene models were supported by full-length cDNAs; ca. 6,000 showed some transcription unit modifications. We also found ca. 580 novel gene models, including 362 newly identified in Bd21. Using the updated transcription start sites, we searched a total of 580 plant cis-motifs in the −3 kb promoter regions and determined a genome-wide Brachypodium promoter architecture. Furthermore, we integrated the Brachypodium full-length cDNAs and updated gene structures with available sequence resources in wheat and barley in a web-accessible database, the RIKEN Brachypodium FL cDNA database. The database represents a “one-stop” information resource for all genomic information in the Pooideae, facilitating functional analysis of genes in this model grass plant and seamless knowledge transfer to the Triticeae crops.

Nitrogen-metabolism related genes in barley - haplotype diversity, linkage mapping and associations with malting and kernel quality parameters

Background

Several studies report about intra-specific trait variation of nitrogen-metabolism related traits, such as N(itrogen)-use efficiency, protein content, N-storage and remobilization in barley and related grass species. The goal of this study was to assess the intra-specific genetic diversity present in primary N-metabolism genes of barley and to investigate the associations of the detected haplotype diversity with malting and kernel quality related traits.

Results

Partial sequences of five genes related to N-metabolism in barley (Hordeum vulgare L.) were obtained, i.e. nitrate reductase 1, glutamine synthetase 2, ferredoxin-dependent glutamate synthase, aspartate aminotransferase and asparaginase. Two to five haplotypes in each gene were discovered in a set of 190 various varieties. The development of 33 SNP markers allowed the genotyping of all these barley varieties consisting of spring and winter types. Furthermore, these markers could be mapped in several doubled haploid populations. Cluster analysis based on haplotypes revealed a more uniform pattern of the spring barleys as compared to the winter barleys. Based on linear model approaches associations to several malting and kernel quality traits including soluble N and protein were identified.

Conclusions

A study was conducted to investigate the presence of sequence variation of several genes related to the primary N-metabolism in barley. The detected diversity could be related to particular phenotypic traits. Specific differences between spring and winter barleys most likely reflect different breeding aims. The developed markers can be used as tool for further genetic studies and marker-assisted selection during breeding of barley.

SNP-revealed genetic diversity in wild emmer wheat correlates with ecological factors

Background

Patterns of genetic diversity between and within natural plant populations and their driving forces are of great interest in evolutionary biology. However, few studies have been performed on the genetic structure and population divergence in wild emmer wheat using a large number of EST-related single nucleotide polymorphism (SNP) markers.

Results

In the present study, twenty-five natural wild emmer wheat populations representing a wide range of ecological conditions in Israel and Turkey were used. Genetic diversity and genetic structure were investigated using over 1,000 SNP markers. A moderate level of genetic diversity was detected due to the biallelic property of SNP markers. Clustering based on Bayesian model showed that grouping pattern is related to the geographical distribution of the wild emmer wheat. However, genetic differentiation between populations was not necessarily dependent on the geographical distances. A total of 33 outlier loci under positive selection were identified using a F_ST-outlier method. Significant correlations between loci and ecogeographical factors were observed.

Conclusions

Natural selection appears to play a major role in generating adaptive structures in wild emmer wheat. SNP markers are appropriate for detecting selectively-channeled adaptive genetic diversity in natural populations of wild emmer wheat. This adaptive genetic diversity is significantly associated with ecological factors.

High mature grain phytase activity in the Triticeae has evolved by duplication followed by neofunctionalization of the purple acid phosphatase phytase (PAPhy) gene

The phytase activity in food and feedstuffs is an important nutritional parameter. Members of the Triticeae tribe accumulate purple acid phosphatase phytases (PAPhy) during grain filling. This accumulation elevates mature grain phytase activities (MGPA) up to levels between ~650 FTU/kg for barley and 6000 FTU/kg for rye. This is notably more than other cereals. For instance, rice, maize, and oat have MGPAs below 100 FTU/kg. The cloning and characterization of the PAPhy gene complement from wheat, barley, rye, einkorn, and Aegilops tauschii is reported here. The Triticeae PAPhy genes generally consist of a set of paralogues, PAPhy_a and PAPhy_b, and have been mapped to Triticeae chromosomes 5 and 3, respectively. The promoters share a conserved core but the PAPhy_a promoter have acquired a novel cis-acting regulatory element for expression during grain filling while the PAPhy_b promoter has maintained the archaic function and drives expression during germination. Brachypodium is the only sequenced Poaceae sharing the PAPhy duplication. As for the Triticeae, the duplication is reflected in a high MGPA of ~4200 FTU/kg in Brachypodium. The sequence conservation of the paralogous loci on Brachypodium chromosomes 1 and 2 does not extend beyond the PAPhy gene. The results indicate that a single-gene segmental duplication may have enabled the evolution of high MGPA by creating functional redundancy of the parent PAPhy gene. This implies that similar MGPA levels may be out of reach in breeding programs for some Poaceae, e.g. maize and rice, whereas Triticeae breeders should focus on PAPhy_a.

Sequence variations of the partially dominant DELLA gene Rht-B1c in wheat and their functional impacts

Rht-B1c, allelic to the DELLA protein-encoding gene Rht-B1a, is a natural mutation documented in common wheat (Triticum aestivum). It confers variation to a number of traits related to cell and plant morphology, seed dormancy, and photosynthesis. The present study was conducted to examine the sequence variations of Rht-B1c and their functional impacts. The results showed that Rht-B1c was partially dominant or co-dominant for plant height, and exhibited an increased dwarfing effect. At the sequence level, Rht-B1c differed from Rht-B1a by one 2kb Veju retrotransposon insertion, three coding region single nucleotide polymorphisms (SNPs), one 197bp insertion, and four SNPs in the 1kb upstream sequence. Haplotype investigations, association analyses, transient expression assays, and expression profiling showed that the Veju insertion was primarily responsible for the extreme dwarfing effect. It was found that the Veju insertion changed processing of the Rht-B1c transcripts and resulted in DELLA motif primary structure disruption. Expression assays showed that Rht-B1c caused reduction of total Rht-1 transcript levels, and up-regulation of GATA-like transcription factors and genes positively regulated by these factors, suggesting that one way in which Rht-1 proteins affect plant growth and development is through GATA-like transcription factor regulation.

Identification of a novel homolog for a calmodulin-binding protein that is upregulated in alloplasmic wheat showing pistillody

Intracellular signaling pathways between the mitochondria and the nucleus are important in both normal and abnormal development in plants. The homeotic transformation of stamens into pistil-like structures (a phenomenon termed pistillody) in cytoplasmic substitution (alloplasmic) lines of bread wheat (Triticum aestivum) has been suggested to be induced by mitochondrial retrograde signaling, one of the forms of intracellular communication. We showed previously that the mitochondrial gene orf260 could alter the expression of nuclear class B MADS-box genes to induce pistillody. To elucidate the interactions between orf260 and nuclear homeotic genes, we performed a microarray analysis to compare gene expression patterns in the young spikes of a pistillody line and a normal line. We identified five genes that showed higher expression levels in the pistillody line. Quantitative expression analysis using real-time PCR indicated that among these five genes, Wheat Calmodulin-Binding Protein 1 (WCBP1) was significantly upregulated in young spikes of the pistillody line. The amino acid sequence of WCBP1 was predicted from the full-length cDNA sequence and found to encode a novel plant calmodulin-binding protein. RT-PCR analysis indicated that WCBP1 was preferentially expressed in young spikes at an early stage and decreased during spike maturation, indicating that it was associated with spikelet/floret development. Furthermore, in situ hybridization analysis suggested that WCBP1 was highly expressed in the pistil-like stamens at early to late developmental stages. These results indicate that WCBP1 plays a role in formation and development of pistil-like stamens induced by mitochondrial retrograde signaling.

Comprehensive functional analyses of expressed sequence tags in common wheat (Triticum aestivum)

About 1 million expressed sequence tag (EST) sequences comprising 125.3 Mb nucleotides were accreted from 51 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including abiotic stresses and pathogen challenges in common wheat (Triticum aestivum). Expressed sequence tags were assembled with stringent parameters after processing with inbuild scripts, resulting in 37,138 contigs and 215,199 singlets. In the assembled sequences, 10.6% presented no matches with existing sequences in public databases. Functional characterization of wheat unigenes by gene ontology annotation, mining transcription factors, full-length cDNA, and miRNA targeting sites were carried out. A bioinformatics strategy was developed to discover single-nucleotide polymorphisms (SNPs) within our large EST resource and reported the SNPs between and within (homoeologous) cultivars. Digital gene expression was performed to find the tissue-specific gene expression, and correspondence analysis was executed to identify common and specific gene expression by selecting four biotic stress-related libraries. The assembly and associated information cater a framework for future investigation in functional genomics.

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

NAD (+) -dependent Formate Dehydrogenase from Plants

NAD(+)-dependent formate dehydrogenase (FDH, EC 1.2.1.2) widely occurs in nature. FDH consists of two identical subunits and contains neither prosthetic groups nor metal ions. This type of FDH was found in different microorganisms (including pathogenic ones), such as bacteria, yeasts, fungi, and plants. As opposed to microbiological FDHs functioning in cytoplasm, plant FDHs localize in mitochondria. Formate dehydrogenase activity was first discovered as early as in 1921 in plant; however, until the past decade FDHs from plants had been considerably less studied than the enzymes from microorganisms. This review summarizes the recent results on studying the physiological role, properties, structure, and protein engineering of plant formate dehydrogenases.

Comprehensive sequence analysis of 24,783 barley full-length cDNAs derived from 12 clone libraries

Full-length cDNA (FLcDNA) libraries consisting of 172,000 clones were constructed from a two-row malting barley cultivar (Hordeum vulgare 'Haruna Nijo') under normal and stressed conditions. After sequencing the clones from both ends and clustering the sequences, a total of 24,783 complete sequences were produced. By removing duplicates between these and publicly available sequences, 22,651 representative sequences were obtained: 17,773 were novel barley FLcDNAs, and 1,699 were barley specific. Highly conserved genes were found in the barley FLcDNA sequences for 721 of 881 rice (Oryza sativa) trait genes with 50% or greater identity. These FLcDNA resources from our Haruna Nijo cDNA libraries and the full-length sequences of representative clones will improve our understanding of the biological functions of genes in barley, which is the cereal crop with the fourth highest production in the world, and will provide a powerful tool for annotating the barley genome sequences that will become available in the near future.

Specific patterns of gene space organisation revealed in wheat by using the combination of barley and wheat genomic resources

Background

Because of its size, allohexaploid nature and high repeat content, the wheat genome has always been perceived as too complex for efficient molecular studies. We recently constructed the first physical map of a wheat chromosome (3B). However gene mapping is still laborious in wheat because of high redundancy between the three homoeologous genomes. In contrast, in the closely related diploid species, barley, numerous gene-based markers have been developed. This study aims at combining the unique genomic resources developed in wheat and barley to decipher the organisation of gene space on wheat chromosome 3B.

Results

Three dimensional pools of the minimal tiling path of wheat chromosome 3B physical map were hybridised to a barley Agilent 15K expression microarray. This led to the fine mapping of 738 barley orthologous genes on wheat chromosome 3B. In addition, comparative analyses revealed that 68% of the genes identified were syntenic between the wheat chromosome 3B and barley chromosome 3 H and 59% between wheat chromosome 3B and rice chromosome 1, together with some wheat-specific rearrangements. Finally, it indicated an increasing gradient of gene density from the centromere to the telomeres positively correlated with the number of genes clustered in islands on wheat chromosome 3B.

Conclusion

Our study shows that novel structural genomics resources now available in wheat and barley can be combined efficiently to overcome specific problems of genetic anchoring of physical contigs in wheat and to perform high-resolution comparative analyses with rice for deciphering the organisation of the wheat gene space.

Genomics and bioinformatics resources for crop improvement

Recent remarkable innovations in platforms for omics-based research and application development provide crucial resources to promote research in model and applied plant species. A combinatorial approach using multiple omics platforms and integration of their outcomes is now an effective strategy for clarifying molecular systems integral to improving plant productivity. Furthermore, promotion of comparative genomics among model and applied plants allows us to grasp the biological properties of each species and to accelerate gene discovery and functional analyses of genes. Bioinformatics platforms and their associated databases are also essential for the effective design of approaches making the best use of genomic resources, including resource integration. We review recent advances in research platforms and resources in plant omics together with related databases and advances in technology.