Molecular basis of genetic improvement for key rice quality traits in Southern China.

Grain qualities including milling quality, appearance quality, eating and cooking quality, and nutritional quality are important indicators in rice breeding. Significant achievements in genetic improvement of rice quality have been made. In this study, we analyzed the variation patterns of 16 traits in 1570 rice varieties and found significant improvements in appearance quality and eating and cooking quality, particularly in hybrid rice. Through genome-wide association study and allelic functional nucleotide polymorphisms analysis of quality trait genes, we found that ALK, FGR1, FLO7, GL7/GW7, GLW7, GS2, GS3, ONAC129, OsGRF8, POW1, WCR1, and Wx were associated with the genetic improvement of rice quality traits in Southern China. Allelic functional nucleotide polymorphisms analysis of 13 important rice quality genes, including fragrance gene fgr, were performed using the polymerase chain reaction amplification refractory mutation system technology. The results showed that Gui516, Gui569, Gui721, Ryousi, Rsimiao, Rbasi, and Yuehui9802 possessed multiple superior alleles. This study elucidates the phenotypic changes and molecular basis of key quality traits of varieties in Southern China. The findings will provide guidance for genetic improvement of rice quality and the development of new varieties.

Genome-wide association study using specific-locus amplified fragment sequencing identifies new genes influencing nitrogen use efficiency in rice landraces.

Nitrogen is essential for crop production. It is a critical macronutrient for plant growth and development. However, excessive application of nitrogen fertilizer is not only a waste of resources but also pollutes the environment. An effective approach to solving this problem is to breed rice varieties with high nitrogen use efficiency (NUE). In this study, we performed a genome-wide association study (GWAS) on 419 rice landraces using 208,993 single nucleotide polymorphisms (SNPs). With the mixed linear model (MLM) in the Tassel software, we identified 834 SNPs associated with root surface area (RSA), root length (RL), root branch number (RBN), root number (RN), plant dry weight (PDW), plant height (PH), root volume (RL), plant fresh weight (PFW), root fractal dimension (RFD), number of root nodes (NRN), and average root diameter (ARD), with a significant level of p < 2.39×10-7. In addition, we found 49 SNPs that were correlated with RL, RBN, RN, PDW, PH, PFW, RFD, and NRN using genome-wide efficient mixed-model association (GEMMA), with a significant level of p < 1×10-6. Additionally, the final results for eight traits associated with 193 significant SNPs by using multi-locus random-SNP-effect mixed linear model (mrMLM) model and 272 significant SNPs associated with 11 traits by using IIIVmrMLM. Within the linkage intervals of significantly associated SNP, we identified eight known related genes to NUE in rice, namely, OsAMT2;3, OsGS1, OsNR2, OsNPF7.4, OsPTR9, OsNRT1.1B, OsNRT2.3, and OsNRT2.2. According to the linkage disequilibrium (LD) decay value of this population, there were 75 candidate genes within the 150-kb regions upstream and downstream of the most significantly associated SNP (Chr5_29804690, Chr5_29956584, and Chr10_17540654). These candidate genes included 22 transposon genes, 25 expressed genes, and 28 putative functional genes. The expression levels of these candidate genes were measured by real-time quantitative PCR (RT-qPCR), and the expression levels of LOC_Os05g51700 and LOC_Os05g51710 in C347 were significantly lower than that in C117; the expression levels of LOC_Os05g51740, LOC_Os05g51780, LOC_Os05g51960, LOC_Os05g51970, and LOC_Os10g33210 were significantly higher in C347 than C117. Among them, LOC_Os10g33210 encodes a peptide transporter, and LOC_Os05g51690 encodes a CCT domain protein and responds to NUE in rice. This study identified new loci related to NUE in rice, providing new genetic resources for the molecular breeding of rice landraces with high NUE.

QTL mapping and identification of candidate genes for cold tolerance at the germination stage in wild rice.

BACKGROUND: Cold damage stress significantly affects rice growth (germination and seedling) and causes serious losses in yield in temperate and high-altitude areas around the globe. OBJECTIVE: This study aimed to explore the cold tolerance (CT) locus of rice and create new cold-tolerant germplasm. We constructed a chromosome segment substitution line (CSSL) with strong CT and fine mapped quantitative trait loci (QTLs) associated with CT by performing the whole-genome resequencing of CSSL with phenotypes under cold treatment. METHODS: A chromosome CSSL, including 271 lines from a cross between the cold-tolerant wild rice Y11 (Oryza rufipogon Griff.) and the cold-sensitive rice variety GH998, was developed to map QTLs conferring CT at the germination stage. The whole-genome resequencing was performed on CSSL for mapping QTLs of associated with CT at the germination stage. RESULTS: A high-density linkage map of the CSSLs was developed using the whole-genome resequencing of 1484 bins. The QTL analysis using 615,466 single-nucleotide polymorphisms (SNPs) led to the identification of 2 QTLs related to germination rate at low-temperature on chromosome 8 (qCTG-8) and chromosome 11 (qCTG-11). The qCTG-8 and qCTG-11 explained 14.55% and 14.31% of the total phenotypic variation, respectively. We narrowed down qCTG-8 and qCTG-11 to 195.5 and 78.83-kb regions, respectively. The expression patterns of important candidate genes in different tissues, and of RNA-sequencing (RNA-seq) in CSSLs, were identified based on gene sequences in qCTG-8 and qCTG-11 cold-induced expression analysis. LOC_Os08g01120 and LOC_Os08g01390 were identified as candidate genes in qCTG-8, and LOC_Os11g32880 was identified as a candidate gene in qCTG-11. CONCLUSIONS: This study demonstrated a general method that could be used to identify useful loci and genes in wild rice and aid in the future cloning of candidate genes of qCTG-8 and qCTG-11. The CSSLs with strong CT were supported for breeding cold-tolerant rice varieties.

Genome-wide selection and introgression of Chinese rice varieties during breeding.

China is the largest rice-producing country, but the genomic landscape of rice diversity has not yet been clarified. In this study, we re-sequence 1070 rice varieties collected from China (400) and other regions in Asia (670). Among the six major rice groups (aus, indica-I, indica-II, aromatic, temperate japonica, and tropical japonica), almost all Chinese varieties belong to the indica-II or temperate japonica group. Most Chinese indica varieties belong to indica-II, which consists of two subgroups developed during different phases of rice breeding. The genomic segments underlying the differences between these subgroups span 36.32 Mb. The Chinese japonica rice varieties fall into the temperate japonica group, consisting of two subgroups based on their geographical distribution. The genomic segments underlying the differences between these subgroups span 27.69 Mb. These differentiated segments in the Chinese indica varieties span 45 genes with nonsynonymous mutations that are closely related to variations in plant height and grain width. Fifty-four genes with nonsynonymous mutations are associated with the differences in heading date between the two Chinese japonica subgroups. These findings provide new insights into rice diversity in China that will facilitate the molecular breeding.

High-Quality Genomes and High-Density Genetic Map Facilitate the Identification of Genes From a Weedy Rice.

Genes have been lost or weakened from cultivated rice during rice domestication and breeding. Weedy rice (Oryza sativa f. spontanea) is usually recognized as the progeny between cultivated rice and wild rice and is also known to harbor an gene pool for rice breeding. Therefore, identifying genes from weedy rice germplasms is an important way to break the bottleneck of rice breeding. To discover genes from weedy rice germplasms, we constructed a genetic map based on w-hole-genome sequencing of a F2 population derived from the cross between LM8 and a cultivated rice variety. We further identified 31 QTLs associated with 12 important agronomic traits and revealed that ORUFILM03g000095 gene may play an important role in grain length regulation and participate in grain formation. To clarify the genomic characteristics from weedy rice germplasms of LM8, we generated a high-quality genome assembly using single-molecule sequencing, Bionano optical mapping, and Hi-C technologies. The genome harbored a total size of 375.8 Mb, a scaffold N50 of 24.1 Mb, and originated approximately 0.32 million years ago (Mya) and was more closely related to Oryza sativa ssp. japonica. and contained 672 unique genes. It is related to the formation of grain shape, heading date and tillering. This study generated a high-quality reference genome of weedy rice and high-density genetic map that would benefit the analysis of genome evolution for related species and suggested an effective way to identify genes related to important agronomic traits for further rice breeding.

OsTTG1, a WD40 repeat gene, regulates anthocyanin biosynthesis in rice.

Anthocyanins play an important role in the growth of plants, and are beneficial to human health. In plants, the MYB-bHLH-WD40 (MBW) complex activates the genes for anthocyanin biosynthesis. However, in rice, the WD40 regulators remain to be conclusively identified. Here, a crucial anthocyanin biosynthesis gene was fine mapped to a 43.4-kb genomic region on chromosome 2, and a WD40 gene OsTTG1 (Oryza sativa TRANSPARENT TESTA GLABRA1) was identified as ideal candidate gene. Subsequently, a homozygous mutant (osttg1) generated by CRISPR/Cas9 showed significantly decreased anthocyanin accumulation in various rice organs. OsTTG1 was highly expressed in various rice tissues after germination, and it was affected by light and temperature. OsTTG1 protein was localized to the nucleus, and can physically interact with Kala4, OsC1, OsDFR and Rc. Furthermore, a total of 59 hub transcription factor genes might affect rice anthocyanin biosynthesis, and LOC_Os01g28680 and LOC_Os02g32430 could have functional redundancy with OsTTG1. Phylogenetic analysis indicated that directional selection has driven the evolutionary divergence of the indica and japonica OsTTG1 alleles. Our results suggest that OsTTG1 is a vital regulator of anthocyanin biosynthesis, and an important gene resource for the genetic engineering of anthocyanin biosynthesis in rice and other plants.

A Combination of Long-Day Suppressor Genes Contributes to the Northward Expansion of Rice.

Growing cultivated rice with a moderate heading date is the key to expanding its cultivation area and maintaining stable yields. The genes that regulate heading date are largely cloned; however, it remains unclear how genetic mutations and their combinations affect the heading date and adaptability of cultivated rice. Here, we report the analysis of genetic variation in eight long-day flowering suppressor genes (Hd1, DTH8, Ghd7, OsCOL4, DTH7, Hd6, Se5, and PhyB) and the phylogenetic relationship of eight genes. Genetic variations in DTH8, Ghd7, Hd1, DTH7, PhyB, and OsCOL4 are correlated with differences in heading date and the correlation between the genetic diversity of Hd6 and Se5 and rice heading data are weak. One group of haplotypes of DTH8, Ghd7, Hd1, DTH7, PhyB, and OsCOL4 are associated with earlier heading dates and appear to have accumulated during the northward expansion of rice cultivation. A minimum of four group A alleles of DTH8, Ghd7, Hd1, DTH7, PhyB, and OsCOL4 are required for the growth of cultivated rice at latitudes above 30°N. This study presents a preliminary investigation of the genetic patterns and adaptation mechanisms of long-day flowering suppressor genes and provides a useful reference for the molecular breeding of rice cultivars for various environments and farming systems.

Genome-wide identification of the peptide transporter family in rice and analysis of the PTR expression modulation in two near-isogenic lines with different nitrogen use efficiency.

BACKGROUND: Nitrogen (N) is a major nutrient element for crop growth. In plants, the members of the peptide transporter (PTR) gene family may involve in nitrate uptake and transport. Here, we identified PTR gene family in rice and analyzed their expression profile in near-isogenic lines. RESULTS: We identified 96, 85 and 78 PTR genes in Nipponbare, R498 and Oryza glaberrima, and the phylogenetic trees were similar in Asian cultivated rice and African cultivated rice. The number of PTR genes was higher in peanut (125) and soybean (127). The 521 PTR genes in rice, maize, sorghum, peanut, soybean and Arabidopsis could be classified into 4 groups, and their distribution was different between monocots and dicots. In Nipponbare genome, the 25 PTR genes were distributed in 5 segmental duplication regions on chromosome 1, 2, 3, 4, 5, 7, 8, 9, and 10. The PTR genes in rice have 0-11 introns and 1-12 exons, and 16 of them have the NPF (NRT1/PTR family) domain. The results of RNA-seq showed that the number of differentially expressed genes (DEGs) between NIL15 and NIL19 at three stages were 928, 1467, and 1586, respectively. Under low N conditions, the number of differentially expressed PTR genes increased significantly. The RNA-seq data was analyzed using WGCNA to predict the potential interaction between genes. We classified the genes with similar expression pattern into one module, and obtained 25 target modules. Among these modules, three modules may be involved in rice N uptake and utilization, especially the brown module, in which hub genes were annotated as protein kinase that may regulate rice N metabolism. CONCLUSIONS: In this study, we comprehensively analyzed the PTR gene family in rice. 96 PTR genes were identified in Nippobare genome and 25 of them were located on five large segmental duplication regions. The Ka/Ks ratio indicated that many PTR genes had undergone positive selection. The RNA-seq results showed that many PTR genes were involved in rice nitrogen use efficiency (NUE), and protein kinases might play an important role in this process. These results provide a fundamental basis to improve the rice NUE via molecular breeding.

Identification of candidate genes for gelatinization temperature, gel consistency and pericarp color by GWAS in rice based on SLAF-sequencing.

Rice is an important cereal in the world. The study of the genetic basis of important agronomic traits in rice landraces and identification of genes will facilitate the breed improvement. Gelatinization temperature (GT), gel consistency (GC) and pericarp color (PC) are important indices of rice cooking and eating quality evaluation and potential nutritional importance, which attract wide attentions in the application of genetic and breeding. To dissect the genetic basis of GT, GC and PC, a total of 419 rice landraces core germplasm collections consisting of 330 indica lines, 78 japonica lines and 11 uncertain varieties were planted, collected, then GT, GC, PC were measured for two years, and sequenced using specific-locus amplified fragment sequencing (SLAF-seq) technology. In this study, 261,385,070 clean reads and 56,768 polymorphic SLAF tags were obtained, which a total of 211,818 single nucleotide polymorphisms (SNPs) were discovered. With 208,993 SNPs meeting the criterion of minor allele frequency (MAF) > 0.05 and integrity> 0.5, the phylogenetic tree and population structure analysis were performed for all 419 rice landraces, and the whole panel mainly separated into six subpopulations based on population structure analysis. Genome-wide association study (GWAS) was carried out for the whole panel, indica subpanel and japonica subpanel with subset SNPs respectively. One quantitative trait locus (QTL) on chromosome 6 for GT was detected in the whole panel and indica subpanel, and one QTL associated with GC was located on chromosome 6 in the whole panel and indica subpanel. For the PC trait, 8 QTLs were detected in the whole panel on chromosome 1, 3, 4, 7, 8, 10 and 11, and 7 QTLs in the indica subpanel on chromosome 3, 4, 7, 8, 10 and 11. For the three traits, no QTL was detected in japonica subpanel, probably because of the polymorphism repartition between the subpanel, or small population size of japonica subpanel. This paper provides new gene resources and insights into the molecular mechanisms of important agricultural trait of rice phenotypic variation and genetic improvement of rice quality variety breeding.

QTL Mapping by Whole Genome Re-sequencing and Analysis of Candidate Genes for Nitrogen Use Efficiency in Rice.

Nitrogen is a major nutritional element in rice production. However, excessive application of nitrogen fertilizer has caused severe environmental pollution. Therefore, development of rice varieties with improved nitrogen use efficiency (NUE) is urgent for sustainable agriculture. In this study, bulked segregant analysis (BSA) combined with whole genome re-sequencing (WGS) technology was applied to finely map quantitative trait loci (QTL) for NUE. A key QTL, designated as qNUE6 was identified on chromosome 6 and further validated by Insertion/Deletion (InDel) marker-based substitutional mapping in recombinants from F2 population (NIL-13B4 × GH998). Forty-four genes were identified in this 266.5-kb region. According to detection and annotation analysis of variation sites, 39 genes with large-effect single-nucleotide polymorphisms (SNPs) and large-effect InDels were selected as candidates and their expression levels were analyzed by qRT-PCR. Significant differences in the expression levels of LOC_Os06g15370 (peptide transporter PTR2) and LOC_Os06g15420 (asparagine synthetase) were observed between two parents (Y11 and GH998). Phylogenetic analysis in Arabidopsis thaliana identified two closely related homologs, AT1G68570 (AtNPF3.1) and AT5G65010 (ASN2), which share 72.3 and 87.5% amino acid similarity with LOC_Os06g15370 and LOC_Os06g15420, respectively. Taken together, our results suggested that qNUE6 is a possible candidate gene for NUE in rice. The fine mapping and candidate gene analysis of qNUE6 provide the basis of molecular breeding for genetic improvement of rice varieties with high NUE, and lay the foundation for further cloning and functional analysis.

Rapid identification of a new gene influencing low amylose content in rice landraces (Oryza sativa L.) using genome-wide association study with specific-locus amplified fragment sequencing.

Starch is the major component of milled rice, and amylose content (AC) affects eating quality. In this study, a genome-wide association study (GWAS) with specific-locus amplified fragment sequencing (SLAF-seq) data was performed for AC on a core collection of 419 rice landraces. Using the compressed mixed linear method based on the Q+K model, we identified a new gene, LAC6 (Chr6: 5.65-6.04 Mb), associated with AC in the low amylose content panel. The LAC6 candidate gene was detected by qRT-PCR in rice panicles. Results indicate that LOC_Os06g11130 was up-regulated, and LOC_Os06g11340 was significantly down-regulated, making it most likely a key candidate gene of LAC6. In conclusion, the findings provide a certain theoretiacal basis of molecular biology for genetic improvement of AC in rice and rice quality variety breeding.