Orthology between genomes of Brachypodium, wheat and rice
© Gupta et al; licensee BioMed Central Ltd. 2009
Received: 09 March 2009
Accepted: 27 May 2009
Published: 27 May 2009
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.
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.
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.
Cereals constitute the most important group of cultivated plants, and are known to have diverged from a common paleopolyploid ancestor ~45–47 million years ago (Mya) . Despite this, a remarkable overall structural and functional similarity exists among different cereal genomes [2, 3], although the size of these genomes differs greatly, ranging from 430 Mb in rice (Oryza sativa) to 16,000 Mb in hexaploid wheat (Triticum aestivum). Due to its small size and availability of whole genome sequence, rice has been used as a model system for a variety of experimental studies including map-based cloning . However, recent studies resolved further the dynamic changes in rice genome sequences, thus questioning the utility of rice as a model crop , and necessitating the need for search of a more efficient model system.
Brachypodium distachyon, a small temperate grass (sub-family Pooideae) has recently emerged as a better model system for the study of temperate grasses. This is particularly, due to several of its desirable biological features and its phylogenetic position [6, 7]. It is postulated that relative to rice genome, Brachypodium genome will exhibit a much higher level of colinearity and synteny to the genomes of temperate cereal crops. In the present study, the available Brachypodium EST contigs (bEST contigs) and supercontigs were utilized to explore further the utility of the Brachypodium genome as a model for carrying out comparative genomics studies in cereals in general, and for wheat genomics in particular. The relationship of Brachypodium genome with wheat and rice genomes has been examined for this purpose, and improved criteria of sequence similarity search were used for more accurate estimation of similarity .
In the present study, EST sequences from Brachypodium were utilized to find out the degree of similarity of Brachypodium genome with EST/genomic sequences of wheat and rice. The orthologous wheat sequences thus identified were also utilized to study the relationship of wheat genome sequences with Brachypodium supercontigs. We have also taken note of the comparisons of chloroplast genomes among eight grass species, which were included in the report on Brachypodium chloroplast genome sequence that was recently worked out .
Orthology between Brachypodium and wheat
As many as 3,818 B. distachyon EST contigs were blasted (BLASTN) against the available wheat EST contigs (containing bin-mapped wESTs) to identify matching wESTs. The analysis revealed that as many as 449 bEST contigs had orthologs in wheat genome.
Analysis of mapped wEST contigs that matched bEST contigs
Distribution of the orthologous loci according to their assignment to wheat chromosomes arranged in two-way classification
Of the above 449 matched wEST contigs (orthologous to bEST contigs), 77 (17.2%) represented unique loci, and the remaining 372 (82.8%) detected multiple loci with 283 (76.1%) having multiple loci on homoeologous chromosomes and 89 (23.9%) having multiple loci on non-homoeologous chromosomes.
Of the 1,154 orthologous loci with known positions on wheat chromosomes, 1,094 (94.8%) loci were known to have earlier been assigned to 159 chromosome bins defined by deletion break points. The remaining 60 (5.2%) loci could be assigned only to individual chromosomes or their arms. A maximum of 386 loci (35.3%) were mapped in the proximal regions (60% of the arm length from centromere; C-0.60) followed by 331 loci (30.3%) mapped to the distal regions (40% terminal arm length; 0.60–1.00). The remaining 377 loci (34.4%) were mapped to the interstitial bins having proximal and distal regions.
Homology between Brachypodium supercontigs and homoeologous groups of wheat
Wheat homoeologous group
bEST mapped contig
Number of Brachypodium supercontigs showing orthology*
Supercontig showing maximum homology
Orthology between bEST contigs and rice genome sequences
Distribution of the orthologous loci on individual rice chromosome
Number of loci
Conserved orthologous sequences among Brachypodium, wheat and rice
Conservation of SSRs among the three genomes
The 183 bEST contig sequences shared by three species (Brachypodium, wheat and rice) were also used for mining SSRs. A total of 100 (54.6%) bEST contigs contained 137 SSRs. As many as 45 of these SSRs showed conservation in wheat and 23 of these SSRs showed conservation in rice. As many as 17 SSRs were conserved across all the three species.
Transferability of conserved orthologous SSRs
List of primers for the 12 conserved SSRs used for wet-lab experiment
Brachypodium contig Id.
Primer sequence 5'-3'
Product size (bp)
Putative uncharac-terized protein
Putative 60S ribos-omal protein L13E
L CCGCACCTCAAGGACTACA R TCGGAGGAGATCTTGGTGAG
Putative uncharac-terized protein
Putative ribosomal protein
Putative uncharac-terized protein
Protein disulfide isomerase 2 precursor
Cytosolic 6-phosphogluconate dehydrogenase
Comparative genomics among grasses initially focused on the analysis of colinearity (gene order) and synteny (gene content) among DNA markers mapped on individual chromosome at a low resolution (10 cM). This led to the identification of 30 rice-independent linkage blocks involved in the constitution of all cereal genomes and allowed identification of a number of rearrangements within individual genomes . However, due to the availability of whole genome sequence of rice, and substantial partial sequences from other cereal genomes, emphasis shifted to a comparison of nucleic acid sequences. In particular, sequences of ~7000 bin-mapped wESTs were aligned with rice genome sequences , allowing improved resolution and discovery of many more rearrangements.
Although rice worked well as a model for all grasses including wheat, and generated useful information, Brachypodium, belonging to subfamily Pooideae (wheat also belongs to Pooideae), is proposed as a better model than rice (subfamily Ehrhartoideae). Recent studies have suggested that relative to rice, Barchypodium is more closely related to wheat and barley and the colinearity between Barchypodium and wheat is better than that between wheat and rice [14, 15]. Chloroplast sequence-based phylogenetic analysis in eight grass species also suggested that Brachypodium is closer to the tribe Triticeae . The possible estimated time of divergence between Brachypodium and Triticeae is also shorter (35 Mya) than that of divergence between wheat and rice (50 Mya)  thus supporting the view that Brachypodium is more closely related with the members of Triticeae.
During the present study, orthologous relationship among bEST contigs, wEST contigs and rice genome sequences was studied using improved criteria of sequence comparison. Observation of higher number of bEST contigs showing orthology with rice genome was mainly attributed to the fact that only a small fraction of wheat genome (0.02%) and almost complete rice genome (95%) were used for sequence comparison with the available Brachypodium EST contigs. If we take into account the proportion of the genome used for comparison, it may be concluded that wheat has higher level of orthology with Brachypodium than with rice.
The mapped loci in different deletion bins of a particular chromosome of wheat matched with same or different supercontigs of Brachypodium. For instance, wheat group 4 chromosomes are highly syntenic to Brachypodium super_1 contig (54.5%) than to other supercontigs, although super_1 contig showed homology with other homoeologous groups also. The mapping information of these Brachypodium supercontigs on individual Brachypodium chromosomes will be useful for developing markers specific to the targeted regions of wheat chromosomes.
It was also observed that although D sub-genome of wheat is smaller in size, the orthologous loci mapped on this sub-genome are no fewer than those mapped on sub-genome B, suggesting closer relationship between Brachypodium and Aegilops tauschii, the donor of the D sub-genome of hexaploid wheat.
The relative abundance of orthologous loci on proximal regions of chromosome arms in wheat is in agreement with the earlier studies in wheat and rice . It seems that higher degree of sequence conservation coincides with the low recombination proximal regions, which is understandable, since high recombination in terminal regions will cause reshuffling of genes during evolution .
The results of the present study indicate that the availability of whole genome sequence of Brachypodium will be of enormous relevance for comparative genomics, gene annotation and evolutionary, structural and functional genomic studies of large genomes of the Triticeae.
Brachypodium, wheat, rice sequence databases
A total of 3,818 Brachypodium EST (bEST) contigs, and a set of 1,015 supercontigs representing 4× coverage of Brachypodium genome, were available in public domain [19, 20]. As many as 3,792 wheat EST (wEST) contigs containing bin-mapped wESTs were available at GrainGenes 2.0  and rice genomic sequences were available at Gramene .
In order to find orthology among Brachypodium, wheat and rice genomes, bEST contigs were blasted against wEST contigs and rice genomic sequences. The pairwise sequence alignment in BLASTN search was improved by using three new parameters . The first parameter, aligned length (AL), corresponds to the sum of the lengths of all the high-scoring segment pairs (HSPs) in a single hit. Second parameter, cumulative identity percentage (CIP) was obtained from the formula, CIP = [Σ Id of HSPs/AL] × 100 and the third parameter, cumulative alignment length percentage (CALP) was calculated as follows: CALP = [AL/QL] × 100, where, QL is the length of query sequence. Last two parameters (CIP and CALP) allow estimation of highest similarity between sequences over the entire length of query sequence. These parameters were applied to all the BLASTN results and values of 60% CIP and 70% CALP were used for identification of orthologs of Brachypodium genomic sequences in wheat (through ESTs) and rice genomes.
Mapping of wheat and rice orthologs
The physical positions of wEST orthologs identified through sequence comparisons were localized to specific bins of wheat chromosomes based on the information about mapped wEST sequences . The rice genomic sequences, which were orthologous to bEST contigs, were also known and were physically localized to specific sites on 12 different rice chromosomes with the help of KaryoView program . The χ2 test for goodness-of-fit was used for testing the random distribution of ortholoci in wheat genome at the level of the three sub-genomes, the seven homoeologous groups, the 21 chromosomes and the 42 chromosome arms. The same was done for the 12 chromosomes of rice.
Assignment of putative function to orthologs
The orthologous sequences belonging to the three genomes (Brachypodium, wheat and rice) were subjected to BLASTX analysis against non-redundant protein database  for assigning putative functions at a cut-off E value of 10-30.
Identification of SSRs in orthologs
The orthologous sequences available in all the three genomes were mined for simple sequence repeats (SSRs) using SSRIT program . The SSRs with a repeat motif of 2–6 nucleotides and a length of ≥ 12 bp were included in the analysis. Primers were designed for the 12 conserved SSRs using PRIMER3 .
Primers for 12 conserved Brachypodium SSRs were synthesized from Invitrogen, USA. PCR was performed separately using the genomic DNA of Brachypodium, wheat and rice in a final volume of 20 μl in an Applied Biosystems 'Veriti Thermal Cycler'. After electrophoresis, polyacrylamide gels were silver stained following Tegelstrom .
This work was supported by the Department of Biotechnology (DBT) and the Department of Science & Technology (DST), Government of India, New Delhi and the Indian National Science Academy (INSA), New Delhi. The support was also received from DST through its FIST-programme and from University Grants Commission (UGC), New Delhi through its SAP-DRS programme. DNA aliquot of B. distachyon was kindly provided by Dr. Azhaguvel Perumal, Texas Agrilife Research Amarillo, TX, USA. Aakash Goyal carried out a part of sequence analysis.
- Paterson AH, Bowers JE, Chapman BA: Ancient polyploidization predating divergence of cereals, and its consequences for comparative genomics. Proc Natl Acad Sci USA. 2004, 101: 9903-9908. 10.1073/pnas.0307901101.PubMed CentralView ArticlePubMedGoogle Scholar
- Gale MD, Devos KM: Plant comparative genetics after 10 years. Science. 1998, 282: 656-659. 10.1126/science.282.5389.656.View ArticlePubMedGoogle Scholar
- Feuillet C, Keller B: Comparative genomics in the grass family: Molecular characterization of grass genome structure and evolution. Annl Bot. 2002, 89: 3-10. 10.1093/aob/mcf008.View ArticleGoogle Scholar
- Xu Y, McCouch SR, Zhang Q: How can we use genomics to improve cereals with rice as a reference genome?. Plant Mol Biol. 2005, 59: 7-26. 10.1007/s11103-004-4681-2.View ArticlePubMedGoogle Scholar
- Gaut BS: Evolutionary dynamics of grass genomes. New Phytol. 2002, 154: 15-28. 10.1046/j.1469-8137.2002.00352.x.View ArticleGoogle Scholar
- Garvin DF: Brachypodium: a new monocot model plant system emerges. J Sci Food Agric. 2007, 87: 1177-1179. 10.1002/jsfa.2868.View ArticleGoogle Scholar
- Huo N, Lazo GR, Vogel JP, You FM, Ma Y, Hayden DM, Coleman-Derr D, Hill TA, Dvorak J, Anderson OD, Luo MC, Gu YQ: The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences. Funct Integr Genomics. 2008, 8: 135-147. 10.1007/s10142-007-0062-7.View ArticlePubMedGoogle Scholar
- Salse J, Bolot S, Throude M, Jouffe V, Piegu B, Quraishi UM, Calcagno T, Cooke R, Delseny M, Feuillet C: Identification and characterization of shared duplications between rice and wheat provide new insight into grass genome evolution. Plant Cell. 2008, 20: 11-24. 10.1105/tpc.107.056309.PubMed CentralView ArticlePubMedGoogle Scholar
- Bortiri E, Coleman-Derr D, Lazo GR, Anderson OD, Gu YQ: The complete chloroplast genome sequence of Brachypodium distachyon: sequence comparison and phylogenetic analysis of eight grass plastomes. BMC Res Notes. 2008, 1: 61-10.1186/1756-0500-1-61.PubMed CentralView ArticlePubMedGoogle Scholar
- Gill BS, Friebe B, Endo TR: Standard karyotype and nomenclature system for description of chromosome bands and structural aberrations in wheat (Triticum aestivum). Genome. 1991, 34: 830-839.View ArticleGoogle Scholar
- Rice Pseudomolecules Information. [http://rice.plantbiology.msu.edu/pseudomolecules/info.shtml]
- Bolot S, Abrouk M, Masood-Quraishi U, Stein N, Messing J, Feuillet C, Salse J: The 'inner circle' of the cereal genome. Curr Opin Plant Biol. 2008, 12: 1-7.View ArticleGoogle Scholar
- Sorrells ME, La Rota M, Bermudez-Kandianis CE, Greene RA, Kantety R, Munkvold JD, Miftahudin , Mahmoud A, Ma X, Gustafson PJ, Qi LL, Echalier B, Gill BS, Matthews DE, Lazo GR, Chao S, Anderson OD, Edwards H, Linkiewicz AM, Dubcovsky J, Akhunov ED, Dvorak J, Zhang D, Nguyen HT, Peng J, Lapitan NLV, Gonzalez-Hernandez JL, Anderson JA, Hossain K, Kalavacharla V, Kianian SF, Choi D-W, Close TJ, Dilbirligi M, Gill KS, Steber C, Walker-Simmons MK, McGuire PE, Qualset CO: Comparative DNA sequence analysis of wheat and rice genomes. Genome Res. 2003, 13: 1818-1827.PubMed CentralPubMedGoogle Scholar
- Huo N, Gu YQ, Lazo GR, Vogel JP, Coleman-Derr D, Luo M-C, Thilmony R, Garvin DF, Anderson OD: Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics. Genome. 2006, 49: 1099-1108. 10.1139/G06-087.View ArticlePubMedGoogle Scholar
- Faris JD, Zhang Z, Fellers JP, Gill BS: Micro-colinearity between rice, Brachypodium and T. monococcum at the wheat domestication locus Q. Funct Integr Genomics. 2008, 8: 149-164. 10.1007/s10142-008-0073-z.View ArticlePubMedGoogle Scholar
- Bossolini E, Wicker T, Knobel PA, Keller B: Comparison of orthologous loci from small grass genomes Brachypodium and rice: implications for wheat genomics and grass genome annotation. Plant J. 2007, 49: 704-717. 10.1111/j.1365-313X.2006.02991.x.View ArticlePubMedGoogle Scholar
- Akhunov ED, Goodyear AW, Geng S, Qi LL, Echalier B, Gill BS, Miftahudin , Gustafson JP, Lazo G, Chao S, Anderson OD, Linkiewicz AM, Dubcovsky J, La Rota M, Sorrells ME, Zhang D, Nguyen HT, Kalavacharla V, Hossain K, Kianian SF, Peng J, Lapitan NLV, Gonzalez-Hernandez JL, Anderson JA, Choi D-W, Close TJ, Dilbirligi M, Gill KS, Walker-Simmons MK, Steber C, McGuire PE, Qualset CQ, Dvorak J: The organization and rate of evolution of wheat genomes are correlated with recombination rates along chromosome arms. Genome Res. 2003, 13: 753-763. 10.1101/gr.808603.PubMed CentralView ArticlePubMedGoogle Scholar
- See DR, Brooks S, Nelson JC, Brown-Guedira G, Friebe B, Gill BS: Gene evolution at the ends of wheat chromosomes. Proc Natl Acad Sci USA. 2006, 103: 4162-4167. 10.1073/pnas.0508942102.PubMed CentralView ArticlePubMedGoogle Scholar
- Genomics and Gene Discovery bEST Resource. [http://brachypodium.pw.usda.gov/]
- Brachypodium BLAST Server. [http://blast.brachybase.org/]
- GrainGenes 2.0 Basic local Alignment Search Tool. [http://wheat.pw.usda.gov/GG2/blast.shtml]
- Gramene Multi Basic local Alignment Search Tool. [http://www.gramene.org/Multi/blastview]
- Bin Mapped wEST Sequences Information. [http://wheat.pw.usda.gov/cgi-bin/westsql/maplocus.cgi]
- Oryza sativa japonica KaryoView. [http://www.gramene.org/Oryza_sativa_japonica/karyoview]
- NCBI Basic Local Alignment Search Tool. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]
- Simple Sequence Repeat Identification Tool. [http://www.gramene.org/db/searches/ssrtool]
- Program to Pick Primers from a DNA Sequence. [http://fokker.wi.mit.edu/cgi-bin/primer3/primer3.cgi]
- Tegelstrom H: Detection of mitochondrial DNA fragments. The Molecular Genetic Analysis of Populations: A Practical Approach. Edited by: Hoelzel AR. 1992, Oxford, IRL Press, 89-114.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.