Deep sequencing transcriptional fingerprinting of rice kernels for dissecting grain quality traits

RAL6 encodes a seed allergenic protein that positively regulates grain weight and seed germination

The rice albumin (RAG) gene family belongs to the Tryp_alpha_amyl family. RAG2, specifically expressed in 14-21 DAP (days after pollination) seeds, regulates grain yield and quality. In this study, we identified another RAG family gene, RAL6, which exhibits specific expression in developing seeds, particularly in 7, 10, and 15 DAP seeds. Employing the CRISPR/Cas9 system, we analyzed functions of RAL6 and found that the ral6 lines (ral6-1, ral6-2, ral6-3, and ral6-4) displayed thinner seeds with significantly decreased 1000-grain weight and grain thickness compared to ZH11. Additionally, the cell width of spikelet cells, total protein and glutelin contents were significantly reduced in ral6. The germination assay and 1% TTC staining revealed a significant decrease in seed vigor among the ral6 lines. The alpha-amylase activity in ral6 mutant seeds was also markedly lower than in ZH11 seeds after 2 days of imbibition. Furthermore, co-expression analysis and GO annotation showed that co-expressed genes were involved in immune response, oligopeptide transport, and the glucan biosynthetic process. Collectively, our findings suggest that RAL6 plays a coordinating role in regulating grain weight and seed germination in rice.

Application of CRISPR/Cas9 system to knock out GluB gene for developing low glutelin rice mutant

The nutritional quality improvement is among the most integral objective for any rice molecular breeding programs. The seed storage proteins (SSPs) have greater role to determine the nutritional quality of any cereal grains. Rice contains relatively balanced amino acid composition and the SSPs are fractioned into albumins (ALB), globulins (GLO), prolamins (PRO) and glutelins (GLU) according to differences in solubility. GLUs are further divided into subfamilies: GluA, GluB, GluC, and GluD depending on resemblance in amino acid. The GLU protein accounts for 60-80% of total protein contents, encoded by 15 genes located on different chromosomes of rice genome. The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) system was employed to knockout Glu-B (LOC-Os02g15070) gene in non-basmati rice PK386 cultivar. The mutant displayed two base pair and three base pair mutation in the targeted regions. The homozygous mutant plant displayed reduction for both in total protein contents and GLU contents whereas, elevation in GLO, ALB and PRO. Moreover, the mutant plant also displayed reduction in physio-chemical properties e.g., total starch, amylose and gel consistency. The agronomic characteristics of both mutant and wild type displayed non-significant differences along with increase in higher percentage of chalkiness in mutant plants. The results obtained from scanning electron microscopy showed the loosely packed starch granules compared to wild type. The gene expression analysis displayed the lower expression of gene at 5 days after flowering (DAF), 10 DAF, 15 DAF and 20 DAF compared to wild type. GUS sub-cellular localization showed the staining in seed which further validated the results obtained from gene expression. Based on these findings it can be concluded Glu-B gene have significant role in controlling GLU contents and can be utilized in breeding programs to enhance the nutritional quality of rice, and may serve as healthy diet for patient allergic with high GLU contents.

Multi-model genome-wide association studies for appearance quality in rice

Improving the quality of the appearance of rice is critical to meet market acceptance. Mining putative quality-related genes has been geared towards the development of effective breeding approaches for rice. In the present study, two SL-GWAS (CMLM and MLM) and three ML-GWAS (FASTmrEMMA, mrMLM, and FASTmrMLM) genome-wide association studies were conducted in a subset of 3K-RGP consisting of 198 rice accessions with 553,831 SNP markers. A total of 594 SNP markers were identified using the mixed linear model method for grain quality traits. Additionally, 70 quantitative trait nucleotides (QTNs) detected by the ML-GWAS models were strongly associated with grain aroma (AR), head rice recovery (HRR, %), and percentage of grains with chalkiness (PGC, %). Finally, 39 QTNs were identified using single- and multi-locus GWAS methods. Among the 39 reliable QTNs, 20 novel QTNs were identified for the above-mentioned three quality-related traits. Based on annotation and previous studies, four functional candidate genes (LOC_Os01g66110, LOC_Os01g66140, LOC_Os07g44910, and LOC_Os02g14120) were found to influence AR, HRR (%), and PGC (%), which could be utilized in rice breeding to improve grain quality traits.

Physicochemical properties and molecular mechanisms of different resistant starch subtypes in rice

Resistant starch (RS) can help prevent diabetes and decrease calorie intake and that from plants are the main source of mankind consumption. Rice is many people's staple food and that with higher RS will help health management. A significantly positive correlation exists between apparent amylose content (AAC) of rice and its RS content. In this study, 72 accessions with moderate or high AAC were selected to explore the regulatory mechanisms and physicochemical properties on different proceeding types of rice RS. RS in raw milled rice (RSm), hot cooked rice (RSc), and retrogradation rice (RSr) showed a wide variation and distinct controlling mechanisms. They were co-regulated by Waxy (Wx), soluble starch synthase (SS) IIb and SSI. Besides that, RSm was also regulated by SSIIa and SSIVb, RSc by granule-bound starch synthase (GBSS) II and RSr by GBSSII and Pullulanase (PUL). Moreover, Wx had significant interactions with SSIIa, SSI, SSIIb and SSIVb on RSm, but only the dominant interactions with SSIIb and SSI on RSc and RSr. Wx was the key factor for the formation of RS, especially the RSc and RSr. The genes had the highest expression at 17 days after flowering and were beneficial for RS formation. The longer the chain length of starch, the higher the RS3 content. RSc and RSr were likely to be contained in medium-size starch granules. The findings favor understanding the biosynthesis of different subtypes of RS.

Genome-Wide Association Study of Seed Morphology Traits in Senegalese Sorghum Cultivars

Sorghum is considered the fifth most important crop in the world. Despite the potential value of Senegalese germplasm for various traits, such as resistance to fungal diseases, there is limited information on the study of sorghum seed morphology. In this study, 162 Senegalese germplasms were evaluated for seed area size, length, width, length-to-width ratio, perimeter, circularity, the distance between the intersection of length & width (IS) and center of gravity (CG), and seed darkness and brightness by scanning and analyzing morphology-related traits with SmartGrain software at the USDA-ARS Plant Science Research Unit. Correlations between seed morphology-related traits and traits associated with anthracnose and head smut resistance were analyzed. Lastly, genome-wide association studies were performed on phenotypic data collected from over 16,000 seeds and 193,727 publicly available single nucleotide polymorphisms (SNPs). Several significant SNPs were found and mapped to the reference sorghum genome to uncover multiple candidate genes potentially associated with seed morphology. The results indicate clear correlations among seed morphology-related traits and potential associations between seed morphology and the defense response of sorghum. GWAS analysis listed candidate genes associated with seed morphologies that can be used for sorghum breeding in the future.

Rice Labeling according to Grain Quality Features Using Laser-Induced Breakdown Spectroscopy

Rice is an important source of nutrition and energy consumed around the world. Thus, quality inspection is crucial for protecting consumers and increasing the rice's value in the productive chain. Currently, methods for rice labeling depending on grain quality features are based on image and/or visual inspection. These methods have shown subjectivity and inefficiency for large-scale analyses. Laser-induced breakdown spectroscopy (LIBS) is an analytical technique showing attractive features due to how quick the analysis can be carried out and its capability of providing spectra that are true fingerprints of the sample's elemental composition. In this work, LIBS performance was evaluated for labeling rice according to grain quality features. The LIBS spectra of samples with their grain quality numerically described as Type 1, 2, and 3 were measured. Several spectral processing methods were evaluated when modeling a k-nearest neighbors (k-NN) classifier. Variable selection was also carried out by principal component analysis (PCA), and then the optimal k-value was selected. The best result was obtained by applying spectrum smoothing followed by normalization by using the first fifteen principal components (PCs) as input variables and k = 9. Under these conditions, the method showed excellent performance, achieving sample classification with 94% overall prediction accuracy. The sensitivities ranged from 90 to 100%, and specificities were in the range of 92-100%. The proposed method has remarkable characteristics, e.g., analytical speed and analysis guided by chemical responses; therefore, the method is not susceptible to subjectivity errors.

Efficient breeding of low glutelin content rice germplasm by simultaneous editing multiple glutelin genes via CRISPR/Cas9

Chronic kidney disease (CKD) and phenylketonuria (PKU) patients need to eat rice with low glutelin content. Therefore, breeding low glutelin content rice varieties with high yield and delicious taste is one of the major goals of rice breeders due to the high demand for the product. In this study, we designed three sgRNAs targeting nine glutelin genes and generated nine T-DNA-free homozygous editing lines with reduced glutelin content compared with the wild-type due to simultaneous mutation(s) in 5-7 glutelin genes. The glutelin content of two lines is even significantly lower than that of the low glutelin content cultivar, LGC-1. Compared to the wild-type, these low glutelin lines showed similar agronomic traits, including yield components and viscosity properties, and can be used as new varieties or parental materials for further breeding.

Identification and Pleiotropic Effect Analysis of GSE5 on Rice Chalkiness and Grain Shape

Chalkiness is one of several major restricting factors for the improvement of rice quality. Although many chalkiness-related quantitative trait loci have been mapped, only a small number of genes have been cloned to date. In this study, the candidate gene GSE5 of a major quantitative trait locus (QTL) for rice chalkiness, qDEC5, was identified by map-based cloning. Phenotyping and haplotype analysis of proActin:GSE5 transgenic line, gse5-cr mutant, and 69 rice varieties further confirmed that GSE5 had the pleiotropic effects and regulated both chalkiness and grain shape. Genetic analysis showed GSE5 was a dominant gene for grain length and a semi-dominant gene for grain width and chalkiness. The DNA interval closely linked to GSE5 was introgressed to Zhenshan 97B (ZB) based on molecular marker-assisted selection, and the improved ZB showed lower chalkiness and longer but smaller grains, which showed that GSE5 played an important role in breeding rice varieties with high yield and good quality. Transcriptomics, proteomics, and qRT-PCR analyses showed that thirty-nine genes associated with carbon and protein metabolism are regulated by GSE5 to affect the formation of chalkiness, including some newly discovered genes, such as OsCESA9, OsHSP70, OsTPS8, OsPFK04, OsSTA1, OsERdj3A, etc. The low-chalkiness lines showed higher amino sugar and nucleotide sugar metabolism at 10 days after pollination (DAP), lower carbohydrate metabolism at 15 DAP, and lower protein metabolism at 10 and 15 DAP. With heat shock at 34/30°C, rice chalkiness increased significantly; OsDjC10 and OsSUS3 were upregulated at 6 and 12 DAP, respectively, and OsGSTL2 was downregulated at 12 DAP. Our results identified the function and pleiotropic effects of qDEC5 dissected its genetic characteristics and the expression profiles of the genes affecting the chalkiness formation, and provided a theoretical basis and application value to harmoniously pursue high yield and good quality in rice production.

Transcriptomic and biochemical investigations support the role of rootstock-scion interaction in grapevine berry quality

Background

In viticulture, rootstock genotype plays a critical role to improve scion physiology, berry quality and to adapt grapevine (Vitis vinifera L.) to different environmental conditions. This study aimed at investigating the effect of two different rootstocks (1103 Paulsen - P - and Mgt 101–14 - M) in comparison with not grafted plants - NGC - on transcriptome (RNA-seq and small RNA-seq) and chemical composition of berry skin in Pinot noir, and exploring the influence of rootstock-scion interaction on grape quality. Berry samples, collected at veraison and maturity, were investigated at transcriptional and biochemical levels to depict the impact of rootstock on berry maturation.

Results

RNA- and miRNA-seq analyses highlighted that, at veraison, the transcriptomes of the berry skin are extremely similar, while variations associated with the different rootstocks become evident at maturity, suggesting a greater diversification at transcriptional level towards the end of the ripening process. In the experimental design, resembling standard agronomic growth conditions, the vines grafted on the two different rootstocks do not show a high degree of diversity. In general, the few genes differentially expressed at veraison were linked to photosynthesis, putatively because of a ripening delay in not grafted vines, while at maturity the differentially expressed genes were mainly involved in the synthesis and transport of phenylpropanoids (e.g. flavonoids), cell wall loosening, and stress response. These results were supported by some differences in berry phenolic composition detected between grafted and not grafted plants, in particular in resveratrol derivatives accumulation.

Conclusions

Transcriptomic and biochemical data demonstrate a stronger impact of 1103 Paulsen rootstock than Mgt 101–14 or not grafted plants on ripening processes related to the secondary metabolite accumulations in berry skin tissue. Interestingly, the MYB14 gene, involved in the feedback regulation of resveratrol biosynthesis was up-regulated in 1103 Paulsen thus supporting a putative greater accumulation of stilbenes in mature berries.

Grapevine comparative early transcriptomic profiling suggests that Flavescence dorée phytoplasma represses plant responses induced by vector feeding in susceptible varieties

Background

Flavescence dorée is the most serious grapevine yellows disease in Europe. It is caused by phytoplasmas which are transmitted from grapevine to grapevine by the leafhopper Scaphoideus titanus. Differences in susceptibility among grapevine varieties suggest the existence of specific genetic features associated with resistance to the phytoplasma and/or possibly with its vector. In this work, RNA-Seq was used to compare early transcriptional changes occurring during the three-trophic interaction between the phytoplasma, its vector and the grapevine, represented by two different cultivars, one very susceptible to the disease and the other scarcely susceptible.

Results

The comparative analysis of the constitutive transcriptomic profiles suggests the existence of passive defense strategies against the insect and/or the phytoplasma in the scarcely-susceptible cultivar. Moreover, the attack by the infective vector on the scarcely-susceptible variety prompted immediate and substantial transcriptomic changes that led to the rapid erection of further active defenses. On the other hand, in the most susceptible variety the response was delayed and mainly consisted of the induction of phytoalexin synthesis. Surprisingly, the jasmonic acid- and ethylene-mediated defense reactions, activated by the susceptible cultivar following FD-free insect feeding, were not detected in the presence of the phytoplasma-infected vector.

Conclusions

The comparison of the transcriptomic response in two grapevine varieties with different levels of susceptibility to Flavescence dorèe highlighted both passive and active defense mechanisms against the vector and/or the pathogen in the scarcely-susceptible variety, as well as the capacity of the phytoplasmas to repress the defense reaction against the insect in the susceptible variety.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5908-6) contains supplementary material, which is available to authorized users.

Genome-wide association study reveals significant genomic regions for improving yield, adaptability of rice under dry direct seeded cultivation condition

BACKGROUND

Puddled transplanted system of rice cultivation despite having several benefits, is a highly labor, water and energy intensive system. In the face of changing climatic conditions, a successful transition from puddled to dry direct seeded rice (DDSR) cultivation system looks must in future. Genome-wide association study was performed for traits including, roots and nutrient uptake (14 traits), plant-morphological (5 traits), lodging-resistance (4 traits) and yield and yield attributing traits (7 traits) with the aim to identify significant marker-trait associations (MTAs) for traits enhancing rice adaptability to dry direct-seeded rice (DDSR) system.

RESULTS

Study identified a total of 37 highly significant MTAs for 20 traits. The false discovery rate (FDR) ranged from 0.264 to 3.69 × 10- 4, 0.0330 to 1.25 × 10- 4, and 0.0534 to 4.60 × 10- 6 in 2015WS, 2016DS and combined analysis, respectively. The percent phenotypic variance (PV) explained by SNPs ranged from 9 to 92%. Among the identified significant MTAs, 15 MTAs associated with the traits including nodal root, root hair length, root length density, stem and culm diameter, plant height and grain yield were reported to be located in the proximity of earlier identified candidate gene. The significant positive correlation of grain-yield with seedling establishment traits, root morphological and nutrient-uptake related traits and grain yield attributing traits pointing towards combining target traits to increase rice yield and adaptability under DDSR. Seven promising progenies with better root morphology, nutrient-uptake and higher grain yield were identified that can further be used in genomics assisted breeding for DDSR varietal development.

CONCLUSIONS

Once validated, the identified MTAs and the SNPs linked with trait of interest could be of direct use in genomic assisted breeding (GAB) to improve grain yield and adaptability of rice under DDSR.

Comparative transcriptome profiling of the response to Pyrenochaeta lycopersici in resistant tomato cultivar Mogeor and its background genotype-susceptible Moneymaker

Breeding for resistance is the most effective tool for controlling the corky root disease of tomato caused by Pyrenochaeta lycopersici. A comparative RNA-Seq-based transcriptomic analysis was conducted at 96 hpi (hours post infection) on two tomato cultivars: resistant Mogeor and its genetic background, and susceptible Moneymaker to investigate the differences in their transcriptomic response and identify the molecular bases of this plant-pathogen interaction. The number of differentially expressed genes (DEGs) identified was much higher in the susceptible than in the resistant genotype; however, the proportion of upregulated genes was higher in Mogeor (70.81%) than in Moneymaker (52.95%). Gene Ontology (GO) analysis enabled identification of 24 terms shared by the two cultivars that were consistent with responses to external stimulus, such as fungal infection. On the other hand, as many as 54 GO were enriched solely in Moneymaker, including terms related to defense response and cell wall metabolism. Our results could support the previous observations in other pathosystems, that susceptibility and resistance have overlapping signaling pathways and responses, suggesting that the P. lycopersici resistance gene pyl might be a recessive allele at a susceptibility locus, for which different candidate genes were identified based on the differences in induction or expression levels, observed between the resistant and susceptible genotype. MapMan analysis highlighted a complex hormone and transcription factors interplay where SA- and JA-induced pathways are modulated in a similar way in both genotypes and thus take part in a common response while the ethylene signaling pathways, induced mainly in susceptible Moneymaker, seem putatively contribute to its susceptibility.

Transcriptome analysis of Oryza sativa in responses to different concentrations of thiocyanate

Effective concentrations of potassium thiocyanate (KSCN) to rice seedlings were experimentally determined using relative growth rate as a sensitive endpoint. Agilent 44-K rice microarray was used to profile the molecular responses of rice seedlings exposed to thiocyanate ion (SCN^-) at three different effective concentrations (EC₁₀, EC₂₀, and EC₅₀). A total of 18,498 known genes were collected from SCN-treated rice microarray analysis. Out of all, 1603, 1882, and 5085 differentially expressed genes (DEGs) were observed at EC₁₀, EC₂₀, and EC₅₀ concentrations, respectively. More upregulated/downregulated DEGs were detected in shoots than in roots after SCN^- exposure. Gene functions and pathway enrichment analysis of DEGs indicated that different effective concentrations of SCN^- resulted in multiple enriched GO categories and KEGG pathways and outcomes were quite tissue-specific. Different regulations and adaptations of gene expression in molecular function (MF), biological process (BP), and cellular components (CC) were observed in rice tissues at different effective concentrations of SCN^-, suggesting their different responsive and adaptive strategies. Information collected here presents a detailed description of SCN-induced alternations of gene expression in rice seedlings and provide valuable information for further searching specific genes participating in transportation, phytotoxic responses, and detoxification of SCN^- in rice seedlings.

Rice grain nutritional traits and their enhancement using relevant genes and QTLs through advanced approaches

BACKGROUND

Rice breeding program needs to focus on development of nutrient dense rice for value addition and helping in reducing malnutrition. Mineral and vitamin deficiency related problems are common in the majority of the population and more specific to developing countries as their staple food is rice.

RESULTS

Genes and QTLs are recently known for the nutritional quality of rice. By comprehensive literature survey and public domain database, we provided a critical review on nutritional aspects like grain protein and amino acid content, vitamins and minerals, glycemic index value, phenolic and flavonoid compounds, phytic acid, zinc and iron content along with QTLs linked to these traits. In addition, achievements through transgenic and advanced genomic approaches have been discussed. The information available on genes and/or QTLs involved in enhancement of micronutrient element and amino acids are summarized with graphical representation.

CONCLUSION

Compatible QTLs/genes may be combined together to design a desirable genotype with superior in multiple grain quality traits. The comprehensive review will be helpful to develop nutrient dense rice cultivars by integrating molecular markers and transgenic assisted breeding approaches with classical breeding.

Genome-Wide Association Study for Traits Related to Plant and Grain Morphology, and Root Architecture in Temperate Rice Accessions

BACKGROUND

In this study we carried out a genome-wide association analysis for plant and grain morphology and root architecture in a unique panel of temperate rice accessions adapted to European pedo-climatic conditions. This is the first study to assess the association of selected phenotypic traits to specific genomic regions in the narrow genetic pool of temperate japonica. A set of 391 rice accessions were GBS-genotyped yielding-after data editing-57000 polymorphic and informative SNPS, among which 54% were in genic regions.

RESULTS

In total, 42 significant genotype-phenotype associations were detected: 21 for plant morphology traits, 11 for grain quality traits, 10 for root architecture traits. The FDR of detected associations ranged from 3 · 10-7 to 0.92 (median: 0.25). In most cases, the significant detected associations co-localised with QTLs and candidate genes controlling the phenotypic variation of single or multiple traits. The most significant associations were those for flag leaf width on chromosome 4 (FDR = 3 · 10-7) and for plant height on chromosome 6 (FDR = 0.011).

CONCLUSIONS

We demonstrate the effectiveness and resolution of the developed platform for high-throughput phenotyping, genotyping and GWAS in detecting major QTLs for relevant traits in rice. We identified strong associations that may be used for selection in temperate irrigated rice breeding: e.g. associations for flag leaf width, plant height, root volume and length, grain length, grain width and their ratio. Our findings pave the way to successfully exploit the narrow genetic pool of European temperate rice and to pinpoint the most relevant genetic components contributing to the adaptability and high yield of this germplasm. The generated data could be of direct use in genomic-assisted breeding strategies.