Whole Genome Sequencing of Elite Rice Cultivars as a Comprehensive Information Resource for Marker Assisted Selection

Current advances in sequencing technologies and bioinformatics revealed the genomic background of rice, a staple food for the poor people, and provided the basis to develop large genomic variation databases for thousands of cultivars. Proper analysis of this massive resource is expected to give novel insights into the structure, function, and evolution of the rice genome, and to aid the development of rice varieties through marker assisted selection or genomic selection. In this work we present sequencing and bioinformatics analyses of 104 rice varieties belonging to the major subspecies of Oryza sativa. We identified repetitive elements and recurrent copy number variation covering about 200 Mbp of the rice genome. Genotyping of over 18 million polymorphic locations within O. sativa allowed us to reconstruct the individual haplotype patterns shaping the genomic background of elite varieties used by farmers throughout the Americas. Based on a reconstruction of the alleles for the gene GBSSI, we could identify novel genetic markers for selection of varieties with high amylose content. We expect that both the analysis methods and the genomic information described here would be of great use for the rice research community and for other groups carrying on similar sequencing efforts in other crops.     

Isozymic variability of traditional rice (Oryza sativa L.) in Africa

Summary Eight hundred and thirtyone traditional varieties of rice, Oryza sativa L., were collected in Africa and analysed for their isozymic variability on 15 enzymatic systems, representing 37 presumed loci. There appears to be a correlation between the type of rice growing and the two groups Indica and Japonica. The degree of genetic diversity is nearly equal in African rice and the Asian one. Alleles due to introgression or mutational events were identified. The results suggest that the evolution of O. sativa is continuous in Africa by means of inter-subspecific or inter-specific crosses.     

Transmission of important chromosomal regions under selection revealed in rice pedigree breeding programs

Genetic analysis across a whole plant genome based on pedigree information offers considerable potential for enhancing genetic gain from plant breeding programs through quantitative trait loci (QTL) mapping and marker-assisted selection. Here, we report its application for graphically genotyping varieties used in Chinese japonica rice (Oryza sativa L.) pedigree breeding programs. We identified 34 important chromosomal regions from the founder parent that are under selection in the breeding programs, and by comparing donor genomic regions that are under selection with QTL locations of agronomic traits, we found that QTL clustered in important genomic regions, in accordance with association analyses of natural populations and other previous studies. The convergence of genomic regions under selection with QTL locations suggests that donor genomic regions harboring key genes/QTL for important agronomic traits have been selected by plant breeders since the 1950s from the founder rice plants. The results provide better understanding of the effects of selection in breeding programs on the traits of rice cultivars. They also provide potentially valuable information for enhancing rice breeding programs through screening candidate parents for targeted molecular markers, improving crop yield potential and identifying suitable genetic material for use in future breeding programs.     

Processes Underpinning Development and Maintenance of Diversity in Rice in West Africa: Evidence from Combining Morphological and Molecular Markers

We assessed the interplay of artificial and natural selection in rice adaptation in low-input farming systems in West Africa. Using 20 morphological traits and 176 molecular markers, 182 farmer varieties of rice (Oryza spp.) from 6 West African countries were characterized. Principal component analysis showed that the four botanical groups (Oryza sativa ssp. indica, O. sativa ssp. japonica, O. glaberrima, and interspecific farmer hybrids) exhibited different patterns of morphological diversity. Regarding O. glaberrima, morphological and molecular data were in greater conformity than for the other botanical groups. A clear difference in morphological features was observed between O. glaberrima rices from the Togo hills and those from the Upper Guinea Coast, and among O. glaberrima rices from the Upper Guinea Coast. For the other three groups such clear patterns were not observed. We argue that this is because genetic diversity is shaped by different environmental and socio-cultural selection pressures. For O. glaberrima, recent socio-cultural selection pressures seemed to restrict genetic diversity while this was not observed for the other botanical groups. We also show that O. glaberrima still plays an important role in the selection practices of farmers and resulting variety development pathways. This is particularly apparent in the case of interspecific farmer hybrids where a relationship was found between pericarp colour, panicle attitude and genetic diversity. Farmer varieties are the product of long and complex trajectories of selection governed by local human agency. In effect, rice varieties have emerged that are adapted to West African farming conditions through genotype × environment × society interactions. The diversity farmers maintain in their rice varieties is understood to be part of a risk-spreading strategy that also facilitates successful and often serendipitous variety innovations. We advocate, therefore, that farmers and farmer varieties should be more effectively involved in crop development.     

An-1 Encodes a Basic Helix-Loop-Helix Protein That Regulates Awn Development,Grain Size,and Grain Number in Rice

Long awns are important for seed dispersal in wild rice (Oryza rufipogon), but are absent in cultivated rice (Oryza sativa). The genetic mechanism involved in loss-of-awn in cultivated rice remains unknown. We report here the molecular cloning of a major quantitative trait locus, An-1, which regulates long awn formation in O. rufipogon. An-1 encodes a basic helix-loop-helix protein, which regulates cell division. The nearly-isogenic line (NIL-An-1) carrying a wild allele An-1 in the genetic background of the awnless indica Guangluai4 produces long awns and longer grains, but significantly fewer grains per panicle compared with Guangluai4. Transgenic studies confirmed that An-1 positively regulates awn elongation, but negatively regulates grain number per panicle. Genetic variations in the An-1 locus were found to be associated with awn loss in cultivated rice. Population genetic analysis of wild and cultivated rice showed a significant reduction in nucleotide diversity of the An-1 locus in rice cultivars, suggesting that the An-1 locus was a major target for artificial selection. Thus, we propose that awn loss was favored and strongly selected by humans, as genetic variations at the An-1 locus that cause awn loss would increase grain numbers and subsequently improve grain yield in cultivated rice.     

Recombinational and mutational analyses of subspecific differentiation in Oryza sativa L.

Recombinational and mutational analyses were performed in an attempt to elucidate the nature of the genic differences and the ancestral relationship between the indica and japonica subspecies of Oryza sativa L. by taking representative varieties from each group. The absence of cytological aberrations at meiosis in indica×japonica hybrids suggest that genetic differences are predominantly genic rather than chromosomal. The pattern of variation induced in the key characters, and the frequency and spectra of mutations obtained in each group, failed to support the hypothesis that indica has given rise to japonica, and suggested that the two subspecies have originated from a common ancestor.     

Genome-Wide Prediction of Traits with Different Genetic Architecture Through Efficient Variable Selection

In genome-based prediction there is considerable uncertainty about the statistical model and method required to maximize prediction accuracy. For traits influenced by a small number of quantitative trait loci (QTL), predictions are expected to benefit from methods performing variable selection [e.g., BayesB or the least absolute shrinkage and selection operator (LASSO)] compared to methods distributing effects across the genome [ridge regression best linear unbiased prediction (RR-BLUP)]. We investigate the assumptions underlying successful variable selection by combining computer simulations with large-scale experimental data sets from rice (Oryza sativa L.), wheat (Triticum aestivum L.), and Arabidopsis thaliana (L.). We demonstrate that variable selection can be successful when the number of phenotyped individuals is much larger than the number of causal mutations contributing to the trait. We show that the sample size required for efficient variable selection increases dramatically with decreasing trait heritabilities and increasing extent of linkage disequilibrium (LD). We contrast and discuss contradictory results from simulation and experimental studies with respect to superiority of variable selection methods over RR-BLUP. Our results demonstrate that due to long-range LD, medium heritabilities, and small sample sizes, superiority of variable selection methods cannot be expected in plant breeding populations even for traits like FRIGIDA gene expression in Arabidopsis and flowering time in rice, assumed to be influenced by a few major QTL. We extend our conclusions to the analysis of whole-genome sequence data and infer upper bounds for the number of causal mutations which can be identified by LASSO. Our results have major impact on the choice of statistical method needed to make credible inferences about genetic architecture and prediction accuracy of complex traits.     

Genetic differentiation among Sri Lankan traditional rice (Oryza sativa) varieties and wild rice species by AFLP markers

Traditional rice varieties are one important component of the biodiversity of Sri Lanka. However, no proper studies have been performed on genetic diversity of the Sri Lankan traditional rice varieties used in our breeding programs. In the present study, the genetic diversity of 46 traditional rice (Oryza sativa) varieties and 5 wild rice species is assesed using Amplified Fragment Length Polymorphism (AFLP) markers. Ten primer combinations generated a total of 784 fragments. Of these, 772 fragments were polymorphic (98.4%). UPGMA analysis based on Jaccard's similarity separated the accessions into four major clusters. A cophenetic correlation with r=0.786 strongly supported the clustering pattern of UPGMA dendrogram. A principal coordinate analysis (PCoA) also confirmed the UPGMA clusters. Accessions referred to the same cluster showed similar morphological characteristics (e.g. height, grain colour etc.) while accessions which are known to be morphologically distinct appeared genetically separated. In addition, the clustering pattern distinctly separated lowland and upland rice varieties. This genetic diversity assessment at the molecular level provides reliable information for selection of germplasm in the development of new rice varieties and in conservation of traditional rice genetic resources.     

Nonfunctional alleles of long‐day suppressor genes independently regulate flowering time

Due to the remarkable adaptability to various environments, rice varieties with diverse flowering times have been domesticated or improved from Oryza rufipogon. Detailed knowledge of the genetic factors controlling flowering time will facilitate understanding the adaptation mechanism in cultivated rice and enable breeders to design appropriate genotypes for distinct preferences. In this study, four genes (Hd1, DTH8, Ghd7 and OsPRR37) in a rice long‐day suppression pathway were collected and sequenced in 154, 74, 69 and 62 varieties of cultivated rice (Oryza sativa) respectively. Under long‐day conditions, varieties with nonfunctional alleles flowered significantly earlier than those with functional alleles. However, the four genes have different genetic effects in the regulation of flowering time: Hd1 and OsPRR37 are major genes that generally regulate rice flowering time for all varieties, while DTH8 and Ghd7 only regulate regional rice varieties. Geographic analysis and network studies suggested that the nonfunctional alleles of these suppression loci with regional adaptability were derived recently and independently. Alleles with regional adaptability should be taken into consideration for genetic improvement. The rich genetic variations in these four genes, which adapt rice to different environments, provide the flexibility needed for breeding rice varieties with diverse flowering times.     

The Birth of a Black Rice Gene and Its Local Spread by Introgression

The origin and spread of novel agronomic traits during crop domestication are complex events in plant evolution. Wild rice (Oryza rufipogon) has red grains due to the accumulation of proanthocyanidins, whereas most cultivated rice (Oryza sativa) varieties have white grains induced by a defective allele in the Rc basic helix-loop-helix (bHLH) gene. Although the events surrounding the origin and spread of black rice traits remain unknown, varieties with black grains due to anthocyanin accumulation are distributed in various locations throughout Asia. Here, we show that the black grain trait originated from ectopic expression of the Kala4 bHLH gene due to rearrangement in the promoter region. Both the Rc and Kala4 genes activate upstream flavonol biosynthesis genes, such as chalcone synthase and dihydroflavonol-4-reductase, and downstream genes, such as leucoanthocyanidin reductase and leucoanthocyanidin dioxygenase, to produce the respective specific pigments. Genome analysis of 21 black rice varieties as well as red- and white-grained landraces demonstrated that black rice arose in tropical japonica and its subsequent spread to the indica subspecies can be attributed to the causal alleles of Kala4. The relatively small size of genomic fragments of tropical japonica origin in some indica varieties indicates that refined introgression must have occurred by natural crossbreeding in the course of evolution of the black trait in rice.     

Rice domestication: histories and mysteries

Domesticated rice (Oryza sativa) is one of the world’s most important food crops, culturally, nutritionally and economically ( Khush 1997 ). Thus, it is no surprise that there is intense curiosity about its genetic and geographical origins, its response to selection under domestication, and the genetic structure of its wild relative, Oryza rufipogon. Studies of Oryza attempting to answer these questions have accompanied each stage of the development of molecular markers, starting with allozymes and continuing to genome sequencing. While many of these studies have been restricted to small sample sizes, in terms of either the number of markers used or the number and distribution of the accessions, costs are now low enough that researchers are including large numbers of molecular markers and accessions. How will these studies relate to previous findings and long‐held assumptions about rice domestication and evolution? If the paper in this issue of Molecular Ecology ( Huang et al. 2012 ) is any indication, there will be some considerable surprises in store. In this study, a geographically and genomically thorough sampling of O. rufipogon and O. sativa revealed two genetically distinct groups of wild rice and also indicated that only one of these groups appears to be related to domesticated rice. While this fits well with previous studies indicating that there are genetic subdivisions within O. rufipogon, it stands in contrast to previous findings that the two major varieties of O. sativa (indica and japonica) were domesticated from two (or more) subpopulations of wild rice.     

Growth promotion and inhibition of the Amazonian wild rice species Oryza grandiglumis to survive flooding

In Asian cultivated rice (Oryza sativa), distinct mechanisms to survive flooding are activated in two groups of varieties. Submergence-tolerant rice varieties possessing the SUBMERGENCE1A (SUB1A) gene display reduced growth during flash floods at the seedling stage and resume growth after the flood recedes, whereas deepwater rice varieties possessing the SNORKEL1 (SK1) and SNORKEL2 (SK2) genes display enhanced growth based on internodal elongation during prolonged submergence at the mature stage. In this study, we investigated the occurrence of these growth responses to submergence in the wild rice species Oryza grandiglumis, which is native to the Amazon floodplains. When subjected to gradual submergence, adult plants of O. grandiglumis accessions showed enhanced internodal elongation with rising water level and their growth response closely resembled that of deepwater varieties of O. sativa with high floating capacity. On the other hand, when subjected to complete submergence, seedlings of O. grandiglumis accessions displayed reduced shoot growth and resumed normal growth after desubmergence, similar to the response of submergence-tolerant varieties of O. sativa. Neither SUB1A nor the SK genes were detected in the O. grandiglumis accessions. These results indicate that the O. grandiglumis accessions are capable of adapting successfully to flooding by activating two contrasting mechanisms as the situation demands and that each mechanism of adaptation to flooding is not mediated by SUB1A or the SK genes.     

High-throughput single nucleotide polymorphism genotyping for breeding applications in rice using the BeadXpress platform

Multiplexed single nucleotide polymorphism (SNP) markers have the potential to increase the speed and cost-effectiveness of genotyping, provided that an optimal SNP density is used for each application. To test the efficiency of multiplexed SNP genotyping for diversity, mapping and breeding applications in rice (Oryza sativa L.), we designed seven GoldenGate VeraCode oligo pool assay (OPA) sets for the Illumina BeadXpress Reader. Validated markers from existing 1536 Illumina SNPs and 44?K Affymetrix SNP chips developed at Cornell University were used to select subsets of informative SNPs for different germplasm groups with even distribution across the genome. A 96-plex OPA was developed for quality control purposes and for assigning a sample into one of the five O. sativa population subgroups. Six 384-plex OPAs were designed for genetic diversity analysis, DNA fingerprinting, and to have evenly-spaced polymorphic markers for quantitative trait locus (QTL) mapping and background selection for crosses between different germplasm pools in rice: Indica/Indica, Indica/Japonica, Japonica/Japonica, Indica/O. rufipogon, and Japonica/O. rufipogon. After testing on a diverse set of rice varieties, two of the SNP sets were re-designed by replacing poor-performing SNPs. Pilot studies were successfully performed for diversity analysis, QTL mapping, marker-assisted backcrossing, and developing specialized genetic stocks, demonstrating that 384-plex SNP genotyping on the BeadXpress platform is a robust and efficient method for marker genotyping in rice.     

Genetic analysis of Indian aromatic and quality rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) germplasm using panels of fluorescently-labeled microsatellite markers

Genetic relationships among Indian aromatic and quality rice (Oryza sativa) germplasm were assessed using 30 fluorescently labeled rice microsatellite markers. The 69 rice genotypes used in this study included 52 Basmati and other scented/quality rice varieties from different parts of India and 17 indica and japonica varieties that served as controls. A total of 235 alleles were detected at the 30 simple sequence repeat (SSR) loci, 62 (26.4%) of which were present only in Basmati and other scented/quality rice germplasm accessions. The number of alleles per locus ranged from 3 to 22, with an average of 7.8, polymorphism information content (PIC) values ranged from 0.2 to 0.9, with an average of 0.6, and the size range between the smallest and the largest allele for a given microsatellite locus varied between 3 bp and 68 bp. Of the 30 SSR markers, 20 could distinguish traditional Basmati rice varieties, and a single panel of eight markers could be used to differentiate the premium traditional Basmati, cross-bred Basmati, and non-Basmati rice varieties having different commercial value in the marketplace. When estimates of inferred ancestry or similarity coefficients were used to cluster varieties, the high-quality Indian aromatic and quality rice genotypes could be distinguished from both indica and japonica cultivars, and crossbred varieties could be distinguished from traditional Basmati rices. The results indicate that Indian aromatic and quality germplasm is genetically distinct from other groups within O. sativa and is the product of a long, independent pattern of evolution. The data also suggest that there is scope for exploiting the genetic diversity of aromatic/quality rice germplasm available in India for national Basmati rice breeding programs.Electronic Supplementary Material Supplementary material is available for this article at .     

Influence of ethnic traditional cultures on genetic diversity of rice landraces under on-farm conservation in southwest China

Background

Crop genetic resources are important components of biodiversity. However, with the large-scale promotion of mono-cropping, genetic diversity has largely been lost. Ex-situ conservation approaches were widely used to protect traditional crop varieties worldwide. However, this method fails to maintain the dynamic evolutionary processes of crop genetic resources in their original habitats, leading to genetic diversity reduction and even loss of the capacity of resistance to new diseases and pests. Therefore, on-farm conservation has been considered a crucial complement to ex-situ conservation. This study aimed at clarifying the genetic diversity differences between ex-situ conservation and on-farm conservation and to exploring the influence of traditional cultures on genetic diversity of rice landraces under on-farm conservation.

Methods

The conservation status of rice landrace varieties, including Indica and Japonica, non-glutinous rice (Oryza sativa) and glutinous rice (Oryza sativa var. glutinosa Matsum), was obtained through ethno-biology investigation method in 12 villages of ethnic groups from Guizhou, Yunnan and Guangxi provinces of China. The genetic diversity between 24 pairs of the same rice landraces from different times were compared using simple sequence repeat (SSR) molecular markers technology. The landrace paris studied were collected in 1980 and maintained ex-situ, while 2014 samples were collected on-farm in southwest of China.

Results

The results showed that many varieties of rice landraces have been preserved on-farm by local farmers for hundreds or thousands of years. The number of alleles (Na), effective number of alleles (Ne), Nei genetic diversity index (He) and Shannon information index (I) of rice landraces were significantly higher by 12.3–30.4 % under on-farm conservation than under ex-situ conservation. Compared with the ex-situ conservation approach, rice landraces under on-farm conservation programs had more alleles and higher genetic diversity. In every site we investigated, ethnic traditional cultures play a positive influence on rice landrace variety diversity and genetic diversity.

Conclusion

Most China’s rice landraces were conserved in the ethnic areas of southwest China. On-farm conservation can effectively promote the allelic variation and increase the genetic diversity of rice landraces over the past 35 years. Moreover, ethnic traditional culture practices are a crucial foundation to increase genetic diversity of rice landraces and implement on-farm conservation.

     

Robustness and Strategies of Adaptation among Farmer Varieties of African Rice (Oryza glaberrima) and Asian Rice (Oryza sativa) across West Africa

This study offers evidence of the robustness of farmer rice varieties (Oryza glaberrima and O. sativa) in West Africa. Our experiments in five West African countries showed that farmer varieties were tolerant of sub-optimal conditions, but employed a range of strategies to cope with stress. Varieties belonging to the species Oryza glaberrima – solely the product of farmer agency – were the most successful in adapting to a range of adverse conditions. Some of the farmer selections from within the indica and japonica subspecies of O. sativa also performed well in a range of conditions, but other farmer selections from within these two subspecies were mainly limited to more specific niches. The results contradict the rather common belief that farmer varieties are only of local value. Farmer varieties should be considered by breeding programmes and used (alongside improved varieties) in dissemination projects for rural food security.     

Molecular Evolution of the TAC1 Gene from Rice （Oryza sativa L.）

Tiller angle is a key feature of the architecture of cultivated rice(Oryza sativa),since it determines planting density and influences rice yield.Our previous work identified Tiller Angle Control 1(TACl) as a major quantitative trait locus that controls rice tiller angle.To further clarify the evolutionary characterization of the TACl gene,we compared a TACl-containing 3164-bp genomic region among 113 cultivated varieties and 48 accessions of wild rice,including 43 accessions of O.rufipogon and five accessions of O.nivara.Only one single nucleotide polymorphism(SNP),a synonymous substitution,was detected in TACl coding regions of the cultivated rice varieties, whereas one synonymous and one nonsynonymous SNP were detected among the TACl coding regions of wild rice accessions.These data indicate that little natural mutation and modification in the TACl coding region occurred within the cultivated rice and its progenitor during evolution.Nucleotide diversities in the TACl gene regions of O.sativa and O.rufipogon of 0.00116 and 0.00112,respectively, further indicate that TACl has been highly conserved during the course of rice domestication.A functional nucleotide polymorphism (FNP) of TACl was only found in the japonica rice group.A neutrality test revealed strong selection,especially in the 3’-flanking region of the TACl coding region containing the FNP in the japonica rice group.However,no selection occurred in the indica and wild-rice groups.A phylogenetic tree derived from TACl sequence analysis suggests that the indica and japonica subspecies arose independently during the domestication of wild rice.     

Genomics-based precision breeding approaches to improve drought tolerance in rice

Rice (Oryza sativa L.), the major staple food crop of the world, faces a severe threat from widespread drought. The development of drought-tolerant rice varieties is considered a feasible option to counteract drought stress. The screening of rice germplasm under drought and its characterization at the morphological, genetic, and molecular levels revealed the existence of genetic variation for drought tolerance within the rice gene pool. The improvements made in managed drought screening and selection for grain yield under drought have significantly contributed to progress in drought breeding programs. The availability of rice genome sequence information, genome-wide molecular markers, and low-cost genotyping platforms now makes it possible to routinely apply marker-assisted breeding approaches to improve grain yield under drought. Grain yield QTLs with a large and consistent effect under drought have been indentified and successfully pyramided in popular rice mega-varieties. Various rice functional genomics resources, databases, tools, and recent advances in “-omics” are facilitating the characterization of genes and pathways involved in drought tolerance, providing the basis for candidate gene identification and allele mining. The transgenic approach is successful in generating drought tolerance in rice under controlled conditions, but field-level testing is necessary. Genomics-assisted drought breeding approaches hold great promise, but a well-planned integration with standardized phenotyping is highly essential to exploit their full potential.     

Genome-wide patterns of nucleotide polymorphism in domesticated rice

Domesticated Asian rice (Oryza sativa) is one of the oldest domesticated crop species in the world, having fed more people than any other plant in human history. We report the patterns of DNA sequence variation in rice and its wild ancestor, O. rufipogon, across 111 randomly chosen gene fragments, and use these to infer the evolutionary dynamics that led to the origins of rice. There is a genome-wide excess of high-frequency derived single nucleotide polymorphisms (SNPs) in O. sativa varieties, a pattern that has not been reported for other crop species. We developed several alternative models to explain contemporary patterns of polymorphisms in rice, including a (i) selectively neutral population bottleneck model, (ii) bottleneck plus migration model, (iii) multiple selective sweeps model, and (iv) bottleneck plus selective sweeps model. We find that a simple bottleneck model, which has been the dominant demographic model for domesticated species, cannot explain the derived nucleotide polymorphism site frequency spectrum in rice. Instead, a bottleneck model that incorporates selective sweeps, or a more complex demographic model that includes subdivision and gene flow, are more plausible explanations for patterns of variation in domesticated rice varieties. If selective sweeps are indeed the explanation for the observed nucleotide data of domesticated rice, it suggests that strong selection can leave its imprint on genome-wide polymorphism patterns, contrary to expectations that selection results only in a local signature of variation.