A meta-QTL analysis highlights genomic hotspots associated with phosphorus use efficiency in rice (Oryza sativa L.).

Phosphorus use efficiency (PUE) is a complex trait, governed by many minor quantitative trait loci (QTLs) with small effects. Advances in molecular marker technology have led to the identification of QTLs underlying PUE. However, their practical use in breeding programs remains challenging due to the unstable effects in different genetic backgrounds and environments, interaction with soil status, and linkage drag. Here, we compiled PUE QTL information from 16 independent studies. A total of 192 QTLs were subjected to meta-QTL (MQTL) analysis and were projected into a high-density SNP consensus map. A total of 60 MQTLs, with significantly reduced number of initial QTLs and confidence intervals (CI), were identified across the rice genome. Candidate gene (CG) mining was carried out for the 38 MQTLs supported by multiple QTLs from at least two independent studies. Genes related to amino and organic acid transport and auxin response were found to be abundant in the MQTLs linked to PUE. CGs were cross validated using a root transcriptome database (RiceXPro) and haplotype analysis. This led to the identification of the eight CGs (OsARF8, OsSPX-MFS3, OsRING141, OsMIOX, HsfC2b, OsFER2, OsWRKY64, and OsYUCCA11) modulating PUE. Potential donors for superior PUE CG haplotypes were identified through haplotype analysis. The distribution of superior haplotypes varied among subspecies being mostly found in indica but were largely scarce in japonica. Our study offers an insight on the complex genetic networks that modulate PUE in rice. The MQTLs, CGs, and superior CG haplotypes identified in our study are useful in the combination of beneficial alleles for PUE in rice.

A Genomic Evaluation of Six Selected Inbred Lines of the Naturalized Plants of Milk Thistle (Silybum marianum L. Gaertn.) in Korea.

Milk thistle (Silybum marianum) belongs to the Asteraceae family and is a medicinal plant native to the Mediterranean Basin. Silymarin in achene is a widely used herbal product for chronic liver disease. There is growing interest in natural medicine using milk thistle in Korea, but the raw material completely relies on imports. Despite its economic importance, phenotypic evaluations of native resources of milk thistle in Korea have not been carried out. In addition, genomic research and molecular marker development are very limited in milk thistle. In this study, we evaluated 220 milk thistle resources consisting of 172 accessions collected from the domestic market, and 48 accessions isolated from 6 accessions distributed by the National Agrobiodiversity Center in Korea. Six plant characteristics (height, seed weight, number of flowers, seed weight per flower, spine length, and color at harvest) were measured, and six samples (M01-M06) were selected to represent the genetic diversity of the population for genomic research. To develop PCR-based and co-dominant insertion/deletion (InDel) markers, we performed genome-wide InDel detection by comparing the whole-genome resequencing data of the six selected accessions with the reference genome sequence (GCA_001541825). As a result, 177 InDel markers with high distinguishability and reproducibility were selected from the 30,845 InDel variants. Unknowingly imported alien plant resources could easily be genetically mixed, and jeopardized seed purity can cause continuous difficulties in the development of high value-added agricultural platforms utilizing natural products. The selected plant materials and 177 validated InDel markers developed via whole-genome resequencing analysis could be valuable resources for breeding, conservation, and ecological studies of natives to Korea, along with acceleration of Silybum marianum industrialization.

Development of a Temperate Climate-Adapted indica Multi-stress Tolerant Rice Variety by Pyramiding Quantitative Trait Loci.

Successful cultivation of rice (Oryza sativa L.) in many Asian countries requires submergence stress tolerance at the germination and early establishment stages. Two quantitative trait loci, Sub1 (conferring submergence tolerance) and AG1 (conferring anaerobic germination), were recently pyramided into a single genetic background, without compromising any desirable agronomic traits, leading to the development of Ciherang-Sub1 + AG1 (CSA). However, little research has been conducted to enhance plant tolerance to abiotic stress (submergence) and biotic stress (rice blast), which occur in a damp climate following flooding. The BC2F5 breeding line was phenotypically characterized using the AvrPi9 isolate. The biotic and abiotic stress tolerance of selected lines was tested under submergence stress and anaerobic germination conditions, and lines tolerant to each stress condition were identified through phenotypic and gene expression analyses. The Ciherang-Sub1 + AG1 + Pi9 (CSA-Pi9) line showed similar agronomic performance to its recurrent parent, CSA, but had significantly reduced chalkiness in field trials conducted in temperate regions. Unexpectedly, the CSA-Pi9 line also showed salinity tolerance. Thus, the breeding line newly developed in this study, CSA-Pi9, functioned under stress conditions, in which Sub1, AG1, and Pi9 play a role and had superior grain quality traits compared to its recurrent parent in temperate regions. We speculate that CSA-Pi9 will enable the establishment of climate-resilient rice cropping systems, particularly in East Asia.

The Effect of Water Level in Rice Cropping System on Phosphorus Uptake Activity of Pup1 in a Pup1+Sub1 Breeding Line.

Pyramiding useful QTLs into an elite variety is a promising strategy to develop tolerant varieties against multiple abiotic stresses. However, some QTLs may not be functionally compatible when they are introgressed into the same variety. Here, we tested the functional compatibility of Pup1 and Sub1, major QTLs for tolerance to phosphorus (P)-deficiency and submergence conditions, respectively. Phenotypic analysis revealed that IR64-Pup1+Sub1 (IPS) plants harboring both Pup1 and Sub1 QTLs show significant tolerance to submerged conditions, similarly to IR64-Sub1, while IPS failed to tolerate P deficiency and mild drought conditions; only IR64-Pup1 showed P deficiency tolerance. In submerged conditions, Sub1A and OsPSTOL1, major genes for Sub1 and Pup1 QTLs, respectively, were expressed at the same levels as in IPS and IR64-Sub1 and in IPS and IR64-Pup1, respectively. On the other hand, in P-non-supplied condition, crown root number, root length, and OsPSTOL1 expression level were significantly lower in IPS compared to those of IR64-Pup1. However, there was no significant difference in P content between IPS and IR64-Pup1. These results imply that Pup1 does not compromise Sub1 function in submerged condition, while Sub1 suppresses Pup1 function in P-non-supplied condition, possibly by regulating the transcript level of Pup1. In conclusion, Pup1 and Sub1 are regarded as functionally compatible under submergence condition but not under P-non-supplied condition. Further study is needed to elucidate the functional incompatibility of Pup1 and Sub1 QTLs in IPS under P-non-supplied condition.

Characterization of the Common Japonica-Originated Genomic Regions in the High-Yielding Varieties Developed from Inter-Subspecific Crosses in Temperate Rice (Oryza sativa L.).

The inter-subspecific crossing between indica and japonica subspecies in rice have been utilized to improve the yield potential of temperate rice. In this study, a comparative study of the genomic regions in the eight high-yielding varieties (HYVs) was conducted with those of the four non-HYVs. The Next-Generation Sequencing (NGS) mapping on the Nipponbare reference genome identified a total of 14 common genomic regions of japonica-originated alleles. Interestingly, the HYVs shared japonica-originated genomic regions on nine chromosomes, although they were developed through different breeding programs. A panel of 94 varieties was classified into four varietal groups with 38 single nucleotide polymorphism (SNP) markers from 38 genes residing in the japonica-originated genomic regions and 16 additional trait-specific SNPs. As expected, the japonica-originated genomic regions were only present in the japonica (JAP) and HYV groups, except for Chr4-1 and Chr4-2. The Wx gene, located within Chr6-1, was present in the HYV and JAP variety groups, while the yield-related genes were conserved as indica alleles in HYVs. The japonica-originated genomic regions and alleles shared by HYVs can be employed in molecular breeding programs to further develop the HYVs in temperate rice.

Mapping and Validation of QTLs for the Amino Acid and Total Protein Content in Brown Rice.

Highly nutritious rice production will be benefited with the improvement of amino acid content (AAC) and protein content (PC). The identification of quantitative trait loci (QTLs) associated with the PC and AAC of rice grains could provide a basis for improving the nutritional value of rice grains. Here, we conducted QTL analyses using recombinant inbred lines from the cross between indica (Milyang 23 or M23) and japonica (Tong 88-7 or T887) rice varieties, afterward employing genotyping-by-sequencing to obtain a high-density genetic map. A total of 17 and 3 QTLs were detected for AAC and PC, respectively. Among them, two QTLs associated with more than 10 AACs, qAAC6.1 and qAAC7.1, were identified for the first time in this study. Each favorable allele that increased the AAC of the two QTLs was derived from M23 and T887, respectively. Allelic combination of qAAC6.1 M23 and qAAC7.1T887 showed significantly higher content of associated amino acids (AAs) than other allelic combinations. Near-isogenic line (NIL) possessing qAAC7.1T887 with M23 genetic background had significantly higher AACs than both parents. These results indicate that the pyramiding of QTLs would be useful in developing brown rice with a high AA and protein content.