Diversification in flowering time due to tandem FT-like gene duplication, generating novel Mendelian factors in wild and cultivated rice

The complex structure of a single Mendelian factor widespread in the Asian cultivated rice ( Oryza sativa ) and its wild progenitor ( Oryza rufipogon ) that caused diverse phenotypes in the timing of flowering under natural field conditions was investigated in near isogenic lines. These near isogenic lines showed differences in flowering time despite all eight accessions collected from tropical regions possessing a recessive gene allelic to the se-pat gene. Fine mapping in two of these near-isogenic lines revealed that cultivated (Patpaku) and wild (W593) accessions had three and two linked quantitative trait loci (QTL) in the candidate regions, respectively, showing that Patpaku and W593 possessed linked QTLs with different effects in addition to the commonly-observed recessive gene ( se-pat ). Molecular dissection suggested that the tandemly duplicated FT-like genes ( Hd3a and RFT1 ) could be the candidate genes for these QTLs. Interestingly, the linked QTLs differed in their epistases, degree of dominance, and genotype × environment interactions. The nucleotide sequences showed that RFT1 has diverged more rapidly than Hd3a during rice evolution, suggesting phenotypic diversification of the two genes. Phylogenetic analysis implied that the se-pat ⁺ alleles might have emerged in different lineages within O. sativa . The present results strongly suggest that nucleotide divergence and shuffling of the linked QTLs by recombination might have created novel Mendelian factors that probably contribute to responding to local environments.     

Levels and patterns of nucleotide variation in domestication QTL regions on rice chromosome 3 suggest lineage-specific selection

Oryza sativa or Asian cultivated rice is one of the major cereal grass species domesticated for human food use during the Neolithic. Domestication of this species from the wild grass Oryza rufipogon was accompanied by changes in several traits, including seed shattering, percent seed set, tillering, grain weight, and flowering time. Quantitative trait locus (QTL) mapping has identified three genomic regions in chromosome 3 that appear to be associated with these traits. We would like to study whether these regions show signatures of selection and whether the same genetic basis underlies the domestication of different rice varieties. Fragments of 88 genes spanning these three genomic regions were sequenced from multiple accessions of two major varietal groups in O. sativa--indica and tropical japonica--as well as the ancestral wild rice species O. rufipogon. In tropical japonica, the levels of nucleotide variation in these three QTL regions are significantly lower compared to genome-wide levels, and coalescent simulations based on a complex demographic model of rice domestication indicate that these patterns are consistent with selection. In contrast, there is no significant reduction in nucleotide diversity in the homologous regions in indica rice. These results suggest that there are differences in the genetic and selective basis for domestication between these two Asian rice varietal groups.     

Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes

Rice is a staple crop that has undergone substantial phenotypic and physiological changes during domestication. Here we resequenced the genomes of 40 cultivated accessions selected from the major groups of rice and 10 accessions of their wild progenitors (Oryza rufipogon and Oryza nivara) to >15 × raw data coverage. We investigated genome-wide variation patterns in rice and obtained 6.5 million high-quality single nucleotide polymorphisms (SNPs) after excluding sites with missing data in any accession. Using these population SNP data, we identified thousands of genes with significantly lower diversity in cultivated but not wild rice, which represent candidate regions selected during domestication. Some of these variants are associated with important biological features, whereas others have yet to be functionally characterized. The molecular markers we have identified should be valuable for breeding and for identifying agronomically important genes in rice.     

Genomic adaptation of flowering‐time genes during the expansion of rice cultivation area

The diversification of flowering time in response to natural environments is critical for the spread of crops to diverse geographic regions. In contrast with recent advances in understanding the molecular basis of photoperiodic flowering in rice (Oryza sativa), little is known about how flowering‐time diversification is structured within rice subspecies. By analyzing genome sequencing data and a set of 429 chromosome segment substitution lines (CSSLs) originating from 10 diverse rice accessions with wide distributions, we revealed diverse effects of allelic variations for common flowering‐time quantitative trait loci in the recipient's background. Although functional variations associated with a few loci corresponded to standing variations among subspecies, the identified functional nucleotide polymorphisms occurred recently after rice subgroup differentiation, indicating that the functional diversity of flowering‐time gene sequences was not particularly associated with phylogenetic relationship between rice subspecies. Intensive analysis of the Hd1 genomic region identified the signature of an early introgression of the Hd1 with key mutation(s) in aus and temperate japonica accessions. Our data suggested that, after such key introgressions, new mutations were selected and accelerated the flowering‐time diversity within subspecies during the expansion of rice cultivation area. This finding may imply that new genome‐wide changes for flowering‐time adaptation are one of the critical determinants for establishing genomic architecture of local rice subgroups. In‐depth analyses of various rice genomes coupling with the genetically confirmed phenotypic changes in a large set of CSSLs enabled us to demonstrate how rice genome dynamics has coordinated with the adaptation of cultivated rice during the expansion of cultivation area.     

Genetic differentiation of wild relatives of rice as assessed by RFLP analysis

To study genetic diversity and relationships of wild relatives of rice, 58 accessions of Oryza rufipogon, Oryza nivara, Oryza sativa f. spontanea and the cultivated Oryza sativa, representing a wide range of their distribution, were analyzed using the restriction fragment length polymorphism (RFLP) technique. All 30-used RFLP probes detected polymorphisms among the Oryza accessions, with an average of 3.8 polymorphic fragments per probe. Considerable genetic diversity was scored among the Oryza accessions with a similarity coefficient ranging from 0.28 to 0.93; but the cluster analysis of the accessions did not show an apparent grouping based on the species classification, instead they were scattered randomly in different groups. Noticeably, the Oryza accessions from the same geographic region, or near-by geographic regions, tended to be clustered in the same groups. The indica rice varieties showed relatively high genetic diversity and were scattered in different groups of their wild relatives, but the japonica varieties showed a relatively low variation and formed an independent group. It is concluded from the molecular analytical result that: (1) the four Oryza taxa have a remarkably close relationship and their independent species status need to be carefully reviewed; (2) geographic isolation has played a significant role in the differentiation of the Oryza accessions; therefore, a wide geographic range needs to be covered in collecting wild rice germplasm for ex situ conservation; and (3) the conventional conclusion of indica rice being directly domesticated from its ancestral wild species, and japonica rice being derived from indica, gains support from our data.     

Multiple introgression events surrounding the Hd1 flowering-time gene in cultivated rice, Oryza sativa L.

Flowering time is a major determinant for the local adaptation of crops. Hd1 is a key flowering-time gene in rice and is orthologous to the Arabidopsis CONSTANS gene. To elucidate the role of Hd1 in selection, we examined the Hd1 alleles of 60 landraces of Asian cultivated rice (Oryza sativa L.) originating from all regions of Asia, which comprised three cultivar groups, indica, japonica, and aus. The identified alleles were classified into four allele groups. The functional Hd1 alleles in allele groups I and II corresponded to indica and japonica, respectively. Non-functional alleles in these groups were not clearly associated with cultivar groups or locations. Allele groups III and IV corresponded to the aus cultivar group. The ancestry of each cultivar group was identified by the coalescent approach for Hd1 molecular evolution using the haplotype patterns of 14 regions over the 1.1 Mb chromosomal region surrounding Hd1 and the pSINE patterns of two loci, 1.4 and 4.4 Mb apart from Hd1. The haplotype patterns clearly revealed that Hd1 allele migration was caused by multiple and complex introgression events between cultivar groups. The Hd1 haplotypes among dozens of accessions of the wild species O. rufipogon were strongly divergent and only two of the haplotype clusters in O. rufipogon were closely related to those in cultivated rice. This strongly suggested that multiple introgression events have played an important role in the shaping and diversification of adaptation in addition to primary selection steps at the beginning of domestication.     

Genetic control of a transition from black to straw-white seed hull in rice domestication

The genetic mechanism involved in a transition from the black-colored seed hull of the ancestral wild rice (Oryza rufipogon and Oryza nivara) to the straw-white seed hull of cultivated rice (Oryza sativa) during grain ripening remains unknown. We report that the black hull of O. rufipogon was controlled by the Black hull4 (Bh4) gene, which was fine-mapped to an 8.8-kb region on rice chromosome 4 using a cross between O. rufipogon W1943 (black hull) and O. sativa indica cv Guangluai 4 (straw-white hull). Bh4 encodes an amino acid transporter. A 22-bp deletion within exon 3 of the bh4 variant disrupted the Bh4 function, leading to the straw-white hull in cultivated rice. Transgenic study indicated that Bh4 could restore the black pigment on hulls in cv Guangluai 4 and Kasalath. Bh4 sequence alignment of all taxa with the outgroup Oryza barthii showed that the wild rice maintained comparable levels of nucleotide diversity that were about 70 times higher than those in the cultivated rice. The results from the maximum likelihood Hudson-Kreitman-Aguade test suggested that the significant reduction in nucleotide diversity in rice cultivars could be caused by artificial selection. We propose that the straw-white hull was selected as an important visual phenotype of nonshattered grains during rice domestication.     

Sequence polymorphisms in wild,weedy, and cultivated rice suggest seed‐shattering locus sh4 played a minor role in Asian rice domestication

The predominant view regarding Asian rice domestication is that the initial origin of nonshattering involved a single gene of large effect, specifically, the sh4 locus via the evolutionary replacement of a dominant allele for shattering with a recessive allele for reduced shattering. Data have accumulated to challenge this hypothesis. Specifically, a few studies have reported occasional seed‐shattering plants from populations of the wild progenitor of cultivated rice (Oryza rufipogon complex) being homozygous for the putative “nonshattering” sh4 alleles. We tested the sh4 hypothesis for the domestication of cultivated rice by obtaining genotypes and phenotypes for a diverse set of samples of wild, weedy, and cultivated rice accessions. The cultivars were fixed for the putative “nonshattering” allele and nonshattering phenotype, but wild rice accessions are highly polymorphic for the putative “nonshattering” allele (frequency ~26%) with shattering phenotype. All weedy rice accessions are the “nonshattering” genotype at the sh4 locus but with shattering phenotype. These data challenge the widely accepted hypothesis that a single nucleotide mutation (“G”/“T”) of the sh4 locus is the major driving force for rice domestication. Instead, we hypothesize that unidentified shattering loci are responsible for the initial domestication of cultivated rice through reduced seed shattering.     

Phylogeny and origin of pearl millet (Pennisetum glaucum [L.] R. Br) as revealed by microsatellite loci

During the last 12,000 years, different cultures around the world have domesticated cereal crops. Several studies investigated the evolutionary history and domestication of cereals such as wheat in the Middle East, rice in Asia or maize in America. The domestication process in Africa has led to the emergence of important cereal crops like pearl millet in Sahelian Africa. In this study, we used 27 microsatellite loci to analyze 84 wild accessions and 355 cultivated accessions originating from the whole pearl millet distribution area in Africa and Asia. We found significantly higher diversity in the wild pearl millet group. The cultivated pearl millet sample possessed 81% of the alleles and 83% of the genetic diversity of the wild pearl millet sample. Using Bayesian approaches, we identified intermediate genotypes between the cultivated and wild groups. We then analyzed the phylogenetic relationship among accessions not showing introgression and found that a monophyletic origin of cultivated pearl millet in West Africa is the most likely scenario supported by our data set.     

Natural Variation of the RICE FLOWERING LOCUS T 1 Contributes to Flowering Time Divergence in Rice

In rice (Oryza sativa L.), there is a diversity in flowering time that is strictly genetically regulated. Some indica cultivars show extremely late flowering under long-day conditions, but little is known about the gene(s) involved. Here, we demonstrate that functional defects in the florigen gene RFT1 are the main cause of late flowering in an indica cultivar, Nona Bokra. Mapping and complementation studies revealed that sequence polymorphisms in the RFT1 regulatory and coding regions are likely to cause late flowering under long-day conditions. We detected polymorphisms in the promoter region that lead to reduced expression levels of RFT1. We also identified an amino acid substitution (E105K) that leads to a functional defect in Nona Bokra RFT1. Sequencing of the RFT1 region in rice accessions from a global collection showed that the E105K mutation is found only in indica, and indicated a strong association between the RFT1 haplotype and extremely late flowering in a functional Hd1 background. Furthermore, SNPs in the regulatory region of RFT1 and the E105K substitution in 1,397 accessions show strong linkage disequilibrium with a flowering time–associated SNP. Although the defective E105K allele of RFT1 (but not of another florigen gene, Hd3a) is found in many cultivars, relative rate tests revealed no evidence for differential rate of evolution of these genes. The ratios of nonsynonymous to synonymous substitutions suggest that the E105K mutation resulting in the defect in RFT1 occurred relatively recently. These findings indicate that natural mutations in RFT1 provide flowering time divergence under long-day conditions.     

水稻长日抑制基因PhyB的多样性及区域适应性初探

Phy 是在长日照条件下抑制水稻开花的关键基因,但目前对水稻PhyB基因的遗传基础还不清楚,研究其分子遗传机制,对于培育光周期适应性广的品种以及扩大水稻种植区域具有重要意义。本研究选择78份亚洲栽培稻(34份籼稻和44份粳稻)及47份野生稻进行测序,对Phy B基因的核苷酸多态性、单倍型进行分析,计算籼稻、粳稻和野生稻的遗传多样性。结果表明,Phy B基因共有28个单倍型,其中有2个高频率的单倍型分别存在于2个栽培稻亚种中。从Network图可以看出栽培稻分为2组(A组和B组),A组栽培稻包括全部的籼稻和4个粳稻个体,B组栽培稻全是粳稻品种。亲缘地理学分析发现,A、B两组栽培稻具有明显不同的地理分布格局,且A组和B组开花时间差异显著,说明Phy B基因的2个高频率单倍型在2个栽培稻亚种中具有区域适应性,Phy B基因在栽培稻中具有明显的驯化信号,随着水稻种植区域的扩大,进化出适应不同地域特有的等位基因,导致开花时间对不同地区的区域适应性及多样性。     

Domestication and geographic origin of Oryza sativa in China: insights from multilocus analysis of nucleotide variation of O. sativa and O. rufipogon

Previous studies have indicated that China is one of the domestication centres of Asian cultivated rice (Oryza sativa), and common wild rice (O. rufipogon) is the progenitor of O. sativa. However, the number of domestication times and the geographic origin of Asian cultivated rice in China are still under debate. In this study, 100 accessions of Asian cultivated rice and 111 accessions of common wild rice in China were selected to examine the relationship between O. sativa and O. rufipogon and thereby infer the domestication and evolution of O. sativa in China through sequence analyses of six gene regions, trnC‐ycf6 in chloroplast genomes, cox3 in mitochondrial genomes and ITS, Ehd1, Waxy, Hd1 in nuclear genomes. The results indicated that the two subspecies of O. sativa (indica and japonica) were domesticated independently from different populations of O. rufipogon with gene flow occurring later from japonica to indica; Southern China was the genetic diversity centre of O. rufipogon, and the Pearl River basin near the Tropic of Cancer was the domestication centre of O. sativa in China.     

Diversity and selective sweep in the OsAMT1;1 genomic region of rice

Background  

Ammonium is one of the major forms in which nitrogen is available for plant growth. OsAMT1;1 is a high-affinity ammonium transporter in rice (Oryza sativa L.), responsible for ammonium uptake at low nitrogen concentration. The expression pattern of the gene has been reported. However, variations in its nucleotides and the evolutionary pathway of its descent from wild progenitors are yet to be elucidated. In this study, nucleotide diversity of the gene OsAMT1;1 and the diversity pattern of seven gene fragments spanning a genomic region approximately 150 kb long surrounding the gene were surveyed by sequencing a panel of 216 rice accessions including both cultivated rice and wild relatives.     

Origins of functional nucleotide polymorphisms in a major quantitative trait locus, qLTG3-1, controlling low-temperature germinability in rice

qLTG3-1 is a major quantitative trait locus (QTL) controlling tolerance to low-temperature at the seed germination stage (termed low-temperature germinability) in rice using a population derived from the cross between Italica Livorno from Italy and Hayamasari from Japan. Map-based cloning identified that qLTG3-1 encodes a protein of unknown function. The molecular identification of this major QTL could make it possible to identify allelic variation and favorable alleles for rice breeding programs. The present study examined the identification of qLTG3-1 alleles and their distribution among 62 landraces of Asian cultivated rice (Oryza sativa L.) collected from 19 different countries, termed the rice core collection. In the coding region, a single non-synonymous substitution and 3 in-frame insertion/deletion polymorphisms (indels) were detected. The almost completely conserved protein alignment of qLTG3-1 was also identified among 5 Oryza species, suggesting that the function of qLTG3-1 is critical for seed germination or for rice growth by pleiotropic effects of the gene. The functional nucleotide polymorphisms (FNPs), a 71-bp deletion found in Hayamasari and an amino acid substitution found in Nipponbare, was identified in varieties from Japan. These alleles with FNPs might be adapted to rice cultivation in specific local conditions. The present results may contribute to the utilization of favorable alleles of qLTG3-1 for the improvement of low-temperature germinability in rice breeding programs.     

Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice

Varying degrees of reduction of genetic diversity in crops relative to their wild progenitors occurred during the process of domestication. Such information, however, has not been available for the Asian cultivated rice (Oryza sativa) despite its importance as a staple food and a model organism. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationship and demographic history of O. sativa and its close relatives (Oryza rufipogon and Oryza nivara), we investigated nucleotide diversity data from 10 unlinked nuclear loci in species-wide samples of these species. The results indicated that O. rufipogon and O. nivara possessed comparable levels of nucleotide variation ((sil) = 0.0077 approximately 0.0095) compared with the relatives of other crops. In contrast, nucleotide diversity of O. sativa was as low as (sil) = 0.0024 and even lower ((sil) = 0.0021 for indica and 0.0011 for japonica), if we consider the 2 subspecies separately. Overall, only 20-10% of the diversity in the wild species was retained in 2 subspecies of the cultivated rice (indica and japonica), respectively. Because statistic tests did not reject the assumption of neutrality for all 10 loci, we further used coalescent to simulate bottlenecks under various lengths and population sizes to better understand the domestication process. Consistent with the dramatic reduction in nucleotide diversity, we detected a severe domestication bottleneck and demonstrated that the sequence diversity currently found in the rice genome could be explained by a founding population of 1,500 individuals if the initial domestication event occurred over a 3,000-year period. Phylogenetic analyses revealed close genetic relationships and ambiguous species boundary of O. rufipogon and O. nivara, providing additional evidence to treat them as 2 ecotypes of a single species. Lowest linkage disequilibrium (LD) was found in the perennial O. rufipogon where the r(2) value dropped to a negligible level within 400 bp, and the highest in the japonica rice where LD extended to the entirely sequenced region ( approximately 900 bp), implying that LD mapping by genome scans may not be feasible in wild rice due to the high density of markers needed.     

Association between nonsynonymous mutations of starch synthase IIa and starch quality in rice (Oryza sativa)

Starch quality is one of the most important agronomic traits in Asian rice, Oryza sativa. Starch synthase IIa (SsIIa) is a major candidate gene for starch quality variation. Within SsIIa, three nonsynonymous mutations in exon 8 have been shown to affect enzyme activity when expressed in Escherichia coli. To search for the variation in SsIIa that is responsible for starch quality variation in rice, we sequenced the SsIIa exon 8 region and measured starch quality as starch disintegration in alkali for 289 accessions of cultivated rice and 57 accessions of its wild ancestor, Oryza rufipogon. A general linear model and nested clade analysis were used to identify the associations between the three nonsynonymous single nucleotide polymorphisms (SNPs) and starch quality. Among the three nonsynonymous SNPs, we found strong evidence of association at one nucleotide site ('SNP 3'), corresponding to a Leu/Phe replacement at codon 781. A second SNP, corresponding to a Val/Met replacement at codon 737, could potentially show an association with increased sample sizes. Variation in SsIIa enzyme activity is associated with the cohesiveness of rice grains when cooked, and our findings are consistent with selection for more cohesive grains during the domestication of tropical japonica rice.     

Evolution and association analysis of Ghd7 in rice

Plant height, heading date, and yield are the main targets for rice genetic improvement. Ghd7 is a pleiotropic gene that controls the aforementioned traits simultaneously. In this study, a rice germplasm collection of 104 accessions (Oryza sativa) and 3 wild rice varieties (O.rufipogon) was used to analyze the evolution and association of Ghd7 with plant height, heading date, and yield. Among the 104 accessions, 76 single nucleotide polymorphisms (SNPs) and six insertions and deletions were found within a 3932-bp DNA fragment of Ghd7. A higher pairwise π and θ in the promoter indicated a highly diversified promoter of Ghd7. Sixteen haplotypes and 8 types of Ghd7 protein were detected. SNP changes between haplotypes indicated that Ghd7 evolved from two distinct ancestral gene pools, and independent domestication processes were detected in indica and japonica varietals respectively. In addition to the previously reported premature stop mutation in the first exon of Ghd7, which caused phenotypic changes of multiple traits, we found another functional C/T mutation (SNP S_555) by structure-based association analysis. SNP S_555 is located in the promoter and was related to plant height probably by altering gene expression. Moreover, another seven SNP mutations in complete linkage were found to be associated with the number of spikelets per panicle, regardless of the photoperiod. These associations provide the potential for flexibility of Ghd7 application in rice breeding programs.     

Leveraging natural diversity: back through the bottleneck

Plant breeders have long recognized the existence of useful genetic variation in the wild ancestors of our domesticated crop species. In cultivated rice (Oryza sativa), crosses between high-yielding elite cultivars and low-yielding wild accessions often give rise to superior offspring, with wild alleles conferring increased performance in the context of the elite cultivar genetic background. Because the breeding value of wild germplasm cannot be determined by examining the performance of wild accessions, a phylogenetic approach is recommended to determine which interspecific combinations are most likely to be useful in a breeding program. As we deepen our understanding of how genetic diversity is partitioned within and between cultivated and wild gene pools of Oryza, breeders will have increased power to make predictions about the most efficient strategies for utilizing wild germplasm for rice improvement.