A new comprehensive annotation of leucine-rich repeat-containing receptors in rice.

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Assessment of the roles of SPO11-2 and SPO11-4 in meiosis in rice using CRISPR/Cas9 mutagenesis.

In Arabidopsis, chromosomal double-strand breaks at meiosis are presumably catalyzed by two distinct SPO11 transesterases, AtSPO11-1 and AtSPO11-2, together with M-TOPVIB. To clarify the roles of the SPO11 paralogs in rice, we used CRISPR/Cas9 mutagenesis to produce null biallelic mutants in OsSPO11-1, OsSPO11-2, and OsSPO11-4. Similar to Osspo11-1, biallelic mutations in the first exon of OsSPO11-2 led to complete panicle sterility. Conversely, all Osspo11-4 biallelic mutants were fertile. To generate segregating Osspo11-2 mutant lines, we developed a strategy based on dual intron targeting. Similar to Osspo11-1, the pollen mother cells of Osspo11-2 progeny plants showed an absence of bivalent formation at metaphase I, aberrant segregation of homologous chromosomes, and formation of non-viable tetrads. In contrast, the chromosome behavior in Osspo11-4 male meiocytes was indistinguishable from that in the wild type. While similar numbers of OsDMC1 foci were revealed by immunostaining in wild-type and Osspo11-4 prophase pollen mother cells (114 and 101, respectively), a surprisingly high number (85) of foci was observed in the sterile Osspo11-2 mutant, indicative of a divergent function between OsSPO11-1 and OsSPO11-2. This study demonstrates that whereas OsSPO11-1 and OsSPO11-2 are the likely orthologs of AtSPO11-1 and AtSPO11-2, OsSPO11-4 has no major role in wild-type rice meiosis.

Homeobox transcription factor OsZHD2 promotes root meristem activity in rice by inducing ethylene biosynthesis.

Root meristem activity is the most critical process influencing root development. Although several factors that regulate meristem activity have been identified in rice, studies on the enhancement of meristem activity in roots are limited. We identified a T-DNA activation tagging line of a zinc-finger homeobox gene, OsZHD2, which has longer seminal and lateral roots due to increased meristem activity. The phenotypes were confirmed in transgenic plants overexpressing OsZHD2. In addition, the overexpressing plants showed enhanced grain yield under low nutrient and paddy field conditions. OsZHD2 was preferentially expressed in the shoot apical meristem and root tips. Transcriptome analyses and quantitative real-time PCR experiments on roots from the activation tagging line and the wild type showed that genes for ethylene biosynthesis were up-regulated in the activation line. Ethylene levels were higher in the activation lines compared with the wild type. ChIP assay results suggested that OsZHD2 induces ethylene biosynthesis by controlling ACS5 directly. Treatment with ACC (1-aminocyclopropane-1-carboxylic acid), an ethylene precursor, induced the expression of the DR5 reporter at the root tip and stele, whereas treatment with an ethylene biosynthesis inhibitor, AVG (aminoethoxyvinylglycine), decreased that expression in both the wild type and the OsZHD2 overexpression line. These observations suggest that OsZHD2 enhances root meristem activity by influencing ethylene biosynthesis and, in turn, auxin.

Engineering meiotic recombination pathways in rice.

In the last 15 years, outstanding progress has been made in understanding the function of meiotic genes in the model dicot and monocot plants Arabidopsis and rice (Oryza sativa L.), respectively. This knowledge allowed to modulate meiotic recombination in Arabidopsis and, more recently, in rice. For instance, the overall frequency of crossovers (COs) has been stimulated 2.3- and 3.2-fold through the inactivation of the rice FANCM and RECQ4 DNA helicases, respectively, two genes involved in the repair of DNA double-strand breaks (DSBs) as noncrossovers (NCOs) of the Class II crossover pathway. Differently, the programmed induction of DSBs and COs at desired sites is currently explored by guiding the SPO11-1 topoisomerase-like transesterase, initiating meiotic recombination in all eukaryotes, to specific target regions of the rice genome. Furthermore, the inactivation of 3 meiosis-specific genes, namely PAIR1, OsREC8 and OsOSD1, in the Mitosis instead of Meiosis (MiMe) mutant turned rice meiosis into mitosis, thereby abolishing recombination and achieving the first component of apomixis, apomeiosis. The successful translation of Arabidopsis results into a crop further allowed the implementation of two breakthrough strategies that triggered parthenogenesis from the MiMe unreduced clonal egg cell and completed the second component of diplosporous apomixis. Here, we review the most recent advances in and future prospects of the manipulation of meiotic recombination in rice and potentially other major crops, all essential for global food security.

Production of low-Cs+ rice plants by inactivation of the K+ transporter OsHAK1 with the CRISPR-Cas system.

The occurrence of radiocesium in food has raised sharp health concerns after nuclear accidents. Despite being present at low concentrations in contaminated soils (below µm), cesium (Cs+ ) can be taken up by crops and transported to their edible parts. This plant capacity to take up Cs+ from low concentrations has notably affected the production of rice (Oryza sativa L.) in Japan after the nuclear accident at Fukushima in 2011. Several strategies have been put into practice to reduce Cs+ content in this crop species such as contaminated soil removal or adaptation of agricultural practices, including dedicated fertilizer management, with limited impact or pernicious side-effects. Conversely, the development of biotechnological approaches aimed at reducing Cs+ accumulation in rice remain challenging. Here, we show that inactivation of the Cs+ -permeable K+ transporter OsHAK1 with the CRISPR-Cas system dramatically reduced Cs+ uptake by rice plants. Cs+ uptake in rice roots and in transformed yeast cells that expressed OsHAK1 displayed very similar kinetics parameters. In rice, Cs+ uptake is dependent on two functional properties of OsHAK1: (i) a poor capacity of this system to discriminate between Cs+ and K+ ; and (ii) a high capacity to transport Cs+ from very low external concentrations that is likely to involve an active transport mechanism. In an experiment with a Fukushima soil highly contaminated with 137 Cs+ , plants lacking OsHAK1 function displayed strikingly reduced levels of 137 Cs+ in roots and shoots. These results open stimulating perspectives to smartly produce safe food in regions contaminated by nuclear accidents.

A plausible mechanism, based upon Short-Root movement, for regulating the number of cortex cell layers in roots.

Formation of specialized cells and tissues at defined times and in specific positions is essential for the development of multicellular organisms. Often this developmental precision is achieved through intercellular signaling networks, which establish patterns of differential gene expression and ultimately the specification of distinct cell fates. Here we address the question of how the Short-root (SHR) proteins from Arabidopsis thaliana (AtSHR), Brachypodium distachyon (BdSHR), and Oryza sativa (OsSHR1 and OsSHR2) function in patterning the root ground tissue. We find that all of the SHR proteins function as mobile signals in A. thaliana and all of the SHR homologs physically interact with the AtSHR binding protein, Scarecow (SCR). Unlike AtSHR, movement of the SHR homologs was not limited to the endodermis. Instead, the SHR proteins moved multiple cell layers and determined the number of cortex, not endodermal, cell layers formed in the root. Our results in A. thaliana are consistent with a mechanism by which the regulated movement of the SHR transcription factor determines the number of cortex cell layers produced in the roots of B. distachyon and O. sativa. These data also provide a new model for ground tissue patterning in A. thaliana in which the ability to form a functional endodermis is spatially limited independently of SHR.

GreenPhylDB v2.0: comparative and functional genomics in plants.

GreenPhylDB is a database designed for comparative and functional genomics based on complete genomes. Version 2 now contains sixteen full genomes of members of the plantae kingdom, ranging from algae to angiosperms, automatically clustered into gene families. Gene families are manually annotated and then analyzed phylogenetically in order to elucidate orthologous and paralogous relationships. The database offers various lists of gene families including plant, phylum and species specific gene families. For each gene cluster or gene family, easy access to gene composition, protein domains, publications, external links and orthologous gene predictions is provided. Web interfaces have been further developed to improve the navigation through information related to gene families. New analysis tools are also available, such as a gene family ontology browser that facilitates exploration. GreenPhylDB is a component of the South Green Bioinformatics Platform (http://southgreen.cirad.fr/) and is accessible at http://greenphyl.cirad.fr. It enables comparative genomics in a broad taxonomy context to enhance the understanding of evolutionary processes and thus tends to speed up gene discovery.

In-depth molecular and phenotypic characterization in a rice insertion line library facilitates gene identification through reverse and forward genetics approaches.

We report here the molecular and phenotypic features of a library of 31,562 insertion lines generated in the model japonica cultivar Nipponbare of rice (Oryza sativa L.), called Oryza Tag Line (OTL). Sixteen thousand eight hundred and fourteen T-DNA and 12,410 Tos17 discrete insertion sites have been characterized in these lines. We estimate that 8686 predicted gene intervals--i.e. one-fourth to one-fifth of the estimated rice nontransposable element gene complement--are interrupted by sequence-indexed T-DNA (6563 genes) and/or Tos17 (2755 genes) inserts. Six hundred and forty-three genes are interrupted by both T-DNA and Tos17 inserts. High quality of the sequence indexation of the T2 seed samples was ascertained by several approaches. Field evaluation under agronomic conditions of 27,832 OTL has revealed that 18.2% exhibit at least one morphophysiological alteration in the T1 progeny plants. Screening 10,000 lines for altered response to inoculation by the fungal pathogen Magnaporthe oryzae allowed to observe 71 lines (0.7%) developing spontaneous lesions simulating disease mutants and 43 lines (0.4%) exhibiting an enhanced disease resistance or susceptibility. We show here that at least 3.5% (four of 114) of these alterations are tagged by the mutagens. The presence of allelic series of sequence-indexed mutations in a gene among OTL that exhibit a convergent phenotype clearly increases the chance of establishing a linkage between alterations and inserts. This convergence approach is illustrated by the identification of the rice ortholog of AtPHO2, the disruption of which causes a lesion-mimic phenotype owing to an over-accumulation of phosphate, in nine lines bearing allelic insertions.

Leucine-rich repeat receptor kinases are sporadically distributed in eukaryotic genomes.

BACKGROUND: Plant leucine-rich repeat receptor-like kinases (LRR-RLKs) are receptor kinases that contain LRRs in their extracellular domain. In the last 15 years, many research groups have demonstrated major roles played by LRR-RLKs in plants during almost all developmental processes throughout the life of the plant and in defense/resistance against a large range of pathogens. Recently, a breakthrough has been made in this field that challenges the dogma of the specificity of plant LRR-RLKs. RESULTS: We analyzed ~1000 complete genomes and show that LRR-RK genes have now been identified in 8 non-plant genomes. We performed an exhaustive phylogenetic analysis of all of these receptors, revealing that all of the LRR-containing receptor subfamilies form lineage-specific clades. Our results suggest that the association of LRRs with RKs appeared independently at least four times in eukaryotic evolutionary history. Moreover, the molecular evolutionary history of the LRR-RKs found in oomycetes is reminiscent of the pattern observed in plants: expansion with amplification/deletion and evolution of the domain organization leading to the functional diversification of members of the gene family. Finally, the expression data suggest that oomycete LRR-RKs may play a role in several stages of the oomycete life cycle. CONCLUSIONS: In view of the key roles that LRR-RLKs play throughout the entire lifetime of plants and plant-environment interactions, the emergence and expansion of this type of receptor in several phyla along the evolution of eukaryotes, and particularly in oomycete genomes, questions their intrinsic functions in mimicry and/or in the coevolution of receptors between hosts and pathogens.

OryGenesDB 2008 update: database interoperability for functional genomics of rice.

OryGenesDB (http://orygenesdb.cirad.fr/index.html) is a database developed for rice reverse genetics. OryGenesDB contains FSTs (flanking sequence tags) of various mutagens and functional genomics data, collected from both international insertion collections and the literature. The current release of OryGenesDB contains 171,000 FSTs, and annotations divided among 10 specific categories, totaling 78 annotation layers. Several additional tools have been added to the main interface; these tools enable the user to retrieve FSTs and design probes to analyze insertion lines. The major innovation of OryGenesDB 2008, besides updating the data and tools, is a new tool, Orylink, which was developed to speed up rice functional genomics by taking advantage of the resources developed in two related databases, Oryza Tag Line and GreenPhylDB. Orylink was designed to field complex queries across these three databases and store both the queries and their results in an intuitive manner. Orylink offers a simple and powerful virtual workbench for functional genomics. Alternatively, the Web services developed for Orylink can be used independently of its Web interface, increasing the interoperability between these different bioinformatics applications.

Efficient Genome Editing in Rice Protoplasts Using CRISPR/CAS9 Construct.

Genome editing technologies, mainly CRISPR/CAS9, are revolutionizing plant biology and breeding. Since the demonstration of its effectiveness in eukaryotic cells, a very large number of derived technologies has emerged. Demonstrating and comparing the effectiveness of all these new technologies in entire plants is a long, tedious, and labor-intensive process that generally involves the production of transgenic plants and their analysis. Protoplasts, plant cells free of their walls, offer a simple, high-throughput system to test the efficiency of these editing technologies in a few weeks' time span. We have developed a routine protocol using protoplasts to test editing technologies in rice. Our protocol allows to test more than 30 constructs in protoplasts prepared from leaf tissues of 100, 9-11-day-old seedlings. CRISPR/CAS9 construct effectiveness can be clearly established within less than a week. We provide here a full protocol, from designing sgRNA to mutation analysis.

Complex regulation of two target genes encoding SPX-MFS proteins by rice miR827 in response to phosphate starvation.

Here we report on the characterization of rice osa-miR827 and its two target genes, OsSPX-MFS1 and OsSPX-MFS2, which encode SPX-MFS proteins predicted to be implicated in phosphate (Pi) sensing or transport. We first show by Northern blot analysis that osa-miR827 is strongly induced by Pi starvation in both shoots and roots. Hybridization of osa-miR827 in situ confirms its strong induction by Pi starvation, with signals concentrated in mesophyll, epidermis and ground tissues of roots. In parallel, we analyzed the responses of the two OsSPX-MFS1 and OsSPX-MFS2 gene targets to Pi starvation. OsSPX-MFS1 mRNA is mainly expressed in shoots under sufficient Pi supply while its expression is reduced on Pi starvation, revealing a direct relationship between induction of osa-miR827 and down-regulation of OsSPX-MFS1. In contrast, OsSPX-MFS2 responds in a diametrically opposed manner to Pi starvation. The accumulation of OsSPX-MFS2 mRNA is dramatically enhanced under Pi starvation, suggesting the involvement of complex regulation of osa-miR827 and its two target genes. We further produced transgenic rice lines overexpressing osa-miR827 and T-DNA knockout mutant lines in which the expression of osa-miR827 is abolished. Compared with wild-type controls, both target mRNAs exhibit similar changes, their expression being reduced and increased in overexpressing and knockout lines, respectively. This suggests that OsSPX-MFS1 and OsSPX-MFS2 are both negatively regulated by osa-miR827 abundance although they respond differently to external Pi conditions. We propose that this is a complex mechanism comprising fine tuning of spatial or temporal regulation of both targets by osa-miR827.

Oryza Tag Line, a phenotypic mutant database for the Genoplante rice insertion line library.

To organize data resulting from the phenotypic characterization of a library of 30,000 T-DNA enhancer trap (ET) insertion lines of rice (Oryza sativa L cv. Nipponbare), we developed the Oryza Tag Line (OTL) database (http://urgi.versailles.inra.fr/OryzaTagLine/). OTL structure facilitates forward genetic search for specific phenotypes, putatively resulting from gene disruption, and/or for GUSA or GFP reporter gene expression patterns, reflecting ET-mediated endogenous gene detection. In the latest version, OTL gathers the detailed morpho-physiological alterations observed during field evaluation and specific screens in a first set of 13,928 lines. Detection of GUS or GFP activity in specific organ/tissues in a subset of the library is also provided. Search in OTL can be achieved through trait ontology category, organ and/or developmental stage, keywords, expression of reporter gene in specific organ/tissue as well as line identification number. OTL now contains the description of 9721 mutant phenotypic traits observed in 2636 lines and 1234 GUS or GFP expression patterns. Each insertion line is documented through a generic passport data including production records, seed stocks and FST information. 8004 and 6101 of the 13,928 lines are characterized by at least one T-DNA and one Tos17 FST, respectively that OTL links to the rice genome browser OryGenesDB.

Origin and Diversity of Plant Receptor-Like Kinases.

Because of their high level of diversity and complex evolutionary histories, most studies on plant receptor-like kinase subfamilies have focused on their kinase domains. With the large amount of genome sequence data available today, particularly on basal land plants and Charophyta, more attention should be paid to primary events that shaped the diversity of the RLK gene family. We thus focus on the motifs and domains found in association with kinase domains to illustrate their origin, organization, and evolutionary dynamics. We discuss when these different domain associations first occurred and how they evolved, based on a literature review complemented by some of our unpublished results.

A fast, efficient and high-throughput procedure involving laser microdissection and RT droplet digital PCR for tissue-specific expression profiling of rice roots.

BACKGROUND: In rice, the cortex and outer tissues play a key role in submergence tolerance. The cortex differentiates into aerenchyma, which are air-containing cavities that allow the flow of oxygen from shoots to roots, whereas exodermis suberification and sclerenchyma lignification limit oxygen loss from the mature parts of roots by forming a barrier to root oxygen loss (ROL). The genes and their networks involved in the cellular identity and differentiation of these tissues remain poorly understood. Identification and characterization of key regulators of aerenchyma and ROL barrier formation require determination of the specific expression profiles of these tissues. RESULTS: We optimized an approach combining laser microdissection (LM) and droplet digital RT-PCR (ddRT-PCR) for high-throughput identification of tissue-specific expression profiles. The developed protocol enables rapid (within 3 days) extraction of high-quality RNA from root tissues with a low contamination rate. We also demonstrated the possibility of extracting RNAs from paraffin blocks stored at 4 °C without any loss of quality. We included a detailed troubleshooting guide that should allow future users to adapt the proposed protocol to other tissues and/or species. We demonstrated that our protocol, which combines LM with ddRT-PCR, can be used as a complementary tool to in situ hybridization for tissue-specific characterization of gene expression even with a low RNA concentration input. We illustrated the efficiency of the proposed approach by validating three of four potential tissue-specific candidate genes detailed in the RiceXpro database. CONCLUSION: The detailed protocol and the critical steps required to optimize its use for other species will democratize tissue-specific transcriptome approaches combining LM with ddRT-PCR for analyses of plants.

Beyond Seek and Destroy: how to Generate Allelic Series Using Genome Editing Tools.

Genome editing tools have greatly facilitated the functional analysis of genes of interest by targeted mutagenesis. Many usable genome editing tools, including different site-specific nucleases and editor databases that allow single-nucleotide polymorphisms (SNPs) to be introduced at a given site, are now available. These tools can be used to generate high allelic diversity at a given locus to facilitate gene function studies, including examining the role of a specific protein domain or a single amino acid. We compared the effects, efficiencies and mutation types generated by our LbCPF1, SpCAS9 and base editor (BECAS9) constructs for the OsCAO1 gene. SpCAS9 and LbCPF1 have similar efficiencies in generating mutations but differ in the types of mutations induced, with the majority of changes being single-nucleotide insertions and short deletions for SpCAS9 and LbCPF1, respectively. The proportions of heterozygotes also differed, representing a majority in our LbCPF1, while with SpCAS9, we obtained a large number of biallelic mutants. Finally, we demonstrated that it is possible to specifically introduce stop codons using the BECAS9 with an acceptable efficiency of approximately 20%. Based on these results, a rational choice among these three alternatives may be made depending on the type of mutation that one wishes to introduce, the three systems being complementary. SpCAS9 remains the best choice to generate KO mutations in primary transformants, while if the desired gene mutation interferes with regeneration or viability, the use of our LbCPF1 construction will be preferred, because it produces mainly heterozygotes. LbCPF1 has been described in other studies as being as effective as SpCAS9 in generating homozygous and biallelic mutations. It will remain to be clarified in the future, whether the different LbCFP1 constructions have different efficiencies and determine the origin of these differences. Finally, if one wishes to specifically introduce stop codons, BECAS9 is a viable and efficient alternative, although it has a lower efficiency than SpCAS9 and LbCPF1 for creating KO mutations.

Phylogenomics of plant genomes: a methodology for genome-wide searches for orthologs in plants.

BACKGROUND: Gene ortholog identification is now a major objective for mining the increasing amount of sequence data generated by complete or partial genome sequencing projects. Comparative and functional genomics urgently need a method for ortholog detection to reduce gene function inference and to aid in the identification of conserved or divergent genetic pathways between several species. As gene functions change during evolution, reconstructing the evolutionary history of genes should be a more accurate way to differentiate orthologs from paralogs. Phylogenomics takes into account phylogenetic information from high-throughput genome annotation and is the most straightforward way to infer orthologs. However, procedures for automatic detection of orthologs are still scarce and suffer from several limitations. RESULTS: We developed a procedure for ortholog prediction between Oryza sativa and Arabidopsis thaliana. Firstly, we established an efficient method to cluster A. thaliana and O. sativa full proteomes into gene families. Then, we developed an optimized phylogenomics pipeline for ortholog inference. We validated the full procedure using test sets of orthologs and paralogs to demonstrate that our method outperforms pairwise methods for ortholog predictions. CONCLUSION: Our procedure achieved a high level of accuracy in predicting ortholog and paralog relationships. Phylogenomic predictions for all validated gene families in both species were easily achieved and we can conclude that our methodology outperforms similarly based methods.

Modulating rice stress tolerance by transcription factors.

Plants are non-mobile organisms and have to adapt to environmental stresses mostly by modulating their growth and development in addition to physiological and biochemical changes. Transcription factors (TFs) regulate genome expression in response to environmental and physiological signals, and some of them switch on plant adaptive developmental and physiological pathways. One TF is encoded by a single gene but regulates the expression of several other genes leading to the activation of complex adaptive mechanisms and hence represents major molecular targets to genetically improve the tolerance of crop plants against different stresses. In this review an updated account of the discovery of TFs involved in biotic and abiotic stress tolerance in the model monocotyledonous plant, rice (Oryza sativa L.) is presented. We illustrate how the elucidation of the function of these TFs can be used to set up genetic engineering strategies and to rationalize molecular breeding using molecular assisted selection towards enhancement of rice tolerance to various stresses. Attempts have also been made to provide information on the molecular mechanisms involved in stress resistance or tolerance processes. We discuss how the comparison of the action of TFs isolated from the dicotyledonous model plant Arabidopsis thaliana in rice and vice versa can contribute to determine whether common or divergent mechanisms underlie stress tolerance in the two plant species. Lastly, we discuss the necessity to discover TFs controlling specifically the root adaptive development which constitutes a major way for the plant to escape to several stresses such as water deficit or mineral nutrient deficiency.

A protocol combining multiphoton microscopy and propidium iodide for deep 3D root meristem imaging in rice: application for the screening and identification of tissue-specific enhancer trap lines.

BACKGROUND: The clear visualization of 3D organization at the cellular level in plant tissues is needed to fully understand plant development processes. Imaging tools allow the visualization of the main fluorophores and in vivo growth monitoring. Confocal microscopy coupled with the use of propidium iodide (PI) counter-staining is one of the most popular tools used to characterize the structure of root meristems in A. thaliana. However, such an approach is relatively ineffective in species with more complex and thicker root systems. RESULTS: We adapted a PI counter-staining protocol to visualize the internal 3D architecture of rice root meristems using multiphoton microscopy. This protocol is simple and compatible with the main fluorophores (CFP, GFP and mCherry). The efficiency and applicability of this protocol were demonstrated by screening a population of 57 enhancer trap lines. We successfully characterized GFP expression in all of the lines and identified 5 lines with tissue-specific expression. CONCLUSIONS: All of these resources are now available for the rice community and represent critical tools for future studies of root development.

Tolerance to mild salinity stress in japonica rice: A genome-wide association mapping study highlights calcium signaling and metabolism genes.

Salinity tolerance is an important quality for European rice grown in river deltas. We evaluated the salinity tolerance of a panel of 235 temperate japonica rice accessions genotyped with 30,000 SNP markers. The panel was exposed to mild salt stress (50 mM NaCl; conductivity of 6 dS m-1) at the seedling stage. Eight different root and shoot growth parameters were measured for both the control and stressed treatments. The Na+ and K+ mass fractions of the stressed plants were measured using atomic absorption spectroscopy. The salt treatment affected plant growth, particularly the shoot parameters. The panel showed a wide range of Na+/K+ ratio and the temperate accessions were distributed over an increasing axis, from the most resistant to the most susceptible checks. We conducted a genome-wide association study on indices of stress response and ion mass fractions in the leaves using a classical mixed model controlling structure and kinship. A total of 27 QTLs validated by sub-sampling were identified. For indices of stress responses, we also used another model that focused on marker × treatment interactions and detected 50 QTLs, three of which were also identified using the classical method. We compared the positions of the significant QTLs to those of approximately 300 genes that play a role in rice salt tolerance. The positions of several QTLs were close to those of genes involved in calcium signaling and metabolism, while other QTLs were close to those of kinases. These results reveal the salinity tolerance of accessions with a temperate japonica background. Although the detected QTLs must be confirmed by other approaches, the number of associations linked to candidate genes involved in calcium-mediated ion homeostasis highlights pathways to explore in priority to understand the salinity tolerance of temperate rice.