The study, published in BMC Plant Biology, began with a large phenotyping effort. The researchers screened 116 rice genotypes in hydroponic systems deprived of iron, scoring how well each line coped with the stress at the seedling stage. From this diversity panel they selected two contrasting parents: RA23, which tolerated iron deficiency with relative ease, and Lalat MAS, which succumbed more readily. Crossing these two lines produced an F2 mapping population, and the extreme iron-deficiency phenotypes among the F2 plants were pooled into tolerant and susceptible bulks for bulked segregant analysis by sequencing, better known as QTL-seq.

QTL-seq is a powerful shortcut for trait mapping. Instead of genotyping every individual plant, researchers sequence the DNA of the extreme bulks and look for regions of the genome where allele frequencies have shifted dramatically between the tolerant and susceptible pools. In this study, two statistical approaches were applied in parallel: the ΔSNP-index method, which measures differences in the proportion of parental alleles across sliding genomic windows, and the G′ statistic, an independent test of allelic differentiation. Both methods converged on the same two genomic regions, quantitative trait loci designated qFe2.1 on chromosome 2 and qFe9.1 on chromosome 9, giving the team high confidence that these intervals genuinely harbor genes influencing iron-deficiency response.

Locating a QTL, however, is only half the battle. Traditional mapping can leave dozens of genes within a candidate interval, and distinguishing the causal gene from bystanders requires an independent line of evidence. The researchers therefore turned to RNA sequencing of the two parental lines grown under iron-deficient conditions, asking which genes within the mapped intervals changed their expression when iron became scarce. This transcriptomic filter narrowed the two QTL regions down to just 14 candidate genes, a manageable list that could be scrutinized one by one.

One gene stood out from the pack with remarkable clarity. OsFe2200, catalogued as Os02g0152200 and LOC_Os02g05810, belongs to the DUF1677 family, a group of proteins whose name, Domain of Unknown Function 1677, signals that their biochemical role has remained obscure. When the researchers compared expression between the two parents, they found that OsFe2200 was expressed roughly 80-fold higher in the iron-deficiency-susceptible parent Lalat MAS than in the tolerant RA23. Quantitative reverse-transcription PCR confirmed the RNA-seq result, validating one of the strongest expression contrasts reported for a candidate gene in this kind of study. Intriguingly, the direction of the effect suggests that susceptibility may be associated with overexpression of this gene, or that the tolerant parent has evolved regulatory mechanisms that keep it quiet under stress.

To understand what kind of protein OsFe2200 encodes, the team turned to structural prediction. Sequence-based searches had yielded little, but structural modeling revealed that the protein adopts a putative Kelch β-propeller-like fold. β-propellers are elegant molecular architectures in which repeated blade-like segments are arranged radially around a central axis, resembling the blades of a propeller. Kelch-repeat propellers in particular are known throughout biology as platforms for protein-protein interactions, scaffolding complexes that range from cytoskeletal regulators to plant hormone signaling components. Further structural comparisons identified β-propeller and Kelch-like proteins as moderate-confidence homologs, leading the authors to describe OsFe2200 as structurally conserved but sequence divergent, a protein whose shape has been preserved over evolution even as its amino acid sequence has drifted beyond easy recognition.

The comparison between the two parents also revealed differences within the protein itself. Pairwise alignment of the OsFe2200 sequences from RA23 and Lalat MAS identified several amino acid substitutions, including a D90V change shared by the tolerant haplotype, along with R236M, N260D, V261M, A276T and K340N. Computational modeling with tools such as Dynamut2 suggested that these substitutions have variable and sometimes opposing effects on local protein flexibility. Notably, three of the mutations, R236M, N260D and V261M, affect a common structural region, while A276T and K340N induce opposite flexibility changes in a nearby region, a pattern consistent with compensatory effects within the combined haplotype. Whether these subtle structural shifts alter protein function remains an open question, but they provide concrete hypotheses for future experiments.

Regulation may matter as much as protein sequence. Scanning the promoter region upstream of OsFe2200, the researchers found iron-responsive motifs, including a canonical E-box (CACGTG) that is conserved in both parents, indicating that differential expression between the lines is unlikely to stem from loss of this core binding site. More interestingly, the tolerant parent RA23 carries a 21-base-pair insertion in its promoter that is absent from Lalat MAS. The authors are careful to note that this insertion has not yet been shown directly to cause the lower expression seen in RA23, but the correlation between the insertion, the promoter architecture and the expression contrast makes it a prime suspect for follow-up functional testing.

To see whether OsFe2200 variation matters beyond the two parents, the team performed haplotype analysis in the Bengal and Assam Aus Panel, a diverse collection of aus rice landraces from eastern South Asia. The analysis showed greater phenotypic variability associated with OsFe2200, and the haplotype carrying the same D90V mutation found in RA23 was superior and significantly different from others for the number of crown roots, an important architectural trait for foraging nutrients in poor soils. The effect on SPAD chlorophyll readings, a proxy for the leaf greening that iron deficiency erodes, was not statistically significant, hinting that the gene’s influence may be tissue-specific or trait-specific rather than a blanket tolerance effect.

The study stops short of proving causation, and the authors are explicit that OsFe2200 is a strong putative candidate rather than a validated tolerance gene. Functional studies, such as knockout mutants, overexpression lines and promoter-swap experiments, will be needed to confirm its role and to dissect whether the 21-bp promoter insertion, the D90V substitution or other variants drive the phenotype. Even so, the work demonstrates the power of integrating QTL-seq with transcriptomics and structural prediction to drag a gene out of functional obscurity and into the spotlight of crop genetics. For rice breeders working on direct-seeded and aerobic systems, where iron deficiency is a growing constraint as water scarcity pushes farmers away from flooded paddies, a confirmed iron-deficiency tolerance gene would be a valuable tool for marker-assisted selection. For biologists, OsFe2200 offers a tantalizing glimpse of how a protein family defined by what we do not know about it may turn out to hold answers to one of agriculture’s most stubborn soil problems.

Subject of Research: Genetic basis of iron-deficiency tolerance in rice seedlings

Article Title: Integrated QTL-seq and transcriptome analyses identify a DUF1677 candidate gene with a Predicted Kelch β-propeller-like fold associated with iron deficiency response in rice

Article References: Panda, S., Banerjee, N., Sinha, A., Fayaz, M., Sowdhamini, R., & Anandan, A. (2026). Integrated QTL-seq and transcriptome analyses identify a DUF1677 candidate gene with a Predicted Kelch β-propeller-like fold associated with iron deficiency response in rice. BMC Plant Biology. https://doi.org/10.1186/s12870-026-10001-1

Image Credits: AI Generated

DOI: 10.1186/s12870-026-10001-1

Keywords: rice, iron deficiency, QTL-seq, transcriptome, DUF1677, Kelch β-propeller, OsFe2200, haplotype analysis, direct-seeded rice, aerobic rice, plant genetics, candidate gene