Multilocus analysis of nucleotide variation and speciation in Oryza officinalis and its close relatives

Nucleotide variation in 10 unlinked nuclear genes was investigated in species-wide samples of Oryza officinalis and its close relatives (Oryza eichingeri and Oryza rhizomatis). Average estimates of nucleotide diversity were the lowest in O. rhizomatis (theta(sil) = 0.0038) and the highest in O. eichingeri (theta(sil) = 0.0057) that is disjunctly distributed in Africa and Sri Lanka. These wild rice species appeared to harbor relatively low levels of nucleotide variation relative to other plant species because the diversity level of O. eichingeri is only 23-46% of those in Zea species and 35% of that in Arabidopsis thaliana. The lower nucleotide diversity in these Oryza species could be best explained by their smaller historic effective population sizes. The speciation model test indicated that O. officinalis and its close relatives might have undergone a process of population contraction since divergence from their ancestor. Incongruent topologies among 10 gene trees, particularly regarding the positions of O. eichingeri and O. rhizomatis accessions might be attributed to lineage sorting arising from ancient polymorphism and hybridization/introgression between the Sri Lankan O. eichingeri and O. rhizomatis. However, the null hypothesis of the isolation model was not rejected for any contrast between taxa, which suggested that no subsequent gene flow shaped the present patterns of nucleotide variation since their divergence and that introgression was not pervasive in this group of species. Our molecular dating provides an approximate divergence time of 0.37 Myr between 2 geographical races of O. eichingeri, much more recent compared with the times of other speciation events in this group (0.63-0.68 Myr). A long-distance dispersal from West Africa to Sri Lanka was more likely to play a role in the disjunct distribution of O. eichingeri.