Meta-analysis of identified genomic regions and candidate genes underlying salinity tolerance in rice (Oryza sativa L.)

Rice output has grown globally, yet abiotic factors are still a key cause for worry. Salinity stress seems to have the more impact on crop production out of all abiotic stresses. Currently one of the most significant challenges in paddy breeding for salinity tolerance with the help of QTLs, is to determine the QTLs having the best chance of improving salinity tolerance with the least amount of background noise from the tolerant parent. Minimizing the size of the QTL confidence interval (CI) is essential in order to primarily include the genes responsible for salinity stress tolerance. By considering that, a genome-wide meta-QTL analysis on 768 QTLs from 35 rice populations published from 2001 to 2022 was conducted to identify consensus regions and the candidate genes underlying those regions responsible for the salinity tolerance, as it reduces the confidence interval (CI) to many folds from the initial QTL studies. In the present investigation, a total of 65 MQTLs were extracted with an average CI reduced from 17.35 to 1.66 cM including the smallest of 0.01 cM. Identification of the MQTLs for individual traits and then classifying the target traits into correlated morphological, physiological and biochemical aspects, resulted in more efficient interpretation of the salinity tolerance, identifying the candidate genes and to understand the salinity tolerance mechanism as a whole. The results of this study have a huge potential to improve the rice genotypes for salinity tolerance with the help of MAS and MABC.


Initial studies of QTL mapping for MQTL analysis
An exhaustive preview of the QTL mapping studies was carried out for the seedling stage and reproductive stage salinity tolerance in rice from 2001 to 2021 (Fig. 1A).Thirty-four studies having thirty-five independent mapping populations with all the information on mapping population, its size, molecular markers, genetic map, LOD score and the phenotypic variance explained (R 2 ) were selected for the MQTL analysis.The studies with the missing information were excluded from the analysis.The mapping populations used in the studies were generated from the cross combinations of 51 susceptible and tolerant parents and the Genome Wide Association Mapping Study (GWAS) included 179 landraces evaluated under saline conditions.The size of the mapping population ranges from 62 to 281 and the molecular markers used for generating the genetic map were ranges from 29 to 56,897.The studies included various types of genotyping platforms i.e., marker-based linkage mapping, Array-based, GWAS and GBS which utilized various molecular markers viz., RFLP, AFLP, InDel, ESTs, SSRs and SNPs.A total of 768 QTLs were extracted for 71 different traits correlated with salinity tolerance at seedling and reproductive stages viz., 95 for salt evaluation score (SEC), 68 for root dry weight (RDW), 51 for shoot length (SL), 50 for shoot K + concentration (SKC), 40 for shoot Na + concentration (SNC) etc.The traits related to salinity tolerance were regrouped into morphological attributes and ionic concentrations (Na + , K + and Na + /K + ratio) for which combined and separate analysis was conducted.Detailed information regarding the parents of the mapping populations, population type, population size, marker used, marker type and the traits under study are represented in Tables 1 and 2.
The successful integration of the Genome-Wide Association Study (GWAS) into the meta-analysis was achieved through a meticulous process involving the conversion of base pair distances into centimorgans (cM), facilitating linkage map preparation.Subsequently, the p-values, serving as indicators of the significance of marker-trait associations, underwent a transformation into logarithm of odds (LOD) scores using the following equation 25 .
Additionally, to ascertain the representative position of each distinct QTL, the mean position was computed based on the flanking markers, offering an accurate characterization of their genomic location within the metaanalysis framework.Through this comprehensive approach, all the essential information necessary for the metaanalysis was successfully amassed.
The confidence interval (95%) for each of the QTL was calculated as the difference between the left and right positions of the QTLs.Wherever the positions were not given it was calculated by the formula CI = X/(N × R 2 ) , where the X is 530 for BC and F 2 , 287 for DH lines and 163 for RILs and ILs, N is the size of the population and R 2 is the phenotypic variance explained 26 .

Preparation of consensus map and projection of QTLs
The consensus map was prepared using the genetic map information from all the 35 maps and the reference maps retrieved from the supplementary tables 17 and 18 of the IRGSP, 2005 (https:// archi ve.grame ne.org), Orjuela et al.  2010 and Cornell SSR 2001 27 .The reference maps are having marker density of 8740 markers with an average map of 1527.22 cM and the average length of each chromosome is 127.27 cM.The consensus genetic map was prepared with the help of the LPmerge (Ver.1.7) software that uses linear programming with ensuring that the order of the markers is preserved as in the linkage maps 28 .All the 768 QTLs were successfully projected on the consensus map as per the information available viz., CIs, LOD score R 2 , peak positions and flanking positions with the help of the software BioMercator V4.2.3 22 .LOD = − log 10 (P)   22 (Fig. 2).

Identification and functional analysis of candidate genes for salinity tolerance
The MQTLs selected for the functional analysis for candidate gene identification based on the criteria of greater than 5 candidate QTLs with an average LOD score of 3 and phenotypic variance of more than 8 per cent.The region of the MQTL for functional annotation is calculated by adding and subtracting the CI on either side of the identified MQTL and the nearest marker to that region is selected.The marker position is retrieved from the GWAS study in base pairs and used as the start and endpoint of the MQTL region 30 .The MQTL region which does not have the flanking SNPs were physically located using the Rice Annotation Project Database (RAPDB) (http:// rapdb.dna.affrc.go.jp), or the sequence of forward and reverse primer is subjected to nucleotide blast in NCBI (https:// blast.ncbi.nlm.nih.gov/ Blast.cgi? PROGR AM= blastn) to detect the sequence in the Nipponbare reference genome.The candidate genes underlying that region was investigated on the genome (IRGSP-1.0)using "Biomart" of the Ensemble plants database and the Gramene (https:// archi ve.grame ne.org/ qtl/).The retrieved GO terms are enriched by removing the redundant terms and summarized with the help of Revigo (http:// revigo.irb.hr/).After identification of candidate genes in MQTL region, FunRiceGenens (https:// funri cegen es.github.io/) and Ensamble Plants (https:// plants.ensem bl.org/ index.html) databases were searched to functionally annotate these genes.All the aforementioned information about the consensus map, QTL density, MQTLs identified for salinity tolerance and the candidate genes were shown in the circus plot (Fig. 3).www.nature.com/scientificreports/

Identification of orthologous genes in other cereals
The candidate genes identified in our study and had the previous references were screened to extract out the orthologous genes in the wheat, maize and sorghum.The candidate gene stable ids were searched against the genomes, Triticum aestivum (IWGSC), Zea mays (Zm-B73-REFERENCE-NAM-5.0) and Sorghum bicolor (Sor-ghum_bicolor_NCBIv3) in "Biomart" of the Ensamble plants database to identify the orthologous genes.The gene stable id, chromosome number, gene start and end were tabulated.The GO term and its definition was extracted from the respective databases of the species for functional analysis and based on the GO terms, the genes responsible for salinity tolerance were selected for searching the previous references for salinity tolerance.

Distribution of QTLs for salinity tolerance
A total of 768 relevant QTLs with the traits that have a significant correlation with salinity tolerance at seedling as well as reproductive stages in rice were selected to discover the associated consensus genomic regions (Supplementary file 1).The QTLs were retrieved from 34 studies having 35 independent mapping populations including RILs (12 case studies), F 2 (10), ILs (9), F 2:3 (2), BC (2) and landraces (1) with population sizes ranging from 61 to 281 reported from 2001 to 2021.The QTLs are identified for the 71 different trait hoods and then they are regrouped into three classes like, traits related to ionic stresses (Na + and K + ), traits associated with morphological attributes and combined analysis of salinity tolerance as a whole.The information of QTLs for all the studied traits and their distribution on the chromosomes are shown in Fig. 1A and B. Highest number of QTLs (148 QTLs) were observed on chromosome 1 followed by chromosome 2 (81 QTLs) while the lowest number of QTLs harbored by chromosome 10 (39 QTLs).Among the traits studied for salinity tolerance, salt evaluation score (SES) having the maximum number of QTLs (95 QTLs) followed by seedling dry weight (SDW, 68 QTLs) and Shoot potassium content (SK, 48 QTLs).The phenotypic variance of the QTLs studied ranged from 1.3% to 81.56% with an average of 13.88% and the LOD score ranged from 1.92 to 32 with an average of 4.43.The confidence interval for the QTLs studied was ranging from 0.3 to 89.4 cM with an average of 17.35 cM (Fig. 1C).

The consensus map of rice
The consensus map prepared from the reference maps and genetic maps of the studies has the marker density of 57,806 with an average of 4817 markers per chromosome and the genetic map of 3041.19 cM long with an average 253.43 cM region per chromosome (Supplementary file 2).The size of the individual chromosome ranges from 153.64 cM (chromosome 11) to 408 cM (chromosome 1) and the number of markers ranges from 3220 (chromosome 10) to 7272 (chromosome 1) (Fig. 1D).The average marker density of the consensus map was 19 markers per cM with the highest on chromosome 11 (30.69 markers per cM) and the lowest on chromosome 12 (14.11markers per cM).

MQTLs identified for salinity tolerance traits
All of the 768 QTLs were projected on the consensus map and from it the QTLs projected outside the chromosomal region were excluded from the further analysis.Also, the QTLs lacking the appropriate marker information were excluded.www.nature.com/scientificreports/MQTLs) while the lowest was on chromosome 7 (3 MQTLs).Correlation between the number of initially identified QTLs and the MQTLs density was low (0.51) suggesting MQTL identification is independent of the prior number and they are based on the CI, LOD and the phenotypic variance explained (Fig. 1B).A higher number of consensus groups does not necessarily imply greater stability of QTL.Instead, the stability of QTL is determined by various factors and should not be solely based on the number of consensus groups.The MQTL1.2 is derived from 61 different populations that have the highest number of candidate QTLs (75), followed by MQTL1.7 and MQTL1.5 (35) which may prove robust and stable in different environments and years.www.nature.com/scientificreports/Highest number of MQTLs was identified for the trait salt evaluation score (SES, 46 MQTLs) with an average of 2.06 SES QTLs per MQTL.Highest number of SES QTLs found on MQTL was 1.2 (9 QTLs) while the lowest for MQTL was 2.3 (1 QTL).For the traits seedling fresh and dry weight (SDW and SFW) the MQTLs, MQTL1.1 (8 QTLs, CI: 0.73 cM), MQTL1.2 (4 QTLs, CI: 0.1 cM), MQTL1.7 (4 QTLs, CI: 0.23 cM), MQTL2.1 (4 QTLs, CI: 0.47 cM), MQTL4.3 (4 QTLs, CI: 1.17 cM) and MQTL6.5 (4 QTLs, 2.13 cM) were promising.MQTL1.7 (9 QTLs, CI: 0.23 cM), MQTL1.2 (6 QTLs, CI: 0.1 cM) and MQTL1.6 (3 QTLs, 1.35 cM) were found associated with the shoot length of the seedling under saline conditions.The root length was positively correlated with the MQTL1.6 (CI: 1.35 cM), MQTL3.1 (CI: 1.22), MQTL6.6 (CI: 1.81) and MQTL9.2 (CI: 0.09 cM).Significant MQTLs for ionic concentrations and their ratio in the shoot were MQTL1.2(12 QTLs), MQTL1.5 (7 QTLs), MQTL2.6 (6 QTLs) and MQTL1.1 (4 QTLs).While for ionic concentrations in the root, MQTL1.2 (6 QTLs), MQTL1.5 (6 QTLs), MQTL4.4 (4 QTLs) and MQTL9.4 (4 QTLs) were responsible with the CI of 0.1 cM, 0.14 cM, 0.16 cM and 5.2 cM, respectively (Fig. 1D).As the tolerance to salinity of the plants against salt stress involves complex interconnecting signaling pathways that ultimately lead to the development of measurable phenotypes.The seedling and reproductive stage salinity traits are the results of these pathways.The MQTLs explained further are based on the mechanisms of the salinity tolerance for which we have identified consensus genomic regions and the candidate genes that ultimately complete the network (Fig. 2).

Identification of candidate genes for MQTLs
Exploring all the 64 MQTL regions for candidate genes resulted in the identification of 6973 candidate genes spanning all the chromosomes (Supplementary file 3).The highest number of candidate genes were identified on chromosome 6 (689 genes) while the lowest was on chromosome 11 (365 genes).The enrichment of the identified candidate genes enabled to filter them into three classes i.e. biological processes, molecular function and cellular components.Among them, 553 candidate genes were identified based on their potential role in salinity tolerance in rice.After screening the genes from the previous references, 56 candidate genes were identified which express in shoot, root and reproductive parts imparting salinity tolerance at seedling and reproductive stages (Table 4).While, the remaining 497 unique genes may have a potential role in maintaining osmotic and ionic balance during salinity stress.In lateral roots, the vascular bundle, the water pore and the basal section of seedling shoots, OsNHX1 or OsNHX5 promoter activity was seen.Salt stress, hyperosmotic stress and ABA increase the expression of OsNHX1, OsNHX2, OsNHX3 and OsNHX5 in rice tissues in distinct ways

MQTLs and candidate genes for ionic stresses
The saltol is the most reviewed major QTL with 43.7% phenotypic variance for Na + /K + homeostasis and is responsible for the salinity tolerance at the seedling stage 43 .This region is identified from the landrace Pokkali showing high tolerance to salinity and also a similar region was identified from the landrace Nona Bokara.In this work, MQTL1.1 on chromosome 1 between 11.23 Mb and 12.15 Mb was discovered which is in the same region as saltol and SKC, with a modest confidence interval of 0.92 Mb and 28 candidate QTLs.Furthermore, compared to the initial identification of that locus, this found area has lowered the confidence interval from 32.7 to 0.73 cM.In previous investigations, the candidate genes Os01g0303600 (OsRFP), Os01g0304100 (OsCCC2), Os01g0307500 (OsHKT1;5, SKC1, OsHKT8) and Os01g0309800 (OsHypB) from this area validated their significance in salinity tolerance.Because the OsHKT and OsSKC genes are responsible for Na + transporters located in the plasmalemma, they move sodium from the xylem to the parenchyma cells, limiting Na + buildup in the shoot 44 .SOS2, a major component of the SOS pathway is required to maintain intracellular Na + and K + homeostasis during the salinity stress and is found in MQTL6.3 on chromosome 6.SOS2 was found to be working for both seedling and reproductive stage salinity tolerance 45 .A potassium transporter, OsHAK present on MQTL3.4 and MQTL9.3 is involved in K + uptake and translocation involved in ion homeostasis 46 .Under salinity stress, OsHAK21 is required to maintain Na + /K + homeostasis and support seed germination and seedling establishment 47 .The NHX genes that are Na + /H + antiporters associated with tissue tolerance are found in MQTL 9.3 with 14 candidate QTLs 48 .Many of the chloride transporters were found to be active as a response to the salinity stress.Among them, OsCCC2 (MQTL1.1)has significance in developmental processes and Cl -ion homeostasis 49 .Higher accumulation of Cl -ions in the cytosol is limited by sCLC-2 (MQTL2.3)by accumulating them into vacuoles 50 .

MQTLs and candidate genes for reproductive stage salinity traits
At reproductive stage salinity tolerance, the major effects of salt are observed in pollen development and maturity which reduces the number of fertile spikelets and ultimately the grain yield.The MQTL4.1, MQTL6.3 and MQTL10.2 on chromosomes 4, 6 and 10 were deduced from 5, 10 and 6 initial QTLs, respectively.They collectively harbored 27 screened genes including OsSCP and OsSPX1 among which the OsSCP is responsible for  www.nature.com/scientificreports/pollen maturation and creation of pollen tube while, OsSPX is responsible for glucose metabolism and sugar transport during the pollen development 35,51 .The MQTL9.4 on chromosome 9 had the members of gene OsJMJ-C having higher expression in flag leaf during salinity stress 52 .The MQTL6.7 was found to be linked with the grain yield per plant under salt stress conditions including nine initial QTLs having the gene MOC1 which is one of the important regulators of tiller number 53 .MQTL6.3 with 10 candidate QTLs was found to be responsible for panicle number, 1000 grain weight, biomass content and grain yield per plant.OsWRKY, PE-1 and OsSOS2 are candidate genes (MQTL6.3)discovered in roots, stems, leaves and leaf sheaths that have been linked to lateral root development, elongation, improved photosynthesis and reproductive stage salinity tolerance 35,54 .

Orthologous genes identified in other cereals
After exploring the Wheat, Maize and Sorghum databases, 48, 23 and 47 orthologues were identified, respectively.Among them, only six genes were having previous references against salinity tolerance and the others are new potential genes that may have a role in salinity tolerance.Wheat gene TraesCS7D02G374400 (chromosome 7D), an orthologue of Os06g0606000 (Rice chromosome 6) is responsible for sensing and signalling to osmotic stress (Yue et al. 2021) 55 .Zm00001eb197200_T003 (chromosome 4), a rice orthologue (Os01g0304100) of maize has a role in maintaining osmotic balance during salinity stress.While in sorghum, an orthologue SORBI_3001G389700 on chromosome 1 (Os03g0320600) was having a role in response to osmotic stress.For potassium ion transmembrane transporter activity, orthologous gene in maize, Zm00001eb016160 on chromosome 1 (Os03g0337500) was found responsible for the transfer of potassium ions (K + ) from one side of a membrane to the other.While in wheat and sorghum, the genes TraesCS4A02G136300 (chromosome 4A, Os03g0337500) and SORBI_3001G379900 (chromosome 1, Os03g0337500) were found to be associated with potassium ion transport 56 .The gene MOC1 for active tillering responsible for active meristem initiation and secondary shoot formation was found in wheat and rice but not in maize and sorghum 57 .The list of orthologous genes, its GO term and GO definition is given in supplementary files 4, 5 and 6.

Discussion
The MQTL analysis identifies the true QTLs from the initial QTL studies that had unevenly distributed and varying genetic regions associated with salinity tolerance.Detection of the MQTLs was found to be independent of the genetic background and with very low environmental effects, as they are identified from a large number of different mapping populations (35 populations) developed from the paired crossing of 51 different genotypes.These populations were taken from the studies (2001-2022) covering 22 years of independent phenotyping, reducing its environmental effects in MQTL analysis.Along with the phenotypic screening, the marker density of the genetic map is also a limiting factor to identify the true QTLs with minimum confidence interval (CI).The consensus map developed in this study is much more informative as it includes the genetic maps from numerous studies with diverse molecular markers including SNPs with the marker density as high as 56,897 (GBS study) 14,26 .As a result, MQTLs were identified with the average confidence interval reduced from 17.35 to 1.66 cM including the smallest of 0.01 cM.To identify MQTLs, association mapping study was also included and compared with the derived MQTLs, which helped to locate the true physical regions (Mb) from the genetic ones (cM).Therefore, the MQTL analysis is the best reliable method to identify the true genomic regions underlying salinity tolerance and can be used very efficiently in marker-assisted backcrossing (MABC).MQTL analysis reduced initial 768 QTLs into 65 MQTLs present on all of the twelve chromosomes, indicating the power of MQTL analysis in narrowing down the genomic regions controlling salinity tolerance.The MQTL1.2 harboured as high as 75 candidate QTLs suggesting the robustness of the method to identify consensus regions.Here, 33 out of 65 MQTLs (CI < 1 cM) were identified with minimum linkage drag qualified as breeders QTLs, have potential use in MABC programmes.The MQTLs analysis ran separately for morphological traits, traits related to ionic concentrations and all the traits combined to ensure that, all the major regions responsible for salinity tolerance were captured.However, assuming that a lot of traits investigated are pleiotropically associated, it is more effective and robust to pool diverse correlated traits reported in the same population 21 .The analysis and interpretation of the target traits as a whole were more potent than analyzing them individually, which allows to identify the hotspots for salinity tolerance for osmotic and ionic stress that gives a large number of phenotypic effects on different traits 14 .In order to transfer the salinity tolerance, one must have knowledge about the region responsible for osmotic stress, ionic stress or a combination of both.Our approach to identify the MQTLs for individual traits and then classifying the target traits into correlated morphological, physiological and biochemical aspects, resulted in more efficient interpretation of the salinity tolerance, identifying the candidate genes and to understand the salinity tolerance mechanism as a whole.
The identification of various sodium, potassium and chloride transporters for ionic stress tolerance are of great importance as they maintain the Na + and Cl content below the damaging level in the cytosol.The identified MQTLs explicitly described the regions responsible for ion homeostasis mechanisms.The MQTL1.1 region (11.23-12.15Mb) responsible for Na + /K + homeostasis was identified as a comparable region to saltol and was narrowed down from 10.3-15.3Mb 24,59 .The important candidate genes OsHKT1,5 (K + transporter), SKC1 (K + transporter) and OsHKT8 (Na + /K + symporter) were all present in the MQTLs, demonstrating that the essential regions were not lost even when the confidence interval was reduced.The Na + /K + symporter or Na + uniporter (OsHKT, MQTL1.1) is present in the plasma membrane and expressed in xylem parenchyma limiting the Na + accumulation in photosynthetic parts of the plant.The identified Na + /H antiporter OsNHX (MQTL9.1)enhances the salt tolerance by compartmentalizing the Na + into the vocules 60 .High-affinity K + transporter (OsHKT, MQTL1.1) and the high-affinity K + uptake (OsHAK, MQTL9.3)mediate the cytosolic spike in potassium ions to regulate the homeostasis in plants under salt stress.OsHAK is also to be responsible for the accumulation of Na + in old leaves.OsCLC (Cl-channel) present on vocule, responsible for Cl sequestration to vocules and expressed only in the roots, nodes, internodes and leaf sheaths were identified in MQTL2.3 50.These are the genes and MQTLs significantly associated with the salinity tolerance and have the potential to use in marker-assisted selection to improve tolerance to salinity.For the reproductive stage salinity tolerance, highly correlated traits for salinity tolerance are pollen viability and stigma receptivity to increased Na + concentration in floral parts 61 .We identified the important genes for pollen and anther development in MQTL4.1(OsSCP) and MQTL6.3 (OsSPX).The MQTL area discovered (MQTL6.3)for grain yield per plant with a CI of 0.54 Mb in the study shows a comparable region from the previous study 30 .
The rice orthologous gene of OsSPX in wheat (TraesCS7D02G374400) located on chromosome 7D is responsible for maintaining osmotic balance under salinity conditions also identified in comparative transcriptomics study 55 .The identified orthologous for OsHAK in the wheat (TraesCS4A02G136300) chromosome 4A, maize (Zm00001eb016160) chromosome 1 and sorghum (SORBI_3001G379900) chromosome 1 were also reported previously as an important potassium transporter of the plasma membrane 56,62 .These potential candidate genes have a crucial role in tolerance to salinity as they have been found in four of the major cereal crops of the Poaceae family and potential to use in MAS to improve the genotypes for salinity tolerance.
To the best of our knowledge, this is the first comprehensive study to report the MQTLs with very small confidence intervals including the major candidate genes identified for the salinity tolerance and also being previously reported.These results have a huge possibility to improve the rice genotypes for salinity tolerance.

Conclusion
In the present investigation, a total of 65 MQTLs were extracted with an average CI reduced from 17.35 to 1.66 cM including the smallest of 0.01 cM.Identification of the MQTLs for individual traits and then classifying the target traits into correlated morphological, physiological and biochemical aspects, resulted in more efficient interpretation of the salinity tolerance, identifying the candidate genes and to understand the salinity tolerance mechanism as a whole.The results of this study have a huge potential to improve the rice genotypes for salinity tolerance with the help of MAS and MABC.

Figure 1 .
Figure 1.(A) QTLs retrieved from 35 independent mapping populations and used for MQTL analysis for salinity tolerance at seedling and reproductive stages (SES: Salt evaluation score, SDW: Shoot dry weight, SL: Shoot length, SK: Shoot K + concentration, SNA: Shoot Na + concentration, SFW: Shoot fresh weight, RL: Root length, SNAK: Shoot Na + /K + concentration, RDW: Root dry weight, CHL: Chlorophyll content, RNA: Root Na + concentration, SDS: Survival days, GP: Germination per cent, IR: Imbibition rate, RK: Root K + concentration, BM: Biomass, PH: Plant height, RNAK: Root Na + /K + concentration, PF: Pollen fertility, RLSDW: Relative shoot dry weight, SH: Shoot height, STE: Spikelet sterility, FNAK: Flag leaf Na + /K + concentration, RTDW: Relative total dry weight, SFD: Shoot fresh/dry weight, TSP: Total spikelet).(B) Chromosome wise QTLs and MQTLs distribution on twelve chromosomes of rice.(C) Frequency distribution of the QTLs based on confidence intervals for seedling and reproductive stage salinity tolerance.(D) Frequency distribution of the MQTLs based on confidence intervals for seedling and reproductive stage salinity tolerance.

Figure 3 .
Figure 3. Visualization of circos plot after MQTL analysis (A) Chromosome number (B) Marker density (C) LOD of the initial QTLs (D) R2 of initial QTLs (E) Confidence interval of initial QTLs (F) Number of candidate QTLs (G) Confidence interval of MQTLs (H) Genes identified in MQTL region (I) Homologous regions of the identified genes.

Table 1 .
Summary of the studies of QTL mapping used for MQTL analysis for salinity tolerance in rice.

Table 2 .
The traits for which The QTLs compiled.

Table 3 .
Summary of the MQTLs detected for salinity tolerance in rice.

Table 4 .
Summary of candidate genes identified from MQTLs for salinity tolerance in rice.