Fifteen Million Years of Evolution in the Oryza Genus Shows Extensive Gene Family Expansion

In analyzing gene families in the whole-genome sequences available for O. sativa （AA）, O. glaberrima （AA）, and O. brachyantha （FF）, we observed large size expansions in the AA genomes compared to FF genomes for the superfamilies F-box and NB-ARC, and five additional families： the Aspartic proteases, BTB/POZ proteins （BTB）, Glutaredoxins, Trypsin a-amylase inhibitor proteins, and Zf-Dof proteins. Their evolutionary dynamic was investigated to understand how and why such important size variations are observed between these closely related species. We show that expansions resulted from both amplification, largely by tandem duplications, and contraction by gene losses. For the F-box and NB-ARC gene families, the genes conserved in all species were under strong purifying selection while expanded orthologous genes were under more relaxed purifying selection. In F-box, NB-ARC, and BTB, the expanded groups were enriched in genes with little evidence of expression, in comparison with conserved groups. We also detected 87 loci under positive selection in the expanded groups. These results show that most of the duplicated copies in the expanded groups evolve neutrally after duplication because of functional redundancy but a fraction of these genes were preserved following neofunctionalization. Hence, the lineage-specific expansions observed between Oryza species were partly driven by directional selection.