Large-scale DNA polymorphism study of Oryza sativa and O. rufipogon reveals the origin and divergence of Asian rice

Reverse genetic approaches for breeding nutrient-rich and climate-resilient cereal and food legume crops

In the last decade, advancements in genomics tools and techniques have led to the discovery of many genes. Most of these genes still need to be characterized for their associated function and therefore, such genes remain underutilized for breeding the next generation of improved crop varieties. The recent developments in different reverse genetic approaches have made it possible to identify the function of genes controlling nutritional, biochemical, and metabolic traits imparting drought, heat, cold, salinity tolerance as well as diseases and insect-pests. This article focuses on reviewing the current status and prospects of using reverse genetic approaches to breed nutrient-rich and climate resilient cereal and food legume crops.

Genetic Diversity, Population Structure and Linkage Disequilibrium Analyses in Tropical Maize Using Genotyping by Sequencing

Several maize breeding programs in India have developed numerous inbred lines but the lines have not been characterized using high-density molecular markers. Here, we studied the molecular diversity, population structure, and linkage disequilibrium (LD) patterns in a panel of 314 tropical normal corn, two sweet corn, and six popcorn inbred lines developed by 17 research centers in India, and 62 normal corn from the International Maize and Wheat Improvement Center (CIMMYT). The 384 inbred lines were genotyped with 60,227 polymorphic single nucleotide polymorphisms (SNPs). Most of the pair-wise relative kinship coefficients (58.5%) were equal or close to 0, which suggests the lack of redundancy in the genomic composition in the majority of inbred lines. Genetic distance among most pairs of lines (98.3%) varied from 0.20 to 0.34 as compared with just 1.7% of the pairs of lines that differed by <0.20, which suggests greater genetic variation even among sister lines. The overall average of 17% heterogeneity was observed in the panel indicated the need for further inbreeding in the high heterogeneous genotypes. The mean nucleotide diversity and frequency of polymorphic sites observed in the panel were 0.28 and 0.02, respectively. The model-based population structure, principal component analysis, and phylogenetic analysis revealed three to six groups with no clear patterns of clustering by centers-wise breeding lines, types of corn, kernel characteristics, maturity, plant height, and ear placement. However, genotypes were grouped partially based on their source germplasm from where they derived.

Multi-trait association study identifies loci associated with tolerance of low phosphorus in Oryza sativa and its wild relatives

We studied variation in adaptive traits and genetic association to understand the low P responses, including the symbiotic association of arbuscular mycorrhizal (AM) fungal colonization in Oryza species (O. sativa, O. nivara, and O. rufipogon). In the present experiment, we performed the phenotypic variability of the morphometric and geometric traits for P deficiency tolerance and conducted the association studies in GLM and MLM methods. A positive association between the geometric trait of the top-view area and root traits suggested the possibility of exploring a non-destructive approach in screening genotypes under low P. The AMOVA revealed a higher proportion of variation among the individuals as they belonged to different species of Oryza and the NM value was 2.0, indicating possible gene flow between populations. A sub-cluster with superior-performing accessions had a higher proportion of landraces (42.85%), and O. rufipogon (33.3%) was differentiated by four Pup1-specific markers. Association mapping identified seven notable markers (RM259, RM297, RM30, RM6966, RM242, RM184, and PAP1) and six potential genotypes (IC459373, Chakhao Aumbi, AC100219, AC100062, Sekri, and Kumbhi Phou), which will be helpful in the marker-assisted breeding to improve rice for P-deprived condition. In addition, total root surface area becomes a single major trait that helps in P uptake under deficit P up to 33% than mycorrhizal colonization. Further, the phenotypic analysis of the morphometric and geometric trait variations and their interactions provides excellent potential for selecting donors for improving P-use efficiency. The identified potential candidate genes and markers offered new insights into our understanding of the molecular and physiological mechanisms driving PUE and improving grain yield under low-P conditions.

Genetic Dissection of Grain Yield Component Traits Under High Nighttime Temperature Stress in a Rice Diversity Panel

To dissect the genetic complexity of rice grain yield (GY) and quality in response to heat stress at the reproductive stage, a diverse panel of 190 rice accessions in the United States Department of Agriculture (USDA) rice mini-core collection (URMC) diversity panel were treated with high nighttime temperature (HNT) stress at the reproductive stage of panicle initiation. The quantifiable yield component response traits were then measured. The traits, panicle length (PL), and number of spikelets per panicle (NSP) were evaluated in subsets of the panel comprising the rice subspecies Oryza sativa ssp. Indica and ssp. Japonica. Under HNT stress, the Japonica ssp. exhibited lower reductions in PL and NSP and a higher level of genetic variation compared with the other subpopulations. Whole genome sequencing identified 6.5 million single nucleotide polymorphisms (SNPs) that were used for the genome-wide association studies (GWASs) of the PL and NSP traits. The GWAS analysis in the Combined, Indica, and Japonica populations under HNT stress identified 83, 60, and 803 highly significant SNPs associated with PL, compared to the 30, 30, and 11 highly significant SNPs associated with NSP. Among these trait-associated SNPs, 140 were coincident with genomic regions previously reported for major GY component quantitative trait loci (QTLs) under heat stress. Using extents of linkage disequilibrium in the rice populations, Venn diagram analysis showed that the highest number of putative candidate genes were identified in the Japonica population, with 20 putative candidate genes being common in the Combined, Indica and Japonica populations. Network analysis of the genes linked to significant SNPs associated with PL and NSP identified modules that were involved in primary and secondary metabolisms. The findings in this study could be useful to understand the pathways/mechanisms involved in rice GY and its components under HNT stress for the acceleration of rice-breeding programs and further functional analysis by molecular geneticists.

A genome-wide association study in Indian wild rice accessions for resistance to the root-knot nematode Meloidogyne graminicola

Rice root-knot nematode (RRKN), Meloidogyne graminicola is one of the major biotic constraints in rice-growing countries of Southeast Asia. Host plant resistance is an environmentally-friendly and cost-effective mean to mitigate RRKN damage to rice. Considering the limited availability of genetic resources in the Asian rice (Oryza sativa) cultivars, exploration of novel sources and genetic basis of RRKN resistance is necessary. We screened 272 diverse wild rice accessions (O. nivara, O. rufipogon, O. sativa f. spontanea) to identify genotypes resistant to RRKN. We dissected the genetic basis of RRKN resistance using a genome-wide association study with SNPs (single nucleotide polymorphism) genotyped by 50K "OsSNPnks" genic Affymetrix chip. Population structure analysis revealed that these accessions were stratified into three major sub-populations. Overall, 40 resistant accessions (nematode gall number and multiplication factor/MF < 2) were identified, with 17 novel SNPs being significantly associated with phenotypic traits such as number of galls, egg masses, eggs/egg mass and MF per plant. SNPs were localized to the quantitative trait loci (QTL) on chromosome 1, 2, 3, 4, 6, 10 and 11 harboring the candidate genes including NBS-LRR, Cf2/Cf5 resistance protein, MYB, bZIP, ARF, SCARECROW and WRKY transcription factors. Expression of these identified genes was significantly (P < 0.01) upregulated in RRKN-infected plants compared to mock-inoculated plants at 7 days after inoculation. The identified SNPs enrich the repository of candidate genes for future marker-assisted breeding program to alleviate the damage of RRKN in rice.

Novel cis-acting regulatory elements in wild Oryza species impart improved rice bran quality by lowering the expression of phospholipase D alpha1 enzyme (OsPLDα1)

Rice bran oil is good quality edible oil, rich in antioxidants and comprised typically of oleic-linoleic type fatty acids. However, presence of a highly lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) increases free fatty acid content in the oil which further leads to stale flavor and rancidity of the oil, making it unfit for human consumption. In this study, we compared the upstream regions of OsPLDα1 orthologs across 34 accessions representing 5 wild Oryza species and 8 cultivars, to uncover sequence variations and identify cis-elements involved in differential transcription of orthologs. Alignment of the upstream sequences to the Nipponbare OsPLDα1 reference sequence revealed the presence of 39 SNPs. Phylogenetic analysis showed that all the selected cultivars and wild species accessions are closely related to the reference except for three accessions of O. rufipogon (IRGC89224, IRGC104425, and IRGC105902). Furthermore, using exon-specific qRT-PCR, OsPLDα1 expression patterns in immature grains indicated significant differences in transcript abundance between the wild species accessions. In comparison to the control, lowest gene expression was observed in IRGC89224 accession (0.20-fold) followed by IRGC105902 (0.26-fold) and IRGC104425 (0.41-fold) accessions. In-silico analysis of the OsPLDα1 promoter revealed that the copy number variations of CGCGBOXAT, GT1CONSENSUS, IBOXCORE, NODCON2GM, OSE2ROOTNODULE, SURECOREATSULTR11, and SORLIP1AT cis-elements play an important role in the transcriptional activities of orthologous genes. Owing to the presence of ARFAT and SEBF elements only in the IRGC89224 accession, which had the lowest gene expression as well, these putative upstream regulatory sequences have been identified as novel cis-elements which may act as repressors in regulating the OsPLDα1 gene expression. The accessions identified with low OsPLDα1 expressions could be further deployed as potential donors of ideal OsPLDα1 allele for transfer of the desired trait into elite rice cultivars.

Combined analysis and miRNA expression profiles of the flowering related genes in common wild rice (oryza rufipogon Griff.)

Common wild rice (Oryza rufipogon Griff.) is the most closely related ancestral species to Asian cultivated rice (Oryza sativa L.). It contains various valuable traits with regard to tolerance to cold, drought and salinity, flowering diversity and many quantitative trait loci with agronomic important traits. Flowering is one of the most important agronomic traits. However, flowering-related transcriptome and how to be regulated by miRNAs have not been estimated in O.rufipogon. To identify how the genes and miRNAs regulating flowering in O.rufipogon. Three O.rufipogon RNA libraries, two vegetative stages (CWRT-V1 and CWRT-V2) and one flowering stage (CWRT-F2) were constructed using leaves tissue and sequenced using Illumina deep sequencing. 27,405, 27,333, 28,979 unique genes were obtained after mapping to the reference genome from CWRT-V1, CWRT-V2 and CWRT-F2, respectively. Then differentially expressed genes (DEGs) were screened and got 1419 unique genes are likely to involve in flower development. Detailed information showed that MADS box and floral meristem identity genes, such as MADS 1, MADS14, Hd1 are involved in common wild rice. Then, combined analysis of miRNA and mRNA expression profiles was performed. Twenty three known miRNA-mRNA pairs and five new candidates were presented an anti-correlationship. Interestingly, 12 miRNAs were negatively correlated with 20 mRNAs encoding flowering-related proteins, indicating that miRNAs regulated target genes to promote flowering in CWRT-F2 group. The results provided here genomic resources for flowering related genes and how these flowering genes were regulated by miRNAs in common wild rice.

Chromosome painting and its applications in cultivated and wild rice

BACKGROUND

The chromosome-specific probe is a fundamental tool of chromosome painting and has been commonly applied in mammalian species. The technology, however, has not been widely applied in plants due to a lack of methodologies for probe development. Identification and labeling of a large number of oligonucleotides (oligos) specific to a single chromosome offers us an opportunity to establish chromosome-specific probes in plants. However, never before has whole chromosome painting been performed in rice.

RESULTS

We developed a pooled chromosome 9-specific probe in rice, which contains 25,000 oligos based on the genome sequence of a japonica rice (Oryza sativa L., AA, 2n = 2× = 24). Chromosome 9 was easily identified in both japonica and indica rice using this chromosome 9-painting probe. The probe was also successfully used to identify and characterize chromosome 9 in additional lines of O. sativa, a translocation line, two new aneuploids associated with chromosome 9 and a wild rice (Oryza eichingeri A. Peter, CC, 2n = 2× = 24).

CONCLUSION

The study reveals that a pool of oligos specific to a chromosome is a useful tool for chromosome painting in rice.

Asian wild rice is a hybrid swarm with extensive gene flow and feralization from domesticated rice

The domestication history of rice remains controversial, with multiple studies reaching different conclusions regarding its origin(s). These studies have generally assumed that populations of living wild rice, O. rufipogon, are descendants of the ancestral population that gave rise to domesticated rice, but relatively little attention has been paid to the origins and history of wild rice itself. Here, we investigate the genetic ancestry of wild rice by analyzing a diverse panel of rice genomes consisting of 203 domesticated and 435 wild rice accessions. We show that most modern wild rice is heavily admixed with domesticated rice through both pollen- and seed-mediated gene flow. In fact, much presumed wild rice may simply represent different stages of feralized domesticated rice. In line with this hypothesis, many presumed wild rice varieties show remnants of the effects of selective sweeps in previously identified domestication genes, as well as evidence of recent selection in flowering genes possibly associated with the feralization process. Furthermore, there is a distinct geographical pattern of gene flow from aus, indica, and japonica varieties into colocated wild rice. We also show that admixture from aus and indica is more recent than gene flow from japonica, possibly consistent with an earlier spread of japonica varieties. We argue that wild rice populations should be considered a hybrid swarm, connected to domesticated rice by continuous and extensive gene flow.

Novel Alleles of Phosphorus-Starvation Tolerance 1 Gene (PSTOL1) from Oryza rufipogon Confers High Phosphorus Uptake Efficiency

Limited phosphorus availability in the soil is one of the major constraints to the growth and productivity of rice across Asian, African and South American countries, where 50% of the rice is grown under rain-fed systems on poor and problematic soils. With an aim to determine novel alleles for enhanced phosphorus uptake efficiency in wild species germplasm of rice Oryza rufipogon, we investigated phosphorus uptake1 (Pup1) locus with 11 previously reported SSR markers and sequence characterized the phosphorus-starvation tolerance 1 (PSTOL1) gene. In the present study, we screened 182 accessions of O. rufipogon along with Vandana as a positive control with SSR markers. From the analysis, it was inferred that all of the O. rufipogon accessions undertaken in this study had an insertion of 90 kb region, including Pup1-K46, a diagnostic marker for PSTOL1, however, it was absent among O. sativa cv. PR114, PR121, and PR122. The complete PSTOL1 gene was also sequenced in 67 representative accessions of O. rufipogon and Vandana as a positive control. From comparative sequence analysis, 53 mutations (52 SNPs and 1 nonsense mutation) were found in the PSTOL1 coding region, of which 28 were missense mutations and 10 corresponded to changes in the amino acid polarity. These 53 mutations correspond to 17 haplotypes, of these 6 were shared and 11 were scored only once. A major shared haplotype was observed among 44 accessions of O. rufipogon along with Vandana and Kasalath. Out of 17 haplotypes, accessions representing 8 haplotypes were grown under the phosphorus-deficient conditions in hydroponics for 60 days. Significant differences were observed in the root length and weight among all the genotypes when grown under phosphorus deficiency conditions as compared to the phosphorus sufficient conditions. The O. rufipogon accession IRGC 106506 from Laos performed significantly better, with 2.5 times higher root weight and phosphorus content as compared to the positive control Vandana. In terms of phosphorus uptake efficiency, the O. rufipogon accessions IRGC 104639, 104712, and 105569 also showed nearly two times higher phosphorus content than Vandana. Thus, these O. rufipogon accessions could be used as the potential donor for improving phosphorus uptake efficiency of elite rice cultivars.

Allele Mining and Selective Patterns of Pi9 Gene in a Set of Rice Landraces from India

Allelic variants of the broad-spectrum blast resistance gene, Pi9 (nucleotide binding site-leucine-rich repeat region) have been analyzed in Indian rice landraces. They were selected from the list of 338 rice landraces phenotyped in the rice blast nursery at central Rainfed Upland Rice Research Station, Hazaribag. Six of them were further selected on the basis of their resistance and susceptible pattern for virulence analysis and selective pattern study of Pi9 gene. The sequence analysis and phylogenetic study illustrated that such sequences are vastly homologous and clustered into two groups. All the blast resistance Pi9 alleles were grouped into one cluster, whereas Pi9 alleles of susceptible landraces formed another cluster even though these landraces have a low level of DNA polymorphisms. A total number of 136 polymorphic sites comprising of transitions, transversions, and insertion and deletions (InDels) were identified in the 2.9 kb sequence of Pi9 alleles. Lower variation in the form of mutations (77) (Transition + Transversion), and InDels (59) were observed in the Pi9 alleles isolated from rice landraces studied. The results showed that the Pi9 alleles of the selected rice landraces were less variable, suggesting that the rice landraces would have been exposed to less number of pathotypes across the country. The positive Tajima's D (0.33580), P > 0.10 (not significant) was observed among the seven rice landraces, which suggests the balancing selection of Pi9 alleles. The value of synonymous substitution (-0.43337) was less than the non-synonymous substitution (0.78808). The greater non-synonymous substitution than the synonymous means that the coding region, mainly the leucine-rich repeat domain was under diversified selection. In this study, the Pi9 gene has been subjected to balancing selection with low nucleotide diversity which is different from the earlier reports, this may be because of the closeness of the rice landraces, cultivated in the same region, and under low pathotype pressure.

QTL Detection and Elite Alleles Mining for Stigma Traits in Oryza sativa by Association Mapping

Stigma traits are very important for hybrid seed production in Oryza sativa, which is a self-pollinated crop; however, the genetic mechanism controlling the traits is poorly understood. In this study, we investigated the phenotypic data of 227 accessions across 2 years and assessed their genotypic variation with 249 simple sequence repeat (SSR) markers. By combining phenotypic and genotypic data, a genome-wide association (GWA) map was generated. Large phenotypic variations in stigma length (STL), stigma brush-shaped part length (SBPL) and stigma non-brush-shaped part length (SNBPL) were found. Significant positive correlations were identified among stigma traits. In total, 2072 alleles were detected among 227 accessions, with an average of 8.3 alleles per SSR locus. GWA mapping detected 6 quantitative trait loci (QTLs) for the STL, 2 QTLs for the SBPL and 7 QTLs for the SNBPL. Eleven, 5, and 12 elite alleles were found for the STL, SBPL, and SNBPL, respectively. Optimal cross designs were predicted for improving the target traits. The detected genetic variation in stigma traits and QTLs provides helpful information for cloning candidate STL genes and breeding rice cultivars with longer STLs in the future.

Rice Varieties in Archaic East Asia: Reduction of Its Diversity from Past to Present Times

The Asian cultivated rice, Oryza sativa, is one of the most important crops feeding more than a third of global population. In spite of the studies for several decades, the origin and domestication history of rice varietal groups, japonica and indica, have not been fully unveiled. Genetic information of ancient rice remains is essential for direct and exclusive insight into the domestication history of rice. We performed ancient DNA analysis of 950- to 2,800-year-old rice remains excavated from Japan and Korea. We found the presence of both japonica- and indica-type varieties in the Yayoi period and the middle ages of Japan and the middle part of Korea Peninsula 2,000 years ago. It is popularly considered that japonica has been exclusively cultivated in northern part of East Asia including Japan and Korea. Our result disclosed unexpectedly wide diversity of rice varieties in archaic East Asia. The present results from ancient rice DNA reveal an exclusive insight for the domestication history of rice which is not provided as far as contemporary rice.

EcoTILLING revealed SNPs in GhSus genes that are associated with fiber- and seed-related traits in upland cotton

Cotton is the most important textile crop in the world due to its cellulose-enriched fibers. Sucrose synthase genes (Sus) play pivotal roles in cotton fiber and seed development. To mine and pyramid more favorable alleles for cotton molecular breeding, single nucleotide polymorphisms (SNPs) of GhSus family genes were investigated across 277 upland cotton accessions by EcoTILLING. As a result, a total of 24 SNPs in the amplified regions of eight GhSus genes were identified. These SNPs were significantly associated with at least one fiber- or seed-related trait measured in Nanjing, Anyang and Kuche in 2007-2009. Four main-effect quantitative trait nucleotides (QTNs) and five epistatic QTNs, with 0.76-3.56% of phenotypic variances explained by each QTN (PVE), were found to be associated with yield-related traits; six epistatic QTNs, with the 0.43-3.48% PVE, were found to be associated with fiber quality-related traits; and one main-effect QTN and one epistatic QTN, with the PVE of 1.96% and 2.53%, were found to be associated with seed oil content and protein content, respectively. Therefore, this study provides new information for molecular breeding in cotton.

Genetic Diversity and Elite Allele Mining for Grain Traits in Rice (Oryza sativa L.) by Association Mapping

Mining elite alleles for grain size and weight is of importance for the improvement of cultivated rice and selection for market demand. In this study, association mapping for grain traits was performed on a selected sample of 628 rice cultivars using 262 SSRs. Grain traits were evaluated by grain length (GL), grain width (GW), grain thickness (GT), grain length to width ratio (GL/GW), and 1000-grain weight (TGW) in 2013 and 2014. Our result showed abundant phenotypic and genetic diversities found in the studied population. In total, 2953 alleles were detected with an average of 11.3 alleles per locus. The population was divided into seven subpopulations and the levels of linkage disequilibrium (LD) ranged from 34 to 84 cM. Genome-wide association mapping detected 10 marker trait association (MTAs) loci for GL, 1MTAs locus for GW, 7 MTAs loci for GT, 3 MTAs loci for GL/GW, and 1 MTAs locus for TGW. Twenty-nine, 2, 10, 5, and 3 elite alleles were found for the GL, GW, GT, GL/GW, and TGW, respectively. Optimal cross designs were predicted for improving the target traits. The accessions containing elite alleles for grain traits mined in this study could be used for breeding rice cultivars and cloning the candidate genes.

A genome-wide association study of a global rice panel reveals resistance in Oryza sativa to root-knot nematodes

The root-knot nematode Meloidogyne graminicola is one of the most serious nematode pests worldwide and represents a major constraint on rice production. While variation in the susceptibility of Asian rice (Oryza sativa) exists, so far no strong and reliable resistance has been reported. Quantitative trait loci for partial resistance have been reported but no underlying genes have been tagged or cloned. Here, 332 accessions of the Rice Diversity Panel 1 were assessed for gall formation, revealing large variation across all subpopulations of rice and higher susceptibility in temperate japonica accessions. Accessions Khao Pahk Maw and LD 24 appeared to be resistant, which was confirmed in large pot experiments where no galls were observed. Detailed observations on these two accessions revealed no nematodes inside the roots 2 days after inoculation and very few females after 17 days (5 in Khao Pahk Maw and <1 in LD 24, in comparison with >100 in the susceptible controls). These two cultivars appear ideal donors for breeding root-knot nematode resistance. A genome-wide association study revealed 11 quantitative trait loci, two of which are close to epistatic loci detected in the Bala x Azucena population. The discussion highlights a small number of candidate genes worth exploring further, in particular many genes with lectin domains and genes on chromosome 11 with homology to the Hordeum Mla locus.

Association Analysis in Rice: From Application to Utilization

Association analysis based on linkage disequilibrium (LD) is an efficient way to dissect complex traits and to identify gene functions in rice. Although association analysis is an effective way to construct fine maps for quantitative traits, there are a few issues which need to be addressed. In this review, we will first summarize type, structure, and LD level of populations used for association analysis of rice, and then discuss the genotyping methods and statistical approaches used for association analysis in rice. Moreover, we will review current shortcomings and benefits of association analysis as well as specific types of future research to overcome these shortcomings. Furthermore, we will analyze the reasons for the underutilization of the results within association analysis in rice breeding.

Review of functional markers for improving cooking, eating, and the nutritional qualities of rice

After yield, quality is one of the most important aspects of rice breeding. Preference for rice quality varies among cultures and regions; therefore, rice breeders have to tailor the quality according to the preferences of local consumers. Rice quality assessment requires routine chemical analysis procedures. The advancement of molecular marker technology has revolutionized the strategy in breeding programs. The availability of rice genome sequences and the use of forward and reverse genetics approaches facilitate gene discovery and the deciphering of gene functions. A well-characterized gene is the basis for the development of functional markers, which play an important role in plant genotyping and, in particular, marker-assisted breeding. In addition, functional markers offer advantages that counteract the limitations of random DNA markers. Some functional markers have been applied in marker-assisted breeding programs and have successfully improved rice quality to meet local consumers' preferences. Although functional markers offer a plethora of advantages over random genetic markers, the development and application of functional markers should be conducted with care. The decreasing cost of sequencing will enable more functional markers for rice quality improvement to be developed, and application of these markers in rice quality breeding programs is highly anticipated.

Nucleotide Sequence Diversity and Linkage Disequilibrium of Four Nuclear Loci in Foxtail Millet (Setaria italica)

Foxtail millet (Setaria italica (L.) Beauv) is one of the earliest domesticated grains, which has been cultivated in northern China by 8,700 years before present (YBP) and across Eurasia by 4,000 YBP. Owing to a small genome and diploid nature, foxtail millet is a tractable model crop for studying functional genomics of millets and bioenergy grasses. In this study, we examined nucleotide sequence diversity, geographic structure, and levels of linkage disequilibrium at four nuclear loci (ADH1, G3PDH, IGS1 and TPI1) in representative samples of 311 landrace accessions across its cultivated range. Higher levels of nucleotide sequence and haplotype diversity were observed in samples from China relative to other sampled regions. Genetic assignment analysis classified the accessions into seven clusters based on nucleotide sequence polymorphisms. Intralocus LD decayed rapidly to half the initial value within ~1.2 kb or less.

Whole-genome analysis revealed the positively selected genes during the differentiation of indica and temperate japonica rice

To investigate the selective pressures acting on the protein-coding genes during the differentiation of indica and japonica, all of the possible orthologous genes between the Nipponbare and 93–11 genomes were identified and compared with each other. Among these genes, 8,530 pairs had identical sequences, and 27,384 pairs shared more than 90% sequence identity. Only 2,678 pairs of genes displaying a Ka/Ks ratio significantly greater than one were revealed, and most of these genes contained only nonsynonymous sites. The genes without synonymous site were further analyzed with the SNP data of 1529 O. sativa and O. rufipogon accessions, and 1068 genes were identified to be under positive selection during the differentiation of indica and temperate japonica. The positively selected genes (PSGs) are unevenly distributed on 12 chromosomes, and the proteins encoded by the PSGs are dominant with binding, transferase and hydrolase activities, and especially enriched in the plant responses to stimuli, biological regulations, and transport processes. Meanwhile, the most PSGs of the known function and/or expression were involved in the regulation of biotic/abiotic stresses. The evidence of pervasive positive selection suggested that many factors drove the differentiation of indica and japonica, which has already started in wild rice but is much lower than in cultivated rice. Lower differentiation and less PSGs revealed between the Or-It and Or-IIIt wild rice groups implied that artificial selection provides greater contribution on the differentiation than natural selection. In addition, the phylogenetic tree constructed with positively selected sites showed that the japonica varieties exhibited more diversity than indica on differentiation, and Or-III of O. rufipogon exhibited more than Or-I.

Extensive sequence variation in rice blast resistance gene Pi54 makes it broad spectrum in nature

Rice blast resistant gene, Pi54 cloned from rice line, Tetep, is effective against diverse isolates of Magnaporthe oryzae. In this study, we prospected the allelic variants of the dominant blast resistance gene from a set of 92 rice lines to determine the nucleotide diversity, pattern of its molecular evolution, phylogenetic relationships and evolutionary dynamics, and to develop allele specific markers. High quality sequences were generated for homologs of Pi54 gene. Using comparative sequence analysis, InDels of variable sizes in all the alleles were observed. Profiling of the selected sites of SNP (Single Nucleotide Polymorphism) and amino acids (N sites ≥ 10) exhibited constant frequency distribution of mutational and substitutional sites between the resistance and susceptible rice lines, respectively. A total of 50 new haplotypes based on the nucleotide polymorphism was also identified. A unique haplotype (H_3) was found to be linked to all the resistant alleles isolated from indica rice lines. Unique leucine zipper and tyrosine sulfation sites were identified in the predicted Pi54 proteins. Selection signals were observed in entire coding sequence of resistance alleles, as compared to LRR domains for susceptible alleles. This is a maiden report of extensive variability of Pi54 alleles in different landraces and cultivated varieties, possibly, attributing broad-spectrum resistance to Magnaporthe oryzae. The sequence variation in two consensus region: 163 and 144 bp were used for the development of allele specific DNA markers. Validated markers can be used for the selection and identification of better allele(s) and their introgression in commercial rice cultivars employing marker assisted selection.

Genetic diversity and association mapping of seed vigor in rice (Oryza sativa L.)

Seed vigor is closely related to direct seeding in rice (Oryza sativa L.). Previous quantitative trait locus (QTL) studies for seed vigor were mainly derived from bi-parental segregating populations and no report from natural populations. In this study, association mapping for seed vigor was performed on a selected sample of 540 rice cultivars (419 from China and 121 from Vietnam). Population structure was estimated on the basis of 262 simple sequence repeat (SSR) markers. Seed vigor was evaluated by root length (RL), shoot length (SL) and shoot dry weight in 2011 and 2012. Abundant phenotypic and genetic diversities were found in the studied population. The population was divided into seven subpopulations, and the levels of linkage disequilibrium (LD) ranged from 10 to 80 cM. We identified 27 marker-trait associations involving 18 SSR markers for three traits. According to phenotypic effects for alleles of the detected QTLs, elite alleles were mined. These elite alleles could be used to design parental combinations and the expected results would be obtained by pyramiding or substituting the elite alleles per QTL (apart from possible epistatic effects). Our results demonstrate that association mapping can complement and enhance previous QTL information for marker-assisted selection and breeding by design.

Archaeological and genetic insights into the origins of domesticated rice

Rice (Oryza sativa) is one of the most important cereal grains in the world today and serves as a staple food source for more than half of the world's population. Research into when, where, and how rice was brought into cultivation and eventually domesticated, along with its development into a staple food source, is thus essential. These questions have been a point of nearly continuous research in both archaeology and genetics, and new information has continually come to light as theory, data acquisition, and analytical techniques have advanced over time. Here, we review the broad history of our scientific understanding of the rice domestication process from both an archaeological and genetic perspective and examine in detail the information that has come to light in both of these fields in the last 10 y. Current findings from genetics and archaeology are consistent with the domestication of O. sativa japonica in the Yangtze River valley of southern China. Interestingly, although it appears rice was cultivated in the area by as early 8000 BP, the key domestication trait of nonshattering was not fixed for another 1,000 y or perhaps longer. Rice was also cultivated in India as early as 5000 BP, but the domesticated indica subspecies currently appears to be a product of the introgression of favorable alleles from japonica. These findings are reshaping our understanding of rice domestication and also have implications for understanding the complex evolutionary process of plant domestication.

The tillering phenotype of the rice plastid terminal oxidase (PTOX) loss-of-function mutant is associated with strigolactone deficiency

The significance of plastid terminal oxidase (PTOX) in phytoene desaturation and chloroplast function has been demonstrated using PTOX-deficient mutants, particularly in Arabidopsis. However, studies on its role in monocots are lacking. Here, we report cloning and characterization of the rice (Oryza sativa) PTOX1 gene. Using Ecotype Targeting Induced Local Lesions IN Genomes (EcoTILLING) and TILLING as forward genetic tools, we identified the causative mutation of an EMS mutant characterized by excessive tillering, semi-dwarfism and leaf variegation that corresponded to the PTOX1 gene. The tillering and semi-dwarf phenotypes of the ptox1 mutant are similar to phenotypes of known strigolactone (SL)-related rice mutants, and both phenotypic traits could be rescued by application of the synthetic SL GR24. The ptox1 mutant accumulated phytoene in white leaf sectors with a corresponding deficiency in β-carotene, consistent with the expected function of PTOX1 in promoting phytoene desaturase activity. There was also no accumulation of the carotenoid-derived SL ent-2'-epi-5-deoxystrigol in root exudates. Elevated concentrations of auxin were detected in the mutant, supporting previous observations that SL interaction with auxin is important in shoot branching control. Our results demonstrate that PTOX1 is required for both carotenoid and SL synthesis resulting in SL-deficient phenotypes in rice.

QTL map meets population genomics: an application to rice

Genes involved in the transition from wild to cultivated crop species should be of great agronomic importance. Population genomic approaches utilizing genome resequencing data have been recently applied for this purpose, although it only reports a large list of candidate genes with no biological information. Here, by resequencing more than 30 genomes altogether of wild rice Oryza rufipogon and cultivated rice O. sativa, we identified a number of regions with clear footprints of selection during the domestication process. We then focused on identifying candidate domestication genes in these regions by utilizing the wealth of QTL information in rice. We were able to identify a number of interesting candidates such as transcription factors that should control key domestication traits such as shattering, awn length, and seed dormancy. Other candidates include those that might have been related to the improvement of grain quality and those that might have been involved in the local adaptation to dry conditions and colder environments. Our study shows that population genomic approaches and QTL mapping information can be used together to identify genes that might be of agronomic importance.

A review of microsatellite markers and their applications in rice breeding programs to improve blast disease resistance

Over the last few decades, the use of molecular markers has played an increasing role in rice breeding and genetics. Of the different types of molecular markers, microsatellites have been utilized most extensively, because they can be readily amplified by PCR and the large amount of allelic variation at each locus. Microsatellites are also known as simple sequence repeats (SSR), and they are typically composed of 1-6 nucleotide repeats. These markers are abundant, distributed throughout the genome and are highly polymorphic compared with other genetic markers, as well as being species-specific and co-dominant. For these reasons, they have become increasingly important genetic markers in rice breeding programs. The evolution of new biotypes of pests and diseases as well as the pressures of climate change pose serious challenges to rice breeders, who would like to increase rice production by introducing resistance to multiple biotic and abiotic stresses. Recent advances in rice genomics have now made it possible to identify and map a number of genes through linkage to existing DNA markers. Among the more noteworthy examples of genes that have been tightly linked to molecular markers in rice are those that confer resistance or tolerance to blast. Therefore, in combination with conventional breeding approaches, marker-assisted selection (MAS) can be used to monitor the presence or lack of these genes in breeding populations. For example, marker-assisted backcross breeding has been used to integrate important genes with significant biological effects into a number of commonly grown rice varieties. The use of cost-effective, finely mapped microsatellite markers and MAS strategies should provide opportunities for breeders to develop high-yield, blast resistance rice cultivars. The aim of this review is to summarize the current knowledge concerning the linkage of microsatellite markers to rice blast resistance genes, as well as to explore the use of MAS in rice breeding programs aimed at improving blast resistance in this species. We also discuss the various advantages, disadvantages and uses of microsatellite markers relative to other molecular marker types.

Molecular diversity in rice blast resistance gene Pi-ta makes it highly effective against dynamic population of Magnaporthe oryzae

Rice blast is one of the important diseases of rice which can be effectively managed by the deployment of resistance genes. Pi-ta is one of the major blast resistant genes effective against pathogen populations in different parts of India. We analysed allelic variants of Pi-ta from 48 rice lines selected after phenotyping of 529 rice landraces across three eco-geographical blast hot spot regions. Besides, Pi-ta orthologue sequences of 220 rice accessions belonging to wild and cultivated species of rice were also included in the study for a better evo-devo perspective of the diversity present in the gene and the selection pressures acting on this locus. We obtained high nucleotide variations (SNPs and insertion-deletions) in the intronic region. We also identified 64 haplotypes based on nucleotide polymorphism in these alleles. Pi-ta orthologues of Indian landraces were scattered in eight major haplotypes indicating its heterogenous nature. We identified a total of 47 different Pi-ta protein variants on the basis of deduced amino acid residues amongst the orthologues. Five unique and novel Pi-ta variants were identified for the first time in rice landraces exhibiting different reaction types against the Magnaporthe oryzae population. A high value of Pi(non/syn) was observed only in the leucine-rich domain of the alleles cloned from Indian landraces, indicating strong selective forces acting on this region. The detailed molecular analysis of the Pi-ta orthologues provides insights to a high degree of inter- and intraspecific relationships amongst the Oryza species. We identified rice landraces possessing the effective alleles of this resistance gene which can be used in future blast resistance breeding programmes.

Origin of Oryza sativa in China inferred by nucleotide polymorphisms of organelle DNA

China is rich of germplasm resources of common wild rice (Oryza rufipogon Griff.) and Asian cultivated rice (O. sativa L.) which consists of two subspecies, indica and japonica. Previous studies have shown that China is one of the domestication centers of O. sativa. However, the geographic origin and the domestication times of O. sativa in China are still under debate. To settle these disputes, six chloroplast loci and four mitochondrial loci were selected to examine the relationships between 50 accessions of Asian cultivated rice and 119 accessions of common wild rice from China based on DNA sequence analysis in the present study. The results indicated that Southern China is the genetic diversity center of O. rufipogon and it might be the primary domestication region of O. sativa. Molecular dating suggested that the two subspecies had diverged 0.1 million years ago, much earlier than the beginning of rice domestication. Genetic differentiations and phylogeography analyses indicated that indica was domesticated from tropical O. rufipogon while japonica was domesticated from O. rufipogon which located in higher latitude. These results provided molecular evidences for the hypotheses of (i) Southern China is the origin center of O. sativa in China and (ii) the two subspecies of O. sativa were domesticated multiple times.

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.

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 (TAC1) as a major quantitative trait locus that controls rice tiller angle. To further clarify the evolutionary characterization of the TAC1 gene, we compared a TAC1-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 TAC1 coding regions of the cultivated rice varieties, whereas one synonymous and one nonsynonymous SNP were detected among the TAC1 coding regions of wild rice accessions. These data indicate that little natural mutation and modification in the TAC1 coding region occurred within the cultivated rice and its progenitor during evolution. Nucleotide diversities in the TAC1 gene regions of O. sativa and O. rufipogon of 0.00116 and 0.00112, respectively, further indicate that TAC1 has been highly conserved during the course of rice domestication. A functional nucleotide polymorphism (FNP) of TAC1 was only found in the japonica rice group. A neutrality test revealed strong selection, especially in the 3'-flanking region of the TAC1 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 TAC1 sequence analysis suggests that the indica and japonica subspecies arose independently during the domestication of wild rice.

SNPs in stress-responsive rice genes: validation, genotyping, functional relevance and population structure

Background

Single nucleotide polymorphism (SNP) validation and large-scale genotyping are required to maximize the use of DNA sequence variation and determine the functional relevance of candidate genes for complex stress tolerance traits through genetic association in rice. We used the bead array platform-based Illumina GoldenGate assay to validate and genotype SNPs in a select set of stress-responsive genes to understand their functional relevance and study the population structure in rice.

Results

Of the 384 putative SNPs assayed, we successfully validated and genotyped 362 (94.3%). Of these 325 (84.6%) showed polymorphism among the 91 rice genotypes examined. Physical distribution, degree of allele sharing, admixtures and introgression, and amino acid replacement of SNPs in 263 abiotic and 62 biotic stress-responsive genes provided clues for identification and targeted mapping of trait-associated genomic regions. We assessed the functional and adaptive significance of validated SNPs in a set of contrasting drought tolerant upland and sensitive lowland rice genotypes by correlating their allelic variation with amino acid sequence alterations in catalytic domains and three-dimensional secondary protein structure encoded by stress-responsive genes. We found a strong genetic association among SNPs in the nine stress-responsive genes with upland and lowland ecological adaptation. Higher nucleotide diversity was observed in indica accessions compared with other rice sub-populations based on different population genetic parameters. The inferred ancestry of 16% among rice genotypes was derived from admixed populations with the maximum between upland aus and wild Oryza species.

Conclusions

SNPs validated in biotic and abiotic stress-responsive rice genes can be used in association analyses to identify candidate genes and develop functional markers for stress tolerance in rice.

A novel blast resistance gene, Pi54rh cloned from wild species of rice, Oryza rhizomatis confers broad spectrum resistance to Magnaporthe oryzae

The dominant rice blast resistance gene, Pi54 confers resistance to Magnaporthe oryzae in different parts of India. In our effort to identify more effective forms of this gene, we isolated an orthologue of Pi54 named as Pi54rh from the blast-resistant wild species of rice, Oryza rhizomatis, using allele mining approach and validated by complementation. The Pi54rh belongs to CC-NBS-LRR family of disease resistance genes with a unique Zinc finger (C(3)H type) domain. The 1,447 bp Pi54rh transcript comprises of 101 bp 5'-UTR, 1,083 bp coding region and 263 bp 3'-UTR, driven by pathogen inducible promoter. We showed the extracellular localization of Pi54rh protein and the presence of glycosylation, myristoylation and phosphorylation sites which implicates its role in signal transduction process. This is in contrast to other blast resistance genes that are predicted to be intracellular NBS-LRR-type resistance proteins. The Pi54rh was found to express constitutively at basal level in the leaves, but upregulates 3.8-fold at 96 h post-inoculation with the pathogen. Functional validation of cloned Pi54rh gene using complementation test showed high degree of resistance to seven isolates of M. oryzae collected from different geographical locations of India. In this study, for the first time, we demonstrated that a rice blast resistance gene Pi54rh cloned from wild species of rice provides broad spectrum resistance to M. oryzae hence can be used in rice improvement breeding programme.

Population structure and genetic diversity in a rice core collection (Oryza sativa L.) investigated with SSR markers

The assessment of genetic diversity and population structure of a core collection would benefit to make use of these germplasm as well as applying them in association mapping. The objective of this study were to (1) examine the population structure of a rice core collection; (2) investigate the genetic diversity within and among subgroups of the rice core collection; (3) identify the extent of linkage disequilibrium (LD) of the rice core collection. A rice core collection consisting of 150 varieties which was established from 2260 varieties of Ting's collection of rice germplasm were genotyped with 274 SSR markers and used in this study. Two distinct subgroups (i.e. SG 1 and SG 2) were detected within the entire population by different statistical methods, which is in accordance with the differentiation of indica and japonica rice. MCLUST analysis might be an alternative method to STRUCTURE for population structure analysis. A percentage of 26% of the total markers could detect the population structure as the whole SSR marker set did with similar precision. Gene diversity and MRD between the two subspecies varied considerably across the genome, which might be used to identify candidate genes for the traits under domestication and artificial selection of indica and japonica rice. The percentage of SSR loci pairs in significant (P<0.05) LD is 46.8% in the entire population and the ratio of linked to unlinked loci pairs in LD is 1.06. Across the entire population as well as the subgroups and sub-subgroups, LD decays with genetic distance, indicating that linkage is one main cause of LD. The results of this study would provide valuable information for association mapping using the rice core collection in future.

Genetic architecture of aluminum tolerance in rice (Oryza sativa) determined through genome-wide association analysis and QTL mapping

Aluminum (Al) toxicity is a primary limitation to crop productivity on acid soils, and rice has been demonstrated to be significantly more Al tolerant than other cereal crops. However, the mechanisms of rice Al tolerance are largely unknown, and no genes underlying natural variation have been reported. We screened 383 diverse rice accessions, conducted a genome-wide association (GWA) study, and conducted QTL mapping in two bi-parental populations using three estimates of Al tolerance based on root growth. Subpopulation structure explained 57% of the phenotypic variation, and the mean Al tolerance in Japonica was twice that of Indica. Forty-eight regions associated with Al tolerance were identified by GWA analysis, most of which were subpopulation-specific. Four of these regions co-localized with a priori candidate genes, and two highly significant regions co-localized with previously identified QTLs. Three regions corresponding to induced Al-sensitive rice mutants (ART1, STAR2, Nrat1) were identified through bi-parental QTL mapping or GWA to be involved in natural variation for Al tolerance. Haplotype analysis around the Nrat1 gene identified susceptible and tolerant haplotypes explaining 40% of the Al tolerance variation within the aus subpopulation, and sequence analysis of Nrat1 identified a trio of non-synonymous mutations predictive of Al sensitivity in our diversity panel. GWA analysis discovered more phenotype-genotype associations and provided higher resolution, but QTL mapping identified critical rare and/or subpopulation-specific alleles not detected by GWA analysis. Mapping using Indica/Japonica populations identified QTLs associated with transgressive variation where alleles from a susceptible aus or indica parent enhanced Al tolerance in a tolerant Japonica background. This work supports the hypothesis that selectively introgressing alleles across subpopulations is an efficient approach for trait enhancement in plant breeding programs and demonstrates the fundamental importance of subpopulation in interpreting and manipulating the genetics of complex traits in rice.

Mapping QTLs for improving grain yield using the USDA rice mini-core collection

Yield is the most important and complex trait for genetic improvement in crops, and marker-assisted selection enhances the improvement efficiency. The USDA rice mini-core collection derived from over 18,000 accessions of global origins is an ideal panel for association mapping. We phenotyped 203 O. sativa accessions for 14 agronomic traits and identified 5 that were highly and significantly correlated with grain yield per plant: plant height, plant weight, tillers, panicle length, and kernels/branch. Genotyping with 155 genome-wide molecular markers demonstrated 5 main cluster groups. Linkage disequilibrium (LD) decayed at least 20 cM and marker pairs with significant LD ranged from 4.64 to 6.06% in four main groups. Model comparisons revealed that different dimensions of principal component analysis affected yield and its correlated traits for mapping accuracy, and kinship did not improve the mapping in this collection. Thirty marker-trait associations were highly significant, 4 for yield, 3 for plant height, 6 for plant weight, 9 for tillers, 5 for panicle length and 3 for kernels/branch. Twenty-one markers contributed to the 30 associations, because 8 markers were co-associated with 2 or more traits. Allelic analysis of OSR13, RM471 and RM7003 for their co-associations with yield traits demonstrated that allele 126 bp of RM471 and 108 bp of RM7003 should receive greater attention, because they had the greatest positive effect on yield traits. Tagging the QTLs responsible for multiple yield traits may simultaneously help dissect the complex yield traits and elevate the efficiency to improve grain yield using marker-assisted selection in rice.

Molecular evidence for a single evolutionary origin of domesticated rice

Asian rice, Oryza sativa, is one of world's oldest and most important crop species. Rice is believed to have been domesticated ∼9,000 y ago, although debate on its origin remains contentious. A single-origin model suggests that two main subspecies of Asian rice, indica and japonica, were domesticated from the wild rice O. rufipogon. In contrast, the multiple independent domestication model proposes that these two major rice types were domesticated separately and in different parts of the species range of wild rice. This latter view has gained much support from the observation of strong genetic differentiation between indica and japonica as well as several phylogenetic studies of rice domestication. We reexamine the evolutionary history of domesticated rice by resequencing 630 gene fragments on chromosomes 8, 10, and 12 from a diverse set of wild and domesticated rice accessions. Using patterns of SNPs, we identify 20 putative selective sweeps on these chromosomes in cultivated rice. Demographic modeling based on these SNP data and a diffusion-based approach provide the strongest support for a single domestication origin of rice. Bayesian phylogenetic analyses implementing the multispecies coalescent and using previously published phylogenetic sequence datasets also point to a single origin of Asian domesticated rice. Finally, we date the origin of domestication at ∼8,200-13,500 y ago, depending on the molecular clock estimate that is used, which is consistent with known archaeological data that suggests rice was first cultivated at around this time in the Yangtze Valley of China.

Association mapping of grain color, phenolic content, flavonoid content and antioxidant capacity in dehulled rice

Phytochemicals such as phenolics and flavonoids in rice grain are antioxidants that are associated with reduced risk of developing chronic diseases including cardiovascular disease, type-2 diabetes and some cancers. Understanding the genetic basis of these traits is necessary for the improvement of nutritional quality by breeding. Association mapping based on linkage disequilibrium has emerged as a powerful strategy for identifying genes or quantitative trait loci (QTL) underlying complex traits in plants. In this study, genome-wide association mapping using models controlling both population structure (Q) and relative kinship (K) were performed to identify the marker loci/QTLs underlying the naturally occurring variations of grain color and nutritional quality traits in 416 rice germplasm accessions including red and black rice. A total of 41 marker loci were identified for all the traits, and it was confirmed that Ra (i.e., Prp-b for purple pericarp) and Rc (brown pericarp and seed coat) genes were main-effect loci for rice grain color and nutritional quality traits. RM228, RM339, fgr (fragrance gene) and RM316 were important markers associated with most of the traits. Association mapping for the traits of the 361 white or non-pigmented rice accessions (i.e., excluding the red and black rice) revealed a total of 11 markers for four color parameters, and one marker (RM346) for phenolic content. Among them, Wx gene locus was identified for the color parameters of lightness (L*), redness (a*) and hue angle (H (o)). Our study suggested that the markers identified in this study can feasibly be used to improve nutritional quality or health benefit properties of rice by marker-assisted selection if the co-segregations of the marker-trait associations are validated in segregating populations.

Ecological divergence in the presence of gene flow in two closely related Oryza species (Oryza rufipogon and O. nivara)

Ecological divergence plays a prominent role in the process of speciation, but how divergence occurs in the face of gene flow is still less clear, and remains controversial among evolutionists. Here we investigated the nucleotide diversity, divergence and gene flow between Oryza nivara and O. rufipogon using sequences of seven chloroplast and nuclear loci. By analysing samples from 26 wild populations across the geographic ranges of the two species, we showed that both species were highly structured and O. rufipogon maintained a higher level of species-wide diversity than O. nivara. Notably, phylogenetic, amova and FST analyses were unable to detect significant nucleotide differentiation between the two species. We estimated that the two species began to diverge at c. 0.16 million years ago. Our coalescent-based simulations strongly rejected the simple isolation model of zero migration between species, but rather provided unambiguous evidence of bidirectional gene flow between species, particularly from O. rufipogon to O. nivara. Our simulations also indicated that gene flow was recurrent during the divergence process rather than arising from secondary contact after allopatric divergence. In conjunction with different morphological and life-history traits and habitat preference in the two species, this study supports the hypothesis that these Oryza species are better treated as ecotypes that diverged quite recently and are still under the process of divergence. Importantly, we demonstrate the ecological divergence between O. rufipogon and O. nivara in the presence of significant gene flow, implying that natural selection plays a primary role in driving the divergence of the two Oryza species.

Towards population genomics of effector-effector target interactions

Pathogen-plant host coevolutionary interactions exert strong natural selection on both organisms, specifically on the genes coding for effectors (pathogens), as well as on those coding for effector targets and R proteins (plant hosts). Natural selection leaves behind DNA sequence signatures on such genes and on linked genomic regions. These signatures can readily be detected by studying the patterns of intraspecies polymorphisms and interspecies divergence of the DNA sequences. Recent developments in DNA sequencing technology have made whole-genome studies on patterns of DNA polymorphisms : divergence possible. This type of analysis, called 'population genomics', appears to be powerful enough to identify novel effector-effector target genes. Here, we provide an overview of the statistical tools used for population genomics and their applications. This is followed by a brief review of evolutionary studies on plant genes involved in host-pathogen interactions. Finally we provide an example from our study on Magnaporthe oryzae.

Genetic diversity and population structure of a diverse set of rice germplasm for association mapping

Germplasm diversity is the mainstay for crop improvement and genetic dissection of complex traits. Understanding genetic diversity, population structure, and the level and distribution of linkage disequilibrium (LD) in target populations is of great importance and a prerequisite for association mapping. In this study, 100 genome-wide simple sequence repeat (SSR) markers were used to assess genetic diversity, population structure, and LD of 416 rice accessions including landraces, cultivars and breeding lines collected mostly in China. A model-based population structure analysis divided the rice materials into seven subpopulations. 63% of the SSR pairs in these accessions were in LD, which was mostly due to an overall population structure, since the number of locus pairs in LD was reduced sharply within each subpopulation, with the SSR pairs in LD ranging from 5.9 to 22.9%. Among those SSR pairs showing significant LD, the intrachromosomal LD had an average of 25-50 cM in different subpopulations. Analysis of the phenotypic diversity of 25 traits showed that the population structure accounted for an average of 22.4% of phenotypic variation. An example association mapping for starch quality traits using both the candidate gene mapping and genome-wide mapping strategies based on the estimated population structure was conducted. Candidate gene mapping confirmed that the Wx and starch synthase IIa (SSIIa) genes could be identified as strongly associated with apparent amylose content (AAC) and pasting temperature (PT), respectively. More importantly, we revealed that the Wx gene was also strongly associated with PT. In addition to the major genes, we found five and seven SSRs were associated with AAC and PT, respectively, some of which have not been detected in previous linkage mapping studies. The results suggested that the population may be useful for the genome-wide marker-trait association mapping. This new association population has the potential to identify quantitative trait loci (QTL) with small effects, which will aid in dissecting complex traits and in exploiting the rich diversity present in rice germplasm.

Application of association mapping to understanding the genetic diversity of plant germplasm resources

Compared to the conventional linkage mapping, linkage disequilibrium (LD)-mapping, using the nonrandom associations of loci in haplotypes, is a powerful high-resolution mapping tool for complex quantitative traits. The recent advances in the development of unbiased association mapping approaches for plant population with their successful applications in dissecting a number of simple to complex traits in many crop species demonstrate a flourish of the approach as a “powerful gene tagging” tool for crops in the plant genomics era of 21st century. The goal of this review is to provide nonexpert readers of crop breeding community with (1) the basic concept, merits, and simple description of existing methodologies for an association mapping with the recent improvements for plant populations, and (2) the details of some of pioneer and recent studies on association mapping in various crop species to demonstrate the feasibility, success, problems, and future perspectives of the efforts in plants. This should be helpful for interested readers of international plant research community as a guideline for the basic understanding, choosing the appropriate methods, and its application.

Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome

Background

The domestication of Asian rice (Oryza sativa) was a complex process punctuated by episodes of introgressive hybridization among and between subpopulations. Deep genetic divergence between the two main varietal groups (Indica and Japonica) suggests domestication from at least two distinct wild populations. However, genetic uniformity surrounding key domestication genes across divergent subpopulations suggests cultural exchange of genetic material among ancient farmers.

Methodology/Principal Findings

In this study, we utilize a novel 1,536 SNP panel genotyped across 395 diverse accessions of O. sativa to study genome-wide patterns of polymorphism, to characterize population structure, and to infer the introgression history of domesticated Asian rice. Our population structure analyses support the existence of five major subpopulations (indica, aus, tropical japonica, temperate japonica and GroupV) consistent with previous analyses. Our introgression analysis shows that most accessions exhibit some degree of admixture, with many individuals within a population sharing the same introgressed segment due to artificial selection. Admixture mapping and association analysis of amylose content and grain length illustrate the potential for dissecting the genetic basis of complex traits in domesticated plant populations.

Conclusions/Significance

Genes in these regions control a myriad of traits including plant stature, blast resistance, and amylose content. These analyses highlight the power of population genomics in agricultural systems to identify functionally important regions of the genome and to decipher the role of human-directed breeding in refashioning the genomes of a domesticated species.

Association analysis of the amino acid contents in rice

The main objective of the present study was to identify simple sequence repeat (SSR) markers associated with the amino acid content of rice (Oryza sativa L.). SSR markers were selected by prescreening for the relationship to amino acid content. Eighty-four rice landrace accessions from Korea were evaluated for 16 kinds of amino acids in brown rice and genotyped with 25 SSR markers. Analysis of population structure revealed four subgroups in the population. Linkage disequilibrium (LD) patterns and distributions are of fundamental importance for genome-wide mapping associations. The mean r(2) value for all intrachromosomal loci pairs was 0.033. LD between linked markers decreased with distance. Marker-trait associations were investigated using the unified mixed-model approach, considering both population structure (Q) and kinship (K). A total of 42 marker-trait associations with amino acids (P < 0.05) were identified using 15 different SSR markers covering three chromosomes and explaining more than 40% of the total variation. These results suggest that association analysis in rice is a viable alternative to quantitative trait loci mapping and should help rice breeders develop strategies for improving rice quality.

Association genetics reveals three novel avirulence genes from the rice blast fungal pathogen Magnaporthe oryzae

To subvert rice (Oryza sativa) host defenses, the devastating ascomycete fungus pathogen Magnaporthe oryzae produces a battery of effector molecules, including some with avirulence (AVR) activity, which are recognized by host resistance (R) proteins resulting in rapid and effective activation of innate immunity. To isolate novel avirulence genes from M. oryzae, we examined DNA polymorphisms of secreted protein genes predicted from the genome sequence of isolate 70-15 and looked for an association with AVR activity. This large-scale study found significantly more presence/absence polymorphisms than nucleotide polymorphisms among 1032 putative secreted protein genes. Nucleotide diversity of M. oryzae among 46 isolates of a worldwide collection was extremely low (theta=8.2x10(-5)), suggestive of recent pathogen dispersal. However, no association between DNA polymorphism and AVR was identified. Therefore, we used genome resequencing of Ina168, an M. oryzae isolate that contains nine AVR genes. Remarkably, a total of 1.68 Mb regions, comprising 316 candidate effector genes, were present in Ina168 but absent in the assembled sequence of isolate 70-15. Association analyses of these 316 genes revealed three novel AVR genes, AVR-Pia, AVR-Pii, and AVR-Pik/km/kp, corresponding to five previously known AVR genes, whose products are recognized inside rice cells possessing the cognate R genes. AVR-Pia and AVR-Pii have evolved by gene gain/loss processes, whereas AVR-Pik/km/kp has evolved by nucleotide substitutions and gene gain/loss.

Contrasting population genetic structure and gene flow between Oryza rufipogon and Oryza nivara

The cross compatible wild relatives of crops have furnished valuable genes for crop improvement. Understanding the genetics of these wild species may enhance their further use in breeding. In this study, sequence variation of the nuclear Lhs1 gene was used to investigate the population genetic structure and gene flow of Oryza rufipogon and O. nivara, two wild species most closely related to O. sativa. The two species diverge markedly in life history and mating system, with O. rufipogon being perennial and outcrossing and O. nivara being annual and predominantly inbreeding. Based on sequence data from 105 plants representing 11 wild populations covering the entire geographic range of these wild species, we detected significantly higher nucleotide variation in O. rufipogon than in O. nivara at both the population and species levels. At the population level the diversity in O. rufipogon (Hd = 0.712; theta (sil) = 0.0017) is 2-3 folds higher than that in O. nivara (Hd = 0.306; theta (sil) = 0.0005). AMOVA partitioning indicated that genetic differentiation among O. nivara populations (78.2%) was much higher than that among O. rufipogon populations (52.3%). The different level of genetic diversity and contrasting population genetic structure between O. rufipogon and O. nivara might be explained by their distinct life histories and mating systems. Our simulation using IM models demonstrated significant gene flow from O. nivara to O. rufipogon, indicating a directional introgression from the annual and selfing species into the perennial and outcrossing species. The ongoing introgression has played an important role in shaping current patterns of genetic diversity of these two wild species.

The evolution of sex-independent transmission ratio distortion involving multiple allelic interactions at a single locus in rice

Transmission ratio distortion (TRD) is frequently observed in inter- and intraspecific hybrids of plants, leading to a violation of Mendelian inheritance. Sex-independent TRD (siTRD) was detected in a hybrid between Asian cultivated rice and its wild ancestor. Here we examined how siTRD caused by an allelic interaction at a specific locus arose in Asian rice species. The siTRD is controlled by the S6 locus via a mechanism in which the S6 allele acts as a gamete eliminator, and both the male and female gametes possessing the opposite allele (S6a) are aborted only in heterozygotes (S6/S6a). Fine mapping revealed that the S6 locus is located near the centromere of chromosome 6. Testcross experiments using near-isogenic lines (NILs) carrying either the S6 or S6a alleles revealed that Asian rice strains frequently harbor an additional allele (S6n) the presence of which, in heterozygotic states (S6/S6n and S6a/S6n), does not result in siTRD. A prominent reduction in the nucleotide diversity of S6 or S6a carriers relative to that of S6n carriers was detected in the chromosomal region. These results suggest that the two incompatible alleles (S6 and S6a) arose independently from S6n and established genetically discontinuous relationships between limited constituents of the Asian rice population.

Molecular evolution of the Pi-ta gene resistant to rice blast in wild rice (Oryza rufipogon)

Rice blast disease resistance to the fungal pathogen Magnaporthe grisea is triggered by a physical interaction between the protein products of the host R (resistance) gene, Pi-ta, and the pathogen Avr (avirulence) gene, AVR-pita. The genotype variation and resistant/susceptible phenotype at the Pi-ta locus of wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. sativa), was surveyed in 36 locations worldwide to study the molecular evolution and functional adaptation of the Pi-ta gene. The low nucleotide polymorphism of the Pi-ta gene of O. rufipogon was similar to that of O. sativa, but greatly differed from what has been reported for other O. rufipogon genes. The haplotypes can be subdivided into two divergent haplogroups named H1 and H2. H1 is derived from H2, with nearly no variation and at a low frequency. H2 is common and is the ancestral form. The leucine-rich repeat (LRR) domain has a high pi(non)/pi(syn) ratio, and the low polymorphism of the Pi-ta gene might have primarily been caused by recurrent selective sweep and constraint by other putative physiological functions. Meanwhile, we provide data to show that the amino acid Ala-918 of H1 in the LRR domain has a close relationship with the resistant phenotype. H1 might have recently arisen during rice domestication and may be associated with the scenario of a blast pathogen-host shift from Italian millet to rice.