Genome-wide haplotype changes produced by artificial selection during modern rice breeding in Japan

Genome-wide association analysis identifies natural allelic variants associated with panicle architecture variation in African rice, Oryza glaberrima Steud

African rice (Oryza glaberrima Steud), a short-day cereal crop closely related to Asian rice (Oryza sativa L.), has been cultivated in Sub-Saharan Africa for ∼ 3,000 years. Although less cultivated globally, it is a valuable genetic resource in creating high-yielding cultivars that are better adapted to diverse biotic and abiotic stresses. While inflorescence architecture, a key trait for rice grain yield improvement, has been extensively studied in Asian rice, the morphological and genetic determinants of this complex trait are less understood in African rice. In this study, using a previously developed association panel of 162 O. glaberrima accessions and new SNP variants characterized through mapping to a new version of the O. glaberrima reference genome, we conducted a genome-wide association study of four major morphological panicle traits. We have found a total of 41 stable genomic regions that are significantly associated with these traits, of which 13 co-localized with previously identified QTLs in O. sativa populations and 28 were unique for this association panel. Additionally, we found a genomic region of interest on chromosome 3 that was associated with the number of spikelets and primary and secondary branches. Within this region was localized the O. sativa ortholog of the PHYTOCHROME B gene (Oglab_006903/OgPHYB). Haplotype analysis revealed the occurrence of natural sequence variants at the OgPHYB locus associated with panicle architecture variation through modulation of the flowering time phenotype, whereas no equivalent alleles were found in O. sativa. The identification in this study of genomic regions specific to O. glaberrima indicates panicle-related intra-specific genetic variation in this species, increasing our understanding of the underlying molecular processes governing panicle architecture. Identified candidate genes and major haplotypes may facilitate the breeding of new African rice cultivars with preferred panicle traits.

Identification of Novel Quantitative Trait Loci for Culm Thickness of Rice Derived from Strong-Culm Landrace in Japan, Omachi

Increasing the lodging resistance of rice through genetic improvement has been an important target in breeding. To further enhance the lodging resistance of high-yielding rice varieties amidst climate change, it is necessary to not only shorten culms but strengthen them as well. A landrace rice variety, Omachi, which was established more than 100 years ago, has the largest culm diameter and bending moment at breaking in the basal internodes among 135 temperate japonica accessions. Using unused alleles in such a landrace is an effective way to strengthen the culm. In this study, we performed quantitative trait locus (QTL) analysis to identify the genetic factors of culm strength of Omachi using recombinant inbred lines (RILs) derived from a cross between Omachi and Koshihikari, a standard variety in Japan. We identified three QTLs for the culm diameter of the 5th internode on chromosomes 3 (qCD3) and 7 (qCD7-1, qCD7-2). Among them, qCD7-2 was verified by QTL analysis using the F₂ population derived from a cross between one of the RILs and Koshihikari. RNA-seq analysis of shoot apex raised 10 candidate genes underlying the region of qCD7-2. The increase in culm strength by accumulating Omachi alleles of qCD3, qCD7-1 and qCD7-2 was 25.0% in 2020. These QTLs for culm diameter pleiotropically increased spikelet number per panicle but did not affect days to heading or culm length. These results suggest that the Omachi alleles of qCD3, qCD7-1 and qCD7-2 are useful for breeding to increase lodging resistance and yield.

CONSTANS Polymorphism Modulates Flowering Time and Maturity in Soybean

CONSTANS (CO) plays a critical role in the photoperiodic flowering pathway. However, the function of soybean CO orthologs and the molecular mechanisms in regulating flowering remain largely unknown. This study characterized the natural variations in CO family genes and their association with flowering time and maturity in soybeans. A total of 21 soybean CO family genes (GmCOLs) were cloned and sequenced in 128 varieties covering 14 known maturity groups (MG 0000-MG X from earliest to latest maturity). Regarding the whole genomic region involving these genes, GmCOL1, GmCOL3, GmCOL8, GmCOL9, GmCOL10, and GmCOL13 were conserved, and the remaining 15 genes showed genetic variation that was brought about by mutation, namely, all single-nucleotide polymorphisms (SNPs) and insertions-deletions (InDels). In addition, a few genes showed some strong linkage disequilibrium. Point mutations were found in 15 GmCOL genes, which can lead to changes in the potential protein structure. Early flowering and maturation were related to eight genes (GmCOL1/3/4/8/13/15/16/19). For flowering and maturation, 11 genes (GmCOL2/5/6/14/20/22/23/24/25/26/28) expressed divergent physiognomy. Haplotype analysis indicated that the haplotypes of GmCOL5-Hap2, GmCOL13-Hap2/3, and GmCOL28-Hap2 were associated with flowering dates and soybean maturity. This study helps address the role of GmCOL family genes in adapting to diverse environments, particularly when it is necessary to regulate soybean flowering dates and maturity.

SSR marker-based genetic diversity analysis and SNP haplotyping of genes associating abiotic and biotic stress tolerance, rice growth and development and yield across 93 rice landraces

BACKGROUND

As rice is the staple food for more than half of the world's population, enhancing grain yield irrespective of the variable climatic conditions is indispensable. Many traditionally cultivated rice landraces are well adapted to severe environmental conditions and have high genetic diversity that could play an important role in crop improvement.

METHODS AND RESULTS

The present study revealed a high level of genetic diversity among the unexploited rice landraces cultivated by the farmers of Kerala. Twelve polymorphic markers detected a total of seventy- seven alleles with an average of 6.416 alleles per locus. Polymorphic Information Content (PIC) value ranged from 0.459 to 0.809, and to differentiate the rice genotypes, RM 242 was found to be the most appropriate marker with a high value of 0.809. The current study indicated that the rice landraces are highly diverse with higher values of the adequate number of alleles, PIC, and Shannon information index. Utilizing these informative SSR markers for future molecular characterization and population genetic studies in rice landraces are advisable. Haplotypes are sets of genomic regions within a chromosome inherited together, and haplotype-based breeding is a promising strategy for designing next-generation rice varieties. Here, haplotype analysis explored 270 haplotype blocks and 775 haplotypes from all the chromosomes of landraces under study. The number of SNPs in each haplotype block ranged from two to 28. Haplotypes of genes related to biotic and abiotic stress tolerance, yield-enhancing, and growth and development in rice landraces were also elucidated in the current study.

CONCLUSIONS

The present investigation revealed the genetic diversity of rice landraces and the haplotype analysis will open the way for genome-wide association studies, QTL identification, and marker-assisted selection in the unexplored rice landraces collected from Kerala.

Translation of continuous artificial selection on phenotype into genotype during rice breeding programs

Understanding genetic diversity among local populations is a primary goal of modern crop breeding programs. Here, we demonstrated the genetic relationships of rice varieties in Hokkaido, Japan, one of the northern limits of rice cultivation around the world. Furthermore, artificial selection during rice breeding programs has been characterized using genome sequences. We utilized 8,565 single nucleotide polymorphisms and insertion/deletion markers distributed across the genome in genotype-by-sequencing for genetic diversity analyses. Phylogenetics, genetic population structure, and principal component analysis showed that a total of 110 varieties were classified into four distinct clusters according to different populations geographically and historically. Furthermore, the genome sequences of 19 rice varieties along with historic representations in Hokkaido, nucleotide diversity and F_ST values in each cluster revealed that artificial selection of elite phenotypes focused on chromosomal regions. These results clearly demonstrated the history of the selections on agronomic traits as genome sequences among current rice varieties from Hokkaido.

Genome-Wide Association Study Reveals Novel Marker-Trait Associations (MTAs) Governing the Localization of Fe and Zn in the Rice Grain

Micronutrient malnutrition due to Fe and Zn, affects around two billion people globally particularly in the developing countries. More than 90% of the Asian population is dependent on rice-based diets, which is low in these micronutrients. In the present study, a set of 192 Indian rice germplasm accessions, grown at two locations, were evaluated for Fe and Zn in brown rice (BR) and milled rice (MR). A significant variation was observed in the rice germplasm for these micronutrients. The grain Fe concentration was in the range of 6.2–23.1 ppm in BR and 0.8–12.3 ppm in MR, while grain Zn concentration was found to be in the range of 11.0–47.0 ppm and 8.2–40.8 ppm in the BR and MR, respectively. Grain Fe exhibited maximum loss upon milling with a mean retention of 24.9% in MR, while Zn showed a greater mean retention of 74.2% in MR. A genome-wide association study (GWAS) was carried out implementing the FarmCPU model to control the population structure and kinship, and resulted in the identification of 29 marker-trait associations (MTAs) with significant associations for traits viz. FeBR (6 MTAs), FeMR (7 MTAs), ZnBR (11 MTAs), and ZnMR (5 MTAs), which could explain the phenotypic variance from 2.1 to as high as 53.3%. The MTAs governing the correlated traits showed co-localization, signifying the possibility of their simultaneous improvement. The robust MTAs identified in the study could be valuable resource for enhancing Fe and Zn concentration in the rice grain and addressing the problem of Fe and Zn malnutrition among rice consumers.

Marker-assisted selection in rice breeding programs in Hokkaido

Rice breeding programs in Hokkaido over the past 100 years have dramatically increased productivity and improved the eating quality of rice. Commercial varieties with high yield and good eating quality, such as Kirara 397, Hoshinoyume, and Nanatsuboshi, have been continuously registered since 1990. Furthermore, varieties with better eating quality using Wx1-1, which reduces amylose content to improve the taste of sticky rice, such as Oborozuki and Yumepirika, were registered in 2006 and 2008, respectively. However, to the best of our knowledge the genomic changes associated with these improvements have not been determined. Better understanding of the relationships between DNA sequences and agricultural traits could facilitate rice breeding programs in Hokkaido. Marker-assisted selection (MAS), which can select the plants with chromosomal regions tagged with DNA markers for desirable traits, is an advanced technology to manage genetic improvements. Here, we summarize the current states of MAS in rice breeding programs in Hokkaido before huge data sets of genome sequences using next-generation sequencing technology come into practical use in rice breeding programs.

High Levels of Variation Within Gene Sequences of Olea europaea L

Gene sequence variation in cultivated olive (Olea europaea L. subsp. europaea var. europaea), the most important oil tree crop of the Mediterranean basin, has been poorly evaluated up to now. A deep sequence analysis of fragments of four genes, OeACP1, OeACP2, OeLUS and OeSUT1, in 90 cultivars, revealed a wide range of polymorphisms along all recognized allele forms and unexpected allele frequencies and genotype combinations. High linkage values among most polymorphisms were recorded within each gene fragment. The great sequence variability corresponded to a low number of alleles and, surprisingly, to a small fraction of genotype combinations. The distribution, frequency, and combination of the different alleles at each locus is possibly due to natural and human pressures, such as selection, ancestrality, or fitness. Phylogenetic analyses of allele sequences showed distant and complex patterns of relationships among cultivated olives, intermixed with other related forms, highlighting an evolutionary connection between olive cultivars and the O. europaea subspecies cuspidata and cerasiformis. This study demonstrates how a detailed and complete sequence analysis of a few gene portions and a thorough genotyping on a representative set of cultivars can clarify important issues related to sequence polymorphisms, reconstructing the phylogeny of alleles, as well as the genotype combinations. The identification of regions representing blocks of recombination could reveal polymorphisms that represent putatively functional markers. Indeed, specific mutations found on the analyzed OeACP1 and OeACP2 fragments seem to be correlated to the fruit weight.

Endosperm enzyme activity is responsible for texture and eating quality of cooked rice grains in Japanese cultivars

Eating quality of cooked rice grains is an important determinant of its market price and consumer acceptance. To comprehensively assess the variation of eating-quality traits in 152 Japanese rice cultivars, we evaluated activities of eight endosperm enzymes related to degradation of starch and cell-wall polysaccharides. Endosperm enzyme activities showed a wide range of variations and were lower in recently developed cultivars than in landraces and old improved cultivars. Activities of most endosperm enzymes correlated significantly with the eating-quality score and surface texture of cooked rice grains. Principal component analysis revealed that rice cultivars with high eating-quality scores had high stickiness of the grain surface and low levels of endosperm enzyme activities. These results suggest that endosperm enzyme activities control texture and eating quality of cooked rice grains in Japanese rice cultivars.

A follow-up study for biomass yield QTLs in rice

The biomass yield (plant weight) of rice fluctuates from year to year. In a previous study, we demonstrated that six quantitative trait loci (QTLs) contribute to the variation in the plant weight of recombinant inbred lines (RILs) of high-yielding Japanese rice cultivars. However, it remains unclear whether the effects of those QTLs are stable over multiple years. Therefore, we evaluated the effect of the alleles on the plant weight of RILs over multiple years, including a change of fertilization level (i.e., in different environments). Even though the biomass yields of all RILs fluctuated among environments, RILs that were selected on the basis of the genotypes of the detected QTLs had a stable rank order of plant weight that corresponded to their genotypes. This multiple-environment experiment reveals the highly significant contribution of both genotypic and environmental variances to the observed variance in plant weight. A marginally significant QTL–environment interaction was detected at only one of the six QTLs, with a subtle contribution. These results support the idea that the biomass yield of rice can be improved through QTL-based allele selection.

Haplotype-based genotyping-by-sequencing in oat genome research

In a de novo genotyping-by-sequencing (GBS) analysis of short, 64-base tag-level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag-level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635-line diversity panel were used to infer chromosome-level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high-resolution genome analysis and genomic selection in oats. A combined genome-wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component-based genome-wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS-derived markers facilitate genome analysis and genomic selection in oat.

Haplotype-based allele mining in the Japan-MAGIC rice population

Multi-parent advanced generation inter-cross (MAGIC) lines have broader genetic variation than bi-parental recombinant inbred lines. Genome-wide association study (GWAS) using high number of DNA polymorphisms such as single-nucleotide polymorphisms (SNPs) is a popular tool for allele mining in MAGIC populations, in which the associations of phenotypes with SNPs are investigated; however, the effects of haplotypes from multiple founders on phenotypes are not considered. Here, we describe an improved method of allele mining using the newly developed Japan-MAGIC (JAM) population, which is derived from eight high-yielding rice cultivars in Japan. To obtain information on the haplotypes in the JAM lines, we predicted the haplotype blocks in the whole chromosomes using 16,345 SNPs identified via genotyping-by-sequencing analysis. Using haplotype-based GWAS, we clearly detected the loci controlling the glutinous endosperm and culm length traits. Information on the alleles of the eight founders, which was based on the effects of mutations revealed by the analysis of next-generation sequencing data, was used to narrow down the candidate genes and reveal the associations between alleles and phenotypes. The haplotype-based allele mining (HAM) proposed in this study is a promising approach to the detection of allelic variation in genes controlling agronomic traits in MAGIC populations.

Genetic dissection of agronomically important traits in closely related temperate japonica rice cultivars

Many quantitative trait loci (QTLs) for agronomically important traits such as grain yield, disease resistance, and stress tolerance of rice (Oryza sativa L.) have been detected by using segregating populations derived from crosses between indica and japonica subspecies or with wild relatives. However, the QTLs involved in the control of natural variation in agronomic traits among closely related cultivars are still unclear. Decoding the whole genome sequences of Nipponbare and other temperate japonica rice cultivars has accelerated the collection of a huge number of single nucleotide polymorphisms (SNPs). These SNPs are good resource for developing polymorphic DNA markers and for detecting QTLs distributed across all rice chromosomes. The temperate japonica rice cultivar Koshihikari has remained the top cultivar for about 40 years since 1979 in Japan. Unraveling the genetic factors in Koshihikari will provide important insights into improving agronomic traits in temperate japonica rice cultivars. Here we describe recent progress in our studies as an example of genetic analysis in closely related cultivars.

Barcode System for Genetic Identification of Soybean [Glycine max (L.) Merrill] Cultivars Using InDel Markers Specific to Dense Variation Blocks

For genetic identification of soybean [Glycine max (L.) Merrill] cultivars, insertions/deletions (InDel) markers have been preferred currently because they are easy to use, co-dominant and relatively abundant. Despite their biological importance, the investigation of InDels with proven quality and reproducibility has been limited. In this study, we described soybean barcode system approach based on InDel makers, each of which is specific to a dense variation block (dVB) with non-random recombination due to many variations. Firstly, 2,274 VBs were mined by analyzing whole genome data in six soybean cultivars (Backun, Sinpaldal 2, Shingi, Daepoong, Hwangkeum, and Williams 82) for transferability to dVB-specific InDel markers. Secondly, 73,327 putative InDels in the dVB regions were identified for the development of soybean barcode system. Among them, 202 dVB-specific InDels from all soybean cultivars were selected by gel electrophoresis, which were converted as 2D barcode types according to comparing amplicon polymorphisms in the five cultivars to the reference cultivar. Finally, the polymorphism of the markers were assessed in 147 soybean cultivars, and the soybean barcode system that allows a clear distinction among soybean cultivars is also detailed. In addition, the changing of the dVBs in a chromosomal level can be quickly identified due to investigation of the reshuffling pattern of the soybean cultivars with 27 maker sets. Especially, a backcross-inbred offspring, "Singang" and a recurrent parent, "Sowon" were identified by using the 27 InDel markers. These results indicate that the soybean barcode system enables not only the minimal use of molecular markers but also comparing the data from different sources due to no need of exploiting allele binning in new varieties.

Single Marker and Haplotype-Based Association Analysis of Semolina and Pasta Colour in Elite Durum Wheat Breeding Lines Using a High-Density Consensus Map

Association mapping is usually performed by testing the correlation between a single marker and phenotypes. However, because patterns of variation within genomes are inherited as blocks, clustering markers into haplotypes for genome-wide scans could be a worthwhile approach to improve statistical power to detect associations. The availability of high-density molecular data allows the possibility to assess the potential of both approaches to identify marker-trait associations in durum wheat. In the present study, we used single marker- and haplotype-based approaches to identify loci associated with semolina and pasta colour in durum wheat, the main objective being to evaluate the potential benefits of haplotype-based analysis for identifying quantitative trait loci. One hundred sixty-nine durum lines were genotyped using the Illumina 90K Infinium iSelect assay, and 12,234 polymorphic single nucleotide polymorphism (SNP) markers were generated and used to assess the population structure and the linkage disequilibrium (LD) patterns. A total of 8,581 SNPs previously localized to a high-density consensus map were clustered into 406 haplotype blocks based on the average LD distance of 5.3 cM. Combining multiple SNPs into haplotype blocks increased the average polymorphism information content (PIC) from 0.27 per SNP to 0.50 per haplotype. The haplotype-based analysis identified 12 loci associated with grain pigment colour traits, including the five loci identified by the single marker-based analysis. Furthermore, the haplotype-based analysis resulted in an increase of the phenotypic variance explained (50.4% on average) and the allelic effect (33.7% on average) when compared to single marker analysis. The presence of multiple allelic combinations within each haplotype locus offers potential for screening the most favorable haplotype series and may facilitate marker-assisted selection of grain pigment colour in durum wheat. These results suggest a benefit of haplotype-based analysis over single marker analysis to detect loci associated with colour traits in durum wheat.

Genetic characterization and population structure of Indian rice cultivars and wild genotypes using core set markers

Genetic diversity among 23 rice genotypes including wild species and cultivars of indica, japonica, aus and aromatic type was investigated using 165 genomewide core set microsatellite (SSR) markers. This genotypic characterization was undertaken to know the genetic similarity among the parental lines to be used in developing a set of chromosome segment substitution lines. In all, 253 alleles were identified using 77 polymorphic SSRs, and polymorphism information content ranged from 0.31 to 0.97 with a mean of 0.79. Cluster analysis grouped the genotypes into three clusters at a genetic similarity of 0.26–0.75. Wild accessions grouped together in cluster-I, indica cultivars formed cluster-II, and aromatic, japonica and aus types came under cluster-III. Principal component analysis also showed similar results. The genotypic data was analyzed using STRUCTURE, and genotypes were grouped into four populations. RM1018 on chromosome 4, RM8009 on chromosome 7, and RM273 on chromosome 12 amplified alleles specific to wild accessions. The information obtained from core set markers would help in selecting diverse parents including wild accessions and for tracking alleles in mapping or breeding populations.

Improvement of Rice Biomass Yield through QTL-Based Selection

Biomass yield of rice (Oryza sativa L.) is an important breeding target, yet it is not easy to improve because the trait is complex and phenotyping is laborious. Using progeny derived from a cross between two high-yielding Japanese cultivars, we evaluated whether quantitative trait locus (QTL)-based selection can improve biomass yield. As a measure of biomass yield, we used plant weight (aboveground parts only), which included grain weight and stem and leaf weight. We measured these and related traits in recombinant inbred lines. Phenotypic values for these traits showed a continuous distribution with transgressive segregation, suggesting that selection can affect plant weight in the progeny. Four significant QTLs were mapped for plant weight, three for grain weight, and five for stem and leaf weight (at α = 0.05); some of them overlapped. Multiple regression analysis showed that about 43% of the phenotypic variance of plant weight was significantly explained (P < 0.0001) by six of the QTLs. From F2 plants derived from the same parental cross as the recombinant inbred lines, we divergently selected lines that carried alleles with positive or negative additive effects at these QTLs, and performed successive selfing. In the resulting F6 lines and parents, plant weight significantly differed among the genotypes (at α = 0.05). These results demonstrate that QTL-based selection is effective in improving rice biomass yield.

Open access resources for genome-wide association mapping in rice

Increasing food production is essential to meet the demands of a growing human population, with its rising income levels and nutritional expectations. To address the demand, plant breeders seek new sources of genetic variation to enhance the productivity, sustainability and resilience of crop varieties. Here we launch a high-resolution, open-access research platform to facilitate genome-wide association mapping in rice, a staple food crop. The platform provides an immortal collection of diverse germplasm, a high-density single-nucleotide polymorphism data set tailored for gene discovery, well-documented analytical strategies, and a suite of bioinformatics resources to facilitate biological interpretation. Using grain length, we demonstrate the power and resolution of our new high-density rice array, the accompanying genotypic data set, and an expanded diversity panel for detecting major and minor effect QTLs and subpopulation-specific alleles, with immediate implications for rice improvement.

Understanding the link between genotype and phenotype can facilitate efforts by breeders to utilize natural variation and develop new crop varieties. Here the authors present a diverse germplasm collection, a high-density genotyping array and a set of bioinformatic tools to enable association mapping in rice.

Quantitative trait loci for rice blast resistance detected in a local rice breeding population by genome-wide association mapping

Plant breeding programs aim to develop cultivars with high adaptability to the specific conditions in a local region. As a result, unique genes and gene combinations have been accumulated in local elite breeding populations during the long history of plant breeding. Genetic analyses on such genes and combinations may be useful for developing new cultivars with more-desirable agronomic traits. Here, we attempted to detect quantitative trait loci (QTL) for rice blast resistance (BR) using a local breeding rice population from Hokkaido, Japan. Using genotyping data on single nucleotide polymorphisms and simple sequence repeat markers distributed throughout the whole genomic region, we detected genetic regions associated with phenotypic variation in BR by a genome-wide association mapping study (GWAS). An additional association analysis using other breeding cultivars verified the effect and inheritance of the associated region. Furthermore, the existence of a gene for BR in the associated region was confirmed by QTL mapping. The results from these studies enabled us to estimate potential of the Hokkaido rice population as a gene pool for improving BR. The results of this study could be useful for developing novel cultivars with vigorous BR in rice breeding programs.

Allele mining and enhanced genetic recombination for rice breeding

Traditional rice varieties harbour a large store of genetic diversity with potential to accelerate rice improvement. For a long time, this diversity maintained in the International Rice Genebank has not been fully used because of a lack of genome information. The publication of the first reference genome of Nipponbare by the International Rice Genome Sequencing Project (IRGSP) marked the beginning of a systematic exploration and use of rice diversity for genetic research and breeding. Since then, the Nipponbare genome has served as the reference for the assembly of many additional genomes. The recently completed 3000 Rice Genomes Project together with the public database (SNP-Seek) provides a new genomic and data resource that enables the identification of useful accessions for breeding. Using disease resistance traits as case studies, we demonstrated the power of allele mining in the 3,000 genomes for extracting accessions from the GeneBank for targeted phenotyping. Although potentially useful landraces can now be identified, their use in breeding is often hindered by unfavourable linkages. Efficient breeding designs are much needed to transfer the useful diversity to breeding. Multi-parent Advanced Generation InterCross (MAGIC) is a breeding design to produce highly recombined populations. The MAGIC approach can be used to generate pre-breeding populations with increased genotypic diversity and reduced linkage drag. Allele mining combined with a multi-parent breeding design can help convert useful diversity into breeding-ready genetic resources.

Genome-wide indel markers shared by diverse Asian rice cultivars compared to Japanese rice cultivar 'Koshihikari'

Insertion-deletion (indel) polymorphisms, such as simple sequence repeats, have been widely used as DNA markers to identify QTLs and genes and to facilitate rice breeding. Recently, next-generation sequencing has produced deep sequences that allow genome-wide detection of indels. These polymorphisms can potentially be used to develop high-accuracy polymerase chain reaction (PCR)-based markers. Here, re-sequencing of 5 indica, 2 aus, and 3 tropical japonica cultivars and Japanese elite cultivar 'Koshihikari' was performed to extract regions containing large indels (10-51 bp) shared by diverse cultivars. To design indel markers for the discrimination of genomic regions between 'Koshihikari' and other diverse cultivars, we subtracted the indel regions detected in 'Koshihikari' from those shared in other cultivars. Two sets of indel markers, KNJ8-indel (shared in eight or more cultivars, including 'Khao Nam Jen' as a representative tropical japonica cultivar) and C5-indel (shared in five to eight cultivars), were established, with 915 and 9,899 indel regions, respectively. Validation of the two marker sets by using 23 diverse cultivars showed a high PCR success rate (≥95%) for 83.3% of the KNJ8-indel markers and 73.9% of the C5-indel markers. The marker sets will therefore be useful for the effective breeding of Japanese rice cultivars.

Metabolome-genome-wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome-wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.

Hybrid breakdown caused by epistasis-based recessive incompatibility in a cross of rice (Oryza sativa L.)

Viability and fertility in organisms depend on epistatic interactions between loci maintained in lineages. Here, we describe reduced fitness of segregants (hybrid breakdown, HB) that emerged in an F2 population derived from a cross between 2 rice (Oryza sativa L.) cultivars, "Tachisugata" (TS) and "Hokuriku 193" (H193), despite both parents and F1s showing normal fitness. Quantitative trait locus (QTL) analyses detected 13 QTLs for 4 morphological traits associated with the HB and 6 associated with principal component scores calculated from values of the morphological traits in the F2 population. Two-way analysis of variance of the putative QTLs identified 4 QTL pairs showing significant epistasis; among them, a pair on chromosomes 1 and 12 made the greatest contribution to HB. The finding was supported by genetic experiments using F3 progeny. HB emerged only when a plant was homozygous for the TS allele at the QTL on chromosome 1 and homozygous for the H193 allele at the QTL on chromosome 12, indicating that each allele behaves as recessive to the other. Our results support the idea that epistasis is an essential part of hybrid fitness.

Variation block-based genomics method for crop plants

BACKGROUND

In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits.

RESULTS

We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color.

CONCLUSIONS

We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

Genetic variation in resistance to blast disease (Pyricularia oryzae Cavara) in Japanese rice (Oryza sativa L.), as determined using a differential system

A total of 324 Japanese rice accessions, including landrace, improved, and weedy types were used to 1) investigate genetic variations in blast resistance to standard differential isolates, and 2) across the genome using polymorphism data on 64 SSR markers. From the polymorphism data, the accessions were classified into two clusters,. Accessions from irrigated lowland areas were included mainly in cluster I, and upland and Indica types were mainly in cluster II. The accessions were classified into three resistance subgroups, A2, B1 and B2, based on the reaction patterns to blast isolates. The accessions in A2 were postulated to have at least two resistance genes Pish and Pik-s, whereas those in B1 had various combinations of the resistance genes Pish, Pia, Pii, Pi3, Pi5(t), and Pik alleles. The B2 accessions were resistant to almost all isolates, and many accessions of cluster II were included, and had Pish, Pia, Pii, Pi3, Pi5(t), certain Pik, Piz and Pita alleles, and unknown genes. The frequencies of accessions of B1 originating in Hokkaido, and those of B2 originating in the Kanto and Tohoku regions were remarkably higher than in the other regions.

Genomic regions involved in yield potential detected by genome-wide association analysis in Japanese high-yielding rice cultivars

BACKGROUND

High-yielding cultivars of rice (Oryza sativa L.) have been developed in Japan from crosses between overseas indica and domestic japonica cultivars. Recently, next-generation sequencing technology and high-throughput genotyping systems have shown many single-nucleotide polymorphisms (SNPs) that are proving useful for detailed analysis of genome composition. These SNPs can be used in genome-wide association studies to detect candidate genome regions associated with economically important traits. In this study, we used a custom SNP set to identify introgressed chromosomal regions in a set of high-yielding Japanese rice cultivars, and we performed an association study to identify genome regions associated with yield.

RESULTS

An informative set of 1152 SNPs was established by screening 14 high-yielding or primary ancestral cultivars for 5760 validated SNPs. Analysis of the population structure of high-yielding cultivars showed three genome types: japonica-type, indica-type and a mixture of the two. SNP allele frequencies showed several regions derived predominantly from one of the two parental genome types. Distinct regions skewed for the presence of parental alleles were observed on chromosomes 1, 2, 7, 8, 11 and 12 (indica) and on chromosomes 1, 2 and 6 (japonica). A possible relationship between these introgressed regions and six yield traits (blast susceptibility, heading date, length of unhusked seeds, number of panicles, surface area of unhusked seeds and 1000-grain weight) was detected in eight genome regions dominated by alleles of one parental origin. Two of these regions were near Ghd7, a heading date locus, and Pi-ta, a blast resistance locus. The allele types (i.e., japonica or indica) of significant SNPs coincided with those previously reported for candidate genes Ghd7 and Pi-ta.

CONCLUSIONS

Introgression breeding is an established strategy for the accumulation of QTLs and genes controlling high yield. Our custom SNP set is an effective tool for the identification of introgressed genome regions from a particular genetic background. This study demonstrates that changes in genome structure occurred during artificial selection for high yield, and provides information on several genomic regions associated with yield performance.

Historical changes in population structure during rice breeding programs in the northern limits of rice cultivation

The rice local population was clearly differentiated into six groups over the 100-year history of rice breeding programs in the northern limit of rice cultivation over the world. Genetic improvements in plant breeding programs in local regions have led to the development of new cultivars with specific agronomic traits under environmental conditions and generated the unique genetic structures of local populations. Understanding historical changes in genome structures and phenotypic characteristics within local populations may be useful for identifying profitable genes and/or genetic resources and the creation of new gene combinations in plant breeding programs. In the present study, historical changes were elucidated in genome structures and phenotypic characteristics during 100-year rice breeding programs in Hokkaido, the northern limit of rice cultivation in the world. We selected 63 rice cultivars to represent the historical diversity of this local population from landraces to the current breeding lines. The results of the phylogenetic analysis demonstrated that these cultivars clearly differentiated into six groups over the history of rice breeding programs. Significant differences among these groups were detected in five of the seven traits, indicating that the differentiation of the Hokkaido rice population into these groups was correlated with these phenotypic changes. These results demonstrated that breeding practices in Hokkaido have created new genetic structures for adaptability to specific environmental conditions and breeding objectives. They also provide a new strategy for rice breeding programs in which such unique genes in local populations in the world can explore the genetic potentials of the local populations.

HapRice, an SNP haplotype database and a web tool for rice

Genome-wide single nucleotide polymorphism (SNP) analysis is a promising tool to examine the genetic diversity of rice populations and genetic traits of scientific and economic importance. Next-generation sequencing technology has accelerated the re-sequencing of diverse rice varieties and the discovery of genome-wide SNPs. Notably, validation of these SNPs by a high-throughput genotyping system, such as an SNP array, could provide a manageable and highly accurate SNP set. To enhance the potential utility of genome-wide SNPs for geneticists and breeders, analysis tools need to be developed. Here, we constructed an SNP haplotype database, which allows visualization of the allele frequency of all SNPs in the genome browser. We calculated the allele frequencies of 3,334 SNPs in 76 accessions from the world rice collection and 3,252 SNPs in 177 Japanese rice accessions; all these SNPs have been validated in our previous studies. The SNP haplotypes were defined by the allele frequency in each cultivar group (aus, indica, tropical japonica and temperate japonica) for the world rice accessions, and in non-irrigated and three irrigated groups (three variety registration periods) for Japanese rice accessions. We also developed web tools for finding polymorphic SNPs between any two rice accessions and for the primer design to develop cleaved amplified polymorphic sequence markers at any SNP. The 'HapRice' database and the web tools can be accessed at http://qtaro.abr.affrc.go.jp/index.html. In addition, we established a core SNP set consisting of 768 SNPs uniformly distributed in the rice genome; this set is of a practically appropriate size for use in rice genetic analysis.

Toward better annotation in plant metabolomics: isolation and structure elucidation of 36 specialized metabolites from Oryza sativa (rice) by using MS/MS and NMR analyses

Metabolomics plays an important role in phytochemical genomics and crop breeding; however, metabolite annotation is a significant bottleneck in metabolomic studies. In particular, in liquid chromatography-mass spectrometry (MS)-based metabolomics, which has become a routine technology for the profiling of plant-specialized metabolites, a substantial number of metabolites detected as MS peaks are still not assigned properly to a single metabolite. Oryza sativa (rice) is one of the most important staple crops in the world. In the present study, we isolated and elucidated the structures of specialized metabolites from rice by using MS/MS and NMR. Thirty-six compounds, including five new flavonoids and eight rare flavonolignan isomers, were isolated from the rice leaves. The MS/MS spectral data of the isolated compounds, with a detailed interpretation of MS fragmentation data, will facilitate metabolite annotation of the related phytochemicals by enriching the public mass spectral data depositories, including the plant-specific MS/MS-based database, ReSpect.

Molecular Genetic Diversity of Major Indian Rice Cultivars over Decadal Periods

Genetic diversity in representative sets of high yielding varieties of rice released in India between 1970 and 2010 was studied at molecular level employing hypervariable microsatellite markers. Of 64 rice SSR primer pairs studied, 52 showed polymorphism, when screened in 100 rice genotypes. A total of 184 alleles was identified averaging 3.63 alleles per locus. Cluster analysis clearly grouped the 100 genotypes into their respective decadal periods i.e., 1970s, 1980s, 1990s and 2000s. The trend of diversity over the decadal periods estimated based on the number of alleles (Na), allelic richness (Rs), Nei’s genetic diversity index (He), observed heterozygosity (Ho) and polymorphism information content (PIC) revealed increase of diversity over the periods in year of releasewise and longevitywise classification of rice varieties. Analysis of molecular variance (AMOVA) suggested more variation in within the decadal periods than among the decades. Pairwise comparison of population differentiation (Fst) among decadal periods showed significant difference between all the pairs except a few. Analysis of trends of appearing and disappearing alleles over decadal periods showed an increase in the appearance of alleles and decrease in disappearance in both the categories of varieties. It was obvious from the present findings, that genetic diversity was progressively on the rise in the varieties released during the decadal periods, between 1970s and 2000s.

Genetical and morphological characterization of cold tolerance at fertilization stage in rice

Cold temperature during the reproductive phase leads to seed sterility, which reduces yield and decreases the grain quality of rice. The fertilization stage, ranging from pollen maturation to the completion of fertilization, is sensitive to unsuitable temperature. Improving cold tolerance at the fertilization stage (CTF) is an important objective of rice breeding program in cold temperature areas. In this study, we characterized fertilization behavior under cold temperature to define the phenotype of CTF and identified quantitative trait loci (QTLs) for CTF. A wide variation in CTF levels has been identified among local cultivars in Hokkaido, which is one of the most northern regions for rice cultivation in the world. Clear varietal differences in pollen germination, and pollen tube elongation due to cold temperature have been observed. These differences may confer a degree of CTF among this population. We conducted QTL analysis for CTF using 120 backcrossed inbred lines derived from a cross between Eikei88223 (vigorous CTF) and Suisei (very weak CTF). Three QTLs for CTF were identified. A clear effect by QTL, qCTF7, for increasing the level of CTF was validated using advanced progeny. These results will facilitate marker-assist selection for desirable QTLs for CTF in rice breeding program.

Genes offering the potential for designing yield-related traits in rice

Breeding of high-yielding rice is crucial for meeting the food demand of the increasing world population. New technologies have facilitated identification of genes involved in quantitative traits, and many genes underpinning quantitative trait loci involved in rice crop yield have been isolated. Meanwhile, various kinds of mutants have been intensively studied, leading to characterization of many genes related to yield traits. A combination of quantitative trait locus analysis and studies of such mutants has made it possible to compile a list of genes available for breeding rice with higher yield.

Genome-wide identification and analysis of miRNA-related single nucleotide polymorphisms (SNPs) in rice

BACKGROUND

MiRNAs are key regulators in the miRNA-mediated regulatory networks. Single nucleotide polymorphisms (SNPs) that occur at miRNA-related regions may cause serious phenotype changes. To gain new insights into the evolution of miRNAs after SNP variation, we performed a genome-wide scan of miRNA-related SNPs, and analyzed their effects on the stability of miRNAs structure and the alteration of target spectrum in rice.

RESULTS

We find that the SNP density in pre-miRNAs is significantly higher than that in the flanking regions, owing to the rapid evolution of a large number of species-specific miRNAs in rice. In contrast, it is obvious that deeply conserved miRNAs are under strong purifying selection during evolution. In most cases, the SNPs in stem regions may result in the miRNA hairpin structures changing from stable to unstable status; And SNPs in mature miRNAs have great potential to have either newly created or disrupted the miRNA-target interactions. However, the total number of gained targets is over 2.5 times greater than that are lost after mutation. Notably, 12 putative domestication-related miRNAs have been identified, where the SNP density is significantly lower.

CONCLUSIONS

The present study provides the first outline of SNP variations occurred in rice pre-miRNAs at the whole genome-wide level. These analyses may deepen our understanding on the effects of SNPs on the evolution of miRNAs in the rice genome.

Roles of the Hd5 gene controlling heading date for adaptation to the northern limits of rice cultivation

During the diversification of cultivated rice after domestication, rice was grown in diverse geographic regions using genetic variations attributed to the combination of alleles in loci for adaptability to various environmental conditions. To elucidate the key gene for adaptation in rice cultivars to the northern limit of rice cultivation, we conducted genetic analyses of heading date using extremely early-heading cultivars. The Hd5 gene controlling heading date (flowering time) generated variations in heading date among cultivars adapted to Hokkaido, where is the northernmost region of Japan and one of the northern limits of rice cultivation in the world. The association of the Hd5 genotype with heading date and genetical analysis clearly showed that the loss-of-function Hd5 has an important role in exhibiting earlier heading among a local population in Hokkaido. Distinct distribution of the loss-of-function Hd5 revealed that this mutation event of the 19-bp deletion occurred in a local landrace Bouzu and that this mutation may have been selected as an early-heading variety in rice breeding programs in Hokkaido in the early 1900s. The loss-of-function Hd5 was then introduced into the rice variety Fanny from France and contributed to its extremely early heading under the presence of functional Ghd7. These results demonstrated that Hd5 plays roles not only in generating early heading in variations of heading date among a local population in Hokkaido, but also in extremely early heading for adaptation to northern limits of rice cultivation.