Comparison of genes among cereals

Molecular evidence for adaptive evolution of drought tolerance in wild cereals

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

Systems biology approach in plant abiotic stresses

Plant abiotic stresses are the major constraint on plant growth and development, causing enormous crop losses across the world. Plants have unique features to defend themselves against these challenging adverse stress conditions. They modulate their phenotypes upon changes in physiological, biochemical, molecular and genetic information, thus making them tolerant against abiotic stresses. It is of paramount importance to determine the stress-tolerant traits of a diverse range of genotypes of plant species and integrate those traits for crop improvement. Stress-tolerant traits can be identified by conducting genome-wide analysis of stress-tolerant genotypes through the highly advanced structural and functional genomics approach. Specifically, whole-genome sequencing, development of molecular markers, genome-wide association studies and comparative analysis of interaction networks between tolerant and susceptible crop varieties grown under stress conditions can greatly facilitate discovery of novel agronomic traits that protect plants against abiotic stresses.

The major threshability genes soft glume (sog) and tenacious glume (Tg), of diploid and polyploid wheat, trace their origin to independent mutations at non-orthologous loci

Threshability is an important crop domestication trait. The wild wheat progenitors have tough glumes enveloping the floret that make spikes difficult to thresh, whereas cultivated wheats have soft glumes and are free-threshing. In hexaploid wheat, the glume tenacity gene Tg along with the major domestication locus Q control threshability. The Q gene was isolated recently and found to be a member of the AP2 class of transcription factors. However, only a few studies have reported on the tough glume trait. Here, we report comparative mapping of the soft glume (sog) gene of diploid Triticum monococcum L. and tenacious glume (Tg) gene of hexaploid T. aestivum L. using chromosome-specific SSR and RFLP markers. The sog gene was flanked by Xgwm71 and Xbcd120 in a 6.8 cM interval on chromosome 2A(m)S of T. monococcum whereas Tg was targeted to a 8.1 cM interval flanked by Xwmc503 and Xfba88 on chromosome 2DS of T. aestivum. Deletion bin mapping of the flanking markers assigned sog close to the centromere on 2AS, whereas Tg was mapped to the most distal region on 2DS. Both 2AS and 2DS maps were colinear ruling out the role of chromosome rearrangements for their non-syntenic positions. Therefore, sog and Tg are not true orthologues suggesting the possibility of a diverse origin.

The problem of the eukaryotic genome size

The current state of knowledge concerning the unsolved problem of the huge interspecific eukaryotic genome size variations not correlating with the species phenotypic complexity (C-value enigma also known as C-value paradox) is reviewed. Characteristic features of eukaryotic genome structure and molecular mechanisms that are the basis of genome size changes are examined in connection with the C-value enigma. It is emphasized that endogenous mutagens, including reactive oxygen species, create a constant nuclear environment where any genome evolves. An original quantitative model and general conception are proposed to explain the C-value enigma. In accordance with the theory, the noncoding sequences of the eukaryotic genome provide genes with global and differential protection against chemical mutagens and (in addition to the anti-mutagenesis and DNA repair systems) form a new, third system that protects eukaryotic genetic information. The joint action of these systems controls the spontaneous mutation rate in coding sequences of the eukaryotic genome. It is hypothesized that the genome size is inversely proportional to functional efficiency of the anti-mutagenesis and/or DNA repair systems in a particular biological species. In this connection, a model of eukaryotic genome evolution is proposed.

Quantitative trait loci (QTL) analysis for rice grain width and fine mapping of an identified QTL allele gw-5 in a recombination hotspot region on chromosome 5

Rice grain width and shape play a crucial role in determining grain quality and yield. The genetic basis of rice grain width was dissected into six additive quantitative trait loci (QTL) and 11 pairs of epistatic QTL using an F(7) recombinant inbred line (RIL) population derived from a single cross between Asominori (japonica) and IR24 (indica). QTL by environment interactions were evaluated in four environments. Chromosome segment substitution lines (CSSLs) harboring the six additive effect QTL were used to evaluate gene action across eight environments. A major, stable QTL, qGW-5, consistently decreased rice grain width in both the Asominori/IR24 RIL and CSSL populations with the genetic background Asominori. By investigating the distorted segregation of phenotypic values of rice grain width and genotypes of molecular markers in BC(4)F(2) and BC(4)F(3) populations, qGW-5 was dissected into a single recessive gene, gw-5, which controlled both grain width and length-width ratio. gw-5 was narrowed down to a 49.7-kb genomic region with high recombination frequencies on chromosome 5 using 6781 BC(4)F(2) individuals and 10 newly developed simple sequence repeat markers. Our results provide a basis for map-based cloning of the gw-5 gene and for marker-aided gene/QTL pyramiding in rice quality breeding.

Analysis of genome organization, composition and microsynteny using 500 kb BAC sequences in chickpea

The small genome size (740 Mb), short life cycle (3 months) and high economic importance as a food crop legume make chickpea (Cicer arietinum L.) an important system for genomics research. Although several genetic linkage maps using various markers and genomic tools have become available, sequencing efforts and their use are limited in chickpea genomic research. In this study, we explored the genome organization of chickpea by sequencing approximately 500 kb from 11 BAC clones (three representing ascochyta blight resistance QTL1 (ABR-QTL1) and eight randomly selected BAC clones). Our analysis revealed that these sequenced chickpea genomic regions have a gene density of one per 9.2 kb, an average gene length of 2,500 bp, an average of 4.7 exons per gene, with an average exon and intron size of 401 and 316 bp, respectively, and approximately 8.6% repetitive elements. Other features analyzed included exon and intron length, number of exons per gene, protein length and %GC content. Although there are reports on high synteny among legume genomes, the microsynteny between the 500 kb chickpea and available Medicago truncatula genomic sequences varied depending on the region analyzed. The GBrowse-based annotation of these BACs is available at http://www.genome.ou.edu/plants_totals.html . We believe that our work provides significant information that supports a chickpea genome sequencing effort in the future.

Comparative sequence analysis of Solanum and Arabidopsis in a hot spot for pathogen resistance on potato chromosome V reveals a patchwork of conserved and rapidly evolving genome segments

Background

Quantitative phenotypic variation of agronomic characters in crop plants is controlled by environmental and genetic factors (quantitative trait loci = QTL). To understand the molecular basis of such QTL, the identification of the underlying genes is of primary interest and DNA sequence analysis of the genomic regions harboring QTL is a prerequisite for that. QTL mapping in potato (Solanum tuberosum) has identified a region on chromosome V tagged by DNA markers GP21 and GP179, which contains a number of important QTL, among others QTL for resistance to late blight caused by the oomycete Phytophthora infestans and to root cyst nematodes.

Results

To obtain genomic sequence for the targeted region on chromosome V, two local BAC (bacterial artificial chromosome) contigs were constructed and sequenced, which corresponded to parts of the homologous chromosomes of the diploid, heterozygous genotype P6/210. Two contiguous sequences of 417,445 and 202,781 base pairs were assembled and annotated. Gene-by-gene co-linearity was disrupted by non-allelic insertions of retrotransposon elements, stretches of diverged intergenic sequences, differences in gene content and gene order. The latter was caused by inversion of a 70 kbp genomic fragment. These features were also found in comparison to orthologous sequence contigs from three homeologous chromosomes of Solanum demissum, a wild tuber bearing species. Functional annotation of the sequence identified 48 putative open reading frames (ORF) in one contig and 22 in the other, with an average of one ORF every 9 kbp. Ten ORFs were classified as resistance-gene-like, 11 as F-box-containing genes, 13 as transposable elements and three as transcription factors. Comparing potato to Arabidopsis thaliana annotated proteins revealed five micro-syntenic blocks of three to seven ORFs with A. thaliana chromosomes 1, 3 and 5.

Conclusion

Comparative sequence analysis revealed highly conserved collinear regions that flank regions showing high variability and tandem duplicated genes. Sequence annotation revealed that the majority of the ORFs were members of multiple gene families. Comparing potato to Arabidopsis thaliana annotated proteins suggested fragmented structural conservation between these distantly related plant species.

Improvement of whole-genome annotation of cereals through comparative analyses

Rice is an important model species for the Poaceae and other monocotyledonous plants. With the availability of a near-complete, finished, and annotated rice genome, we performed genome level comparisons between rice and all plant species in which large genomic or transcriptomic data sets are available to determine the utility of cross-species sequence for structural and functional annotation of the rice genome. Through comparative analyses with four plant genome sequence data sets and transcript assemblies from 185 plant species, we were able to confirm and improve the structural annotation of the rice genome. Support for 38,109 (89.3%) of the total 42,653 nontransposable element-related genes in the rice genome in the form of a rice expressed sequence tag, full-length cDNA, or plant homolog from our comparative analyses could be found. Although the majority of the putative homologs were obtained from Poaceae species, putative homologs were identified in dicotyledonous angiosperms, gymnosperms, and other plants such as algae, moss, and fern. A set of rice genes (7669) lacking a putative homolog was identified which may be lineage-specific genes that evolved after speciation and have a role in species diversity. Improvements to the current rice gene structural annotation could be identified from our comparative alignments and we were able to identify 487 genes which were mostly likely missed in the current rice genome annotation and another 500 genes for structural annotation review. We were able to demonstrate the utility of cross-species comparative alignments in the identification of noncoding sequences and in confirmation of gene nesting in rice.

The TIGR Rice Genome Annotation Resource: improvements and new features

In The Institute for Genomic Research Rice Genome Annotation project (), we have continued to update the rice genome sequence with new data and improve the quality of the annotation. In our current release of annotation (Release 4.0; January 12, 2006), we have identified 42 653 non-transposable element-related genes encoding 49 472 gene models as a result of the detection of alternative splicing. We have refined our identification methods for transposable element-related genes resulting in 13 237 genes that are related to transposable elements. Through incorporation of multiple transcript and proteomic expression data sets, we have been able to annotate 24 799 genes (31 739 gene models), representing ∼50% of the total gene models, as expressed in the rice genome. All structural and functional annotation is viewable through our Rice Genome Browser which currently supports 59 tracks. Enhanced data access is available through web interfaces, FTP downloads and a Data Extractor tool developed in order to support discrete dataset downloads.

Molecular plant breeding: achievements in green biotechnology and future perspectives

Since one decade ago, transgenic crop plants are globally grown; in 2004, it was estimated to cover a total of 81 Mio ha in 17 countries. At present, four plant species (soybean, maize, cotton and rapeseed) dominate with two traits (herbicide tolerance and insect resistance). The traits on which research concentrates and the constructs which might come next onto the market are outlined. The procedure on how to clone such genes of interest, e.g. via map-based cloning, and some other helpful approaches of green biotechnology, like high throughput techniques and functional markers, are summarised, and a rough calculation about the market value of transgenic crops in US dollars is quoted.

Computational identification and systematic analysis of the ACR gene family in Oryza sativa

Based on sequence similarity search and domain detection, nine ACT domain repeat protein-coding genes (the "ACR" genes) in rice were identified, which were mainly distributed on the chromosomes 2, 3, 4, and 8. An InterPro database search indicated that four copies of the ACT domain linearly occupied the entire polypeptide. The first three ACT domains were linked by two different sequences. However, the fourth ACT domain was extremely close to ACT3. Gene structure comparisons showed large differences in exon numbers, from three to eight, among members of the rice ACR gene family. In addition, it appeared that gene duplication might be operative when the compositions of exons and introns were analyzed. Phylogenetic analysis divided the ACR gene family into five distinct groups, and this division was generally according to the expression patterns of the ACR genes. The Arabidopsis and rice ACR proteins were clustered across together, suggesting that these ACR genes might originate from an ancient common ancestor. Notably, the identification of orthologues and paralogues would be useful for rice gene functional annotation.

Uneven distribution of expressed sequence tag loci on maize pachytene chromosomes

Examining the relationships among DNA sequence, meiotic recombination, and chromosome structure at a genome-wide scale has been difficult because only a few markers connect genetic linkage maps with physical maps. Here, we have positioned 1195 genetically mapped expressed sequence tag (EST) markers onto the 10 pachytene chromosomes of maize by using a newly developed resource, the RN-cM map. The RN-cM map charts the distribution of crossing over in the form of recombination nodules (RNs) along synaptonemal complexes (SCs, pachytene chromosomes) and allows genetic cM distances to be converted into physical micrometer distances on chromosomes. When this conversion is made, most of the EST markers used in the study are located distally on the chromosomes in euchromatin. ESTs are significantly clustered on chromosomes, even when only euchromatic chromosomal segments are considered. Gene density and recombination rate (as measured by EST and RN frequencies, respectively) are strongly correlated. However, crossover frequencies for telomeric intervals are much higher than was expected from their EST frequencies. For pachytene chromosomes, EST density is about fourfold higher in euchromatin compared with heterochromatin, while DNA density is 1.4 times higher in heterochromatin than in euchromatin. Based on DNA density values and the fraction of pachytene chromosome length that is euchromatic, we estimate that approximately 1500 Mbp of the maize genome is in euchromatin. This overview of the organization of the maize genome will be useful in examining genome and chromosome evolution in plants.

An analytical study of the number of steady states in gene regulatory networks

Herein, I proposed a model for gene networks and studied the steady states in the dynamics both by numerical and analytical methods. In this model, mRNA and protein levels change continuously with time; each gene alters transcriptional regulation depending on the concentration of transcription factors. The dynamical behavior of continuous model is quite complex and different from that of the discrete model, the Boolean network. Large portion of steady states of this model can be classified into three types. The rough structure of gene interactions, which corresponds to Boolean function, is sufficient to predict the expression level of each gene in these types of steady states. I also determined the expected numbers of two major types of steady states observed in a randomly generated gene network. The results obtained from these formulae contradict previously accepted belief. The results are that neither gene number nor connectivity between genes increases the expected number of steady states in a random gene network. The number of steady states is very small. The number of self-regulatory genes, however, effectively increases the number of steady states in a network. These results imply that increases in gene number may not be the direct driving force for the evolution of a variety of different cell types within organisms. Instead, the number of self-regulatory genes may significantly increase cellular diversity.

Quantitative trait loci and the study of plant domestication

Plant domestication ranks as one of the most important developments in human history, giving human populations the potential to harness unprecedented quantities of the earth's resources. But domestication has also played a more subtle historical role as the foundation of the modern study of evolution and adaptation. Until recently, however, researchers interested in domestication were limited to studying phenotypic changes or the genetics of simple Mendelian traits, when often the characters of most interest--fruit size, yield, height, flowering time, etc.--are quantitative in nature. The goals of this paper are to review some of the recent work on the quantitative genetics of plant domestication, identify some of the common trends found in this literature, and offer some novel interpretations of the data that is currently available.

Enhancing gene targeting efficiency in higher plants: rice is on the move

Meeting the challenge of routine gene targeting (GT) in higher plants is of crucial interest to researchers and plant breeders who are currently in need of a powerful tool to specifically modify a given locus in a genome. Higher plants have long been considered the last lineage resistant to targeting technology. However, a recent report described an efficient method of T-DNA-mediated targeted disruption of a non-selectable locus in rice [Terada et al., Nat Biotechnol 20: 1030-1034 (2002)]. Though this study was an obvious breakthrough, further improvement of GT frequencies may derive from a better understanding of the natural mechanisms that control homologous recombination (HR) processes. In this review, we will focus on what is known about HR and the factors which may hamper the development of routine GT by HR in higher plants. We will also present the current strategies envisaged to overcome these limitations, such as expression of recombination proteins and refinements in the design of the transformation vector.

Development and mapping of EST-derived simple sequence repeat markers for hexaploid wheat

Expressed sequence tags (ESTs) are a valuable source of molecular markers. To enhance the resolution of an existing linkage map and to identify putative functional polymorphic gene loci in hexaploid wheat (Triticum aestivum L.), over 260,000 ESTs from 5 different grass species were analyzed and 5418 SSR-containing sequences were identified. Using sequence similarity analysis, 156 cross-species superclusters and 138 singletons were used to develop primer pairs, which were then tested on the genomic DNA of barley (Hordeum vulgare), maize (Zea mays), rice (Oryza sativa), and wheat. Three-hundred sixty-eight primer pairs produced PCR amplicons from at least one species and 227 primer pairs amplified DNA from two or more species. EST-SSR sequences containing dinucleotide motifs were significantly more polymorphic (74%) than those containing trinucleotides (56%), and polymorphism was similar for markers in both coding and 5' untranslated (UTR) regions. Out of 112 EST-SSR markers, 90 identified 149 loci that were integrated into a reference wheat genetic map. These loci were distributed on 19 of the 21 wheat chromosomes and were clustered in the distal chromosomal regions. Multiple-loci were detected by 39% of the primer pairs. Of the 90 mapped ESTs, putative functions for 22 were identified using BLASTX queries. In addition, 80 EST-SSR markers (104 loci) were located to chromosomes using nullisomic-tetrasomic lines. The enhanced map from this study provides a basis for comparative mapping using orthologous and PCR-based markers and for identification of expressed genes possibly affecting important traits in wheat.

Comparative mapping in the Pinaceae

A comparative genetic map was constructed between two important genera of the family Pinaceae. Ten homologous linkage groups in loblolly pine (Pinus taeda L.) and Douglas fir (Pseudotsuga menziesii [Mirb.] Franco) were identified using orthologous expressed sequence tag polymorphism (ESTP) and restriction fragment length polymorphism (RFLP) markers. The comparative mapping revealed extensive synteny and colinearity between genomes of the Pinaceae, consistent with the hypothesis of conservative chromosomal evolution in this important plant family. This study reports the first comparative map in forest trees at the family taxonomic level and establishes a framework for comparative genomics in Pinaceae.

EST derived SSR markers for comparative mapping in wheat and rice

Structural and functional relationships between the genomes of hexaploid wheat ( Triticum aestivum L.) (2n=6x=42) and rice (Oryza sativa L.) (2n=2x=24) were evaluated using linkage maps supplemented with simple sequence repeat (SSR) loci obtained from publicly available expressed sequence tags (ESTs). EST-SSR markers were developed using two main strategies to design primers for each gene: (1) primer design for multiple species based on supercluster analysis, and (2) species-specific primer design. Amplification was more consistent using the species-specific primer design for each gene. Forty-four percent of the primers designed specifically for wheat sequences were successful in amplifying DNA from both species. Existing genetic linkage maps were enhanced for the wheat and rice genomes using orthologous loci amplified with 58 EST-SSR markers obtained from both wheat and rice ESTs. The PCR-based anchor loci identified by these EST-SSR markers support previous patterns of conservation between wheat and rice genomes; however, there was a high frequency of interrupted colinearity. In addition, multiple loci amplified by these primers made the comparative analysis more difficult. Enhanced comparative maps of wheat and rice provide a useful tool for interpreting and transferring molecular, genetic, and breeding information between these two important species. These EST-SSR markers are particularly useful for constructing comparative framework maps for different species, because they amplify closely related genes to provide anchor points across species.

Structural and functional analysis of rice genome

Rice is an excellent system for plant genomics as it represents a modest size genome of 430 Mb. It feeds more than half the population of the world. Draft sequences of the rice genome, derived by whole-genome shotgun approach at relatively low coverage (4-6 X), were published and the International Rice Genome Sequencing Project (IRGSP) declared high quality (>10 X), genetically anchored, phase 2 level sequence in 2002. In addition, phase 3 level finished sequence of chromosomes 1, 4 and 10 (out of 12 chromosomes of rice) has already been reported by scientists from IRGSP consortium. Various estimates of genes in rice place the number at >50,000. Already, over 28,000 full-length cDNAs have been sequenced, most of which map to genetically anchored genome sequence. Such information is very useful in revealing novel features of macro- and micro-level synteny of rice genome with other cereals. Microarray analysis is unraveling the identity of rice genes expressing in temporal and spatial manner and should help target candidate genes useful for improving traits of agronomic importance. Simultaneously, functional analysis of rice genome has been initiated by marker-based characterization of useful genes and employing functional knock-outs created by mutation or gene tagging. Integration of this enormous information is expected to catalyze tremendous activity on basic and applied aspects of rice genomics.

Physical maps and recombination frequency of six rice chromosomes

We constructed physical maps of rice chromosomes 1, 2, and 6-9 with P1-derived artificial chromosome (PAC) and bacterial artificial chromosome (BAC) clones. These maps, with only 20 gaps, cover more than 97% of the predicted length of the six chromosomes. We submitted a total of 193 Mbp of non-overlapping sequences to public databases. We analyzed the DNA sequences of 1316 genetic markers and six centromere-specific repeats to facilitate characterization of chromosomal recombination frequency and of the genomic composition and structure of the centromeric regions. We found marked changes in the relative recombination rate along the length of each chromosome. Chromosomal recombination at the centromere core and surrounding regions on the six chromosomes was completely suppressed. These regions have a total physical length of about 23 Mbp, corresponding to 11.4% of the entire size of the six chromosomes. Chromosome 6 has the longest quiescent region, with about 5.6 Mbp, followed by chromosome 8, with quiescent region about half this size. Repetitive sequences accounted for at least 40% of the total genomic sequence on the partly sequenced centromeric region of chromosome 1. Rice CentO satellite DNA is arrayed in clusters and is closely associated with the presence of Centromeric Retrotransposon of Rice (CRR)- and RIce RetroElement 7 (RIRE7)-like retroelement sequences. We also detected relatively small coldspot regions outside the centromeric region; their repetitive content and gene density were similar to those of regions with normal recombination rates. Sequence analysis of these regions suggests that either the amount or the organization patterns of repetitive sequences may play a role in the inactivation of recombination.

Functional markers in plants

Different approaches (including association studies) have recently been adopted for the functional characterization of allelic variation in plants and to identify sequence motifs affecting phenotypic variation. We propose the term 'functional markers' for DNA markers derived from such functionally characterized sequence motifs. Functional markers are superior to random DNA markers such as RFLPs, SSRs and AFLPs owing to complete linkage with trait locus alleles. We outline the definition, development, application and prospects of functional markers.