Abstract
Rice (Oryza sativa) is one of the key staple crops, providing food for nearly half of the world's population. The past twenty years have seen significant advances in understanding Oryza species through genome sequencing efforts. However, the stability of Oryza genomes during their divergence has not been well characterized. Here, by performing gene collinearity and comparative genomics analysis, we selected ten Oryza species and three other Poaceae species to check their genome stability, with Leersia perrieri as the reference. Intra- and intergenomic analysis showed a ~30% difference in homologous block numbers and a 35.7% difference in collinear gene numbers per block, indicating that Oryza genomes have undergone extensive DNA permutations. Notably, we found that Oryza chromosomes with smaller ordinal numbers have often preserved larger percentages of genes, while those with bigger numbers have undergone more gene losses. This unique observation may be explained by elevated gene losses incurred by illegitimate or homoeologous recombination between homoeologous chromosomes produced by the grass-common tetraploidization (GCT) ~100 million years ago (Mya), e.g., Chro. 11 and 12. However, the lowered gene loss rates in Chro. 1-3 could be explained by earlier restriction of illegitimate recombination after the GCT due to there often being (larger) neo-chromosomes produced by the fusion of ancestral chromosomes. The enriched NBS-LRR (nucleotide-binding site and leucine-rich repeat) genes in chromosomes 11 and 12 are another explanation for the above observation. Further evidence was obtained from other Poaceae plants. Moreover, we revealed around twice as many differences in tandem genes and their densities among Oryza plants, further showing their divergent levels of genome stability. The present efforts may contribute to the understanding of the stability of the Oryza genome and its formation, evolution, and functional innovation.