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.

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.

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.

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.

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.

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.

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.

Sub-cellular markers highlight intracellular dynamics of membrane proteins in response to abiotic treatments in rice

BACKGROUND

Cell biology approach using membrane protein markers tagged with fluorescent proteins highlights the dynamic behaviour of plant cell membranes, not only in the standard but also in changing environmental conditions. In the past, this strategy has been extensively developed in plant models such as Arabidopsis.

RESULTS

Here, we generated a set of transgenic lines expressing membrane protein markers to extend this approach to rice, one of the most cultivated crop in the world and an emerging plant model. Lines expressing individually eight membrane protein markers including five aquaporins (OsPIP1;1, OsPIP2;4, OsPIP2;5, OsTIP1;1, OsTIP2;2) and three endosomal trafficking proteins (OsRab5a, OsGAP1, OsSCAMP1) were obtained. Importantly, we challenged in roots the aquaporin-expressing transgenic lines upon salt and osmotic stress and uncovered a highly dynamic behaviour of cell membrane.

CONCLUSION

We have uncovered the relocalization and dynamics of plasma membrane aquaporins upon salt and osmotic stresses in rice. Importantly, our data support a model where relocalization of OsPIPs is concomitant with their high cycling dynamics.

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.

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 ¹³⁷ Cs⁺ , plants lacking OsHAK1 function displayed strikingly reduced levels of ¹³⁷ Cs⁺ in roots and shoots. These results open stimulating perspectives to smartly produce safe food in regions contaminated by nuclear accidents.

Corrigendum: New Insights on Leucine-Rich Repeats Receptor-Like Kinase Orthologous Relationships in Angiosperms

[This corrects the article on p. 381 in vol. 8, PMID: 28424707.].

Dynamic Regulation of Auxin Response during Rice Development Revealed by Newly Established Hormone Biosensor Markers

The hormone auxin is critical for many plant developmental processes. Unlike the model eudicot plant Arabidopsis (Arabidopsis thaliana), auxin distribution and signaling in rice tissues has not been systematically investigated due to the absence of suitable auxin response reporters. In this study we observed the conservation of auxin signaling components between Arabidopsis and model monocot crop rice (Oryza sativa), and generated complementary types of auxin biosensor constructs, one derived from the Aux/IAA-based biosensor DII-VENUS but constitutively driven by maize ubiquitin-1 promoter, and the other termed DR5-VENUS in which a synthetic auxin-responsive promoter (DR5rev ) was used to drive expression of the yellow fluorescent protein (YFP). Using the obtained transgenic lines, we observed that during the vegetative development, accumulation of DR5-VENUS signal was at young and mature leaves, tiller buds and stem base. Notably, abundant DR5-VENUS signals were observed in the cytoplasm of cortex cells surrounding lateral root primordia (LRP) in rice. In addition, auxin maxima and dynamic re-localization were seen at the initiation sites of inflorescence and spikelet primordia including branch meristems (BMs), female and male organs. The comparison of these observations among Arabidopsis, rice and maize suggests the unique role of auxin in regulating rice lateral root emergence and reproduction. Moreover, protein localization of auxin transporters PIN1 homologs and GFP tagged OsAUX1 overlapped with DR5-VENUS during spikelet development, helping validate these auxin response reporters are reliable markers in rice. This work firstly reveals the direct correspondence between auxin distribution and rice reproductive and root development at tissue and cellular level, and provides high-resolution auxin tools to probe fundamental developmental processes in rice and to establish links between auxin, development and agronomical traits like yield or root architecture.

New Insights on Leucine-Rich Repeats Receptor-Like Kinase Orthologous Relationships in Angiosperms

Leucine-Rich Repeats Receptor-Like Kinase (LRR-RLK) genes represent a large and complex gene family in plants, mainly involved in development and stress responses. These receptors are composed of an LRR-containing extracellular domain (ECD), a transmembrane domain (TM) and an intracellular kinase domain (KD). To provide new perspectives on functional analyses of these genes in model and non-model plant species, we performed a phylogenetic analysis on 8,360 LRR-RLK receptors in 31 angiosperm genomes (8 monocots and 23 dicots). We identified 101 orthologous groups (OGs) of genes being conserved among almost all monocot and dicot species analyzed. We observed that more than 10% of these OGs are absent in the Brassicaceae species studied. We show that the ECD structural features are not always conserved among orthologs, suggesting that functions may have diverged in some OG sets. Moreover, we looked at targets of positive selection footprints in 12 pairs of OGs and noticed that depending on the subgroups, positive selection occurred more frequently either in the ECDs or in the KDs.

Immunoprofiling of Rice Root Cortex Reveals Two Cortical Subdomains

The formation and differentiation of aerenchyma, i.e., air-containing cavities that are critical for flooding tolerance, take place exclusively in the cortex. The understanding of development and differentiation of the cortex is thus an important issue; however, studies on this tissue are limited, partly because of the lack of available molecular tools. We screened a commercially available library of cell wall antibodies to identify markers of cortical tissue in rice roots. Out of the 174 antibodies screened, eight were cortex-specific. Our analysis revealed that two types of cortical tissues are present in rice root seedlings. We named these cell layers "inner" and "outer" based on their location relative to the stele. We then used the antibodies to clarify cell identity in lateral roots. Without these markers, previous studies could not distinguish between the cortex and sclerenchyma in small lateral roots. By immunostaining lateral root sections, we showed that the internal ground tissue in small lateral roots has outer cortical identity.

The roots of future rice harvests

Rice production faces the challenge to be enhanced by 50% by year 2030 to meet the growth of the population in rice-eating countries. Whereas yield of cereal crops tend to reach plateaus and a yield is likely to be deeply affected by climate instability and resource scarcity in the coming decades, building rice cultivars harboring root systems that can maintain performance by capturing water and nutrient resources unevenly distributed is a major breeding target. Taking advantage of gathering a community of rice root biologists in a Global Rice Science Partnership workshop held in Montpellier, France, we present here the recent progresses accomplished in this area and focal points where an international network of laboratories should direct their efforts.

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.

Surfing along the root ground tissue gene network

Organization of tissues in Arabidopsis thaliana root is made of, from outside in, epidermis, cortex, middle cortex, endodermis, pericycle and vascular tissues. Cortex, middle cortex and endodermis form the ground tissue (GT) system. Functional and molecular characterization of GT patterning mutants' properties has greatly increased our understanding of fundamental processes of plant root development. These studies have demonstrated GT is an elegant model that can be used to study how different cell types and cell fates are specified. This review analyzes GT mutants to provide a detailed account of the molecular network that regulates GT formation in A. thaliana. The most recent results indicate an unexpectedly complex network of transcription factors, epigenetic and hormonal controls that play crucial roles in GT development. Major differences exist between GT formation in dicots and monocots, particularly in the model plant rice, opening the way for evo-devo of GT formation in angiosperm. In rice, adaptation to submergence relies on a multilayered cortex. Moreover, variation in the number of cortex cell layers is also observed between the five root types. A mechanism of control for cortical cell number should then exist in rice and it remains to be determined if any of the Arabidopsis thaliana identified GT network members are also involved in this process in rice. Alternatively, a totally different network may have been invented. However, first available results suggest functional conservation in rice of at least two transcription factors, SHORT ROOT (SHR) and SCARECROW (SCR), involved in ground tissue formation in Arabidopsis.

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.

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.

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.

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.

Genetic control of root development in rice, the model cereal

Cereals possess a fibrous root system that is mainly composed of crown roots that emerge postembryonically from the nodes of the stem. Because the root system is not directly accessible and consequently difficult to study, it remains a target for breeders to improve the ability of plants to exploit the mineral and water resources of the soil. Breeding for root architecture necessitates identifying the genetic determinants of root development. This research is now underway in cereals, particularly in rice, the monocot model species. In this review, we examine recent data identifying genes that govern root development in cereals, such as ARL1/CRL1 in rice and RTCS in maize which encodes a conserved lateral organ boundary domain transcription factor involved in crown root initiation and development in response to auxin. Finally, we discuss the detection and validation of root development quantitative trait loci.

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.

GreenPhylDB: a database for plant comparative genomics

GreenPhylDB (http://greenphyl.cirad.fr) is a comprehensive platform designed to facilitate comparative functional genomics in Oryza sativa and Arabidopsis thaliana genomes. The main functions of GreenPhylDB are to assign O. sativa and A. thaliana sequences to gene families using a semi-automatic clustering procedure and to create ‘orthologous’ groups using a phylogenomic approach. To date, GreenPhylDB comprises the most complete list of plant gene families, which have been manually curated (6421 families). GreenPhylDB also contains all of the phylogenomic relationships computed for 4375 families. A total of 492 TAIR, 1903 InterPro and 981 KEGG families and subfamilies were manually curated using the clusters created with the TribeMCL software. GreenPhylDB integrates information from several other databases including UniProt, KEGG, InterPro, TAIR and TIGR. Several entry points can be used to display phylogenomic relationships for A. thaliana or O. sativa sequences, using TAIR, TIGR gene ID, family name, InterPro, gene alias, UniProt or protein/nucleic sequence. Finally, a powerful phylogenomics tool, GreenPhyl Ortholog Search Tool (GOST), was incorporated into GreenPhylDB to predict orthologous relationships between O. sativa/A. thaliana protein(s) and sequences from other plant species.

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.

Model-assisted physiological analysis of Phyllo, a rice architectural mutant

Studies of phenotype of knockout mutants can provide new insights into physiological, phenological and architectural feedbacks in the plant system. Phyllo, a mutant of Nippon Bare rice (Oryza sativa L.) producing small leaves in rapid succession, was isolated during multiplication of a T-DNA insertion library. Phyllo phenotype was compared with the wild type (WT) during vegetative development in hydroponics culture using a wide range of physiological and biometric measurements. These were integrated with the help of the functional-structural model EcoMeristem, explicitly designed to study interactions between morphogenesis and carbon assimilation. Although the phenotype of the mutant was caused by a single recessive gene, it differed in many ways from the WT, suggesting a pleiotropic effect of this mutation. Phyllochron was 25 (1-4 leaf stage) to 38% (>>4 leaf stage) shorter but showed normal transition from juvenile to adult phase after leaf 4. Leaf size also increased steadily with leaf position as in WT. The mutant had reduced leaf blade length : width and blade : sheath length ratios, particularly during the transition from heterotrophic to autotrophic growth. During the same period, root : shoot dry weight ratio was significantly diminished. Specific leaf area (SLA) was strongly increased in the mutant but showed normal descending patterns with leaf position. Probably related to high SLA, the mutant had much lower light-saturated leaf photosynthetic rates and lower radiation use efficiency (RUE) than the WT. Leaf extension rates were strongly reduced in absolute terms but were high in relative terms (normalised by final leaf length). The application of the EcoMeristem model to these data indicated that the mutant was severely deficient in assimilate, resulting from low RUE and high organ initiation rate causing high assimilate demand. This was particularly pronounced during the heterotrophic-autotrophic transition, probably causing shorter leaf blades relative to sheaths, as well as a temporary reduction of assimilate partitioning to roots. The model accurately simulated the mutant's high leaf mortality and absence of tillering. The simulated assimilate shortage was supported by observed reductions in starch storage in sheaths. Soluble sugar concentrations differed between mutant and WT in roots but not in shoots. Specifically, the hexose : sucrose ratio was 50% lower in the roots of the mutant, possibly indicating low invertase activity. Furthermore, two OsCIN genes coding for cell wall invertases were not expressed in roots, and others were expressed weakly. This was interpreted as natural silencing via sugar signalling. In summary, the authors attributed the majority of observed allometric and metabolic modifications in the mutant to an extreme assimilate shortage caused by hastened shoot organogenesis and inefficient leaf morphology.

OryGenesDB: a database for rice reverse genetics

Insertional mutant databases containing Flanking Sequence Tags (FSTs) are becoming key resources for plant functional genomics. We have developed OryGenesDB (), a database dedicated to rice reverse genetics. Insertion mutants of rice genes are catalogued by Flanking Sequence Tag (FST) information that can be readily accessed by this database. Our database presently contains 44166 FSTs generated by most of the rice insertional mutagenesis projects. The OryGenesDB genome browser is based on the powerful Generic Genome Browser (GGB) developed in the framework of the Generic Model Organism Project (GMOD). The main interface of our web site displays search and analysis interfaces to look for insertions in any candidate gene of interest. Several starting points can be used to exhaustively retrieve the insertions positions and associated genomic information using blast, keywords or gene name search. The toolbox integrated in our database also includes an ‘anchoring’ option that allows immediate mapping and visualization of up to 50 nucleic acid sequences in the rice Genome Browser of OryGenesDB. As a first step toward plant comparative genomics, we have linked the rice and Arabidopsis whole genome using all the predicted pairs of orthologs by best BLAST mutual hit (BBMH) connectors.

High throughput T-DNA insertion mutagenesis in rice: a first step towards in silico reverse genetics

A library of 29,482 T-DNA enhancer trap lines has been generated in rice cv. Nipponbare. The regions flanking the T-DNA left border from the first 12,707 primary transformants were systematically isolated by adapter anchor PCR and sequenced. A survey of the 7480 genomic sequences larger than 30 bp (average length 250 bp), representing 56.4% of the total readable sequences and matching the rice bacterial artificial chromosome/phage artificial chromosome (BAC/PAC) sequences assembled in pseudomolecules allowed the assigning of 6645 (88.8%) T-DNA insertion sites to at least one position in the rice genome of cv. Nipponbare. T-DNA insertions appear to be rather randomly distributed over the 12 rice chromosomes, with a slightly higher insertion frequency in chromosomes 1, 2, 3 and 6. The distribution of 723 independent T-DNA insertions along the chromosome 1 pseudomolecule did not differ significantly from that of the predicted coding sequences in exhibiting a lower insertion density around the centromere region and a higher density in the subtelomeric regions where the gene density is higher. Further establishment of density graphs of T-DNA inserts along the recently released 12 rice pseudomolecules confirmed this non-uniform chromosome distribution. T-DNA appeared less prone to hot spots and cold spots of integration when compared with those revealed by a concurrent assignment of the Tos17 retrotransposon flanking sequences deposited in the National Center for Biotechnology Information (NCBI). T-DNA inserts rarely integrated into repetitive sequences. Based on the predicted gene annotation of chromosome 1, preferential insertion within the first 250 bp from the putative ATG start codon has been observed. Using 4 kb of sequences surrounding the insertion points, 62% of the sequences showed significant similarity to gene encoding known proteins (E-value < 1.00 e(-05)). To illustrate the in silico reverse genetic approach, identification of 83 T-DNA insertions within genes coding for transcription factors (TF) is presented. Based both on the estimated number of members of several large TF gene families (e.g. Myb, WRKY, HD-ZIP, Zinc-finger) and on the frequency of insertions in chromosome 1 predicted genes, we could extrapolate that 7-10% of the rice gene complement is already tagged by T-DNA insertion in the 6116 independent transformant population. This large resource is of high significance while assisting studies unravelling gene function in rice and cereals, notably through in silico reverse genetics.

Molecular and genetic characterization of a non-climacteric phenotype in melon reveals two loci conferring altered ethylene response in fruit

Fruit ripening and abscission are associated with an ethylene burst in several melon (Cucumis melo) genotypes. In cantaloupe as in other climacteric fruit, exogenous ethylene can prematurely induce abscission, ethylene production, and ripening. Melon genotypes without fruit abscission or without ethylene burst also exist and are, therefore, non-climacteric. In the nonabscising melon fruit PI 161375, exogenous ethylene failed to stimulate abscission, loss of firmness, ethylene production, and expression of all target genes tested. However, the PI 161375 etiolated seedlings displayed the usual ethylene-induced triple response. Genetic analysis on a population of recombinant cantaloupe Charentais x PI 161375 inbred lines in segregation for fruit abscission and ethylene production indicated that both characters are controlled by two independent loci, abscission layer (Al)-3 and Al-4. The non-climacteric phenotype in fruit tissues is attributable to ethylene insensitivity conferred by the recessive allelic forms from PI 161375. Five candidate genes (two ACO, two ACS, and ERS) that were localized on the melon genetic map did not exhibit colocalization with Al-3 or Al-4.

A reference map of Cucumis melo based on two recombinant inbred line populations

A composite genetic melon map was generated based on two recombinant inbred line (RI) populations. By analyzing the segregation of 346 AFLPs, 113 IMAs and phenotypic characters on a RI population of 163 individuals derived from the cross Védrantais x PI 161375, a first map was constructed. About 20% of the molecular markers were skewed, and the residual heterozygosity was estimated at 4.43% which was not significantly different from the theoretical value of 4.2%. The genome distribution of molecular markers among the 12 linkage groups was not different from a random distribution with the exception of linkage group XII which was found significantly less populated. The genome distributions of IMAs and AFLPs were complementary. AFLPs were found mainly in the middle of each linkage group and sometimes clustered, whereas IMAs were found mainly at the end. A total of 318 molecular markers, mainly AFLP and IMA markers, were mapped on 63 RIs of the second population, Védrantais x PI 414723. Comparison of the maps enables one to conclude that AFLPs and IMAs of like molecular size, amplified with the same primer combination, correspond to the same genetic locus. Both maps were joined through 116 common markers comprising 106 comigrating AFLPs/IMAs, plus five SSRs and five phenotypic markers. The integrated melon map contained 668 loci issuing from the segregation of 1,093 molecular markers in the two RI populations. The composite map spanned 1,654 cM on 12 linkage groups which is the haploid number of chromosomes in melon. Thirty two known-function probes, i.e. known-function genes (9) and morphological traits (23), were included in this map. In addition, the composite map was anchored to previously published maps through SSRs, RFLPs and phenotypic characters.

Genetic control of fruit shape acts prior to anthesis in melon ( Cucumis melo L.)

Genetic control of fruit shape in Cucumis melo was studied using QTL analysis in two Recombinant Inbred (RI) populations consisting of 163 and 63 individuals, respectively, obtained by crossing the same round-fruited parent with two different elongated-fruit lines. Fruit shape is mainly explained by fruit length in these two populations. Most QTLs for fruit shape and ovary shape detected were found to co-segregate, thus demonstrating early control of fruit shape during ovary development. A high level of correlation between fruit shape and ovary shape was also found in 14 unrelated genetic lines, a finding which suggests that control of fruit shape by gene(s) active early in the ovary is a general feature in C. melo. Two major flower genes, a ( monoecious) and p ( pentamerous), were shown to have major effects on fruit shape. Major tightly linked QTLs for fruit and ovary shape were found close to the a and p genes, probably reflecting their pleiotropic effect on fruit shape. Moreover, one of the two QTLs detected in the Védrantais x PI 414723 population was also found in the Védrantais x PI 161375 population. Variation of fruit shape in melon could be due to variations having quantitative effects on a large set of genes that are probably involved in ovary development.