Sequencing quality and library assessment
Whole-genome sequencing of the traditional rice landrace Iluppai Poo Samba generated high-quality paired-end sequencing data suitable for downstream analysis. Sequencing quality assessment demonstrated consistently high base-calling accuracy across both technical replicates. Quality metrics indicated that more than 97% of bases achieved Q20 scores, while over 94% of bases exceeded Q30 quality thresholds, confirming the generation of high-confidence sequencing reads. The average GC content of the sequencing libraries was approximately 45%, which is consistent with the expected genomic GC composition of
Oryza sativa. Quality score distribution profiles further demonstrated stable sequencing performance across read lengths, indicating minimal decline in base quality across sequencing cycles.
Library quality assessment using the Agilent TapeStation platform confirmed successful library preparation and appropriate fragment size distribution for Illumina sequencing. Average library fragment sizes ranged from approximately 433 to 446 bp, indicating uniform fragmentation and efficient adapter ligation during library preparation. The high sequencing quality and optimal library profiles collectively supported the reliability of downstream genome alignment and variant detection analyses.
Genome-wide SNP distribution and genomic variation patterns
Genome-wide variant analysis revealed extensive sequence variation across all rice chromosomes in the Iluppai Poo Samba genome. SNP density analysis using non-overlapping 100 kb genomic windows demonstrated heterogeneous distribution of polymorphisms throughout the genome, with several chromosomal intervals exhibiting elevated SNP density relative to the reference genome (Fig 1). Distinct SNP-enriched hotspot regions and low-variation intervals were observed across multiple chromosomes, indicating uneven genome-wide polymorphism patterns.
The observed variation in SNP density suggests differential accumulation of sequence variation across chromosomal regions, potentially reflecting differences in evolutionary conservation, recombination frequency, selection pressure and genomic organization
(Huang et al., 2011). Similar heterogeneous SNP distribution patterns have been reported in other rice genomic studies and are often associated with adaptive divergence and cultivar-specific variation
(Trung et al., 2017; Wang et al., 2009). The widespread polymorphism detected in Iluppai Poo Samba indicates substantial genomic differentiation from the
Oryza sativa japonica reference genome and highlights the unique genetic diversity retained within this traditional landrace. The heterogeneous distribution of SNPs across the genome suggest that different chromosomal regions have experienced varying levels of genetic diversification and selective pressure during the evolution of Iluppai Poo Samba. Polymorphism-rich genomic regions may harbor genes associated with environmental adaptation, stress tolerance and other agronomically important traits, whereas relatively conserved regions are likely enriched for genes under stronger functional constraint. These genomic variation patterns provide valuable targets for future trait mapping, functional validation and marker-assisted breeding in traditional rice germplasm
(Huang et al., 2011; Wang et al., 2009).
Chromosome-wise variation profiling
Chromosome-wise SNP distribution analysis revealed considerable variation in polymorphism abundance across individual rice chromosomes. Several chromosomes exhibited broader high-density polymorphic regions, whereas others contained localized SNP-enriched intervals interspersed with relatively conserved regions (Fig 2). In particular, chromosomes CP132240.1, CP132244.1, CP132245.1 and CP132246.1 displayed extensive high-density polymorphic regions, while chromosomes CP132239.1 and CP132243.1 exhibited more localized SNP enrichment patterns.
The uneven distribution of SNPs across chromosomes may reflect differences in recombination dynamics, mutation rates, selective constraints and historical evolutionary processes. Regions exhibiting elevated SNP abundance may represent genomic intervals associated with increased genetic diversity and adaptive variation, whereas regions with lower SNP density may correspond to more conserved genomic segments. These chromosome-specific variation landscapes provide useful targets for future investigations of agronomically important traits, including stress adaptation, grain quality and aroma-associated characteristics in traditional rice germplasm. Chromosome-specific differences in SNP abundance may also reflect historical recombination events and cultivar-specific evolutionary processes that have shaped the genetic architecture of Iluppai Poo Samba. The identification of polymorphism-rich chromosomal intervals provides a valuable framework for future quantitative trait locus (QTL) mapping and genome-wide association studies aimed at identifying loci controlling stress tolerance, grain quality and other economically important traits
(Huang et al., 2011; 2010).
Functional classification of genomic variants
Functional annotation of genome-wide variants revealed extensive sequence variation distributed across coding and non-coding genomic regions of the Iluppai Poo Samba genome. A substantial proportion of identified variants were localized within intergenic and intronic regions, indicating widespread regulatory and non-coding genomic variation across the genome. In addition, large numbers of synonymous and non-synonymous SNPs were identified across all chromosomes, suggesting substantial functional variation within protein-coding genes. High-impact variants including stop-gain and stop-loss mutations were also detected across multiple chromosomal regions, potentially contributing to altered gene function and phenotypic diversification. The presence of extensive coding and regulatory-region polymorphisms indicates that genomic variation within Iluppai Poo Samba extends beyond neutral variation and may influence multiple biological pathways associated with adaptation and agronomic performance.
Transition mutations were consistently more abundant than transversion mutations across all chromosomes, with Ts/Tv ratios ranging from 2.24 to 2.50. These values are consistent with mutation patterns commonly observed in plant genomes and support the reliability and biological consistency of the identified variants. Furthermore, the substantial number of non-synonymous variants identified in coding regions suggests the presence of functional allelic diversity that may influence protein function and contribute to phenotypic variation associated with stress responses, development and environmental adaptation. Although experimental validation is required to confirm their biological effects, these variants represent promising candidate loci for future functional genomics, marker development and rice improvement programs
(Depristo et al., 2011; Wakeley, 1996). Detailed summary statistics of genomic variants identified throughout the Iluppai Poo Samba genome are presented in Table 1.
Candidate genes harboring functional variants
Functional annotation identified variants in biologically important genes involved in stress adaptation, ion transport, signaling pathways and transcriptional regulation. Variants were identified in HKT transporter family genes including HKT1;1, HKT1;3, HKT2;1 and HKT2;4, which are known to regulate sodium transport and ion homeostasis during salinity stress in rice. Genetic variation within these transporter-associated loci may influence salinity tolerance and environmental adaptation by altering ion transport efficiency. Although functional validation is required, these variants represent promising candidate alleles for future studies aimed at improving stress resilience in rice cultivars
(Liu et al., 2022).
In addition, variants were detected in kinase-related genes including CPK7, CPK13 and CPK19, which function as key regulators of cellular signaling pathways involved in plant responses to abiotic and biotic stresses. Variants were also identified in KIN17, a conserved DNA/RNA-binding protein implicated in genome maintenance and cellular homeostasis, suggesting additional functional diversity within genes associated with cellular regulation. Additionally, variants were identified in transcription factor and regulatory genes including HSFA2D, GATA15 and BZIP39. These transcription factors regulate stress-responsive gene expression, developmental processes and metabolic pathways, suggesting that functional polymorphisms within these genes may contribute to cultivar-specific adaptation and phenotypic diversity.
Furthermore, variants were detected in genes associated with developmental regulation and cell wall biosynthesis, including IRX9L and IRX14, indicating additional sources of functional genomic diversity that may influence plant growth and structural development. Collectively, these findings suggest that Iluppai Poo Samba harbors a diverse repertoire of functional genetic variation extending beyond neutral sequence polymorphisms. The identified candidate genes provide promising targets for future functional validation, marker-assisted breeding and the genetic improvement of stress resilience and grain quality in rice.