Integration of QTL Mapping and Whole Genome Sequencing Identifies Candidate Genes for Alkalinity Tolerance in Rice (Oryza sativa)

Soil alkalinity is an important stressor that impairs crop growth and development, resulting in reduced crop productivity. Unlike salinity stress, research efforts to understand the mechanism of plant adaptation to alkaline stress is limited in rice, a major staple food for the world population. We evaluated a population of 193 recombinant inbred lines (RIL) developed from a cross between Cocodrie and N22 under alkaline stress at the seedling stage. Using a linkage map consisting of 4849 SNP markers, 42 additive QTLs were identified. There were seven genomic regions where two or more QTLs for multiple traits colocalized. Three important QTL clusters were targeted, and several candidate genes were identified based on high impact variants using whole genome sequences (WGS) of both parents and differential expression in response to alkalinity stress. These genes included two expressed protein genes, the glucan endo-1,3-beta-glucosidase precursor, F-box domain-containing proteins, double-stranded RNA-binding motif-containing protein, aquaporin protein, receptor kinase-like protein, semialdehyde hydrogenase, and NAD-binding domain-containing protein genes. Tolerance to alkaline stress in Cocodrie was most likely due to the low Na+/K+ ratio resulting from reduced accumulation of Na+ ions and higher accumulation of K+ in roots and shoots. Our study demonstrated the utility of integrating QTL mapping with WGS to identify the candidate genes in the QTL regions. The QTLs and candidate genes originating from the tolerant parent Cocodrie should be targeted for introgression to improve alkalinity tolerance in rice and to elucidate the molecular basis of alkali tolerance.