Studies on Ancient Rice—Where Botanists, Agronomists, Archeologists, Linguists, and Ethnologists Meet

Taiwan’s aboriginal peoples are thought to be related to ancestral Austronesian-speaking peoples. Currently, Taiwan has 14 officially acknowledged aboriginal tribes. The major crops currently farmed in aboriginal areas are rice (Oryza sativa) and foxtail millet (Setaria italica). Archeologists recently excavated the remains of several early cultures in Taiwan. The most plentiful plant remains were carbonated rice and foxtail millet grains. The earliest ¹⁴C date of these excavation sites is ～5,000 bp. These settlements may be those of the earliest ancestral Austronesian speakers in Taiwan. Rice domestication is a complex story. In this study, we identified the functional nucleotide polymorphisms of 16 domestication-related genes using 60 landraces collected from aboriginal Taiwanese villages about 100 years ago. We also screened the phenotypes of these landraces. By integrating pheno- and genotypic data, together with data from archeologists and linguists, we may be able to better understand the history of rice cultivation in Taiwan and nearby areas.     

The Race Structure of the Rice Blast Pathogen Across Southern and Northeastern China

Background

Rice blast, caused by the ascomycete Magnaporthe oryzae (Mo), imposes a major constraint on rice productivity. Managing the disease through the deployment of host resistance requires a close understanding of race structure of the pathogen population.

Results

The host/pathogen interaction between isolates sampled from four Mo populations collected across the rice-producing regions of China was tested using two established panels of differential cultivars. The clearest picture was obtained from the Chinese cultivar panel, for which the frequency of the various races, the race diversity index, the specific race isolate frequency, and the frequency of the three predominant races gave a consistent result, from which it was concluded that the pathogen population present in the southern production region was more diverse than that in the northeastern region. The four blast resistance genes Pi1, Pik, Pik-m, and Piz all still remain effective in the southern China rice production area, as does Pi1 in the northeastern region. The effectiveness of Pita, Pik-p, Piz, and Pib is restricted to single provinces. The distinctive resistance profile shown by the Chinese differential cultivar set implied the presence of at least five as yet unidentified blast resistance genes.

Conclusions

The Chinese differential cultivar set proved to be more informative than the Japanese one for characterizing the race structure of the rice blast pathogen in China. A number of well characterized host resistance genes, in addition to some as yet uncharacterized ones, remain effective across the major rice production regions in China.

     

Higher yields of hybrid rice do not depend on nitrogen fertilization under moderate to high soil fertility conditions

Background

Increasing rice yield with fewer external inputs is critical to ensuring food security, reducing environmental costs, and improving returns. Use of hybrid rice has expanded greatly in China due to its higher yield potential. Meanwhile, large and increasing amounts of nitrogen (N) fertilizers have been used for expanding rice production in China. It is not clear to what extent the success of hybrid rice in China is associated with N fertilizer inputs.

Findings

We observed that the higher grain yield with N fertilizer in hybrid rice was driven more by a higher yield without N fertilizer than by increases in grain yield with N fertilizer under moderate to high soil fertility conditions.

Conclusions

Our results suggest that greater application of N fertilizers is not needed to benefit from hybrid rice production under moderate to high soil fertility conditions, and that improving and maintaining soil fertility should be a focus for sustaining hybrid rice production. Moreover, our study also indicates that zero-N testing may be a potentially useful tool to develop hybrid rice with high yield and without requirement of greater external N inputs under moderate to high soil fertility conditions.

     

<Emphasis Type="Italic">OsLAP6/OsPKS1</Emphasis>, an orthologue of <Emphasis Type="Italic">Arabidopsis PKSA/LAP6</Emphasis>, is critical for proper pollen exine formation

Background

Male fertility is crucial for rice yield, and the improvement of rice yield requires hybrid production that depends on male sterile lines. Although recent studies have revealed several important genes in male reproductive development, our understanding of the mechanisms of rice pollen development remains unclear.

Results

We identified a rice mutant oslap6 with complete male sterile phenotype caused by defects in pollen exine formation. By using the MutMap method, we found that a single nucleotide polymorphism (SNP) variation located in the second exon of OsLAP6/OsPKS1 was responsible for the mutant phenotype. OsLAP6/OsPKS1 is an orthologous gene of Arabidopsis PKSA/LAP6, which functions in sporopollenin metabolism. Several other loss-of-function mutants of OsLAP6/OsPKS1 generated by the CRISPR/Cas9 genomic editing tool also exhibited the same phenotype of male sterility. Our cellular analysis suggested that OsLAP6/OsPKS1 might regulate pollen exine formation by affecting bacula elongation. Expression examination indicated that OsLAP6/OsPKS1 is specifically expressed in tapetum, and its product is localized to the endoplasmic reticulum (ER). Protein sequence analysis indicated that OsLAP6/OsPKS1 is conserved in land plants.

Conclusions

OsLAP6/OsPKS1 is a critical molecular switch for rice male fertility by participating in a conserved sporopollenin precursor biosynthetic pathway in land plants. Manipulation of OsLAP6/OsPKS1 has potential for application in hybrid rice breeding.

     

Updating the elite rice variety Kongyu 131 by improving the <Emphasis Type="Italic">Gn1a</Emphasis> locus

Background

Kongyu 131 is an elite japonica rice variety of Heilongjiang Province, China. It has the characteristics of early maturity, superior quality, high yield, cold tolerance and wide adaptability. However, there is potential to improve the yield of Kongyu 131 because of the relatively few grains per panicle compared with other varieties. Hence, we rebuilt the genome of Kongyu 131 by replacing the GRAIN NUMBER1a (Gn1a) locus with a high-yielding allele from a big panicle indica rice variety, GKBR. High-resolution melting (HRM) analysis was used for single nucleotide polymorphism (SNP) genotyping.

Results

Quantitative trait locus (QTL) analysis of the BC₃F₂ population showed that the introgressed segment carrying the Gn1a allele of GKBR significantly increased the branch number and grain number per panicle. Using 5 SNP markers designed against the sequence within and around Gn1a, the introgressed chromosome segment was shortened to approximately 430 Kb to minimize the linkage drag by screening recombinants in the target region. Genomic components of the new Kongyu 131 were detected using 220 SNP markers evenly distributed across 12 chromosomes, suggesting that the recovery ratio of the recurrent parent genome (RRPG) was 99.89%. Compared with Kongyu 131, the yield per plant of the new Kongyu 131 increased by 8.3% and 11.9% at Changchun and Jiamusi, respectively.

Conclusions

To achieve the high yield potential of Kongyu 131, a minute chromosome fragment carrying the favorable Gn1a allele from the donor parent was introgressed into the genome of Kongyu 131, which resulted in a larger panicle and subsequent yield increase in the new Kongyu 131. These results indicate the feasibility of improving an undesirable trait of an elite variety by replacing only a small chromosome segment carrying a favorable allele.

     

Divergent evolution of rice blast resistance <Emphasis Type="Italic">Pi54</Emphasis> locus in the genus <Emphasis Type="Italic">Oryza</Emphasis>

Background

The rice blast resistance gene Pi54 was cloned from Oryza sativa ssp. indica cv. Tetep, which conferred broad-spectrum resistance against Magnaporthe oryzae. Pi54 allelic variants have been identified in not only domesticates but also wild rice species, but the majority of japonica and some indica cultivars lost the function.

Results

We here found that Pi54 (Os11g0639100) and its homolog Os11g0640600 (named as #11) were closely located on a 25 kbp region in japonica cv. Sasanishiki compared to a 99 kbp region in japonica cv. Nipponbare. Sasanishiki lost at least six genes containing one other R-gene cluster (Os11g0639600, Os11g0640000, and Os11g0640300). Eight AA-genome species including five wild rice species were classified into either Nipponbare or Sasanishiki type. The BB-genome wild rice species O. punctata was Sasanishiki type. The FF-genome wild rice species O. brachyantha (the basal lineage of Oryza) was neither, because Pi54 was absent and the orientation of the R-gene cluster was reversed in comparison with Nipponbare-type species. The phylogenetic analysis showed that #11gene of O. brachyantha was on the root of both Pi54 and #11 alleles. All Nipponbare-type Pi54 alleles were specifically disrupted by 143 and 37/44?bp insertions compared to Tetep and Sasanishiki type. In addition, Pi54 of japonica cv. Sasanishiki lost nucleotide-binding site and leucine-rich repeat (NBS–LRR) domains owing to additional mutations.

Conclusions

These results suggest that Pi54 might be derived from a tandem duplication of the ancestor #11 gene in progenitor FF-genome species. Two divergent structures of Pi54 locus caused by a mobile unit containing the nearby R-gene cluster could be developed before domestication. This study provides a potential genetic resource of rice breeding for blast resistance in modern cultivars sustainability.

     

Identification and validation of a novel major QTL for harvest index in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Background

Harvest index (HI) in rice is defined as the ratio of grain yield (GY) to biomass (BM). Although it has been demonstrated that HI is significantly related to yield and is considered as one of the most important traits in high-yielding rice breeding, HI-based high-yielding rice breeding is difficult due to its polygenic nature and insufficient knowledge on the genetic basis of HI. Therefore, searching for rice varieties with high HI and mapping genes associated with high HI can facilitate marker-assisted breeding for high HI in rice.

Results

Yuexiangzhan, a popular indica cultivar with good reputation of high HI was crossed with Shengbasimiao, an indica cultivar with lower HI to develop a recombinant inbred line population, and QTL mapping for HI and its component traits was conducted. In total, five QTLs for HI, three QTLs for GY, and six QTLs for BM were detected in two-year experiments. Among the three GY QTLs, one co-located with the HI QTL on chromosome 8, while the other two co-located with the two tightly-linked BM QTLs on chromosome 3. The co-located QTLs in each of the chromosomal regions produced additive effects in the same direction. Particularly, the HI QTL on chromosome 8, qHI-8, could be detected across two years and explained 42.8% and 44.5% of the phenotypic variation, respectively. The existence of qHI-8 was confirmed by the evaluation of the near isogenic lines derived from a residual heterozygous line, and this QTL was delimitated to a 1070 kb interval by substitution mapping.

Conclusion

In the present study, the detected GY QTLs overlapped with both HI QTL and BM QTL, suggesting a positive relationship between GY and HI or BM, respectively. With an understanding of the genetic basis for grain yield, harvest index and biomass, it is possible to achieve higher yield through enhancing HI and BM by pyramiding the favorable alleles for the two traits via marker-assisted selection (MAS). As qHI-8 has a large phenotypic effect on HI and expresses stably in different environments, it provides a promising target for further genetic characterization of HI and MAS of high HI in rice breeding.

     

Marker assisted pyramiding of <Emphasis Type="Italic">Bph6</Emphasis> and <Emphasis Type="Italic">Bph9</Emphasis> into elite restorer line 93–11 and development of functional marker for <Emphasis Type="Italic">Bph9</Emphasis>

Background

The brown planthopper (BPH) has become the most destructive and a serious threat to the rice production in Asia. Breeding the resistant varieties with improved host resistance is the most effective and ecosystem-friendly strategy of BPH biological management. As host resistance was always broken down by the presence of the upgrading BPH biotype, the more resistant varieties with novel resistance genes or pyramiding known identified BPH resistance genes would be needed urgently for higher resistant level and more durability of resistance.

Results

Here, we developed near isogenic lines of Bph9 (NIL-Bph9) by backcrossing elite cultivar 93–11 with Pokkali (harboring Bph9) using marker-assisted selection (MAS). Subsequently, we pyramided Bph6 and Bph9 in 93–11 genetic background through MAS. The resulting Bph6 and Bph9 pyramided line LuoYang69 had stronger antixenotic and antibiosis effects on BPH and exhibited significantly enhanced resistance to BPH than near isogenic lines NIL-Bph6 and NIL-Bph9. LuoYang69 derived hybrids, harboring heterozygous Bph6 and Bph9 genes, also conferred high level of resistance to BPH. Furthermore, LuoYang69 did not affect the elite agronomic traits and rice grain quality of 93–11. The current study also developed functional markers for Bph9. Using functional dominant marker, we screened and evaluated worldwide accessions of rice germplasm. Of the 673 varieties tested, 8 cultivars were identified to harbor functional Bph9 gene.

Conclusion

The development of Bph6 and Bph9 pyramided line LuoYang69 provides valuable resource to develop hybrid rice with highly and durable BPH resistance. The development of functional markers will promote MAS of Bph9. The identified Bph9 containing cultivars can be used as new sources for BPH resistance breeding programs.

     

Large-scale deployment of a rice 6 K SNP array for genetics and breeding applications

Background

Fixed arrays of single nucleotide polymorphism (SNP) markers have advantages over reduced representation sequencing in their ease of data analysis, consistently higher call rates, and rapid turnaround times. A 6 K SNP array represents a cost-benefit “sweet spot” for routine genetics and breeding applications in rice. Selection of informative SNPs across species and subpopulations during chip design is essential to obtain useful polymorphism rates for target germplasm groups. This paper summarizes results from large-scale deployment of an Illumina 6 K SNP array for rice.

Results

Design of the Illumina Infinium 6 K SNP chip for rice, referred to as the Cornell_6K_Array_Infinium_Rice (C6AIR), includes 4429 SNPs from re-sequencing data and 1571 SNP markers from previous BeadXpress 384-SNP sets, selected based on polymorphism rate and allele frequency within and between target germplasm groups. Of the 6000 attempted bead types, 5274 passed Illumina’s production quality control. The C6AIR was widely deployed at the International Rice Research Institute (IRRI) for genetic diversity analysis, QTL mapping, and tracking introgressions and was intensively used at Cornell University for QTL analysis and developing libraries of interspecific chromosome segment substitution lines (CSSLs) between O. sativa and diverse accessions of O. rufipogon or O. meridionalis. Collectively, the array was used to genotype over 40,000 rice samples. A set of 4606 SNP markers was used to provide high quality data for O. sativa germplasm, while a slightly expanded set of 4940 SNPs was used for O. sativa X O. rufipogon populations. Biparental polymorphism rates were generally between 1900 and 2500 well-distributed SNP markers for indica x japonica or interspecific populations and between 1300 and 1500 markers for crosses within indica, while polymorphism rates were lower for pairwise crosses within U.S. tropical japonica germplasm. Recently, a second-generation array containing ~7000 SNP markers, referred to as the C7AIR, was designed by removing poor-performing SNPs from the C6AIR and adding markers selected to increase the utility of the array for elite tropical japonica material.

Conclusions

The C6AIR has been successfully used to generate rapid and high-quality genotype data for diverse genetics and breeding applications in rice, and provides the basis for an optimized design in the C7AIR.

     

Genome-wide association mapping for root cone angle in rice

Background

Plant root systems play a major role in anchoring and in water and nutrient uptake from the soil. The root cone angle is an important parameter of the root system architecture because, combined with root depth, it helps to determine the volume of soil explored by the plant. Two genes, DRO1 and SOR1, and several QTLs for root cone angle have been discovered in the last 5 years.

Results

To find other QTLs linked to root cone angle, a genome-wide association mapping study was conducted on two panels of 162 indica and 169 japonica rice accessions genotyped with two sets of SNP markers (genotyping-by-sequencing set with approximately 16,000 markers and high-density-rice-array set with approximately 300,000 markers). The root cone angle of all accessions was measured using a screen protractor on images taken after 1 month of plant growth in the Rhizoscope phenotyping system. The distribution of the root cone angle in the indica panel was Gaussian, but several accessions of the japonica panel (all the bulus from Indonesia and three temperate japonicas from Nepal or India) appeared as outliers with a very wide root cone angle. The data were submitted to association mapping using a mixed model with control of structure and kinship. A total of 15 QTLs for the indica panel and 40 QTLs for the japonica panel were detected. Genes underlying these QTLs (+/?50 kb from the significant markers) were analyzed. We focused our analysis on auxin-related genes, kinases, and genes involved in root developmental processes and identified 8 particularly interesting genes.

Conclusions

The present study identifies new sources of wide root cone angle in rice, proposes ways to bypass some drawbacks of association mapping to further understand the genetics of the trait and identifies candidate genes deserving further investigation.

     

Genome-wide association study and candidate gene analysis of rice cadmium accumulation in grain in a diverse rice collection

Background

Cadmium (Cd) accumulation in rice followed by transfer to the food chain causes severe health problems in humans. Breeding of low Cd accumulation varieties is one of the most economical ways to solve the problem. However, information on the identity of rice germplasm with low Cd accumulation is limited, particularly in indica, and the genetic basis of Cd accumulation in rice is not well understood.

Results

Screening of 312 diverse rice accessions revealed that the grain Cd concentrations of these rice accessions ranged from 0.12 to 1.23?mg/kg, with 24 accessions less than 0.20?mg/kg. Three of the 24 accessions belong to indica. Japonica accumulated significantly less Cd than indica (p < 0.001), while tropical japonica accumulated significantly less Cd than temperate japonica (p < 0.01). GWAS in all accessions identified 14 QTLs for Cd accumulation, with 7 identified in indica and 7 identified in japonica subpopulations. No common QTL was identified between indica and japonica. The previously identified genes (OsHMA3, OsNRAMP1, and OsNRAMP5) from japonica were colocalized with QTLs identified in japonica instead of indica. Expression analysis of OsNRAMP2, the candidate gene of the novel QTL (qCd3–2) identified in the present study, demonstrated that OsNRAMP2 was mainly induced in the shoots of high Cd accumulation accessions after Cd treatment. Four amino acid differences were found in the open reading frame of OsNRAMP2 between high and low Cd accumulation accessions. The allele from low Cd accumulation accessions significantly increased the Cd sensitivity and accumulation in yeast. Subcellular localization analysis demonstrated OsNRAMP2 expressed in the tonoplast of rice protoplast.

Conclusion

The results suggest that grain Cd concentrations are significantly different among subgroups, with Cd concentrations decreasing from indica to temperate japonica to tropical japonica. However, considerable variations exist within subgroups. The fact that no common QTL was identified between indica and japonica implies that there is a different genetic basis for determining Cd accumulation between indica and japonica, or that some QTLs for Cd accumulation in rice are subspecies-specific. Through further integrated analysis, it is speculated that OsNRAMP2 could be a novel functional gene associated with Cd accumulation in rice.