Discovery and mapping of genomic regions governing economically important traits of Basmati rice

Development of miRNA-SSR and target-SSR markers from yield-associate genes and their applicability in the assessment of genetic diversity and association mapping in rice (Oryza sativa L.)

The gene-derived functional markers are considered effective to use in marker-assisted breeding and genetic diversity analysis. As of now, no functional markers have been identified from miRNAs regulating yield traits. The miRNAs play a key role as regulators in controlling the candidate genes involved in grain yield improvement in rice. In this study, 13 miRNA-SSR and their target gene SSR markers were mined from 29 yield-responsive miRNA along with their 29 target genes in rice. The validation of these markers showed that four miRNA-SSRs and one target gene SSR markers had shown polymorphism among 120 diverse rice genotypes. The PIC values ranged from 0.25 (OsARF18-SSR) to 0.72 (miR408-SSR, miR172b-SSR, and miR396f-SSR) with an average value of 0.57. These polymorphic markers grouped 120 rice genotypes into 3 main clusters based on the levels of high genetic diversity. These markers also showed significant association with key yield traits. Among all, miR172b-SSR showed a strong association with plant height in two seasons. This investigation suggests that this new class of molecular markers has great potential in the characterization of rice germplasm by genetic diversity and population structure and in marker-assisted breeding for the development of high-yielding varieties.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-024-01462-z.

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.

Genetic improvement for blast resistance in high-yielding cold-tolerant rice (Oryza sativa L.) cultivar Himalaya 741 by marker-assisted backcross breeding

Blast disease and cold stress are two major yield-limiting factors for rice under temperate climates. Marker-assisted backcross breeding approach (MABB) was employed for the improvement of blast resistance in a popular cold-tolerant variety 'Himalaya741' by introgressing a broad-spectrum resistance locus Pi9 from a Basmati donor PB1637. A combined use of phenotypic selection and marker-based genotypic selection ensured speedy reconstitution of the recurrent parent genome (RPG) in backcross progenies; RPG recovery in most of the progenies was > 96% with three progenies namely, HPU-1-33, -38 and -49 showing complete recovery of recurrent parent genome. Notwithstanding a very higher recovery rate of RPG in introgression lines, the lines still inherited a large linkage block > 13.3 Mb with Pi9 from the donor line PB1637. The donor chromosome segments co-inherited with Pi9 gene, however, did not have any adverse effect on the agronomic performance of the Pi9 introgression lines. Of the eight genetically superior Pi9 introgression lines identified, two exhibited resemblance to Himalaya 741 for most of the agronomic traits in addition to having superior grain length and tiller number. The introgression line HPU-1-81 displayed 44% yield superiority over recurrent parent, primarily due to improvement in yield-contributing traits, namely, tiller number, panicle length, thousand-seed-weight and grain length. All the Pi9 introgression lines displayed a high level of resistance comparable to PB1637 against two highly virulent blast races, which collectively displayed compatibility to 15 different major resistance genes. The introgression lines also possessed reproductive stage cold tolerance similar to recurrent parent under prevailing cold stress conditions. The agronomically superior Pi9 introgression lines developed herein are expected to provide a comparable or better substitute to blast susceptible variety Himalaya 741 for extenuating losses due to cold stress and blast disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03244-w.

Use of Mathematical Methods for the Biosafety Assessment of Agricultural Crops

In Russia and around the world, there are important questions regarding the potential threats to national and biological safety created by genetic technologies and the need to improve or introduce new, justified, and adequate measures for their control, regulation, and prevention. The article shows that a significant volume of the global market is occupied by five major transgenic crops, and producers are ready to switch to crops with an edited genome that has been approved in the United States, Argentina, and other countries. We propose a qualitatively new approach to the risk assessment of edited plants, "Safe Design," and we have also developed an extremely important, fundamentally new approach to the development of methods that combine next-generation sequencing (NGS) and Bioinformatics for the assessment of the crop import biosafety. The proposed mathematical approach provides a detailed analysis of the possible insertions of DNA fragments into the genome of edited crops and a clarification of their biological significance. The developed method can be used in the rapid screening of plants for the presence of potentially dangerous genes, viral sequences, and nonspecific promoter sequences.

Evaluation of Whole-Genome Sequence, Genetic Diversity, and Agronomic Traits of Basmati Rice (Oryza sativa L.)

Basmati is considered a unique varietal group of rice (Oryza sativa L.) because of its aroma and superior grain quality. Previous genetic analyses of rice showed that most of the Basmati varieties are classified into the aromatic group. Despite various efforts, genomic relationship of Basmati rice with other varietal groups and genomic variation in Basmati rice are yet to be understood. In the present study, we resequenced the whole genome of three traditional Basmati varieties at a coverage of more than 25X using Illumina HiSeq2500 and mapped the obtained sequences to the reference genome sequences of Nipponbare (japonica rice), Kasalath (aus rice), and Zhenshan 97 (indica rice). Comparison of these sequences revealed common single nucleotide polymorphisms (SNPs) in the genic regions of three Basmati varieties. Analysis of these SNPs revealed that Basmati varieties showed fewer sequence variations compared with the aus group than with the japonica and indica groups. Gene ontology (GO) enrichment analysis indicated that SNPs were present in genes with various biological, molecular, and cellular functions. Additionally, functional annotation of the Basmati mutated gene cluster shared by Nipponbare, Kasalath, and Zhenshan 97 was found to be associated with the metabolic process involved in the cellular aromatic compound, suggesting that aroma is an important specific genomic feature of Basmati varieties. Furthermore, 30 traditional Basmati varieties were classified into three different groups, aromatic (22 varieties), aus (four varieties), and indica (four varieties), based on genome-wide SNPs. All 22 aromatic Basmati varieties harbored the fragrant-inducing Badh2 allele. We also performed comparative analysis of 13 key agronomic and grain quality traits of Basmati rice and other rice varieties. Three traits including length-to-width ratio of grain (L/W ratio), panicle length (PL), and amylose content (AC) showed significant (P < 0.05 and P < 0.01) differences between the aromatic and indica/aus groups. Comparative analysis of genome structure, based on genome sequence variation and GO analysis, revealed that the Basmati genome was derived mostly from the aus and japonica groups. Overall, whole-genome sequence data and genetic diversity information obtained in this study will serve as an important resource for molecular breeding and genetic analysis of Basmati varieties.

Quantitative trait loci mapping of panicle traits in rice

In this study 90 individuals of recombinant inbred lines (RILs) were developed by crossing subspecies of japonica rice cultivar, ‘Nagdong’ and an indica type cultivar, ‘Cheongcheong’. These individuals were used to identify the quantitative trait loci of panicle traits using SSR markers. A genetic linkage map was constructed using one hundred fifty four simple sequence repeat (SSR) primers covering distance of 1973.6 cM of the whole genome with mean distance of 13.9 cM among markers. QTLs were mapped using composite interval mapping method, nineteen QTLs were recognized for the panicle traits on chromosomes 4, 5, 6, 8, 10, 11, 12 with individual QTL explained 8.8% to 37.9% of phenotypic variation. Two pleiotropic effects loci were found on chromosomes 4 and 6. These QTLs affecting leaf traits, panicle traits and panicle branch traits would be beneficial to high-yield rice improvement.

Regional Association Analysis of MetaQTLs Delineates Candidate Grain Size Genes in Rice

Molecular mapping studies which aim to identify genetic basis of diverse agronomic traits are vital for marker-assisted crop improvement. Numerous Quantitative Trait Loci (QTLs) mapped in rice span long genomic intervals with hundreds to thousands of genes, which limits their utilization for marker-assisted genetic enhancement of rice. Although potent, fine mapping of QTLs is challenging task as it requires screening of large number of segregants to identify suitable recombination events. Association mapping offers much higher resolution as compared to QTL mapping, but detects considerable number of spurious QTLs. Therefore, combined use of QTL and association mapping strategies can provide advantages associated with both these methods. In the current study, we utilized meta-analysis approach to identify metaQTLs associated with grain size/weight in diverse Indian indica and aromatic rice accessions. Subsequently, attempt has been made to narrow-down identified grain size/weight metaQTLs through individual SNP- as well as haplotype-based regional association analysis. The study identified six different metaQTL regions, three of which were successfully revalidated, and substantially scaled-down along with GS3 QTL interval (positive control) by regional association analysis. Consequently, two potential candidate genes within two reduced metaQTLs were identified based on their differential expression profiles in different tissues/stages of rice accessions during seed development. The developed strategy has broader practical utility for rapid delineation of candidate genes and natural alleles underlying QTLs associated with complex agronomic traits in rice as well as major crop plants enriched with useful genetic and genomic information.

Genome-wide generation and use of informative intron-spanning and intron-length polymorphism markers for high-throughput genetic analysis in rice

We developed genome-wide 84634 ISM (intron-spanning marker) and 16510 InDel-fragment length polymorphism-based ILP (intron-length polymorphism) markers from genes physically mapped on 12 rice chromosomes. These genic markers revealed much higher amplification-efficiency (80%) and polymorphic-potential (66%) among rice accessions even by a cost-effective agarose gel-based assay. A wider level of functional molecular diversity (17-79%) and well-defined precise admixed genetic structure was assayed by 3052 genome-wide markers in a structured population of indica, japonica, aromatic and wild rice. Six major grain weight QTLs (11.9-21.6% phenotypic variation explained) were mapped on five rice chromosomes of a high-density (inter-marker distance: 0.98 cM) genetic linkage map (IR 64 x Sonasal) anchored with 2785 known/candidate gene-derived ISM and ILP markers. The designing of multiple ISM and ILP markers (2 to 4 markers/gene) in an individual gene will broaden the user-preference to select suitable primer combination for efficient assaying of functional allelic variation/diversity and realistic estimation of differential gene expression profiles among rice accessions. The genomic information generated in our study is made publicly accessible through a user-friendly web-resource, "Oryza ISM-ILP marker" database. The known/candidate gene-derived ISM and ILP markers can be enormously deployed to identify functionally relevant trait-associated molecular tags by optimal-resource expenses, leading towards genomics-assisted crop improvement in rice.

An Efficient Strategy Combining SSR Markers- and Advanced QTL-seq-driven QTL Mapping Unravels Candidate Genes Regulating Grain Weight in Rice

Development and use of genome-wide informative simple sequence repeat (SSR) markers and novel integrated genomic strategies are vital to drive genomics-assisted breeding applications and for efficient dissection of quantitative trait loci (QTLs) underlying complex traits in rice. The present study developed 6244 genome-wide informative SSR markers exhibiting in silico fragment length polymorphism based on repeat-unit variations among genomic sequences of 11 indica, japonica, aus, and wild rice accessions. These markers were mapped on diverse coding and non-coding sequence components of known cloned/candidate genes annotated from 12 chromosomes and revealed a much higher amplification (97%) and polymorphic potential (88%) along with wider genetic/functional diversity level (16-74% with a mean 53%) especially among accessions belonging to indica cultivar group, suggesting their utility in large-scale genomics-assisted breeding applications in rice. A high-density 3791 SSR markers-anchored genetic linkage map (IR 64 × Sonasal) spanning 2060 cM total map-length with an average inter-marker distance of 0.54 cM was generated. This reference genetic map identified six major genomic regions harboring robust QTLs (31% combined phenotypic variation explained with a 5.7-8.7 LOD) governing grain weight on six rice chromosomes. One strong grain weight major QTL region (OsqGW5.1) was narrowed-down by integrating traditional QTL mapping with high-resolution QTL region-specific integrated SSR and single nucleotide polymorphism markers-based QTL-seq analysis and differential expression profiling. This led us to delineate two natural allelic variants in two known cis-regulatory elements (RAV1AAT and CARGCW8GAT) of glycosyl hydrolase and serine carboxypeptidase genes exhibiting pronounced seed-specific differential regulation in low (Sonasal) and high (IR 64) grain weight mapping parental accessions. Our genome-wide SSR marker resource (polymorphic within/between diverse cultivar groups) and integrated genomic strategy can efficiently scan functionally relevant potential molecular tags (markers, candidate genes and alleles) regulating complex agronomic traits (grain weight) and expedite marker-assisted genetic enhancement in rice.