- Research article
- Open Access
The sequence of rice chromosomes 11 and 12, rich in disease resistance genes and recent gene duplications
BMC Biologyvolume 3, Article number: 20 (2005)
Rice is an important staple food and, with the smallest cereal genome, serves as a reference species for studies on the evolution of cereals and other grasses. Therefore, decoding its entire genome will be a prerequisite for applied and basic research on this species and all other cereals.
We have determined and analyzed the complete sequences of two of its chromosomes, 11 and 12, which total 55.9 Mb (14.3% of the entire genome length), based on a set of overlapping clones. A total of 5,993 non-transposable element related genes are present on these chromosomes. Among them are 289 disease resistance-like and 28 defense-response genes, a higher proportion of these categories than on any other rice chromosome. A three-Mb segment on both chromosomes resulted from a duplication 7.7 million years ago (mya), the most recent large-scale duplication in the rice genome. Paralogous gene copies within this segmental duplication can be aligned with genomic assemblies from sorghum and maize. Although these gene copies are preserved on both chromosomes, their expression patterns have diverged. When the gene order of rice chromosomes 11 and 12 was compared to wheat gene loci, significant synteny between these orthologous regions was detected, illustrating the presence of conserved genes alternating with recently evolved genes.
Because the resistance and defense response genes, enriched on these chromosomes relative to the whole genome, also occur in clusters, they provide a preferred target for breeding durable disease resistance in rice and the isolation of their allelic variants. The recent duplication of a large chromosomal segment coupled with the high density of disease resistance gene clusters makes this the most recently evolved part of the rice genome. Based on syntenic alignments of these chromosomes, rice chromosome 11 and 12 do not appear to have resulted from a single whole-genome duplication event as previously suggested.
Rice (Oryza sativa) is a major staple food and is consumed by nearly half the world's population. It accounts for more than 21% of global human per capita energy and 15% of per capita protein. In the past few decades, although rice production has doubled due to the introduction of high yielding varieties/hybrids and improved cultivation practices, it is still insufficient to cope with ever-increasing global demand, which is expected to increase at the rate of about 1% per annum . At the same time, inappropriate natural resource use, along with biotic and abiotic stress pressure, is casting a shadow on rice production . Access to the rice genome sequence will enable identification of genes responsible for traits and alleles that will be essential to meet the growing demands of food production in the coming years. Towards this end, whole-genome shotgun-based draft sequences of the indica and japonica subspecies of rice were reported previously [3–5] while the International Rice Genome Sequencing Project (IRGSP), using a clone-by-clone approach, focused on generating a high quality, finished sequence of the rice genome . Indeed, access to the rice genome has served as a catalyst for investigations on comparative genomics, functional genomics, map-based gene cloning and molecular breeding in rice [7, 8].
Rice is one of the cereals in the Poaceae family, which collectively provides the largest source of calories for human consumption. Within the cereals are larger genome species such as maize, wheat, millet and sorghum . Synteny in the cereals has been reported previously using molecular markers  and recently at a higher resolution using sequences available from the rice, sorghum, maize, and wheat genomes [3, 11, 14]. Thus, access to a complete, high-resolution rice genome sequence will facilitate research on other cereals with larger, partially sequenced genomes.
One of the emerging features of plant genomes appears to be the recent generation of gene copies that have diverged in regulation and function . A common pathway for such a mechanism in plants is whole-genome duplication (WGD), which occurred for instance in maize as little as 5 million years ago . Other mechanisms involve tandem gene amplifications and segmental duplications. In rice, chromosomes 11 and 12 provide such an example. In addition, several genes of agronomic importance, such as blast, bacterial blight, virus and insect resistance, photoperiod-sensitive male-sterility as well as salt tolerance, have been mapped on these two chromosomes . Here, we report the in-depth analysis of chromosomes 11 and 12 of a japonica cultivar of rice with genetic lengths estimated to be 118.6 cM and 110.1 cM, respectively . We have generated high-quality finished sequences for these two chromosomes, annotated the chromosomes for genes and other features, and used these data to examine novel features such as the organization of disease resistance genes. We also asked whether chromosomes 11 and 12 are likely to have resulted from a WGD event by examining duplication events between the two chromosomes through alignments of mapped genes of wheat along the rice chromosomes.
Results and discussion
General features of chromosomes 11 and 12
A total of 255 and 269 BAC/PAC clones were sequenced from chromosomes 11 and 12, respectively. Chromosome 11 is slightly longer than chromosome 12, with 28.4 Mb and 27.5 Mb, respectively, of sequence identified for them (Table 1, Fig. 1). Excluding the telomeres, a few physical gaps remain on the two pseudomolecules; six and one on chromosomes 11 and 12, respectively (Accession Numbers: chromosome 11: DP000010, chromosome 12: DP000011). This was due to a difficulty in obtaining clones that span these regions, especially the centromeres and telomeres, which present technical challenges in sequencing highly repetitive tracts of the genome.
Of the 4,286 and 4,169 genes identified on chromosomes 11 and 12, respectively, 3,148 and 2,845 could be annotated as non-transposable element related (non-TE; Fig. 1). A similar percentage of genes on the two chromosomes could be assigned a putative function or annotated as encoding an expressed protein, leaving a similar percentage annotated as encoding a hypothetical protein (Table 1). With respect to domain composition, there was a striking difference in representation of some Pfam domains between the chromosomes, suggesting the presence of genes coding for different suites of proteins on the two chromosomes. As shown in Tables 2 and 3, there was a large enrichment of proteins on chromosome 11 (166 proteins) containing a leucine-rich repeat (LRR) domain, which is involved in protein-protein recognition and is a hallmark of disease resistance genes (see below). Besides the LRR-domain, there are two other domains common to disease resistance genes, protein kinase and NB-ARC domains, which were also enriched on chromosome 11 relative to chromosome 12; 106 versus 74 proteins containing protein kinase domains and 102 versus 49 proteins with NB-ARC domains on chromosome 11 versus 12, respectively. Based on alignments by a rice transcript, 1,235 (39%) and 1,221 (43%) of non-TE-related genes on chromosomes 11 and 12, respectively, were active (Additional figure 10 [see Additional file 1]). Of these, 952 and 980 non-TE-related genes from chromosomes 11 and 12, respectively, could be aligned with a FL-cDNA. In addition, 152 non-TE-related genes aligned with ESTs from other monocots on each chromosome. For 26 genes from chromosome 11 and 23 from chromosome 12, only non-monocot evidence was available for expression (Additional figure 10 [see Additional file 1]). A ready set of mutants is available for 800 and 845 genes on chromosomes 11 and 12, respectively, as evidenced by a flanking sequence tag (FST) in or within 500 bp of the transcription unit.
A similar fraction of chromosomes 11 and 12 could be classified as repetitive (29.5% and 31.6%, respectively (Additional tables 5 and 6 [see Additional file 1]), which with the exception of miniature inverted repeat transposable elements (MITEs) were enriched in the centromeric and pericentromeric regions of the two chromosomes (Fig. 1). As noted previously , MITEs, which constitute 18.6 and 16.2% of the repetitive sequences on chromosomes 11 and 12, respectively, were enriched in the euchromatic arms of the two chromosomes (Fig. 1).
A similar number of non-TE-related proteins with putative homologs in model species were present in chromosomes 11 and 12 (Fig. 2). The highest number of putative homologs was seen in Arabidopsis, with 60.9% and 59.6% of the proteins from chromosomes 11 and 12, respectively, having a potential homolog in Arabidopsis (E value cutoff 10-5). Similar ranges of potential homologs were seen between the two chromosomes and the other model organisms, ranging from 7–30.6% in chromosome 11 and 10–33.6% in chromosome 12.
A high percentage of rice genes could be aligned with genomic assemblies from two other cereals, maize and sorghum. Of the 2,405 maize assemblies with a best hit on rice chromosome 11 or 12, 1,879 overlap with non-TE-related genes while 160 overlap with TE-related genes (Additional table 7 [see Additional file 1]). The remaining 366 assemblies align to intergenic regions, which might include genes or TEs that are degenerated and not recognized by our bioinformatics methods. Interestingly, for sorghum, the occurrence of alignments was similar to maize except that there was a higher incidence of alignment to TE or intergenic regions, which could be due the fact that the sorghum dataset includes only methyl-filtrated and not Co t-enriched sequences [20–22]. The 1,469 non-TE-related genes that aligned to both maize and sorghum genomic assemblies suggest that these genes might predate the divergence of the Panicoideae, estimated at 50 mya .
Disease resistance genes in rice chromosomes 11 and 12
Disease resistance genes (R-like genes) conferring resistance to viral, bacterial, fungal and nematode pathogens have been grouped into five classes on the basis of their encoded protein products . For instance, a genome-wide study identified such typical protein domains in the R-like genes of Arabidopsis . Major classes of conserved domains include leucine zipper (LZ), coiled coil (CC), nucleotide binding site (NBS), leucine-rich repeat (LRR), protein kinase, Toll-IL-IR homology region (TIR) and trans-membrane (TM) as well as miscellaneous R-genes. Apart from the genes containing these typical domains, one maize disease resistance gene, HM1, encodes a reductase that detoxifies HC-toxin . We analyzed the sequences of chromosomes 11 and 12 to identify the type and distribution of the disease resistance as well as downstream defense response genes (Fig. 3). We identified 201 R-like gene models in chromosome 11 (Additional table 8 [see Additional file 1]), which is 4.5% of the total number of gene models predicted for this chromosome. Of these, 73 (36.3%) have homology to the NBS-LRR class of R-like genes and 38 show homology to the LRR-TM-like genes. The long arm of rice chromosome 11 (11L) contains almost twice (132 genes) the number of R-like genes compared to the short arm (68 genes). We also identified 17 downstream defense response genes including glucanases, chitinases and thaumatin-like proteins. Most of these R-like genes and defense response-like genes are present in large clusters of tandem arrays indicating their origin by duplication from a few ancestral genes (Fig. 4). Large clusters of genes (31, 35 and 37 members) are present between positions 30–40 cM (between 5,658,160 – 7,920,402 bp), 80–90 cM (between 20,001,952 – 22,632,644 bp) and 110–119 cM (between 26,213,144 – 28,180,239 bp) of chromosome 11, respectively. Nucleotide positions are based on the pseudomolecules of TIGR release 3 . A detailed depiction of duplicated genes in the interval 112–119 cM on chromosome 11 is shown in Fig. 5 and a cladogram of these genes is shown in Additional figure 11 [see Additional file 1], confirming that the genes within a cluster generally belong to the same clade. A tandem array of 12 Xa21-like genes is present on the short arm of chromosome 11 at position 19 cM (between 3,516,492 – 3,665,074 bp), while the actual Xa21 gene is present on the long arm of chromosome 11 . A large cluster of 14 defense response genes, 12 of which are chitinases, is present in tandem at 116.2 cM (between 28,056,455 – 28,122,601 bp). Several other R-like genes are also arranged in similar but small clusters.
On chromosome 12, 88 gene models showed homology to R-like genes (Fig. 3 Additional table 8 [see Additional file 1]), which is 2.0% of the total number of gene models predicted for this chromosome. However, 50 of these 88 gene models belonged to the miscellaneous category including viral resistance, Verticillium wilt resistance, BLB resistance genes and those containing LRR motif but without NBS, CC or LZ motifs. On chromosome 12, only 18 gene models showed homology to the NBS-LRR category. Although the R-like genes and defense response genes were distributed throughout the chromosome, they were present in clusters (Fig. 4). A large cluster of 23 miscellaneous resistance genes was found between positions 40–50 cM (between 5,858,249 and 10,009,727 bp) and a cluster of seven defense response genes was found at position 107.4 cM (between 26,803,729 and 26,870,016 bp). The previously described blast resistance gene Pita  was present near the centromere of chromosome 12. The total number of R-like genes in chromosome 12 was less than half of the number in chromosome 11, even though the size and total number of genes in these two chromosomes are similar. This can be attributed to the enrichment of NBS-LRR, LZ-NBS-LRR and LRR-TM-like genes in chromosome 11, which was 3–4 times higher than chromosome 12 (Fig. 3). The difference in the number of R-like genes between the two chromosomes also has an impact on the extent of tandem gene arrays. Chromosome 11 has a total of 924 genes (29%) and chromosome 12 has 684 genes (24%) that are duplicated at least once within a short distance (Fig. 6).
A previous report on the analysis of R-like genes in the entire rice genome revealed that most of the R-like genes (24.98%) are present on chromosome 11 . Another report found more than 25% of the physically mapped NBS-encoding genes on chromosome 11 and claimed that approximately 20% of the NBS-LRR genes in the Nipponbare genome were predicted to be pseudogenes . Comparison of the total number of genes containing the LRR domain (837) in the finished genome sequence  with the number of such genes on chromosomes 11 (166 genes) and 12 (60 genes) in the present study (Tables 2 and 3) shows a similar percentage (27%) of R-like genes to be on these two chromosomes. However, the previous total estimate of 536 R-like genes is significantly lower than the one from the finished genome (837), indicating that previous draft sequences did not provide sufficient information for capturing all candidate genes for disease resistance and, in addition, that gene finding methods are not directly comparable. Complete sequence information is also important for accurate sequence alignments of R-like genes with cloned disease resistance genes from rice (Pita, Xa21, Pib) Arabidopsis (RPM1), tomato, (Cf2, Cf9), barley (Mla1) and wheat (Lr10). The R-like gene and defense response gene hot spots in these two chromosomes will be invaluable for future mapping and cloning of disease resistance genes from naturally occurring disease resistant rice lines. However, manual re-annotation of these R-like genes, as done in case of Arabidopsis , and functional validation of the candidate resistance genes, would be important for drawing practical benefits from this information.
Recent duplication on chromosomes 11 and 12
Chromosomes 11 and 12 harbor duplicated regions at the distal ends of their short arms as determined by physical and genetic mapping . Based on previous versions of the rice genome sequence, it became evident that segmental duplications occur throughout the genome [5, 7, 33]. Here, we focused on segmental duplications between chromosome 11 and 12 only because the finished sequences of these chromosomes permit us to conduct an analysis of the organization of their genes in a detail that was not possible before. At ≥50% coverage and ≥80% identity and excluding TE-related gene models, 546 gene models were identified as duplicated. A subset of these, 350 from chromosome 11 and 352 from chromosome 12, was found to be unique between these chromosomes, excluding repeated gene models. Based on the chromosomal locations of these unique genes, the maximum extent of duplication was confined within the first 3 Mb of both chromosomes (Fig. 7). Possibly, there could be a difference between japonica and indica varieties, as previous comparisons identified a second duplication between chromosomes 11 and 12 in the indica  but not the japonica  species. Here, most of the duplicated gene models (98%) were found in the same orientation on both the chromosomes (Fig. 8a,b). The 3 Mb region was further analyzed to find the number of duplicated genes at ≥30, ≥40, ≥50, ≥60, ≥70, ≥80 and ≥90% coverage. The numbers of duplicated genes at ≥30% coverage for chromosomes 11 and 12 were 287 and 304, while at ≥90% coverage, the numbers of duplicated genes were 109 and 113, respectively (Additional table 9 [see Additional file 1]). Using alignment to ESTs and FL-cDNAs, a total of 132 (chromosome 11) and 145 (chromosome 12) gene models were expressed; of these, a total of 90 expressed gene models were common to both chromosomes. Interestingly, 42 and 55 gene models were expressed but their homoeologous copies on the respective chromosome were not expressed [Additional table 10; see Additional file 1]. Although it is possible that all ESTs or cDNAs are not represented in the database, it may be that variation in expression reflects diverged expression of one of the homoeologous copies. Divergence of regulation of gene expression after gene amplification has also been reported in maize [15, 34].
To determine the time of the segmental duplication, the full-length coding regions of a subset of homoeologous genes that had been genetically mapped and were distributed over the entire length of the duplication were aligned and pair-wise analyzed for their nucleotide substitution rates (Table 4). By applying a codon likelihood model, synonymous (Ks) and non-synonymous (Ka) substitution rates were compared. The values for each gene pair varied 18-fold for Ks and 556-fold for Ka, indicating that divergence rates differed significantly among genes. We therefore applied the χ2 homogeneity test to the estimates of divergence time for the rice homoeologs . The χ2 value was highly significant (χ2 = 191.2, P < 0.001), implying that these linked genes diverged within the same time interval. Based on the speciation of rice about 50 mya , we calculate the average divergence time for the segmental duplication to be 7.7 mya. We explain the difference from previously published estimates of 25 and 21 mya [3, 5] by the sequence accuracy of finished versus whole-genome shotgun sequences .
The number of unique duplicated gene models for the remaining portion (~24–25 Mb) of chromosomes 11 and 12 was 114 (32.6%) and 102 (29%), respectively, with a high degree of rearrangement. Within this region, there were 163 and 189 gene models at ≥30% coverage, while at ≥90% coverage, there were only 31 genes for both the chromosomes (Additional table 10 [see Additional file 1]). Interestingly, only 18 (15.8%) and 14 (13.7%) genes (chromosome 11 and 12, respectively) were found to be expressed based on homology to an EST or cDNA, in contrast to 56 and 58% in the first 3 Mb region of both the chromosomes. Out of these, nine expressed genes were common to both chromosomes. Also, there were nine and five expressed genes of chromosomes 11 and 12 that had non-expressed counterparts. This analysis supports the view that the first 3 Mb duplication of chromosome 11 and 12 is of recent origin. However, the loss of order as well as the low percentage of expressed genes, on each chromosome or in duplicate gene pairs, does not support a duplication of the remaining regions of chromosomes 11 and 12, unless it is very ancient in origin and signatures of duplication have largely disappeared.
To investigate further whether rice chromosomes 11 and 12 were related to each other and arose by a WGD event, we compared each gene model predicted from these two chromosomes against 586,577 wheat ESTs. Comparison of bin-mapped wheat ESTs with earlier versions of rice genome sequence provided a glimpse of the structural similarities between these two important cereal genomes [11, 13]. However, the present study involves comparative analysis of the finished version of rice chromosomes 11 and 12, which allowed a higher resolution analysis with wheat and a comparative analysis between rice chromosomes 11 and 12. A gene-by-gene comparison of sequence homology of all the predicted rice gene models with the wheat ESTs revealed that 1,588 (35.8%) and 2,220 (51%) of the gene models in chromosomes 11 and 12, respectively, have significant homology with the wheat ESTs at a bit cut-off score of 100. Of these, 416 gene models (26.2%) from chromosome 11 and 552 gene models (24.86%) from chromosome 12 showed significant similarity with bin-mapped wheat EST contigs, a number that is much larger than the levels reported earlier for chromosomes 11 and 12 . A complete list of the matching genes along with their annotated functions is provided in Additional tables 11 and 12 [see Additional file 1].
Out of the 416 gene models from rice chromosome 11, 338 (81.3%) mapped to single wheat homoeologous groups, with the maximum number (35.2%) located on group 4 chromosomes of wheat. The distribution of rice genes that mapped to the wheat chromosome groups 1, 2, 3, 5, 6 and 7 was 9.2%, 9.5%, 12.4%, 15.4%, 8.9% and 9.5%, respectively (Fig. 9 Additional figures 12 and 13 [see Additional file 1]). The majority of rice chromosome 11 gene models mapped to group 4 chromosomes of wheat indicating their common origin. Many of these genes were clustered in the distal region of the long arm of wheat chromosome 4A, but the same gene models mapped to the short arms of wheat chromosomes 4B and 4D, reflecting significant rearrangements and a dynamic state of gene organization in homoeologous chromosomes of wheat and rice chromosome 11. Similar to chromosome 11, a high percentage of chromosome 12 gene models (441, 79.9%) mapped to single wheat homoeologous groups. Of these, the maximum percentage (30.6%) mapped to group 5 chromosomes of wheat whereas the remainder were distributed almost uniformly on the other six wheat homoeologous groups (Fig. 9 Additional figures 12 and 13 [see Additional file 1]).
A comparative distribution of rice gene homologs from chromosomes 11 and 12 to the seven wheat homoeologous groups clearly indicates that the two rice chromosomes have different origins and, apart from the recent 3 Mb duplication, there is not much in common between the two chromosomes. This contradicts earlier observations that the whole of the rice chromosomes 11 and 12 may have evolved as a result of chromosome duplication via polyploidy about 70 mya, i.e. before the divergence of cereals . Furthermore, these results are also consistent with the alignment of the genetic maps of cereal genomes , where rice aligns with two homoeologous regions of maize.
With the completion of the sequences of rice chromosomes 11 and 12, we were able to identify 289 R-like and 28 defense response-like genes, accurately date a 3 Mb recent duplication between the two chromosomes, and show significant synteny between these two chromosomes and wheat chromosome groups 4 and 5, respectively. Rice chromosome 11 has several large clusters of fast evolving disease resistance and defense response genes that have originated by the process of tandem duplication and subsequent divergence under the selective pressure of rice pathogens. This sequence and annotation will be an essential resource for the engineering of rice for tolerance to biotic and abiotic stresses to accommodate growing production needs.
Sequencing and bioinformatic methods
Bacterial artificial chromosome and P1-derived artificial chromosome (PAC) clones from chromosomes 11 and 12 were sequenced to 8–10× coverage using shotgun DNA sequencing [35, 36] and standard high-throughput methods . Gaps and ambiguities were resolved using re-sequencing, alternative chemistries, PCR and primer walking . Pseudomolecules or virtual contigs for chromosomes 11 and 12 were constructed as described previously . The pseudomolecules and annotation are available in GenBank under the accession numbers DP000010 and DP000011. Genes were identified using the ab initio gene finder FGENESH  and gene model structure was improved with EST and full-length cDNA (FL-cDNA) evidence using the Program to Assemble Spliced Alignments . Genes were annotated for function as described . Transposable element-related (TE) genes were identified using either TBLASTN similarity to known repetitive elements in the TIGR Oryza Repeat Database  or the presence of repetitive element-related Pfam domains . Flanking sequence tags (FSTs) from various insertion mutagenesis projects were aligned with the two pseudomolecules using flast  with a cutoff of 95% identity over 80% of the FST length. Transfer RNAs were identified using tRNA-ScanSE . Organellar insertions were determined using BLASTN with the rice mitochondrion and chloroplast genome sequences using a cutoff of 95% identity (Fig. 1). Sequenced genetic markers (4,619 sequences [18, 43]) were searched against the pseudomolecules with flast; a cutoff criterion of ≥95% identity over ≥90% length of the marker was used to align the markers to the chromosome sequence. Repetitive sequences were identified, classified, and quantitated on the pseudomolecules using RepeatMasker  with the TIGR Oryza Repeat Database  and a default RepeatMasker cut-off score of 225.
Non-TE-related genes with rice transcript evidence were identified by searching chromosomes 11 and 12 against the TIGR rice gene index (Release 16 ). To identify genes with FL-cDNA support, we searched 33,678 FL-cDNAs available from the KOME database . Genes with an alignment of ≥95% identity and a minimum of 50% length of the alignment covered by a gene index sequence were considered supported by a rice EST or FL-cDNA. Seven monocot gene indices  were also used to demonstrate support of expressed genes (wheat release 9.0, maize release 15.0, barley release 9.0, sorghum release 8.0, rye release 3.0, onion release 1.0 and sugarcane release 2.1). Cutoff criteria used with the non-rice monocot gene indices were ≥70% identity over 80% of the length of the gene index sequence and a minimum of 50% of the gene model coverage by the gene index sequence. For the analysis of the type and distribution of the disease resistance and defense response genes, coding sequences from the chromosomes 11 and 12 gene models were used in a BLASTX search with the nr database of NCBI  and the top hits were extracted in Excel files. The BLASTX output was then searched manually using auto filters against 24 different keywords/phrases known to represent R-like and defense response genes, and categorized into five main classes as follows: (i) NBS-LRR (matching with NBS-LRR, but not with LZ-NBS-LRR and LRR, CC-NBS-LRR, Pib, Pita, I2C, Rp1-d8, T10RGA, LR10, Mla1 and rust resistance), (ii) LZ-NBS-LRR (matching with LZ-NBS-LRR, but not with NBS-LRR, CC-NBS-LRR, LRR and RPM1), (iii) LRR-TM (matching with Xa21, serine/threonine kinases and Cf2/Cf5 resistance), (iv) miscellaneous category (matching with disease resistance, viral resistance, Yr10, Verticillium wilt resistance, LRR, but not with NBS-LRR, CC-NBS-LRR, LZ-NBS-LRR and bacterial leaf blight resistance) and (v) defense response genes (matching with glucanases, chitinases and thaumatin).
Maize genomic assemblies (243,807 total [Release 4.0, Feb. 23, 2004]) were downloaded from TIGR [20, 49]. Sorghum genomic reads (593,969) were downloaded from NCBI and processed using the Lucy program  to remove vector and low quality sequences. The remaining "good" reads (504,458) were assembled into 163,908 clusters as described . The assemblies from sorghum and maize were searched against the TIGR rice pseudomolecules [Release 3.0, January 2005] using the BLASTZ program . Homologs were identified in other model organisms by searching the non-TE-related proteins from chromosomes 11 and 12 against the E. coli, Synechocystis, yeast, Drosophila, C. elegans, human and Arabidopsis predicted proteomes using BLASTP.
Chromosome 11–12 duplication analysis was performed using the predicted gene models. Chromosome 12 sequences were used as query against a database of chromosome 11 sequences using MegaBLAST . Only non-TE gene models at ≥50% coverage were used for further analysis. To compare the expression pattern of duplicated genes, a homology search was performed against the NCBI EST and KOME cDNA databases using the unique duplicated non-TE gene models at ≥90% coverage and identity. Wheat-rice synteny was determined using BLASTN of the chromosome 11 and 12 gene model sequences against Triticum aestivum sequences obtained from dbEST. Search parameters were as described previously . Matching gene models were compared with wheat contigs containing bin-mapped ESTs  and plotted on the 21 wheat chromosomes.
bacterial artificial chromosome
flanking sequence tag
million years ago
nucleotide binding site
Khush GS: Origin, dispersal, cultivation and variation of rice. Plant Mol Biol. 1997, 35 (1–2): 25-34. 10.1023/A:1005810616885.
Peng S, Huang J, Sheehy JE, Laza RC, Visperas RM, Zhong X, Centeno GS, Khush GS, Cassman KG: Rice yields decline with higher night temperature from global warming. Proc Natl Acad Sci U S A. 2004, 101 (27): 9971-9975. 10.1073/pnas.0403720101.
Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, et al: A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002, 296 (5565): 92-100. 10.1126/science.1068275.
Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, et al: A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002, 296 (5565): 79-92. 10.1126/science.1068037.
Yu J, Wang J, Lin W, Li S, Li H, Zhou J, Ni P, Dong W, Hu S, Zeng C, et al: The Genomes of Oryza sativa: a history of duplications. PLoS Biol. 2005, 3 (2): e38-10.1371/journal.pbio.0030038.
International-Rice-Genome-Sequencing-Project: The map-based sequence of the rice genome. Nature. 2005, 436 (7052): 793-800.
Paterson AH, Bowers JE, Chapman BA: Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics. Proc Natl Acad Sci U S A. 2004, 101 (26): 9903-9908. 10.1073/pnas.0307901101.
Sasaki T, Matsumoto T, Antonio BA, Nagamura Y: From mapping to sequencing, post-sequencing and beyond. Plant Cell Physiol. 2005, 46 (1): 3-13. 10.1093/pcp/pci503.
Arumuganathan K, Earle ED: Nuclear DNA content of some important plant species. Plant Mol Biol Reporter. 1991, 9: 208-218.
Gale MD, Devos KM: Comparative genetics in the grasses. Proc Natl Acad Sci U S A. 1998, 95 (5): 1971-1974. 10.1073/pnas.95.5.1971.
Sorrells ME, La Rota M, Bermudez-Kandianis CE, Greene RA, Kantety R, Munkvold JD, Miftahudin , Mahmoud A, Ma X, Gustafson PJ, et al: Comparative DNA sequence analysis of wheat and rice genomes. Genome Res. 2003, 13 (8): 1818-1827.
Singh NK, Raghuvanshi S, Srivastava SK, Gaur A, Pal AK, Dalal V, Singh A, Ghazi IA, Bhargav A, Yadav M, et al: Sequence analysis of the long arm of rice chromosome 11 for rice-wheat synteny. Funct Integr Genomics. 2004, 4 (2): 102-117. 10.1007/s10142-004-0109-y.
La Rota M, Sorrells ME: Comparative DNA sequence analysis of mapped wheat ESTs reveals the complexity of genome relationships between rice and wheat. Funct Integr Genomics. 2004, 4 (1): 34-46. 10.1007/s10142-003-0098-2.
Lai J, Ma J, Swigonova Z, Ramakrishna W, Linton E, Llaca V, Tanyolac B, Park YJ, Jeong OY, Bennetzen JL, et al: Gene loss and movement in the maize genome. Genome Res. 2004, 14 (10A): 1924-1931. 10.1101/gr.2701104.
Swigonova Z, Bennetzen JL, Messing J: Structure and evolution of the r/b chromosomal regions in rice, maize and sorghum. Genetics. 2005, 169 (2): 891-906. 10.1534/genetics.104.034629.
Swigonova Z, Lai J, Ma J, Ramakrishna W, Llaca V, Bennetzen JL, Messing J: Close split of sorghum and maize genome progenitors. Genome Res. 2004, 14 (10A): 1916-1923. 10.1101/gr.2332504.
Khush GS, Brar DS: Rice genetics from Mendel to functional genomics. Rice Genetics. Edited by: Kush GS, Brar DS, Hardy B. 2001, Los Banos, Philippines: IRRI, IV: 3-25.
Harushima Y, Yano M, Shomura A, Sato M, Shimano T, Kuboki Y, Yamamoto T, Lin SY, Antonio BA, Parco A, et al: A high-density rice genetic linkage map with 2275 markers using a single F2 population. Genetics. 1998, 148 (1): 479-494.
The-Rice-Chromosome-10-Sequencing-Consortium: In-depth view of structure, activity, and evolution of rice chromosome 10. Science. 2003, 300 (5625): 1566-1569. 10.1126/science.1083523.
Whitelaw CA, Barbazuk WB, Pertea G, Chan AP, Cheung F, Lee Y, Zheng L, van Heeringen S, Karamycheva S, Bennetzen JL, et al: Enrichment of gene-coding sequences in maize by genome filtration. Science. 2003, 302 (5653): 2118-2120. 10.1126/science.1090047.
Palmer LE, Rabinowicz PD, O'Shaughnessy AL, Balija VS, Nascimento LU, Dike S, de la Bastide M, Martienssen RA, McCombie WR: Maize genome sequencing by methylation filtration. Science. 2003, 302 (5653): 2115-2117. 10.1126/science.1091265.
Bedell JA, Budiman MA, Nunberg A, Citek RW, Robbins D, Jones J, Flick E, Rholfing T, Fries J, Bradford K, et al: Sorghum genome sequencing by methylation filtration. PLoS Biol. 2005, 3 (1): e13-10.1371/journal.pbio.0030013.
Gaut BS: Evolutionary dynamics of grass genomes. New Phytologist. 2002, 154: 15-28. 10.1046/j.1469-8137.2002.00352.x.
Hulbert SH, Webb CA, Smith SM, Sun Q: Resistance gene complexes: evolution and ugtilization. Annu Rev Phytopathol. 2001, 39: 285-312. 10.1146/annurev.phyto.39.1.285.
Meyers BC, Kozik A, Griego A, Kuang H, Michelmore RW: Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis. Plant Cell. 2003, 15 (4): 809-834. 10.1105/tpc.009308.
Johal GS, Briggs SP: Reductase activity encoded by the HM1 disease resistance gene in maize. Science. 1992, 258 (5084): 985-987.
Yuan Q, Ouyang S, Wang A, Zhu W, Maiti R, Lin H, Hamilton J, Haas B, Sultana R, Cheung F, et al: The institute for genomic research Osa1 rice genome annotation database. Plant Physiol. 2005, 138 (1): 18-26. 10.1104/pp.104.059063.
Song WY, Wang GL, Chen LL, Kim HS, Pi LY, Holsten T, Gardner J, Wang B, Zhai WX, Zhu LH, et al: A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science. 1995, 270 (5243): 1804-1806.
Bryan GT, Wu KS, Farrall L, Jia Y, Hershey HP, McAdams SA, Faulk KN, Donaldson GK, Tarchini R, Valent B: tA single amino acid difference distinguishes resistant and susceptible alleles of the rice blast resistance gene Pi-ta. Plant Cell. 2000, 12 (11): 2033-2046. 10.1105/tpc.12.11.2033.
Zhou T, Wang Y, Chen JQ, Araki H, Jing Z, Jiang K, Shen J, Tian D: Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes. Mol Genet Genomics. 2004, 271 (4): 402-415. 10.1007/s00438-004-0990-z.
Monosi B, Wisser RJ, Pennill L, Hulbert SH: Full-genome analysis of resistance gene homologues in rice. Theor Appl Genet. 2004, 109 (7): 1434-1447. 10.1007/s00122-004-1758-x.
Wu J, Kurata N, Tanoue H, Shimokawa T, Umehara Y, Yano M, Sasaki T: Physical mapping of duplicated genomic regions of two chromosome ends in rice. Genetics. 1998, 150 (4): 1595-1603.
Guyot R, Keller B: Ancestral genome duplication in rice. Genome. 2004, 47 (3): 610-614. 10.1139/g04-016.
Song R, Llaca V, Linton E, Messing J: Sequence, regulation, and evolution of the maize 22-kD alpha zein gene family. Genome Res. 2001, 11 (11): 1817-1825.
Messing J, Crea R, Seeburg PH: A system for shotgun DNA sequencing. Nucleic Acids Res. 1981, 9 (2): 309-321.
Vieira J, Messing J: The pUC plasmids, an M13mp7-derived system for insertion mutagenesis and sequencing with synthetic universal primers. Gene. 1982, 19 (3): 259-268. 10.1016/0378-1119(82)90015-4.
Salamov AA, Solovyev VV: Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000, 10 (4): 516-522. 10.1101/gr.10.4.516.
Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK, Hannick LI, Maiti R, Ronning CM, Rusch DB, Town CD, et al: Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003, 31 (19): 5654-5666. 10.1093/nar/gkg770.
Ouyang S, Buell CR: The TIGR Plant Repeat Databases: a collective resource for the identification of repetitive sequences in plants. Nucleic Acids Res. 2004, 32 (Database): D360-363. 10.1093/nar/gkh099.
Bateman A, Birney E, Cerruti L, Durbin R, Etwiller L, Eddy SR, Griffiths-Jones S, Howe KL, Marshall M, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res. 2002, 30 (1): 276-280. 10.1093/nar/30.1.276.
Yuan Q, Liang F, Hsiao J, Zismann V, Benito MI, Quackenbush J, Wing R, Buell R: Anchoring of rice BAC clones to the rice genetic map in silico. Nucleic Acids Res. 2000, 28 (18): 3636-3641. 10.1093/nar/28.18.3636.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25 (5): 955-964. 10.1093/nar/25.5.955.
Causse MA, Fulton TM, Cho YG, Ahn SN, Chunwongse J, Wu K, Xiao J, Yu Z, Ronald PC, Harrington SE, et al: Saturated molecular map of the rice genome based on an interspecific backcross population. Genetics. 1994, 138 (4): 1251-1274.
Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J: The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucleic Acids Res. 2001, 29 (1): 159-164. 10.1093/nar/29.1.159.
Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H, et al: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301 (5631): 376-379. 10.1126/science.1081288.
Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17 (12): 1093-1104. 10.1093/bioinformatics/17.12.1093.
Schwartz S, Kent WJ, Smit A, Zhang Z, Baertsch R, Hardison RC, Haussler D, Miller W: Human-mouse alignments with BLASTZ. Genome Res. 2003, 13 (1): 103-107. 10.1101/gr.809403.
Zhang Z, Schwartz S, Wagner L, Miller W: A greedy algorithm for aligning DNA sequences. Comput Biol. 2000, 7: 203-214. 10.1089/10665270050081478.
Genoscope: Nathalie Choisne1, Nadia Demange1, Gisela Orjeda1, Sylvie Samain1, Angélique D'Hont2, Laurence Cattolico1, Eric Pelletier1, Arnaud Couloux1, Béatrice Segurens1, Patrick Wincker1, Claude Scarpelli1, Jean Weissenbach1, Marcel Salanoubat1, Francis Quétier1; Indian Initiative for Rice Genome Sequencing (IIRGS), Indian Agricultural Research Institute: Nagendra K Singh3, Trilochan Mohapatra3, Tilak R Sharma3, Kishor Gaikwad3, Archana Singh3, Vivek Dalal3, Subodh K Srivastava3, Anupam Dixit3, Ajit K Pal3, Irfan A Ghazi3, Mahavir Yadav3, Awadhesh Pandit3, Ashutosh Bhargava3, K Sureshbabu3, Rekha Dixit3, Harvinder Singh3, Suresh C Swain3, Sumita Pal3, M Ragiba3, Pradeep K Singh3, Vibha Singhal3, Sangeeta D Mendiratta3, Kamlesh Batra3; Indian Initiative for Rice Genome Sequencing (IIRGS), University of Delhi South Campus: Saurabh Raghuvanshi4, Amitabh Mohanty4, Arvind K. Bharti4,8, Anupama Gaur4, Vikrant Gupta4, Dibyendu Kumar4, Ravi Vydianathan4, Shubha Vij4, Anita Kapur4, Parul Khurana4, Sulabha Sharma4, Paramjit Khurana4, Jitendra P. Khurana4, Akhilesh K. Tyagi4†; The Institute for Genomic Research (TIGR): Qiaoping Yuan5, Shu Ouyang5, Jia Liu5, Wei Zhu5, Aihui Wang5, Haining Lin5, John Hamilton5, Brian Haas5, Jennifer Wortman5, Kristine M. Jones5, Mary Kim5, Larry Overton5, Tamara Tsitrin5, Douglas Fadrosh5, Jayati Bera5, Bruce Weaver5, Shaohua Jin5, Shivani Johri5, Matt Reardon5, Hue Vuong5, Luke Tallon5, Susan Van Aken5, Matthew Lewis5, Teresa Utterback5, Tamara Feldblyum5, Victoria Zismann5, Stacey Iobst5, Joseph Hsiao5, Aymeric R. de Vazeille5, Steven L. Salzberg5, Owen White5, Claire Fraser5, and C. Robin Buell5; Arizona Genomics Institute: Yeisoo Yu6, Teri Rambo6, Jennifer Currie6, Kristi Collura6, Hye Ran Kim6, Diana Stum6, Wenming Wang6, Dave Kudrna6, Christopher Mueller6, Rod A. Wing6; Cold Spring Harbor Laboratory: Melissa Kramer7, Lori Spiegel7, Lidia Nascimento7, Raymond Preston7, Theresa Zutavern7, Lance Palmer7, Andrew O'Shaughnessy7, Sujit Dike7, W. Richard McCombie7; Plant Genome Initiative at Rutgers: Joachim Messing8, Charles Du8,11 Galina Fuks8, Eric Linton8,12 Bahattin Tanyolac8,13; Rice Genome Research Program: Jianzhong Wu9, Nobukazu Namiki9, Hiroshi Mizuno9, Takashi Matsumoto9, Takuji Sasaki9; University of Wisconsin: Weiwei Jin10, Jiming Jiang10
Consortia members' affiliation:
1Genoscope, 2 rue Gaston Crémieux CP5706, 91057 Evry-Cedex (France); 2UMR PIA, Cirad-Amis, TA40-03 avenue Agropolis, 34398 Montpellier Cedex 05 (France); 3National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute, New Delhi-110 012, India; 4Department of Plant Molecular Biology, University of Delhi South Campus, New Delhi-110021, India; 5The Institute for Genomic Research, 9712 Medical Center Dr, Rockville, MD 20850; 6Arizona Genomics Institute, The University of Arizona, Tucson, AZ 85750; 7Cold Spring Harbor Laboratory, Cold Spring Harbor, NY 11724; 8Plant Genome Initiative at Rutgers, Waksman Institute, Rutgers University, Piscataway, New Jersey 08854; 9Rice Genome Research Program, National Institute of Agrobiological Sciences, 1–2 Kannondai, 2-Chome, Tsukuba, Ibaraki 305-8602, Japan; 10University of Wisconsin, Department of Horticulture, Madison, WI 53706, USA.
Consortia members' present addresses:
11Montclair State University, Montclair, New Jersey 07043; 12Plant Biology Labs, Michigan State University, East Lansing, MI 48824; 13Department of Bioengineering, Ege University, Izmir 35100, Turkey.
Consortia members' sponsors:
The IIRGS acknowledges the Department of Biotechnology, Government of India, for financial assistance, and the Indian Council of Agricultural Research, New Delhi, for support. Work on rice chromosome 11 at TIGR was made possible by grants from the U. S. Department of Agriculture Cooperative State Research, Education, and Extension Service (C.R.B. 99-35317-8275, 2003-35317-13173), the National Science Foundation (C.R.B. DBI998282; DBI0321538), and the U. S. Department of Energy (C.R.B. DE-FG02-99ER20357). We acknowledge the assistance of the TIGR Sequencing Facility, the TIGR Informatics Department, the TIGR IT Group, and the J. Craig Venter Joint Technology Center. Funding for work on rice chromosome 11 was provided by a grant from the U. S. Department of Agriculture Cooperative State Research, Education, and Extension Service to R.A.W and W.R.M. (#2002-35317-12414). R.A.W. acknowledges the assistance of the AGI Sequencing, Physical Mapping, BAC/EST Resource Centers and the Arizona Genomics Computational Laboratory. W.R.M. acknowledges the CSHL sequencing and informatics groups. Work at the PGIR was supported by Rutgers, The State University of New Jersey. We thank the Japanese Ministry of Agriculture, Forestry and Fishery for genetic markers, which were utilized to identify seed BACs for sequencing.
Consortia members are listed in order of Mb of DNA sequences produced. CRB (TIGR), NKS (IIRGS-IARI), AT (IIRGS-UDSC), and JM (PGIR) conducted the coordination of data analysis and manuscript conception.
At GENOSCOPE, EP, AC, JW, MS, and FQ provided data analysis or intellectual input while NC, ND, GO, SS, AD, LC, BS, PW and CS contributed to sequence/assembly aspects of the project.
At IIRGS-IARI, AD, IAG contributed in physical mapping and identification of BACs; KG, AS, MY, RD, SP, MR, PKS and VS in shotgun cloning, template preparation; SS, AB, AP, KS, HS, SCS and SDM in sequencing, VD and AKP in sequence assembly and annotation, KB in wheat-rice synteny. TRS, TM and NKS contributed in all activities and provided intellectual inputs in analysis.
At IIRGS-UDSC, SR, AM, AKB, AG, VG, DK, VR, SV, AK, PK, SS contributed to sequence, assembly and data analysis. PK, JPK, AKT provided intellectual input for strategic plan and participated in manuscript preparation. AKT is also Coordinator of IIRGS.
At TIGR, QY, SO, JL, WZ, AW, HL, JH, BH, JW, SLS, OW, CF, and CRB provided data analysis or intellectual input while KMJ, MK, LO, TT, DF, JB, BW, SJ, MR, HV, LT, SVA, ML, TU, TF, VZ, SI, JH and ARV contributed to sequence/assembly aspects of the project.
At AGI, WW, DK and CM contributed to shotgun library construction and clone management. YY, TR, JC, KC, HRK, DS and RAW contributed to shotgun sequencing, sequence assembly and finishing aspects of this project.
At CSHL, MK and LS provided sequence assembly and finishing, LN, RP and TZ production sequencing, LP, AO and SD contributed to data analysis and annotation, and WRM directed the research.
At PGIR, EL contributed to clone management, EL and BT to shotgun library construction and sequence/assembly, CD, GF and JM to data analysis, and JM to writing the manuscript; JM is the corresponding author.
At RGP, JW, HM and TM contributed to physical map construction, Fiber-FISH analysis, and physical map extension, NM carried out data analysis and TS directed the research.
At UW, WJ and JJ contributed to gap analysis using fiber-FISH.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.