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.

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.

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 and evaluation of reference genes for quantitative real-time PCR analysis in Passiflora edulis under stem rot condition.

Passion fruit (Passiflora edulis), an important tropical and subtropical fruit, has a high edible and medicinal value. Stem rot disease is one of the most important diseases of passion fruit. An effective way for control and prevention of this disease is to identify the genes associated with resistance to this disease. Quantitative real-time PCR (RT-qPCR) has mainly been widely applied to detect gene expression because of its simplicity, fastness, low cost and high sensitivity. One of the requirements for RT-qPCR is the availability of suitable reference genes for normalization of gene expression. However, currently, no Passiflora edulis reference genes have been identified andthus it has hindered the gene expression studies in this plant. The present study aimed to address this issue. We analyzed sixteen candidate reference genes, including nine common (GAPDH, UBQ, ACT1, ACT2, EF-1α-1, EF-1α-2, TUA, NADP, and GBP) and seven novel genes (C13615, C24590, C27182, C10445, C21209, C22199, and C22526), in different tissues (stem, leaf, flower and fruit) of two accessions under stem rot condition. We calculated the expression stability in twenty-four samples using the ΔCt, GeNorm, NormFinder, BestKeeper and RefFinder. The results showed that both C21209 and EF-1α-2 were sufficient to normalize gene expression under stem rot, whereas the commonly used reference genes, GAPDH and UBQ, were the least stable ones. The expression patterns of PeUFC under stem rot condition normalized by stable and unstable reference genes indicated the suitability of using the optimal reference genes. To our knowledge, this is the first systematic study of reference genes in Passiflora edulis, which identified a number of reliable reference genes suitable for gene expression studies in Passiflora edulis by RT-qPCR.

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.

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.

Assessment of Rice Sheath Blight Resistance Including Associations with Plant Architecture, as Revealed by Genome-Wide Association Studies.

BACKGROUND: Sheath blight (ShB) disease caused by Rhizoctonia solani Kühn, is one of the most economically damaging rice (Oryza sativa L.) diseases worldwide. There are no known major resistance genes, leaving only partial resistance from small-effect QTL to deploy for cultivar improvement. Many ShB-QTL are associated with plant architectural traits detrimental to yield, including tall plants, late maturity, or open canopy from few or procumbent tillers, which confound detection of physiological resistance. RESULTS: To identify QTL for ShB resistance, 417 accessions from the Rice Diversity Panel 1 (RDP1), developed for association mapping studies, were evaluated for ShB resistance, plant height and days to heading in inoculated field plots in Arkansas, USA (AR) and Nanning, China (NC). Inoculated greenhouse-grown plants were used to evaluate ShB using a seedling-stage method to eliminate effects from height or maturity, and tiller (TN) and panicle number (PN) per plant. Potted plants were used to evaluate the RDP1 for TN and PN. Genome-wide association (GWA) mapping with over 3.4 million SNPs identified 21 targeted SNP markers associated with ShB which tagged 18 ShB-QTL not associated with undesirable plant architecture traits. Ten SNPs were associated with ShB among accessions of the Indica subspecies, ten among Japonica subspecies accessions, and one among all RDP1 accessions. Across the 18 ShB QTL, only qShB4-1 was not previously reported in biparental mapping studies and qShB9 was not reported in the GWA ShB studies. All 14 PN QTL overlapped with TN QTL, with 15 total TN QTL identified. Allele effects at the five TN QTL co-located with ShB QTL indicated that increased TN does not inevitably increase disease development; in fact, for four ShB QTL that overlapped TN QTL, the alleles increasing resistance were associated with increased TN and PN, suggesting a desirable coupling of alleles at linked genes. CONCLUSIONS: Nineteen accessions identified as containing the most SNP alleles associated with ShB resistance for each subpopulation were resistant in both AR and NC field trials. Rice breeders can utilize these accessions and SNPs to develop cultivars with enhanced ShB resistance along with increased TN and PN for improved yield potential.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.).

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F2 population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways' biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.

Comparative transcriptomic analysis reveals the cold acclimation during chilling stress in sensitive and resistant passion fruit (Passiflora edulis) cultivars.

Chilling stress (CS) is an important limiting factor for the growth and development of passion fruit (Passiflora edulis) in winter in South China. However, little is known about how the passion fruit responds and adapts to CS. In this study, we performed transcriptome sequencing of cold-susceptible cultivar Huangjinguo (HJG) and cold-tolerant cultivar Tainong 1 (TN1) under normal temperature (NT) and CS conditions, and a total of 47,353 unigenes were obtained by seven databases. Using differentially expressed unigenes (DEGs) analysis, 3,248 and 4,340 DEGs were identified at two stages, respectively. The Gene Ontology (GO) enrichment analysis showed that the DEGs were mainly related to phosphorylation, membrane protein, and catalytic activity. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, the unigenes of plant-pathogen interaction, plant hormone signal transduction and fatty acid metabolism were enriched. Then, the 12,471 filtered unigenes were clustered into eight co-expression modules, and two modules were correlated with CS. In this two modules, 32 hub unigenes were obtained. Furthermore, the unigenes related to CS were validated using quantitative real-time PCR (RT-qPCR). This work showed that the expression levels of CS-related unigenes were very different in two passion fruit cultivars. The results provide information for the development of passion fruit with increased chilling tolerance.

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.