Genome-wide discovery of DNA polymorphisms by whole genome sequencing differentiates weedy and cultivated rice

Characteristics analysis of Early Responsive to Dehydration genes in Arabidopsis thaliana (AtERD)

Early Responsive to Dehydration (ERD) genes are rapidly induced in response to various biotic and abiotic stresses, such as bacteria, drought, light, temperature and high salt in Arabidopsis thaliana. Sixteen ERD of Arabidopsis thaliana (AtERD) genes have been previously identified. The lengths of the coding region of the genes are 504-2838 bp. They encode 137-745 amino acids. In this study, the AtERD genes structure and promoter are analyzed through bioinformatics, and a overall function is summarized and a systematic signal pathway involving AtERD genes is mapped. AtERD9, AtERD11 and AtERD13 have the GST domain. AtERD10 and AtERD14 have the Dehyd domain. The promoters regions contain 32 light responsive elements, 23 ABA responsive elements, 5 drought responsive elements, 5 meristem expression related elements and 132 core promoter elements. The study provides a theoretical guidance for subsequent studies of AtERD genes.

Genome-wide discovery of genetic variations between rice cultivars with contrasting drought stress response and their potential functional relevance

Drought stress is a serious threat to rice productivity. Investigating genetic variations between drought-tolerant (DT) and drought-sensitive (DS) rice cultivars may decipher the candidate genes/regulatory regions involved in drought stress tolerance/response. In this study, whole-genome resequencing data of four DS and five DT rice cultivars were analyzed. We identified a total of approximately 4.8 million single nucleotide polymorphisms (SNPs) and 0.54 million insertions/deletions (InDels). The genetic variations (162,638 SNPs and 17,217 InDels) differentiating DS and DT rice cultivars were found to be unevenly distributed throughout the rice genome; however, they were more frequent near the transcription start and stop sites than in the genic regions. The cis-regulatory motifs representing the binding sites of stress-related transcription factors (MYB, HB, bZIP, ERF, ARR, and AREB) harboring the SNPs/InDels in the promoter regions of a few differentially expressed genes (DEGs) were identified. Importantly, many of these DEGs were located within the drought-associated quantitative trait loci. Overall, this study provides a valuable large-scale genotyping resource and facilitates the discovery of candidate genes associated with drought stress tolerance in rice.

Discovery of DNA polymorphisms via genome-resequencing and development of molecular markers between two barley cultivars

Genome resequencing uncovers genome-wide DNA polymorphisms that are useful for the development of high-density InDel markers between two barley cultivars. Discovering genomic variations and developing genetic markers are crucial for genetics studies and molecular breeding in cereal crops. Although InDels (insertions and deletions) have become popular because of their abundance and ease of detection, discovery of genome-wide DNA polymorphisms and development of InDel markers in barley have lagged behind other cereal crops such as rice, maize and wheat. In this study, we re-sequenced two barley cultivars, Golden Promise (GP, a classic British spring barley variety) and Hua30 (a Chinese spring barley variety), and mapped clean reads to the reference Morex genome, and identified in total 13,933,145 single nucleotide polymorphisms (SNPs) and 1,240,456 InDels for GP with Morex, 11,297,100 SNPs and 781,687 InDels for Hua30 with Morex, and 13,742,399 SNPs and 1,191,597 InDels for GP with Hua30. We further characterized distinct types, chromosomal distribution patterns, genome location, functional effect, and other features of these DNA polymorphisms. Additionally, we revealed the functional relevance of these identified SNPs/InDels regarding different flowering times between Hua30 and GP within 17 flowering time genes. Furthermore, we developed a series of InDel markers and validated them experimentally in 43 barley core accessions, respectively. Finally, we rebuilt population structure and phylogenetic tree of these 43 barley core accessions. Collectively, all of these genetic resources will facilitate not only the basic research but also applied research in barley.

Exploring Genomic Variations in Nematode-Resistant Mutant Rice Lines

Rice (Oryza sativa) production is seriously affected by the root-knot nematode Meloidogyne graminicola, which has emerged as a menace in upland and irrigated rice cultivation systems. Previously, activation tagging in rice was utilized to identify candidate gene(s) conferring resistance against M. graminicola. T-DNA insertional mutants were developed in a rice landrace (acc. JBT 36/14), and four mutant lines showed nematode resistance. Whole-genome sequencing of JBT 36/14 was done along with the four nematode resistance mutant lines to identify the structural genetic variations that might be contributing to M. graminicola resistance. Sequencing on Illumina NovaSeq 6000 platform identified 482,234 genetic variations in JBT 36/14 including 448,989 SNPs and 33,245 InDels compared to reference indica genome. In addition, 293,238-553,648 unique SNPs and 32,395-65,572 unique InDels were found in the four mutant lines compared to their JBT 36/14 background, of which 93,224 SNPs and 8,170 InDels were common between all the mutant lines. Functional annotation of genes containing these structural variations showed that the majority of them were involved in metabolism and growth. Trait analysis revealed that most of these genes were involved in morphological traits, physiological traits and stress resistance. Additionally, several families of transcription factors, such as FAR1, bHLH, and NAC, and putative susceptibility (S) genes, showed the presence of SNPs and InDels. Our results indicate that subject to further genetic validations, these structural genetic variations may be involved in conferring nematode resistance to the rice mutant lines.

Genome-wide SNP and InDel analysis of three Philippine mango species inferred from whole-genome sequencing

Background

The Philippines is among the top 10 major exporters of mango worldwide. However, genomic studies of Philippine mangoes remain largely unexplored and lacking. Here, we sequenced the whole genome of the three Philippine mango species, namely, Mangifera odorata (Huani), Mangifera altissima (Paho), and Mangifera indica “Carabao” variety using Illumina HiSeq 2500, to identify and analyze their genome-wide variants (SNPs and InDels).

Results

The high confidence variants were identified by successfully mapping 93–95% of the quality-filtered reads to the Alphonso and Tommy Atkins mango reference genomes. Using these two currently available mango genomes, most variants were observed in M. odorata (4,353,063 and 4,277,287), followed by M. altissima (3,392,763 and 3,449,917), and lastly, M. indica Carabao (2,755,267 and 2,852,480). Approximately 50, 46, and 38% of the variants were unique in the three Philippine mango genomes. The analysis of variant effects and functional annotation across the three mango species revealed 56,982 variants with high-impact effects mapped onto 37,746 genes, of which 25% were found to be novel. The affected mango genes include those with potential economic importance such as 6945 genes for defense/resistance/immune response, 323 genes for fruit development, and 338 genes for anthocyanin production.

Conclusions

To date, this is the first sequencing effort to comprehensively analyze genome-wide variants essential for the development of genome-wide markers specific to these mango species native to the Philippines. This study provides an important genomic resource that can be used for the genetic improvement of mangoes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00326-3.

Copper oxide (CuO) nanoparticles affect yield, nutritional quality, and auxin associated gene expression in weedy and cultivated rice (Oryza sativa L.) grains

Weedy rice grows competitively with cultivated rice and significantly diminishes rice grain production worldwide. The different effects of Cu-based nanomaterials on the production of weedy and cultivated rice, especially the grain qualities are not known. Grains were collected from weedy and cultivated rice grown for four months in field soil amended with nanoscale CuO (nCuO), bulk CuO (bCuO), and copper sulfate (CuSO₄) at 0, 75, 150, 300, and 600 mg Cu/kg soil. Cu translocation, essential element accumulation, yield, sugar, starch, protein content, and the expression of auxin associated genes in grains were determined. The grains of weedy and cultivated rice were differentially impacted by CuO-based compounds. At ≥300 mg/kg, nCuO and bCuO treated rice had no grain production. Treatment at 75 mg/kg significantly decreased grain yield as compared to control with the order: bCuO (by 88.7%) > CuSO₄ (by 47.2%) ~ nCuO (by 38.3% only in cultivated rice); at the same dose, the Cu grain content was: nCuO ~ CuSO₄ > bCuO > control. In weedy grains, K, Mg, Zn, and Ca contents were decreased by 75 and 150 mg/kg nCuO by up to 47.4%, 34.3%, 37.6%, and 60.0%, but no such decreases were noted in cultivated rice, and Fe content was increased by up to 88.6%, and 53.2%. In rice spikes, nCuO increased Mg, Ca, Fe, and Zn levels by up to 118.1%, 202.6%, 133.8%, and 103.9%, respectively. Nanoscale CuO at 75 and 150 mg/kg upregulated the transcription of an auxin associated gene by 5.22- and 1.38-fold, respectively, in grains of weedy and cultivated rice. The biodistribution of Cu-based compounds in harvested grain was determined by two-photon microscopy. These findings demonstrate a cultivar-specific and concentration-dependent response of rice to nCuO. A potential use of nCuO at 75 and 150 mg/kg in cultivar-dependent delivery system was suggested based on enhanced grain nutritional quality, although the yield was compromised. This knowledge, at the physiological and molecular level, provides valuable information for the future use of Cu-based nanomaterials in sustainable agriculture.

Whole-Genome Sequencing and RNA-Seq Reveal Differences in Genetic Mechanism for Flowering Response between Weedy Rice and Cultivated Rice

Flowering is a key agronomic trait that influences adaptation and productivity. Previous studies have indicated the genetic complexity associated with the flowering response in a photoinsensitive weedy rice accession PSRR-1 despite the presence of a photosensitive allele of a key flowering gene Hd1. In this study, we used whole-genome and RNA sequencing data from both cultivated and weedy rice to add further insights. The de novo assembly of unaligned sequences predicted 225 genes, in which 45 were specific to PSRR-1, including two genes associated with flowering. Comparison of the variants in PSRR-1 with the 3K rice genome (RG) dataset identified unique variants within the heading date QTLs. Analyses of the RNA-Seq result under both short-day (SD) and long-day (LD) conditions revealed that many differentially expressed genes (DEGs) colocalized with the flowering QTLs, and some DEGs such as Hd1, OsMADS56, Hd3a, and RFT1 had unique variants in PSRR-1. Ehd1, Hd1, OsMADS15, and OsMADS56 showed different alternate splicing (AS) events between genotypes and day length conditions. OsMADS56 was expressed in PSRR-1 but not in Cypress under both LD and SD conditions. Based on variations in both sequence and expression, the unique flowering response in PSRR-1 may be due to the high-impact variants of flowering genes, and OsMADS56 is proposed as a key regulator for its day-neutral flowering response.

Field Screen and Genotyping of Phaseolus vulgaris against Two Begomoviruses in Georgia, USA

Simple Summary

Snap bean (Phaseolus vulgaris) production and quality have been negatively impacted by two whitefly-transmitted begomoviruses: cucurbit leaf crumple virus (CuLCrV) and sida golden mosaic Florida virus (SiGMFV), which often appear as a mixed infection in Georgia. However, there is no information available in terms of resistance to these two viruses in commercial cultivars/genotypes. Hence, commercially available snap bean varieties/genotypes (n = 84 in 2018; n = 80 in 2019; most of the genotypes were common in both years (with a few exceptions) were screened in two field seasons of 2018 and 2019. We also included two commonly grown Lima bean (Phaseolus lunatus) varieties in our field screening. As a result of this screening, we identified twenty Phaseolus genotypes with high-to-moderate levels of resistance and twenty-one genotypes with high levels of susceptibility. While there were differences among the Phaseolus spp. in severity of viral symptoms, suggesting differential susceptibility to viruses (CuLCrV and SiGMFV) and potential field resistance, the resistance mechanism is yet to be characterized. However, based on the greenhouse evaluation with two genotypes-each (susceptible vs. resistant) exposed to viruliferous whiteflies infected with CuLCrV and SiGMFV, we observed that the susceptible genotypes accumulated higher copy numbers of both viruses and displayed severe crumple severity compared to the resistant genotypes, indicating that resistant might potentially be against the virus complex than against the whiteflies. Adult whitefly counts differed among the Phaseolus spp. in both the years, indicating variability in host preference. We further sequenced 82 genotypes (80 snap bean and two Lima bean) to unravel the variations within the genomes. Genome sequencing followed by bioinformatic analyses revealed a considerable number of sequence variants, single nucleotide polymorphisms (SNPs), and insertions and deletions (InDels) in the genomes. Considering the variations in disease response and the underlying variations in the sequenced genomes, it can be speculated that some of the phenotypic variations (against CuLCrV and SiGMFV) could be due to a high level of genomic variation in the host. Future genome-wide association studies with the identified genomic variants may shed some light on this.

Abstract

The production and quality of Phaseolusvulgaris (snap bean) have been negatively impacted by leaf crumple disease caused by two whitefly-transmitted begomoviruses: cucurbit leaf crumple virus (CuLCrV) and sida golden mosaic Florida virus (SiGMFV), which often appear as a mixed infection in Georgia. Host resistance is the most economical management strategy against whitefly-transmitted viruses. Currently, information is not available with respect to resistance to these two viruses in commercial cultivars. In two field seasons (2018 and 2019), we screened Phaseolus spp. genotypes (n = 84 in 2018; n = 80 in 2019; most of the genotypes were common in both years with a few exceptions) for resistance against CuLCrV and/or SiGMFV. We also included two commonly grown Lima bean (Phaseolus lunatus) varieties in our field screening. Twenty Phaseolus spp. genotypes with high to moderate-levels of resistance (disease severity ranging from 5%–50%) to CuLCrV and/or SiGMFV were identified. Twenty-one Phaseolus spp. genotypes were found to be highly susceptible with a disease severity of ≥66%. Furthermore, based on the greenhouse evaluation with two genotypes-each (two susceptible and two resistant; identified in field screen) exposed to viruliferous whiteflies infected with CuLCrV and SiGMFV, we observed that the susceptible genotypes accumulated higher copy numbers of both viruses and displayed severe crumple severity compared to the resistant genotypes, indicating that resistance might potentially be against the virus complex rather than against the whiteflies. Adult whitefly counts differed significantly among Phaseolus genotypes in both years. The whole genome of these Phaseolus spp. [snap bean (n = 82); Lima bean (n = 2)] genotypes was sequenced and genetic variability among them was identified. Over 900 giga-base (Gb) of filtered data were generated and >88% of the resulting data were mapped to the reference genome, and SNP and Indel variants in Phaseolus spp. genotypes were obtained. A total of 645,729 SNPs and 68,713 Indels, including 30,169 insertions and 38,543 deletions, were identified, which were distributed in 11 chromosomes with chromosome 02 harboring the maximum number of variants. This phenotypic and genotypic information will be helpful in genome-wide association studies that will aid in identifying the genetic basis of resistance to these begomoviruses in Phaseolus spp.

Whole genome sequence analysis of rice genotypes with contrasting response to salinity stress

Salinity is a major abiotic constraint for rice farming. Abundant natural variability exists in rice germplasm for salt tolerance traits. Since few studies focused on the genome level variation in rice genotypes with contrasting response to salt stress, genomic resequencing in diverse genetic materials is needed to elucidate the molecular basis of salt tolerance mechanisms. The whole genome sequences of two salt tolerant (Pokkali and Nona Bokra) and three salt sensitive (Bengal, Cocodrie, and IR64) rice genotypes were analyzed. A total of 413 million reads were generated with a mean genome coverage of 93% and mean sequencing depth of 18X. Analysis of the DNA polymorphisms revealed that 2347 nonsynonymous SNPs and 51 frameshift mutations could differentiate the salt tolerant from the salt sensitive genotypes. The integration of genome-wide polymorphism information with the QTL mapping and expression profiling data led to identification of 396 differentially expressed genes with large effect variants in the coding regions. These genes were involved in multiple salt tolerance mechanisms, such as ion transport, oxidative stress tolerance, signal transduction, and transcriptional regulation. The genome-wide DNA polymorphisms and the promising candidate genes identified in this study represent a valuable resource for molecular breeding of salt tolerant rice varieties.

A Novel Mutation of the NARROW LEAF 1 Gene Adversely Affects Plant Architecture in Rice (Oryza sativa L.)

Plant architecture is critical for enhancing the adaptability and productivity of crop plants. Mutants with an altered plant architecture allow researchers to elucidate the genetic network and the underlying mechanisms. In this study, we characterized a novel nal1 rice mutant with short height, small panicle, and narrow and thick deep green leaves that was identified from a cross between a rice cultivar and a weedy rice accession. Bulked segregant analysis coupled with genome re-sequencing and cosegregation analysis revealed that the overall mutant phenotype was caused by a 1395-bp deletion spanning over the last two exons including the transcriptional end site of the nal1 gene. This deletion resulted in chimeric transcripts involving nal1 and the adjacent gene, which were validated by a reference-guided assembly of transcripts followed by PCR amplification. A comparative transcriptome analysis of the mutant and the wild-type rice revealed 263 differentially expressed genes involved in cell division, cell expansion, photosynthesis, reproduction, and gibberellin (GA) and brassinosteroids (BR) signaling pathways, suggesting the important regulatory role of nal1. Our study indicated that nal1 controls plant architecture through the regulation of genes involved in the photosynthetic apparatus, cell cycle, and GA and BR signaling pathways.

Comparative Transcriptomics of Rice Genotypes with Contrasting Responses to Nitrogen Stress Reveals Genes Influencing Nitrogen Uptake through the Regulation of Root Architecture

The indiscriminate use of nitrogenous fertilizers continues unabated for commercial crop production, resulting in air and water pollution. The development of rice varieties with enhanced nitrogen use efficiency (NUE) will require a thorough understanding of the molecular basis of a plant’s response to low nitrogen (N) availability. The global expression profiles of root tissues collected from low and high N treatments at different time points in two rice genotypes, Pokkali and Bengal, with contrasting responses to N stress and contrasting root architectures were examined. Overall, the number of differentially expressed genes (DEGs) in Pokkali (indica) was higher than in Bengal (japonica) during low N and early N recovery treatments. Most low N DEGs in both genotypes were downregulated whereas early N recovery DEGs were upregulated. Of these, 148 Pokkali-specific DEGs might contribute to Pokkali’s advantage under N stress. These DEGs included transcription factors and transporters and were involved in stress responses, growth and development, regulation, and metabolism. Many DEGs are co-localized with quantitative trait loci (QTL) related to root growth and development, chlorate-resistance, and NUE. Our findings suggest that the superior growth performance of Pokkali under low N conditions could be due to the genetic differences in a diverse set of genes influencing N uptake through the regulation of root architecture.

Whole-genome sequencing reveals uniqueness of black-hulled and straw-hulled weedy rice genomes

Both SH and BHA weedy rice genotypes evolved independently and have distinct genomic composition. Different genetic mechanisms may be responsible for their competitiveness and adaptation to diverse environmental conditions. Two major types of weedy rice are recognized in the USA based on morphology: straw-hull (SH) and black-hull awned (BHA) weedy rice. We performed whole-genome resequencing of a SH weedy rice 'PSRR-1', a BHA weedy rice 'BHA1115', and a japonica cultivar 'Cypress' to delineate genome-wide differences and their relevance to genetics and evolution of weedy attributes. The high-quality reads were uniformly distributed with 82-88% genome coverage. The number of genotype-specific SNPs and InDels was highest in Cypress, followed by BHA1115 and PSRR-1. However, more genes were affected in BHA1115 compared with other two genotypes which is evident from the number of high-impact SNPs and InDels. Haplotype analysis of selected genes involved in domestication, adaptation, and agronomic performance not only differentiated SH from BHA weedy rice and supported evolution of weedy rice through de-domestication, but also validated the function of several genes such as qAn-1, qAn-2, Bh4, Rc, SD1, OsLG1, and OsC1. Several candidate genes were identified for previously reported seed dormancy and seed shattering QTLs. The SH and BHA weedy rice have distinct genomic composition, and the BHA weedy rice likely diverged earlier than SH weedy rice. The accumulation of plant development, reproduction, and defense-related genes in weedy rice possibly helped them to compete, survive, and spread under a wide range of environmental conditions by employing novel and diverse mechanisms. The genomic resources will be useful for both weed management and rice improvement by exploring the molecular basis of key agronomic, adaptive, and domestication attributes.

Parallel reaction monitoring revealed tolerance to drought proteins in weedy rice (Oryza sativa f. spontanea)

Drought is a complicated abiotic stress factor with severe effects on rice growth and production. Weedy rice is a valuable genetic resource that possesses a strong capacity for drought tolerance, cold tolerance, and salt tolerance, and is an excellent material for studying rice tolerance. Here, according to comprehensive tolerance to drought index D, accession WR16 was selected based on strong drought tolerance among 133 studied weedy red rice germplasms. WR16 was compared with Oryza sativa ssp. Japonica. cv. IAPAR-9, a reference genotype originating from Brazil. In addition, accession WR24 was classified as moderately tolerant to drought accessions. Transcriptomic and proteomic analyses were combined to identify 38 co-upregulated proteins related to drought tolerance, and targeted parallel reaction monitoring (PRM) was used to precisely quantify and verify nine proteins in the complex backgrounds. Result showed that six proteins were significantly (Fisher's exact P value < 0.05) related to drought tolerance in accessions WR16 and WR24. Among them, OS09T0478300-01, OS09T0530300-01, and OS01T0800500-01 formed a combined defense system to respond to drought stress in weedy rice. Results of these studies provide comprehensive information for precisely identifying and verifying tolerance to drought proteins and lay a solid theoretical foundation for research on drought tolerance mechanisms.

Genome-wide discovery of DNA variants in cucumber somaclonal lines

The emergence of somaclonal variability in in vitro cultures is undesirable during micropropagation, but this phenomenon may be a source of genetic variability sought by breeders. The main factors that affect the appearance of variability are known, but the exact mechanism has not yet been determined. In this paper, we used next-generation sequencing and comparative genomics to study changes in the genomes of cucumber lines resulting from in vitro regeneration and somaclonal mutation in comparison to a reference, the highly inbred B10 line. The total number of obtained polymorphisms differed between the three somaclonal lines S1, S2 and S3, with 8369, 7591 and 44510, respectively. Polymorphisms occurred most frequently in non-coding regions and were mainly SNPs. High-impact changes accounted for 1%-3% of all polymorphisms and most often caused an open reading frame shift. Functional analysis of genes affected by high impact variants showed that they were related to transport, biosynthetic processes, nucleotide-containing compounds and cellular protein modification processes. The obtained results indicated significant factors affecting somaclonal variability and the appearance of changes in the genome, and demonstrated a lack of dependence between phenotype and the number of genomic polymorphisms.

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.

Comprehensive Analysis of the Chitinase Gene Family in Cucumber (Cucumis sativus L.): From Gene Identification and Evolution to Expression in Response to Fusarium oxysporum

Chitinases, a subgroup of pathogenesis-related proteins, are responsible for catalyzing the hydrolysis of chitin. Accumulating reports indicate that chitinases play a key role in plant defense against chitin-containing pathogens and are therefore good targets for defense response studies. Here, we undertook an integrated bioinformatic and expression analysis of the cucumber chitinases gene family to identify its role in defense against Fusarium oxysporum f. sp. cucumerinum. A total of 28 putative chitinase genes were identified in the cucumber genome and classified into five classes based on their conserved catalytic and binding domains. The expansion of the chitinase gene family was due mainly to tandem duplication events. The expression pattern of chitinase genes was organ-specific and 14 genes were differentially expressed in response to F. oxysporum challenge of fusarium wilt-susceptible and resistant lines. Furthermore, a class I chitinase, CsChi23, was constitutively expressed at high levels in the resistant line and may play a crucial role in building a basal defense and activating a rapid immune response against F. oxysporum. Whole-genome re-sequencing of both lines provided clues for the diverse expression patterns observed. Collectively, these results provide useful genetic resource and offer insights into the role of chitinases in cucumber-F. oxysporum interaction.

Comparative Genomic Analysis of Rice with Contrasting Photosynthesis and Grain Production under Salt Stress

Unfavourable environmental conditions, including soil salinity, lead to decreased rice (Oryza sativa L.) productivity, especially at the reproductive stage. In this study, we examined 30 rice varieties, which revealed significant differences in the photosynthetic performance responses under salt stress conditions during the reproductive stage, which ultimately affected yield components after recovery. In rice with a correlation between net photosynthetic rate (P_N) and intercellular CO₂ concentration (C_i) under salt stress, P_N was found to be negatively correlated with filled grain number after recovery. Applying stringent criteria, we identified 130,317 SNPs and 15,396 InDels between two “high-yield rice” varieties and two “low-yield rice” varieties with contrasting photosynthesis and grain yield characteristics. A total of 2089 genes containing high- and moderate-impact SNPs or InDels were evaluated by gene ontology (GO) enrichment analysis, resulting in over-represented terms in the apoptotic process and kinase activity. Among these genes, 262 were highly expressed in reproductive tissues, and most were annotated as receptor-like protein kinases. These findings highlight the importance of variations in signaling components in the genome and these loci can serve as potential genes in rice breeding to produce a variety with salt avoidance that leads to increased yield in saline soil.