Rice Yellow Mottle Virus resistance by genome editing of the Oryza sativa L. ssp. japonica nucleoporin gene OsCPR5.1 but not OsCPR5.2.

Rice yellow mottle virus (RYMV) causes one of the most devastating rice diseases in Africa. Management of RYMV is challenging. Genetic resistance provides the most effective and environment-friendly control. The recessive resistance locus rymv2 (OsCPR5.1) had been identified in African rice (Oryza glaberrima), however, introgression into Oryza sativa ssp. japonica and indica remains challenging due to crossing barriers. Here, we evaluated whether CRISPR/Cas9 genome editing of the two rice nucleoporin paralogs OsCPR5.1 (RYMV2) and OsCPR5.2 can be used to introduce RYMV resistance into the japonica variety Kitaake. Both paralogs had been shown to complement the defects of the Arabidopsis atcpr5 mutant, indicating partial redundancy. Despite striking sequence and structural similarities between the two paralogs, only oscpr5.1 loss-of-function mutants were fully resistant, while loss-of-function oscpr5.2 mutants remained susceptible, intimating that OsCPR5.1 plays a specific role in RYMV susceptibility. Notably, edited lines with short in-frame deletions or replacements in the N-terminal domain (predicted to be unstructured) of OsCPR5.1 were hypersusceptible to RYMV. In contrast to mutations in the single Arabidopsis AtCPR5 gene, which caused severely dwarfed plants, oscpr5.1 and oscpr5.2 single and double knockout mutants showed neither substantial growth defects nor symptoms indicative lesion mimic phenotypes, possibly reflecting functional differentiation. The specific editing of OsCPR5.1, while maintaining OsCPR5.2 activity, provides a promising strategy for generating RYMV-resistance in elite Oryza sativa lines as well as for effective stacking with other RYMV resistance genes or other traits.

Genome editing of an African elite rice variety confers resistance against endemic and emerging Xanthomonas oryzae pv. oryzae strains.

Bacterial leaf blight (BB) of rice, caused by Xanthomonas oryzae pv. oryzae (Xoo), threatens global food security and the livelihood of small-scale rice producers. Analyses of Xoo collections from Asia, Africa and the Americas demonstrated complete continental segregation, despite robust global rice trade. Here, we report unprecedented BB outbreaks in Tanzania. The causative strains, unlike endemic African Xoo, carry Asian-type TAL effectors targeting the sucrose transporter SWEET11a and iTALes suppressing Xa1. Phylogenomics clustered these strains with Xoo from Southern-China. African rice varieties do not carry effective resistance. To protect African rice production against this emerging threat, we developed a hybrid CRISPR-Cas9/Cpf1 system to edit all known TALe-binding elements in three SWEET promoters of the East African elite variety Komboka. The edited lines show broad-spectrum resistance against Asian and African strains of Xoo, including strains recently discovered in Tanzania. The strategy could help to protect global rice crops from BB pandemics.

The Complete Genome Resource of Xanthomonas oryzae pv. oryzae CIX2779 Includes the First Sequence of a Plasmid for an African Representative of This Rice Pathogen.

The bacterial plant pathogen Xanthomonas oryzae pv. oryzae is responsible for the foliar rice bacterial blight disease. Genetically contrasted, continent-specific, sublineages of this species can cause important damages to rice production both in Asia and Africa. We report on the genome of the CIX2779 strain of this pathogen, previously named NAI1 and originating from Niger. Oxford Nanopore long reads assembly and Illumina short reads polishing produced a genome sequence composed of a 4,725,792-bp circular chromosome and a 39,798-bp-long circular plasmid designated pCIX2779_1. The chromosome structure and base-level sequence are highly related to reference strains of African X. oryzae pv. oryzae and encode identical transcription activator-like effectors for virulence. Importantly, our in silico analysis strongly indicates that pCIX2779_1 is a genuine conjugative plasmid, the first indigenous one sequenced from an African strain of the X. oryzae species. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

The Rice ILI2 Locus Is a Bidirectional Target of the African Xanthomonas oryzae pv. oryzae Major Transcription Activator-like Effector TalC but Does Not Contribute to Disease Susceptibility.

Xanthomonas oryzae pv. oryzae (Xoo) strains that cause bacterial leaf blight (BLB) limit rice (Oryza sativa) production and require breeding more resistant varieties. Transcription activator-like effectors (TALEs) activate transcription to promote leaf colonization by binding to specific plant host DNA sequences termed effector binding elements (EBEs). Xoo major TALEs universally target susceptibility genes of the SWEET transporter family. TALE-unresponsive alleles of clade III OsSWEET susceptibility gene promoter created with genome editing confer broad resistance on Asian Xoo strains. African Xoo strains rely primarily on the major TALE TalC, which targets OsSWEET14. Although the virulence of a talC mutant strain is severely impaired, abrogating OsSWEET14 induction with genome editing does not confer equivalent resistance on African Xoo. To address this contradiction, we postulated the existence of a TalC target susceptibility gene redundant with OsSWEET14. Bioinformatics analysis identified a rice locus named ATAC composed of the INCREASED LEAF INCLINATION 2 (ILI2) gene and a putative lncRNA that are shown to be bidirectionally upregulated in a TalC-dependent fashion. Gain-of-function approaches with designer TALEs inducing ATAC sequences did not complement the virulence of a Xoo strain defective for SWEET gene activation. While editing the TalC EBE at the ATAC loci compromised TalC-mediated induction, multiplex edited lines with mutations at the OsSWEET14 and ATAC loci remained essentially susceptible to African Xoo strains. Overall, this work indicates that ATAC is a probable TalC off-target locus but nonetheless documents the first example of divergent transcription activation by a native TALE during infection.

An atypical class of non-coding small RNAs is produced in rice leaves upon bacterial infection.

Non-coding small RNAs (sRNA) act as mediators of gene silencing and regulate plant growth, development and stress responses. Early insights into plant sRNAs established a role in antiviral defense and they are now extensively studied across plant-microbe interactions. Here, sRNA sequencing discovered a class of sRNA in rice (Oryza sativa) specifically associated with foliar diseases caused by Xanthomonas oryzae bacteria. Xanthomonas-induced small RNAs (xisRNAs) loci were distinctively upregulated in response to diverse virulent strains at an early stage of infection producing a single duplex of 20-22 nt sRNAs. xisRNAs production was dependent on the Type III secretion system, a major bacterial virulence factor for host colonization. xisRNA loci overlap with annotated transcripts sequences, with about half of them encoding protein kinase domain proteins. A number of the corresponding rice cis-genes have documented functions in immune signaling and xisRNA loci predominantly coincide with the coding sequence of a conserved kinase motif. xisRNAs exhibit features of small interfering RNAs and their biosynthesis depend on canonical components OsDCL1 and OsHEN1. xisRNA induction possibly mediates post-transcriptional gene silencing but they do not broadly suppress cis-genes expression on the basis of mRNA-seq data. Overall, our results identify a group of unusual sRNAs with a potential role in plant-microbe interactions.

Improved bacterial leaf blight disease resistance in the major elite Vietnamese rice cultivar TBR225 via editing of the OsSWEET14 promoter.

TBR225 is one of the most popular commercial rice varieties in Northern Vietnam. However, this variety is highly susceptible to bacterial leaf blight (BLB), a disease caused by Xanthomonas oryzae pv. oryzae (Xoo) which can lead to important yield losses. OsSWEET14 belongs to the SWEET gene family that encodes sugar transporters. Together with other Clade III members, it behaves as a susceptibility (S) gene whose induction by Asian Xoo Transcription-Activator-Like Effectors (TALEs) is absolutely necessary for disease. In this study, we sought to introduce BLB resistance in the TBR225 elite variety. First, two Vietnamese Xoo strains were shown to up-regulate OsSWEET14 upon TBR225 infection. To investigate if this induction is connected with disease susceptibility, nine TBR225 mutant lines with mutations in the AvrXa7, PthXo3 or TalF TALEs DNA target sequences of the OsSWEET14 promoter were obtained using the CRISPR/Cas9 editing system. Genotyping analysis of T0 and T1 individuals showed that mutations were stably inherited. None of the examined agronomic traits of three transgene-free T2 edited lines were significantly different from those of wild-type TBR225. Importantly, one of these T2 lines, harboring the largest homozygous 6-bp deletion, displayed decreased OsSWEET14 expression as well as a significantly reduced susceptibility to a Vietnamese Xoo strains and complete resistance to another one. Our findings indicate that CRISPR/Cas9 editing conferred an improved BLB resistance to a Vietnamese commercial elite rice variety.

Broad-spectrum resistance to bacterial blight in rice using genome editing.

Bacterial blight of rice is an important disease in Asia and Africa. The pathogen, Xanthomonas oryzae pv. oryzae (Xoo), secretes one or more of six known transcription-activator-like effectors (TALes) that bind specific promoter sequences and induce, at minimum, one of the three host sucrose transporter genes SWEET11, SWEET13 and SWEET14, the expression of which is required for disease susceptibility. We used CRISPR-Cas9-mediated genome editing to introduce mutations in all three SWEET gene promoters. Editing was further informed by sequence analyses of TALe genes in 63 Xoo strains, which revealed multiple TALe variants for SWEET13 alleles. Mutations were also created in SWEET14, which is also targeted by two TALes from an African Xoo lineage. A total of five promoter mutations were simultaneously introduced into the rice line Kitaake and the elite mega varieties IR64 and Ciherang-Sub1. Paddy trials showed that genome-edited SWEET promoters endow rice lines with robust, broad-spectrum resistance.

Functional and Genome Sequence-Driven Characterization of tal Effector Gene Repertoires Reveals Novel Variants With Altered Specificities in Closely Related Malian Xanthomonas oryzae pv. oryzae Strains.

Rice bacterial leaf blight (BLB) is caused by Xanthomonas oryzae pv. oryzae (Xoo) which injects Transcription Activator-Like Effectors (TALEs) into the host cell to modulate the expression of target disease susceptibility genes. Xoo major-virulence TALEs universally target susceptibility genes of the SWEET sugar transporter family. TALE-unresponsive alleles of OsSWEET genes have been identified in the rice germplasm or created by genome editing and confer resistance to BLB. In recent years, BLB has become one of the major biotic constraints to rice cultivation in Mali. To inform the deployment of alternative sources of resistance in this country, rice lines carrying alleles of OsSWEET14 unresponsive to either TalF (formerly Tal5) or TalC, two important TALEs previously identified in West African Xoo, were challenged with a panel of strains recently isolated in Mali and were found to remain susceptible to these isolates. The characterization of TALE repertoires revealed that talF and talC specific molecular markers were simultaneously present in all surveyed Malian strains, suggesting that the corresponding TALEs are broadly deployed by Malian Xoo to redundantly target the OsSWEET14 gene promoter. Consistent with this, the capacity of most Malian Xoo to induce OsSWEET14 was unaffected by either talC- or talF-unresponsive alleles of this gene. Long-read sequencing and assembly of eight Malian Xoo genomes confirmed the widespread occurrence of active TalF and TalC variants and provided a detailed insight into the diversity of TALE repertoires. All sequenced strains shared nine evolutionary related tal effector genes. Notably, a new TalF variant that is unable to induce OsSWEET14 was identified. Furthermore, two distinct TalB variants were shown to have lost the ability to simultaneously induce two susceptibility genes as previously reported for the founding members of this group from strains MAI1 and BAI3. Yet, both new TalB variants retained the ability to induce one or the other of the two susceptibility genes. These results reveal molecular and functional differences in tal repertoires and will be important for the sustainable deployment of broad-spectrum and durable resistance to BLB in West Africa.

Targeted promoter editing for rice resistance to Xanthomonas oryzae pv. oryzae reveals differential activities for SWEET14-inducing TAL effectors.

As a key virulence strategy to cause bacterial leaf blight, Xanthomonas oryzae pv. oryzae (Xoo) injects into the plant cell DNA-binding proteins called transcription activator-like effectors (TALEs) that bind to effector-binding elements (EBEs) in a sequence-specific manner, resulting in host gene induction. TALEs AvrXa7, PthXo3, TalC and Tal5, found in geographically distant Xoo strains, all target OsSWEET14, thus considered as a pivotal TALE target acting as major susceptibility factor during rice-Xoo interactions. Here, we report the generation of an allele library of the OsSWEET14 promoter through stable expression of TALE-nuclease (TALEN) constructs in rice. The susceptibility level of lines carrying mutations in AvrXa7, Tal5 or TalC EBEs was assessed. Plants edited in AvrXa7 or Tal5 EBEs were resistant to bacterial strains relying on the corresponding TALE. Surprisingly, although indels within TalC EBE prevented OsSWEET14 induction in response to BAI3 wild-type bacteria relying on TalC, loss of TalC responsiveness failed to confer resistance to this strain. The TalC EBE mutant line was, however, resistant to a strain expressing an artificial SWEET14-inducing TALE whose EBE was also edited in this line. This work offers the first set of alleles edited in TalC EBE and uncovers a distinct, broader range of activities for TalC compared to AvrXa7 or Tal5. We propose the existence of additional targets for TalC beyond SWEET14, suggesting that TALE-mediated plant susceptibility may result from induction of several, genetically redundant, host susceptibility genes by a single effector.

The Casuarina NIN gene is transcriptionally activated throughout Frankia root infection as well as in response to bacterial diffusible signals.

Root nodule symbioses (RNS) allow plants to acquire atmospheric nitrogen by establishing an intimate relationship with either rhizobia, the symbionts of legumes or Frankia in the case of actinorhizal plants. In legumes, NIN (Nodule INception) genes encode key transcription factors involved in nodulation. Here we report the characterization of CgNIN, a NIN gene from the actinorhizal tree Casuarina glauca using both phylogenetic analysis and transgenic plants expressing either ProCgNIN::reporter gene fusions or CgNIN RNAi constructs. We have found that CgNIN belongs to the same phylogenetic group as other symbiotic NIN genes and CgNIN is able to complement a legume nin mutant for the early steps of nodule development. CgNIN expression is correlated with infection by Frankia, including preinfection stages in developing root hairs, and is induced by culture supernatants. Knockdown mutants were impaired for nodulation and early root hair deformation responses were severely affected. However, no mycorrhizal phenotype was observed and no induction of CgNIN expression was detected in mycorrhizas. Our results indicate that elements specifically required for nodulation include NIN and possibly related gene networks derived from the nitrate signalling pathways.

Transcriptome Changes in Hirschfeldia incana in Response to Lead Exposure.

Hirschfeldia incana, a pseudometallophyte belonging to the Brassicaceae family and widespread in the Mediterranean region, was selected for its ability to grow on soils contaminated by lead (Pb). The global comparison of gene expression using microarrays between a plant susceptible to Pb (Arabidopsis thaliana) and a Pb tolerant plant (H. incana) enabled the identification of a set of specific genes expressed in response to lead exposure. Three groups of genes were particularly over-represented by the Pb exposure in the biological processes categorized as photosynthesis, cell wall, and metal handling. Each of these gene groups was shown to be directly involved in tolerance or in protection mechanisms to the phytotoxicity associated with Pb. Among these genes, we demonstrated that MT2b, a metallothionein gene, was involved in lead accumulation, confirming the important role of metallothioneins in the accumulation and the distribution of Pb in leaves. On the other hand, several genes involved in biosynthesis of ABA were shown to be up-regulated in the roots and shoots of H. incana treated with Pb, suggesting that ABA-mediated signaling is a possible mechanism in response to Pb treatment in H. incana. This latest finding is an important research direction for future studies.

Lead tolerance and accumulation in Hirschfeldia incana, a Mediterranean Brassicaceae from metalliferous mine spoils.

Lead is a heavy metal of particular concern with respect to environmental quality and health. The lack of plant species that accumulate and tolerate Pb is a limiting factor to understand the molecular mechanisms involved in Pb tolerance. In this study we identified Hirschfeldia incana, a Brassicaceae collected from metalliferous mine spoils in Morocco, as a Pb accumulator plant. H. incana exhibited high Pb accumulation in mine soils and in hydroponic cultures. Major Pb accumulation occurred in the roots and a part of Pb translocated from the roots to the shoots, even to the siliques. These findings demonstrated that H. incana is a Pb accumulator species. The expression of several candidate genes after Pb-exposure was measured by quantitative PCR and two of them, HiHMA4 and HiMT2a, coding respectively for a P1B-type ATPase and a metallothionein, were particularly induced by Pb-exposure in both roots and leaves. The functional characterization of HiHMA4 and HiMT2a was achieved using Arabidopsis T-DNA insertional mutants. Pb content and primary root growth analysis confirmed the role of these two genes in Pb tolerance and accumulation. H. incana could be considered as a good experimental model to identify genes involved in lead tolerance and accumulation in plants.

Transcriptomics of actinorhizal symbioses reveals homologs of the whole common symbiotic signaling cascade.

Comparative transcriptomics of two actinorhizal symbiotic plants, Casuarina glauca and Alnus glutinosa, was used to gain insight into their symbiotic programs triggered following contact with the nitrogen-fixing actinobacterium Frankia. Approximately 14,000 unigenes were recovered in roots and 3-week-old nodules of each of the two species. A transcriptomic array was designed to monitor changes in expression levels between roots and nodules, enabling the identification of up- and down-regulated genes as well as root- and nodule-specific genes. The expression levels of several genes emblematic of symbiosis were confirmed by quantitative polymerase chain reaction. As expected, several genes related to carbon and nitrogen exchange, defense against pathogens, or stress resistance were strongly regulated. Furthermore, homolog genes of the common and nodule-specific signaling pathways known in legumes were identified in the two actinorhizal symbiotic plants. The conservation of the host plant signaling pathway is all the more surprising in light of the lack of canonical nod genes in the genomes of its bacterial symbiont, Frankia. The evolutionary pattern emerging from these studies reinforces the hypothesis of a common genetic ancestor of the Fabid (Eurosid I) nodulating clade with a genetic predisposition for nodulation.

Activation of the isoflavonoid pathway in actinorhizal symbioses.

We investigated the involvement of flavonoids in the actinorhizal nodulation process resulting from the interaction between the tropical tree Casuarina glauca Sieb. ex Spreng. and the actinomycete Frankia. Eight C. glauca genes involved in flavonoid biosynthesis: chalcone synthase (CHS), chalcone isomerase (CHI), isoflavone reductase (IFR), flavonoid-3-hydroxylase (F3H), flavonoid 3'-hydroxylase (F3'H), flavonoid 3',5' hydroxylase (F3'5'H), dihydroflavonol 4-reductase (DFR) and flavonol synthase (FLS), were identified from a unigene database and gene expression patterns were monitored by quantitative real-time PCR (qRT-PCR) during the nodulation time course. Results showed that FLS and F3'5'H transcripts accumulated in mature nodules whereas CHI and IFR transcripts accumulated preferentially early after inoculation with Frankia. Comparison of IFR and CHI expression in inoculated plants and in control plants cultivated with or without nitrogen confirmed that early expression of IFR is specifically linked to symbiosis. Taken together, these data suggest for the first time that isoflavonoids are implicated in actinorhizal nodulation.

Post-transcriptional gene silencing in the root system of the actinorhizal tree Allocasuarina verticillata.

In recent years, RNA interference has been exploited as a tool for investigating gene function in plants. We tested the potential of double-stranded RNA interference technology for silencing a transgene in the actinorhizal tree Allocasuarina verticillata. The approach was undertaken using stably transformed shoots expressing the beta-glucuronidase (GUS) gene under the control of the constitutive promoter 35S; the shoots were further transformed with the Agrobacterium rhizogenes A4RS containing hairpin RNA (hpRNA) directed toward the GUS gene, and driven by the 35S promoter. The silencing and control vectors contained the reporter gene of the green fluorescent protein (GFP), thus allowing a screening of GUS-silenced composite plantlets for autofluorescence. With this rapid procedure, histochemical data established that the reporter gene was strongly silenced in both fluorescent roots and actinorhizal nodules. Fluorometric data further established that the level of GUS silencing was usually greater than 90% in the hairy roots containing the hairpin GUS sequences. We found that the silencing process of the reporter gene did not spread to the aerial part of the composite A. verticillata plants. Real-time quantitative polymerase chain reaction showed that GUS mRNAs were substantially reduced in roots and, thereby, confirmed the knock-down of the GUS transgene in the GFP(+) hairy roots. The approach described here will provide a versatile tool for the rapid assessment of symbiotically related host genes in actinorhizal plants of the Casuarinaceae family.

SymRK defines a common genetic basis for plant root endosymbioses with arbuscular mycorrhiza fungi, rhizobia, and Frankiabacteria.

Root endosymbioses vitally contribute to plant nutrition and fitness worldwide. Nitrogen-fixing root nodulation, confined to four plant orders, encompasses two distinct types of associations, the interaction of legumes (Fabales) with rhizobia bacteria and actinorhizal symbioses, where the bacterial symbionts are actinomycetes of the genus Frankia. Although several genetic components of the host-symbiont interaction have been identified in legumes, the genetic basis of actinorhiza formation is unknown. Here, we show that the receptor-like kinase gene SymRK, which is required for nodulation in legumes, is also necessary for actinorhiza formation in the tree Casuarina glauca. This indicates that both types of nodulation symbiosis share genetic components. Like several other legume genes involved in the interaction with rhizobia, SymRK is also required for the interaction with arbuscular mycorrhiza (AM) fungi. We show that SymRK is involved in AM formation in C. glauca as well and can restore both nodulation and AM symbioses in a Lotus japonicus symrk mutant. Taken together, our results demonstrate that SymRK functions as a vital component of the genetic basis for both plant-fungal and plant-bacterial endosymbioses and is conserved between legumes and actinorhiza-forming Fagales.

A Role for auxin during actinorhizal symbioses formation?

The symbiotic interaction between the soil bacteria Frankia and actinorhizal plants leads to the formation of nitrogen-fixing nodules resembling modified lateral roots. Little is known about the signals exchanged between the two partners during the establishment of these endosymbioses. However, a role for plant hormones has been suggested.Recently, we studied the role of auxin influx activity during actinorhizal symbioses. An inhibitor of auxin influx was shown to perturb nodule formation. Moreover we identified a functional auxin influx carrier that is produced specifically in Frankia-infected cells. These results together with previous data showing auxin production by Frankia lead us to propose a model of auxin action during the symbiotic infection process.

Functional analysis of the metallothionein gene cgMT1 isolated from the actinorhizal tree Casuarina glauca.

cgMT1 is a metallothionein (MT)-like gene that was isolated from a cDNA library of young nitrogen-fixing nodules resulting from the symbiotic interaction between Frankia spp. and the actinorhizal tree Casuarina glauca. cgMT1 is highly transcribed in the lateral roots and nitrogen-fixing cells of actinorhizal nodules; it encodes a class I type 1 MT. To obtain insight into the function of cgMT1, we studied factors regulating the expression of the MT promoter region (PcgMT1) using a beta-glucuronidase (gus) fusion approach in transgenic plants of Arabidopsis thaliana. We found that copper, zinc, and cadmium ions had no significant effect on the regulation of PcgMT1-gus expression whereas wounding and H2O2 treatments led to an increase in reporter gene activity in transgenic leaves. Strong PcgMT1-gus expression also was observed when transgenic plants were inoculated with a virulent strain of the bacterial pathogen Xanthomonas campestris pv. campestris. Transgenic Arabidopsis plants expressing cgMT1 under the control of the constitutive 35S promoter were characterized by reduced accumulation of H2O2 when leaves were wounded and by increased susceptibility to the bacterial pathogen X. campestris. These results suggest that cgMT1 could play a role during the oxidative response linked to biotic and abiotic stresses.

Auxin influx activity is associated with Frankia infection during actinorhizal nodule formation in Casuarina glauca.

Plants from the Casuarinaceae family enter symbiosis with the actinomycete Frankia leading to the formation of nitrogen-fixing root nodules. We observed that application of the auxin influx inhibitor 1-naphtoxyacetic acid perturbs actinorhizal nodule formation. This suggests a potential role for auxin influx carriers in the infection process. We therefore isolated and characterized homologs of the auxin influx carrier (AUX1-LAX) genes in Casuarina glauca. Two members of this family were found to share high levels of deduced protein sequence identity with Arabidopsis (Arabidopsis thaliana) AUX-LAX proteins. Complementation of the Arabidopsis aux1 mutant revealed that one of them is functionally equivalent to AUX1 and was named CgAUX1. The spatial and temporal expression pattern of CgAUX1 promoter:beta-glucuronidase reporter was analyzed in Casuarinaceae. We observed that CgAUX1 was expressed in plant cells infected by Frankia throughout the course of actinorhizal nodule formation. Our data suggest that auxin plays an important role during plant cell infection in actinorhizal symbioses.

Expressed sequence-tag analysis in Casuarina glauca actinorhizal nodule and root.

The present study aimed to identify and assess the frequency and tissue specificity of plant genes in the actinorhizal Casuarina glauca-Frankia symbiosis through expressed sequence tag (EST) analysis. Using a custom analysis pipeline for raw sequences of C. glauca uninfected roots and nodules, we obtained an EST databank web interface. Gene expression was studied in nodules vs roots using comparative quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). From roots and nodules, 2028 ESTs were created and clustered in 242 contigs and 1429 singletons, giving a total of 1616 unique genes. Half the nodule transcripts showed no similarity to previously identified genes. Genes of primary metabolism, protein synthesis, cell division and defence were highly represented in the nodule library. Differential expression was observed between roots and nodules for several genes linked to primary metabolism and flavonoid biosynthesis. This comparative EST-based study provides the first picture of the set of genes expressed during actinorhizal symbiosis. We consider our database to be a flexible tool that can be used for the management of EST data from other actinorhizal symbioses.