Population genetic structure of Oryza rufipogon and Oryza nivara: implications for the origin of O. nivara

Genetic architecture of ecological divergence between Oryza rufipogon and Oryza nivara

Ecological divergence due to habitat difference plays a prominent role in the formation of new species, but the genetic architecture during ecological speciation and the mechanism underlying phenotypic divergence remain less understood. Two wild ancestors of rice (Oryza rufipogon and Oryza nivara) are a progenitor-derivative species pair with ecological divergence and provide a unique system for studying ecological adaptation/speciation. Here, we constructed a high-resolution linkage map and conducted a quantitative trait locus (QTL) analysis of 19 phenotypic traits using an F2 population generated from a cross between the two Oryza species. We identified 113 QTLs associated with interspecific divergence of 16 quantitative traits, with effect sizes ranging from 1.61% to 34.1% in terms of the percentage of variation explained (PVE). The distribution of effect sizes of QTLs followed a negative exponential, suggesting that a few genes of large effect and many genes of small effect were responsible for the phenotypic divergence. We observed 18 clusters of QTLs (QTL hotspots) on 11 chromosomes, significantly more than that expected by chance, demonstrating the importance of coinheritance of loci/genes in ecological adaptation/speciation. Analysis of effect direction and v-test statistics revealed that interspecific differentiation of most traits was driven by divergent natural selection, supporting the argument that ecological adaptation/speciation would proceed rapidly under coordinated selection on multiple traits. Our findings provide new insights into the understanding of genetic architecture of ecological adaptation and speciation in plants and help effective manipulation of specific genes or gene cluster in rice breeding.

Multiple domestications of Asian rice

The origin of domesticated Asian rice (Oryza sativa L.) has been controversial for more than half a century. The debates have focused on two leading hypotheses: a single domestication event in China or multiple domestication events in geographically separate areas. These two hypotheses differ in their predicted history of genes/alleles selected during domestication. Here we amassed a dataset of 1,578 resequenced genomes, including an expanded sample of wild rice from throughout its geographic range. We identified 993 selected genes that generated phylogenetic trees on which japonica and indica formed a monophyletic group, suggesting that the domestication alleles of these genes originated only once in either japonica or indica. Importantly, the domestication alleles of most selected genes (~80%) stemmed from wild rice in China, but the domestication alleles of a substantial minority of selected genes (~20%) originated from wild rice in South and Southeast Asia, demonstrating separate domestication events of Asian rice.

Phenotypic Variation and the Impact of Admixture in the Oryza rufipogon Species Complex (ORSC)

Crop wild relatives represent valuable reservoirs of variation for breeding, but their populations are threatened in natural habitats, are sparsely represented in genebanks, and most are poorly characterized. The focus of this study is the Oryza rufipogon species complex (ORSC), wild progenitor of Asian rice (Oryza sativa L.). The ORSC comprises perennial, annual and intermediate forms which were historically designated as O. rufipogon, O. nivara, and O. sativa f. spontanea (or Oryza spp., an annual form of mixed O. rufipogon/O. nivara and O. sativa ancestry), respectively, based on non-standardized morphological, geographical, and/or ecologically-based species definitions and boundaries. Here, a collection of 240 diverse ORSC accessions, characterized by genotyping-by-sequencing (113,739 SNPs), was phenotyped for 44 traits associated with plant, panicle, and seed morphology in the screenhouse at the International Rice Research Institute, Philippines. These traits included heritable phenotypes often recorded as characterization data by genebanks. Over 100 of these ORSC accessions were also phenotyped in the greenhouse for 18 traits in Stuttgart, Arkansas, and 16 traits in Ithaca, New York, United States. We implemented a Bayesian Gaussian mixture model to infer accession groups from a subset of these phenotypic data and ascertained three phenotype-based group assignments. We used concordance between the genotypic subpopulations and these phenotype-based groups to identify a suite of phenotypic traits that could reliably differentiate the ORSC populations, whether measured in tropical or temperate regions. The traits provide insight into plant morphology, life history (perenniality versus annuality) and mating habit (self- versus cross-pollinated), and are largely consistent with genebank species designations. One phenotypic group contains predominantly O. rufipogon accessions characterized as perennial and largely out-crossing and one contains predominantly O. nivara accessions characterized as annual and largely inbreeding. From these groups, 42 "core" O. rufipogon and 25 "core" O. nivara accessions were identified for domestication studies. The third group, comprising 20% of our collection, has the most accessions identified as Oryza spp. (51.2%) and levels of O. sativa admixture accounting for more than 50% of the genome. This third group is potentially useful as a "pre-breeding" pool for breeders attempting to incorporate novel variation into elite breeding lines.

Divergence in flowering time is a major component contributing to reproductive isolation between two wild rice species (Oryza rufipogon and O. nivara)

It is of critical importance for our understanding of speciation process to determine the forms of reproductive isolation and their relative importance in species divergence. Oryza nivara and O. rufipogon are direct ancestors of Asian cultivated rice and a progenitor-daughter species pair. Investigating the reproductive isolation between them provides insights into plant speciation and helps understanding of the rice domestication. Here, we quantitatively measured the major components of reproductive isolation between the two species based on common garden and crossing experiments for three pairs of sympatric populations in Nepal, Cambodia and Laos. We revealed significant differences in the flowering times between species pairs, with O. nivara flowering much earlier than O. rufipogon. A very weak reduction in seed set but no reduction in F1 viability and fertility were detected for the crosses between species relative to those within species. Moreover, we detected asymmetrical compatibility between species and found that emasculation significantly decreased pollination success in O. nivara but not in O. rufipogon. Our study demonstrates that the divergence between O. nivara and O. rufipogon is maintained almost entirely by the difference in flowering times and suggests that differential flowering times contribute to both habitat preferences and reproductive isolation between species.

A chromosome-level genome assembly of the wild rice Oryza rufipogon facilitates tracing the origins of Asian cultivated rice

Oryza rufipogon Griff. is a wild progenitor of the Asian cultivated rice Oryza sativa. To better understand the genomic diversity of the wild rice, high-quality reference genomes of O. rufipogon populations are needed, which also facilitate utilization of the wild genetic resources in rice breeding. In this study, we generated a chromosome-level genome assembly of O. rufipogon using a combination of short-read sequencing, single-molecule sequencing, BioNano and Hi-C platforms. The genome sequence (399.8 Mb) was assembled into 46 scaffolds on the 12 chromosomes, with contig N50 and scaffold N50 of 13.2 Mb and 20.3 Mb, respectively. The genome contains 36,520 protein-coding genes, and 49.37% of the genome consists of repetitive elements. The genome has strong synteny with those of the O. sativa subspecies indica and japonica, but containing some large structural variations. Evolutionary analysis unveiled the polyphyletic origins of O. sativa, in which the japonica and indica genome formations involved different divergent O. rufipogon (including O. nivara) lineages, accompanied by introgression of genomic regions between japonica and indica. This high-quality reference genome provides insight on the genome evolution of the wild rice and the origins of the O. sativa subspecies, and valuable information for basic research and rice breeding.

Novel allelic variation in the Phospholipase D alpha1 gene (OsPLDα1) of wild Oryza species implies to its low expression in rice bran

Rice bran, a by-product after milling, is a rich source of phytonutrients like oryzanols, tocopherols, tocotrienols, phytosterols, and dietary fibers. Moreover, exceptional properties of the rice bran oil make it unparalleled to other vegetable oils. However, a lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) causes rancidity and ‘stale flavor’ in the oil, and thus limits the rice bran usage for human consumption. To improve the rice bran quality, sequence based allele mining at OsPLDα1 locus (3.6 Kb) was performed across 48 accessions representing 11 wild Oryza species, 8 accessions of African cultivated rice, and 7 Oryza sativa cultivars. From comparative sequence analysis, 216 SNPs and 30 InDels were detected at the OsPLDα1 locus. Phylogenetic analysis revealed 20 OsPLDα1 cDNA variants which further translated into 12 protein variants. The O. officinalis protein variant, when compared to Nipponbare, showed maximum variability comprising 22 amino acid substitutions and absence of two peptides and two β-sheets. Further, expression profiling indicated significant differences in transcript abundance within as well as between the OsPLDα1 variants. Also, a new OsPLDα1 transcript variant having third exon missing in it, Os01t0172400-06, has been revealed. An O. officinalis accession (IRGC101152) had lowest gene expression which suggests the presence of novel allele, named as OsPLDα1-1a (GenBank accession no. MF966931). The identified novel allele could be further deployed in the breeding programs to overcome rice bran rancidity in elite cultivars.

Identification of Anther Length QTL and Construction of Chromosome Segment Substitution Lines of Oryza longistaminata

Life histories and breeding systems strongly affect the genetic diversity of seed plants, but the genetic architectures that promote outcrossing in Oryza longistaminata, a perennial wild species in Africa, are not understood. We conducted a genetic analysis of the anther length of O. longistaminata accession W1508 using advanced backcross quantitative trait locus (QTL) analysis and chromosomal segment substitution lines (CSSLs) in the genetic background of O. sativa Taichung 65 (T65), with simple sequence repeat markers. QTL analysis of the BC₃F₁ population (n = 100) revealed that four main QTL regions on chromosomes 3, 5, and 6 were associated to anther length. We selected a minimum set of BC₃F₂ plants for the development of CSSLs to cover as much of the W1508 genome as possible. The additional minor QTLs were suggested in the regional QTL analysis, using 21 to 24 plants in each of the selected BC₃F₂ population. The main QTLs found on chromosomes 3, 5, and 6 were validated and designated qATL3, qATL5, qATL6.1, and qATL6.2, as novel QTLs identified in O. longistaminata in the mapping populations of 94, 88, 70, and 95 BC₃F₄ plants. qATL3, qATL5, and qATL6.1 likely contributed to anther length by cell elongation, whereas qATL6.2 likely contributed by cell multiplication. The QTLs were confirmed again in an evaluation of the W1508ILs. In several chromosome segment substitution lines without the four validated QTLs, the anthers were also longer than those of T65, suggesting that other QTLs also increase anther length in W1508. The cloning and diversity analyses of genes conferring anther length QTLs promotes utilization of the genetic resources of wild species, and the understanding of haplotype evolution on the differentiation of annuality and perenniality in the genus Oryza.

Signatures of differential selection in chloroplast genome between japonica and indica

BACKGROUND

The domestication process of Asian rice (Oryza sativa L.) is complicated. It's well established that Oryza rufipogon is the ancestor of Asian rice, although the number of domestication events still controversial. Recently, numerous types of studies based on rice nuclear genome have been conducted, but the results are quite different. Chloroplasts (cp) are also part of the rice genome and have a conserved cyclic structure that is valuable for plant genetics and evolutionary studies. Therefore, we conducted chloroplast-based studies, aiming to provide more evidence for the domestication of Asian rice.

RESULTS

A total of 1389 variants were detected from the chloroplast genomes of 412 accessions obtained through the world. Oryza sativa L. ssp. japonica exhibited slightly less diversity (π) than Oryza sativa L. indica and wild rice. The fixation index values (F_ST) revealed that indica and japonica exhibited farther genetic distances compared with wild rice. Across cp genome, Tajima's D test demonstrated that different selection sites occurred in Asian rice. Principal component analyses (PCA) and multidimensional scaling (MDS) clearly classify the Asian rice into different groups. Furthermore, introgression patterns identified that indica and japonica shared no introgression events in cp level, and phylogenetic studies showed cultivated rice were well separated from different type of wild rice.

CONCLUSIONS

Here, we focus on the domestication of Asian rice (indica and japonica). Diversity and phylogenetic analyses revealed some selection characteristics in the chloroplast genome that potentially occurred in different Asian rice during the domestication. The results shown that Asian rice had been domesticated at least twice. In additional, japonica may experience a strong positive selection or bottleneck event during the domestication.

Parallel Speciation of Wild Rice Associated with Habitat Shifts

The occurrence of parallel speciation strongly implies the action of natural selection. However, it is unclear how general a phenomena parallel speciation is since it was only shown in a small number of animal species. In particular, the adaptive process and mechanisms underlying the process of parallel speciation remain elusive. Here, we used an integrative approach incorporating population genomics, common garden, and crossing experiments to investigate parallel speciation of the wild rice species Oryza nivara from O. rufipogon. We demonstrated that O. nivara originated multiple times from different O. rufipogon populations and revealed that different O. nivara populations have evolved similar phenotypes under divergent selection, a reflection of recurrent local adaptation of ancient O. rufipogon populations to dry habitats. Almost completed premating isolation was detected between O. nivara and O. rufipogon in the absence of any postmating barriers between and within these species. These results suggest that flowering time is a “magic” trait that contributes to both local adaptation and reproductive isolation in the origin of wild rice species. Our study thus demonstrates a convincing case of parallel ecological speciation as a consequence of adaptation to new environments.

An overview on reproductive isolation in Oryza sativa complex

Reproductive isolation is generally regarded as the essence of the speciation process. Studying closely related species is convenient for understanding the genetic basis of reproductive isolation. Therefore, the present review is restricted to the species and subspecies of the Oryza sativa complex, which includes the two domestic rice cultivars and six wild species. Although closely related, these rice species are separated from each other by a range reproductive barriers. This review presents a comprehensive understanding of the forces that shaped the formation of reproductive barriers among and between the species of the O. sativa complex. We suggest the possibility that domestication and artificial breeding in these rice species can lead to the early stages of speciation. Understanding the evolutionary and molecular mechanisms underlying reproductive isolation in rice will increase our knowledge in speciation and would also offer practical significance for the implementation of crop improvement strategies.

Population genetics and evolutionary history of the wild rice species Oryza rufipogon and O. nivara in Sri Lanka

Genetic diversity and population genetic structure of the wild rice species Oryza rufipogon and O. nivara in Sri Lanka were studied using 33 microsatellite markers. A total of 315 individuals of 11 natural populations collected from the wet, intermediate, and dry zones of the country were used in the study. We found a moderate to high level of genetic diversity at the population level, with the polymorphic loci (P) ranging from 60.6% to 100% (average 81.8%) and the expected heterozygosity (H E) varying from 0.294 to 0.481 (average 0.369). A significant genetic differentiation between species and strong genetic structure within species were also observed. Based on species distribution modeling, we detected the dynamics of the preferred habitats for the two species in Sri Lanka and demonstrated that both O. rufipogon and O. nivara populations have expanded substantially since the last internal glacial. In addition, we showed that the geographical distribution of the two species corresponded to the climate zones and identified a few of key environmental variables that contribute to the distribution of the two species, implying the potential mechanism for ecological adaptation of these two species in Sri Lanka. These studies provided important insights into the population genetics and evolution of these wild species in Sri Lanka and are of great significance to the in situ conservation and utilization of these wild resources in genetic improvement of rice.

Morphological and molecular dissection of wild rices from eastern India suggests distinct speciation between O. rufipogon and O. nivara populations

The inter relationships between the two progenitors is interesting as both wild relatives are known to be the great untapped gene reservoirs. The debate continues on granting a separate species status to Oryza nivara. The present study was conducted on populations of Oryza rufipogon and Oryza nivara from Eastern India employing morphological and molecular characteristics. The cluster analysis of the data on morphological traits could clearly classify the two wild forms into two separate discrete groups without any overlaps i.e. lack of intermediate forms, suggesting the non-sympatric existence of the wild forms. Amplification of hyper variable regions of the genome could reveal 144 alleles suggesting high genetic diversity values (average He = 0.566). Moreover, with 42.37% of uncommon alleles between the two wild relatives, the molecular variance analysis (AMOVA) could detect only 21% of total variation (p < 0.001) among them and rest 59% was within them. The population structure analysis clearly classified these two wild populations into two distinct sub-populations (K = 2) without any overlaps i.e. lack of intermediate forms, suggesting the non-sympatric existence of the wild forms. Clear differentiation into two distinct groups indicates that O. rufipogon and O. nivara could be treated as two different species.

Morphological and Molecular Data Reveal Three Distinct Populations of Indian Wild Rice Oryza rufipogon Griff. Species Complex

Wild relatives of crops possess adaptive mutations for agronomically important traits, which could play significant role in crop improvement for sustainable agriculture. However, global climate change and human activities pose serious threats to the natural habitats leading to erosion of genetic diversity of wild rice populations. The purpose of this study was to explore and characterize India's huge untapped wild rice diversity in Oryza rufipogon Griff. species complex from a wide range of ecological niches. We made strategic expeditions around diversity hot spots in 64 districts of nine different agro-climatic zones of the country and collected 418 wild rice accessions. Significant variation was observed among the accessions for 46 morphological descriptors, allowing classification into O. nivara, O. rufipogon, and O. sativa f. spontanea morpho-taxonomic groups. Genome-specific pSINE1 markers confirmed all the accessions having AA genome, which were further classified using ecotype-specific pSINE1 markers into annual, perennial, intermediate, and an unknown type. Principal component analysis revealed continuous variation for the morphological traits in each ecotype group. Genetic diversity analysis based on multi-allelic SSR markers clustered these accessions into three major groups and analysis of molecular variance for nine agro-climatic zones showed that 68% of the genetic variation was inherent amongst individuals while only 11% of the variation separated the zones, though there was significant correlation between genetic and spatial distances of the accessions. Model based population structure analysis using genome wide bi-allelic SNP markers revealed three sub-populations designated 'Pro-Indica,' 'Pro-Aus,' and 'Mid-Gangetic,' which showed poor correspondence with the morpho-taxonomic classification or pSINE1 ecotypes. There was Pan-India distribution of the 'Pro-Indica' and 'Pro-Aus' sub-populations across agro-climatic zones, indicating a more fundamental grouping based on the ancestry closely related to 'Indica' and 'Aus' groups of rice cultivars. The Pro-Indica population has substantial presence in the Eastern Himalayan Region and Lower Gangetic Plains, whereas 'Pro-Aus' sub-population was predominant in the Upper Gangetic Plains, Western Himalayan Region, Gujarat Plains and Hills, and Western Coastal Plains. In contrast 'Mid-Gangetic' population was largely concentrated in the Mid Gangetic Plains. The information presented here will be useful in the utilization of wild rice resources for varietal improvement.

Origin of rice (Oryza sativa L.) domestication genes

A number of genes that contribute to the domestication traits of cultivated rice have been identified. These include Sh4, Rc, PROG1 and LABA1, which are associated with non-shattering rachis, white pericarp, erect growth and barbless awns, respectively. The mutations giving rise to the "domestication alleles" of these genes are either invariable in cultivated rice, or have variability that is strictly associated with the phenotypic trait. This observation forms the basis to those current rice domestication models that envisage a single origin for the domesticated phenotype. Such models assume that the domestication alleles are absent or rare in wild rice, emerged under cultivation and spread across all rice groups by introgressive hybridization. We examined whole-genome sequencing datasets for wild and cultivated rice to test the former two assumptions. We found that the rc and laba1 alleles occur in wild rice with broad geographical distribution, and reach frequencies as high as 13 and 15%, respectively. These results are in agreement with previous observations of the prog1 and sh4 domestication alleles in wild populations. We also show that the diversity of the genomic regions surrounding the rc, laba1, prog1 and sh4 alleles in wild accessions is greater than that in cultivated rice, suggesting that these alleles emerged prior to domestication. Our findings indicate that the possibility that independent rice groups obtained identical domestication alleles directly from the wild population needs to be considered.

The Rice Paradox: Multiple Origins but Single Domestication in Asian Rice

The origin of domesticated Asian rice (Oryza sativa) has been a contentious topic, with conflicting evidence for either single or multiple domestication of this key crop species. We examined the evolutionary history of domesticated rice by analyzing de novo assembled genomes from domesticated rice and its wild progenitors. Our results indicate multiple origins, where each domesticated rice subpopulation (japonica, indica, and aus) arose separately from progenitor O. rufipogon and/or O. nivara. Coalescence-based modeling of demographic parameters estimate that the first domesticated rice population to split off from O. rufipogon was O. sativa ssp. japonica, occurring at ∼13.1-24.1 ka, which is an order of magnitude older then the earliest archeological date of domestication. This date is consistent, however, with the expansion of O. rufipogon populations after the Last Glacial Maximum ∼18 ka and archeological evidence for early wild rice management in China. We also show that there is significant gene flow from japonica to both indica (∼17%) and aus (∼15%), which led to the transfer of domestication alleles from early-domesticated japonica to proto-indica and proto-aus populations. Our results provide support for a model in which different rice subspecies had separate origins, but that de novo domestication occurred only once, in O. sativa ssp. japonica, and introgressive hybridization from early japonica to proto-indica and proto-aus led to domesticated indica and aus rice.

WRKY transcription factor genes in wild rice Oryza nivara

The WRKY transcription factor family is one of the largest gene families involved in plant development and stress response. Although many WRKY genes have been studied in cultivated rice (Oryza sativa), the WRKY genes in the wild rice species Oryza nivara, the direct progenitor of O. sativa, have not been studied. O. nivara shows abundant genetic diversity and elite drought and disease resistance features. Herein, a total of 97 O. nivara WRKY (OnWRKY) genes were identified. RNA-sequencing demonstrates that OnWRKY genes were generally expressed at higher levels in the roots of 30-day-old plants. Bioinformatic analyses suggest that most of OnWRKY genes could be induced by salicylic acid, abscisic acid, and drought. Abundant potential MAPK phosphorylation sites in OnWRKYs suggest that activities of most OnWRKYs can be regulated by phosphorylation. Phylogenetic analyses of OnWRKYs support a novel hypothesis that ancient group IIc OnWRKYs were the original ancestors of only some group IIc and group III WRKYs. The analyses also offer strong support that group IIc OnWRKYs containing the HVE sequence in their zinc finger motifs were derived from group Ia WRKYs. This study provides a solid foundation for the study of the evolution and functions of WRKY genes in O. nivara.