Whole genome sequencing and analysis of Swarna, a widely cultivated indica rice variety with low glycemic index

Endosperm structure and Glycemic Index of Japonica Italian rice varieties

Introduction

Given that rice serves as a crucial staple food for a significant portion of the global population and with the increasing number of individuals being diagnosed with diabetes, a primary objective in genetic improvement is to identify and cultivate low Glycemic Index (GI) varieties. This must be done while ensuring the preservation of grain quality.

Methods

25 Italian rice genotypes were characterized calculating their GI "in vivo" and, together with other 29 Italian and non-Italian genotypes they were studied to evaluate the grain inner structure through Field Emission Scanning Electron Microscopy (FESEM) technique. Using an ad-hoc developed algorithm, morphological features were extracted from the FESEM images, to be then inspected by means of multivariate data analysis methods.

Results and Discussion

Large variability was observed in GI values (49 to 92 with respect to glucose), as well as in endosperm morphological features. According to the percentage of porosity is possible to distinguish approximately among rice varieties having a crystalline grain (< 1.7%), those intended for the preparation of risotto (> 5%), and a third group having intermediate characteristics. Waxy rice varieties were not united by a certain porosity level, but they shared a low starch granules eccentricity. With reference to morphological features, rice varieties with low GI (<55) seem to be characterized by large starch granules and low porosity values. Our data testify the wide variability of Italian rice cultivation giving interesting information for future breeding programs, finding that the structure of the endosperm can be regarded as a specific characteristic of each variety.

Whole-Genome Sequencing of KMR3 and Oryza rufipogon-Derived Introgression Line IL50-13 (Chinsurah Nona 2/Gosaba 6) Identifies Candidate Genes for High Yield and Salinity Tolerance in Rice

The genomes of an elite rice restorer line KMR3 (salinity-sensitive) and its salinity-tolerant introgression line IL50-13, a popular variety of coastal West Bengal, India, were sequenced. High-quality paired-end reads were obtained for KMR3 (147.6 million) and IL50-13 (131.4 million) with a sequencing coverage of 30X-39X. Scaffolds generated from the pre-assembled contigs of each sequenced genome were mapped separately onto the reference genome of Oryza sativa ssp. japonica cultivar Nipponbare to identify genomic variants in terms of SNPs and InDels. The SNPs and InDels identified for KMR3 and IL50-13 were then compared with each other to identify polymorphic SNPs and InDels unique and common to both the genomes. Functional enrichment analysis of the protein-coding genes with unique InDels identified GO terms involved in protein modification, ubiquitination, deubiquitination, peroxidase activity, and antioxidant activity in IL50-13. Linoleic acid metabolism, circadian rhythm, and alpha-linolenic acid metabolism pathways were enriched in IL50-13. These GO terms and pathways are involved in reducing oxidative damage, thus suggesting their role in stress responses. Sequence analysis of QTL markers or genes known to be associated with grain yield and salinity tolerance showed polymorphism in 20 genes, out of which nine were not previously reported. These candidate genes encoded Nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4 (NB-ARC) domain-containing protein, cyclase, receptor-like kinase, topoisomerase II-associated protein PAT1 domain-containing protein, ion channel regulatory protein, UNC-93 domain-containing protein, subunit A of the heteromeric ATP-citrate lyase, and three conserved hypothetical genes. Polymorphism was observed in the coding, intron, and untranslated regions of the genes on chromosomes 1, 2, 4, 7, 11, and 12. Genes showing polymorphism between the two genomes were considered as sequence-based new candidates derived from Oryza rufipogon for conferring high yield and salinity tolerance in IL50-13 for further functional studies.

Genetic diversity and population structure of rice (Oryza sativa L.) landraces from Kerala, India analyzed through genotyping-by-sequencing

Researchers stand at the vanguard of advancement and application of next-generation sequencing technology for developing dominant strategies for the sustainable management of genetically diverse crops. We attempt to fill the existing research lacuna in the molecular characterization of potent rice landraces in Kerala. Genotyping-by-sequencing (GBS) was performed on 96 Kerala rice accessions to identify single-nucleotide polymorphisms (SNPs), to examine the genetic diversity, population structure, and to delineate linkage disequilibrium (LD) pattern. GBS identified 5856 high-quality SNPs. The structure analysis indicated three subpopulations with the highest probability for population clustering with significant genetic differentiation, confirmed by principal component analysis. The genome-wide LD decay distance was 772 kb, at which the r² dropped to half its maximum value. The analysis of genetic properties of the identified SNP panel with an average polymorphism information content (PIC) value of 0.22 and a minor allele frequency (MAF) > 0.1 unveiled their efficacy in genome-wide association studies (GWAS). High F_ST (0.266) and low Nm (0.692) portray a strong genetic differentiation among the rice landraces, complementing the genetic structuring observed in the studied population. Slow LD decay in the rice landraces reflects their self-pollinating behavior and the indirect selection of desired traits by domestication. Moreover, the high LD entails only a minimum number of SNP markers for detecting marker-trait association. The diverse germplasm utilized in this study can be further utilized to disclose genetic variants associated with phenotypic traits and define signatures of selection via GWAS and selective sweep, respectively.

Gene Expression Dynamics in Rice Peduncles at the Heading Stage

Improving grain yield in the staple food crop rice has been long sought goal of plant biotechnology. One of the traits with significant impact on rice breeding programs is peduncle elongation at the time of heading failing which leads to significant reduction in grain yield due to incomplete panicle exsertion. To decipher transcriptional dynamics and molecular players underlying peduncle elongation, we performed RNA sequencing analysis of elongating and non-elongating peduncles in two Indian cultivars, Swarna and Pokkali, at the time of heading. Along with genes associated with cell division and cell wall biosynthesis, we observed significant enrichment of genes associated with auxins, gibberellins, and brassinosteroid biosynthesis/signaling in the elongating peduncles before heading in both the genotypes. Similarly, genes associated with carbohydrate metabolism and mobilization, abiotic stress response along with cytokinin, abscisic acid, jasmonic acid, and ethylene biosynthesis/signaling were enriched in non-elongating peduncles post heading. Significant enrichment of genes belonging to key transcription factor families highlights their specialized roles in peduncle elongation and grain filling before and after heading, respectively. A comparison with anther/pollen development-related genes provided 76 candidates with overlapping roles in anther/pollen development and peduncle elongation. Some of these are important for carbohydrate remobilization to the developing grains. These can be engineered to combat with incomplete panicle exsertion in male sterile lines and manipulate carbohydrate dynamics in grasses. Overall, this study provides baseline information about potential target genes for engineering peduncle elongation with implications on plant height, biomass composition and grain yields in rice.

Nano-scale zero valent iron modulates Fe/Cd transporters and immobilizes soil Cd for production of Cd free rice

In recent decades, nanoscale zero valent iron (nZVI) has been found to be a promising approach for heavy metal remediation. This study is the first report highlighting the role of nZVI to ameliorate Cadmium (Cd) stress in rice along with its effects in expressions of transporter genes, agronomic parameters and grain nutrient status. Initially, 3 concentration of Cd (10, 50, 250 μM) and nZVI (50, 100, 200 mg L^-1) were selected. PCA analysis based on growth parameters, photosynthetic pigment contents and lipid peroxidation rate confirmed that 100 mg L^-1 nZVI was most suitable for remediation of 10 μM Cd. It was evident that, nZVI can alleviate Cd-induced toxic effects by enhancing antioxidant defense mechanisms and other physiological processes in plants. nZVI treated rice seedlings also showed upregulation of phytochelatins which aided in Cd chelation within vacuoles. Study of root morphology with scanning electron microscopy and ROS imaging with confocal microscopy confirmed that nZVI could alleviate oxidative stress due to Cd uptake. In nZVI treated rice seedlings, gene expressions of iron (Fe) transporters (like, IRT1,IRT2,YSL2,YSL15) which are responsible for both Fe and Cd uptake were significantly down-regulated whereas, OsVIT1 and OsCAX4 genes were over expressed which lead to sequestration of Cd in vacuoles. Cd localization assay with leadmium proved that Cd translocation was reduced with nZVI treatment. To further validate our findings a pot experiment was carried out where it was found that nZVI could immobilize Cd in soil prevented accumulation of Cd in rice grains in addition to improving yield.

Genome- wide structural and functional variant discovery of rice landraces using genotyping by sequencing

Rice landraces are vital genetic resources for agronomic and quality traits but the undeniable collection of Kerala landraces remains poorly delineated. To effectively conserve, manage, and use these resources, understanding the genomic structure of germplasm is essential. Genotyping by sequencing (GBS) enables identification of an immense number of single nucleotide polymorphism (SNP) and insertion deletion (InDel) from 96 rice germplasm. In the present study, a total of 16.9 × 10⁷ reads were generated, and among that 16.3 × 10⁷ reads were mapped to the indica reference genome. Exploring GBS data unfolded a wide genomic variations including 82,59,639 SNPs and 1,07,140 Indels. Both neighbor-joining tree and principal coordinate analysis with InDel markers revealed the selected germplasm in this study as highly diverse in structure. We assembled unmapped reads which were further employed for gene ontology analysis. These unmapped sequences that are generally expelled from subsequent studies of GBS data analysis may exist as an unexplored resort for several novel significant biological findings. The discovery of SNPs from the haplotyping results of GS3 and GIF1 genes provided insight into marker- assisted selection based on grain size and yield and can be utilized for rice yield improvement. To our knowledge, this is the first report on structural variation analysis using the GBS platform in rice landraces collected from Kerala. Genomic information from this study endows with valuable resources for perceptive rice landrace structure and can also facilitate sequencing-based molecular breeding.

Genome-wide analysis of polymorphisms identified domestication-associated long low-diversity region carrying important rice grain size/weight quantitative trait loci

Rice grain size and weight are major determinants of grain quality and yield and so have been under rigorous selection since domestication. However, the genetic basis for contrasting grain size/weight trait among Indian germplasms and their association with domestication-driven evolution is not well understood. In this study, two long (LGG) and two short grain (SGG) genotypes were resequenced. LGG (LGR and PB 1121) differentiated from SGG (Sonasal and Bindli) by 504 439 single nucleotide polymorphisms (SNPs) and 78 166 insertion-and-deletion polymorphisms. The LRK gene cluster was different and a truncation mutation in the LRK8 kinase domain was associated with LGG. Phylogeny with 3000 diverse rice accessions revealed that the four sequenced genotypes belonged to the japonica group and were at the edge of the clades indicating them to be the potential source of genetic diversity available in Indian rice germplasm. Six SNPs were significantly associated with grain size/weight and the top four of these could be validated in mapping a population, suggesting this study as a valuable resource for high-throughput genotyping. A contiguous long low-diversity region (LDR) of approximately 6 Mb carrying a major grain weight quantitative trait loci (harbouring OsTOR gene) was identified on Chromosome 5. This LDR was identified as an evolutionary important site with significant positive selection and multiple selection sweeps, and showed association with many domestication-related traits, including grain size/weight. The aus population retained more allelic variations in the LDR than the japonica and indica populations, suggesting it to be one of the divergence loci. All the data and analyses can be accessed from the RiceSzWtBase database.

In Silico Identification of QTL-Based Polymorphic Genes as Salt-Responsive Potential Candidates through Mapping with Two Reference Genomes in Rice

Recent advances in next generation sequencing have created opportunities to directly identify genetic loci and candidate genes for abiotic stress responses in plants. With the objective of identifying candidate genes within the previously identified QTL-hotspots, the whole genomes of two divergent cultivars for salt responses, namely At 354 and Bg 352, were re-sequenced using Illumina Hiseq 2500 100PE platform and mapped to Nipponbare and R498 genomes. The sequencing results revealed approximately 2.4 million SNPs and 0.2 million InDels with reference to Nipponbare while 1.3 million and 0.07 million with reference to R498 in two parents. In total, 32,914 genes were reported across all rice chromosomes of this study. Gene mining within QTL hotspots revealed 1236 genes, out of which 106 genes were related to abiotic stress. In addition, 27 abiotic stress-related genes were identified in non-QTL regions. Altogether, 32 genes were identified as potential genes containing polymorphic non-synonymous SNPs or InDels between two parents. Out of 10 genes detected with InDels, tolerant haplotypes of Os01g0581400, Os10g0107000, Os11g0655900, Os12g0622500, and Os12g0624200 were found in the known salinity tolerant donor varieties. Our findings on different haplotypes would be useful in developing resilient rice varieties for abiotic stress by haplotype-based breeding studies.

Discovery of Functional SNPs via Genome-Wide Exploration of Malaysian Pigmented Rice Varieties

Recently, rice breeding program has shown increased interests on the pigmented rice varieties due to their benefits to human health. However, the genetic variation of pigmented rice varieties is still scarce and remains unexplored. Hence, we performed genome-wide SNP analysis from the genome resequencing of four Malaysian pigmented rice varieties, representing two black and two red rice varieties. The genome of four pigmented varieties was mapped against Nipponbare reference genome sequences, and 1.9 million SNPs were discovered. Of these, 622 SNPs with polymorphic sites were identified in 258 protein-coding genes related to metabolism, stress response, and transporter. Comparative analysis of 622 SNPs with polymorphic sites against six rice SNP datasets from the Ensembl Plants variation database was performed, and 70 SNPs were identified as novel SNPs. Analysis of SNPs in the flavonoid biosynthetic genes revealed 40 nonsynonymous SNPs, which has potential as molecular markers for rice seed colour identification. The highlighted SNPs in this study show effort in producing valuable genomic resources for application in the rice breeding program, towards the genetic improvement of new and improved pigmented rice varieties.

Rice Genomics: over the Past Two Decades and into the Future

Domestic rice (Oryza sativa L.) is one of the most important cereal crops, feeding a large number of worldwide populations. Along with various high-throughput genome sequencing projects, rice genomics has been making great headway toward direct field applications of basic research advances in understanding the molecular mechanisms of agronomical traits and utilizing diverse germplasm resources. Here, we briefly review its achievements over the past two decades and present the potential for its bright future.

Cas9/sgRNA-based genome editing and other reverse genetic approaches for functional genomic studies in rice

One of the important and direct ways of investigating the function of a gene is to characterize the phenotypic consequences associated with loss or gain-of-function of the corresponding gene. These mutagenesis strategies have been successfully deployed in Arabidopsis, and subsequently extended to crop species including rice. Researchers have made vast advancements in the area of rice genomics and functional genomics, as it is a diploid plant with a relatively smaller genome size unlike other cereals. The advent of rice genome research and the annotation of high-quality genome sequencing along with the developments in databases and computer searches have enabled the functional characterization of unknown genes in rice. Further, with the improvements in the efficiency of regeneration and transformation protocols, it has now become feasible to produce sizable mutant populations in indica rice varieties also. In this review, various mutagenesis methods, the current status of the mutant resources, limitations and strengths of insertional mutagenesis approaches and also results obtained with suitable screens for stress tolerance in rice are discussed. In addition, targeted genome editing using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) or Cas9/single-guide RNA system and its potential applications in generating transgene-free rice plants through genome engineering as an efficient alternative to classical transgenic technology are also discussed.

Assembly of an early-matured japonica (Geng) rice genome, Suijing18, based on PacBio and Illumina sequencing

The early-matured japonica (Geng) rice variety, Suijing18 (SJ18), carries multiple elite traits including durable blast resistance, good grain quality, and high yield. Using PacBio SMRT technology, we produced over 25 Gb of long-read sequencing raw data from SJ18 with a coverage of 62×. Using Illumina paired-end whole-genome shotgun sequencing technology, we generated 59 Gb of short-read sequencing data from SJ18 (23.6 Gb from a 200 bp library with a coverage of 59× and 35.4 Gb from an 800 bp library with a coverage of 88×). With these data, we assembled a single SJ18 genome and then generated a set of annotation data. These data sets can be used to test new programs for variation deep mining, and will provide new insights into the genome structure, function, and evolution of SJ18, and will provide essential support for biological research in general.

Whole-Genome Characteristics and Polymorphic Analysis of Vietnamese Rice Landraces as a Comprehensive Information Resource for Marker-Assisted Selection

Next generation sequencing technologies have provided numerous opportunities for application in the study of whole plant genomes. In this study, we present the sequencing and bioinformatic analyses of five typical rice landraces including three indica and two japonica with potential blast resistance. A total of 688.4 million 100 bp paired-end reads have yielded approximately 30-fold coverage to compare with the Nipponbare reference genome. Among them, a small number of reads were mapped to both chromosomes and organellar genomes. Over two million and eight hundred thousand single nucleotide polymorphisms (SNPs) and insertions and deletions (InDels) in indica and japonica lines have been determined, which potentially have significant impacts on multiple transcripts of genes. SNP deserts, contiguous SNP-low regions, were found on chromosomes 1, 4, and 5 of all genomes of rice examined. Based on the distribution of SNPs per 100 kilobase pairs, the phylogenetic relationships among the landraces have been constructed. This is the first step towards revealing several salient features of rice genomes in Vietnam and providing significant information resources to further marker-assisted selection (MAS) in rice breeding programs.

Rice Improvement Through Genome-Based Functional Analysis and Molecular Breeding in India

Rice is one of the main pillars of food security in India. Its improvement for higher yield in sustainable agriculture system is also vital to provide energy and nutritional needs of growing world population, expected to reach more than 9 billion by 2050. The high quality genome sequence of rice has provided a rich resource to mine information about diversity of genes and alleles which can contribute to improvement of useful agronomic traits. Defining the function of each gene and regulatory element of rice remains a challenge for the rice community in the coming years. Subsequent to participation in IRGSP, India has continued to contribute in the areas of diversity analysis, transcriptomics, functional genomics, marker development, QTL mapping and molecular breeding, through national and multi-national research programs. These efforts have helped generate resources for rice improvement, some of which have already been deployed to mitigate loss due to environmental stress and pathogens. With renewed efforts, Indian researchers are making new strides, along with the international scientific community, in both basic research and realization of its translational impact.

Genome-wide DNA polymorphisms in Kavuni, a traditional rice cultivar with nutritional and therapeutic properties

Although rice genome was sequenced in the year 2002, efforts in resequencing the large number of available accessions, landraces, traditional cultivars, and improved varieties of this important food crop are limited. We have initiated resequencing of the traditional cultivars from India. Kavuni is an important traditional rice cultivar from South India that attracts premium price for its nutritional and therapeutic properties. Whole-genome sequencing of Kavuni using Illumina platform and SNPs analysis using Nipponbare reference genome identified 1 150 711 SNPs of which 377 381 SNPs were located in the genic regions. Non-synonymous SNPs (62 708) were distributed in 19 251 genes, and their number varied between 1 and 115 per gene. Large-effect DNA polymorphisms (7769) were present in 3475 genes. Pathway mapping of these polymorphisms revealed the involvement of genes related to carbohydrate metabolism, translation, protein-folding, and cell death. Analysis of the starch biosynthesis related genes revealed that the granule-bound starch synthase I gene had T/G SNPs at the first intron/exon junction and a two-nucleotide combination, which were reported to favour high amylose content and low glycemic index. The present study provided a valuable genomics resource to study the rice varieties with nutritional and medicinal properties.

Functional cereals for production in new and variable climates

Adaptation of cereal crops to variable or changing climates requires that essential quality attributes are maintained to deliver food that will be acceptable to human consumers. Advances in cereal genomics are delivering insights into the molecular basis of nutritional and functional quality traits in cereals and defining new genetic resources. Understanding the influence of the environment on expression of these traits will support the retention of these essential functional properties during climate adaptation. New cereals for use as whole grain or ground to flour for other food products may be based upon the traditional species such as rice and wheat currently used in these food applications but may also include new options exploiting genomics tools to allow accelerated domestication of new species.