Machine Learning Reveals Spatiotemporal Genome Evolution in Asian Rice Domestication

Domestication is anthropogenic evolution that fulfills mankinds critical food demand. As such, elucidating the molecular mechanisms behind this process promotes the development of future new food resources including crops. With the aim of understanding the long-term domestication process of Asian rice and by employing the Oryza sativa subspecies (indica and japonica) as an Asian rice domestication model, we scrutinized past genomic introgressions between them as traces of domestication. Here we show the genome-wide introgressive region (IR) map of Asian rice, by utilizing 4,587 accession genotypes with a stable outgroup species, particularly at the finest resolution through a machine learning-aided method. The IR map revealed that 14.2% of the rice genome consists of IRs, including both wide IRs (recent) and narrow IRs (ancient). This introgressive landscape with their time calibration indicates that introgression events happened in multiple genomic regions over multiple periods. From the correspondence between our wide IRs and the so-called selective sweep regions, we provide a definitive answer to a long-standing controversy over the evolutionary origin of Asian rice domestication, single or multiple origins: It heavily depends upon which regions you pay attention to, implying that wider genomic regions represent immediate short history of Asian rice domestication as a likely support to the single origin, while its ancient history is interspersed in narrower traces throughout the genome as a possible support to the multiple origin.