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

Abstract

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.

1. Introduction

Bacterial leaf blight (BLB) is a widespread rice disease caused by Xanthomonas oryzae pv. oryzae (Xoo) [1]. The significant economic impact of this disease on rice cultivation [2] has prompted breeding commercial varieties for enhanced resistance traits by conventional introgression of resistance loci from the rice germplasm [3,4]. Recently, pioneering work [5] and the development of straightforward genome editing tools have, however, sparked a wave of reports on the engineering of BLB resistance in elite varieties [6,7,8,9,10,11]. In this strategy, the modified rice genomic sequences correspond to effector binding elements (EBEs), the DNA sequence recognized by a translocated bacterial virulence protein of the transcription activator-like effectors (TALEs) family. Editing a relevant EBE thus renders the altered genotype refractory to TALE-mediated BLB susceptibility gene induction and confers disease resistance.

Upon injection into the plant host cell via the molecular syringe of the bacterial type III secretion system (T3SS), TALEs enter the nucleus and promote the transcription of their target locus. TALEs were initially thought to induce transcription in a single direction by activating the synthesis of a protein-coding mRNA molecule on the same strand and 40–60 bp downstream of their target EBE [12,13,14]. However, this does not seem to be universal because the transcription of sequences on the opposite strand and upstream of the EBE was recently reported to occur in a couple of specific TALE–EBE combinations [15,16]. EBE sequence recognition is mediated by the central domain of the effector that wraps around the major groove of the genomic DNA and makes specific contacts with bases of this sequence. These interactions involve repeat variable diresidues (RVDs) corresponding to positions 12 and 13 of otherwise-conserved ~34 amino acid repeats that make up the central domain. In a colinear manner, amino acid 13 from successive repeats along the TALE central domain sequence recognizes a cognate nucleobase along the EBE strand. Subsequent to its discovery [17,18], the RVD–nucleotide association code has been relatively well characterized and prominent RVDs, such as NI, NG, HD, and NN, were shown to recognize A, T, C, and G/A, respectively [19]. As a general rule and with the exception of RipTALs [20], EBEs start with a T that is not recognized by a canonical repeat but rather by a structurally related helical bundle or a cryptic repeat located before the central domain [21,22,23]. Because this feature contributes to DNA recognition, it is ordinarily referred to as “repeat 0” and the corresponding location of the recognition site as “position 0.”

With a few notable exceptions, including the transcription elongation factor OsTFIIAγ1 and the b-ZIP transcription factor OsTFX1, which were shown to be induced by TALE effectors from Philippine Xoo strains [24] or the African-Xoo-specific target OsERF#123 [25], in rice, susceptibility genes exploited by Xoo major TALEs, which were thus edited for resistance engineering, belong to clade III of the OsSWEET gene family [26]. How their rice membranous sugar transporter products favor parasitism is still an unsettled issue, but they may support bacterial metabolism by releasing sucrose into the apoplast [27]. Virtually all pathogenic Xoo strains harbor at least one TALE inducing one of three OsSWEET clade III homologs. For instance, PthXo1 from the Philippine strain PXO99A [28] targets OsSWEET11, while PthXo2 variants from Asian Xoo target OsSWEET13 [8,29]. Probably as a result of convergent evolution, overlapping or distinct EBEs have been identified on the OsSWEET14 promoter. The cognate TALEs AvrXa7 and PthXo3 originate from Asian Xoo strains [30]. Others, namely TalC and TalF, have been described in African Xoo strains [26,31].

Both functional characterization and long read sequencing of the African Xoo TALE repertoires revealed that TalC is highly conserved and found in all analyzed strains. A large share of African strains also possesses the redundant TalF TALE [8]. Some strains, such as BAI3 [25] and MAI68 [32], contain a TalF variant that is inactive on the OsSWEET14 promoter and, therefore, rely solely on TalC for OsSWEET14 induction.

The genetic divide observed in phylogenetic studies between Xoo strains of the Asian and African lineages [33,34] also appears in the prospects of engineering BLB resistance with modified, TALE-unresponsive, OsSWEET promoter alleles. While this strategy possesses an immense potential for the broad control of Asian Xoo strains, as previously demonstrated [8,9], the situation appears less straightforward regarding African Xoo bacteria. This is because the dependency of African Xoo strains on the SWEET susceptibility component for host colonization differs from that of the Asian strains analyzed so far. Initial clade III OsSWEET promoter editing approaches to engineer resistance to African strains were unsuccessful [35]: While OsSWEET14 induction upon African Xoo BAI3 inoculations was abrogated in TalC-EBE-edited rice lines, those lines, in two different Japonica rice backgrounds, remained essentially susceptible to this Burkinabe strain. A subsequent study obtained equivalent results when challenging Kitaake lines with the Malian MAI68 strain [32]. More recent findings align with this conclusion and provide additional insight. For instance, a single OsSWEET14 knockout event in Kitaake background exhibits moderately reduced susceptibility (about half the lesion length of wild-type plants) to AXO1947, BAI3, or MAI1 African strains [7]. In contrast, an OsSWEET14 knockout mutant in rice cultivar Zhonghua 11 appears to completely abolish lesions caused by AXO1947 [11], suggesting that other polymorphic loci may epistatically influence the disease outcome depending on the rice genotype. The susceptibility of quintuple-EBE-edited lines from Oliva et al. [8] is again modestly decreased when challenged with the AXO1947 strain. However, a double OsSWEET14/OsSWEET13 knockout line in Kitaake is completely resistant to the AXO1947 African strain [7], which represents evidence that those clade III genes are redundant. The expression of other OsSWEET genes from clade III is not modified upon BAI3 inoculation in either wild-type or edited lines, suggesting that the loss of OsSWEET14 induction is not compensated by transcriptional upregulation of redundant homologous genes. Thus, in contrast to Asian Xoo, African Xoo has the capacity to infect rice in the absence of clade III SWEET susceptibility gene induction [35]. These experimental data were unexpected because they are inconsistent with the dramatic phenotype of BAI3-1-1, a BAI3 derivative with a defective talC gene that has been previously shown to be unable to cause significant BLB symptoms and efficiently colonize susceptible rice leaves, leading to the conclusion that this effector is a major virulence TALE [26,31]. To reconcile these observations, we previously postulated the existence of an unidentified susceptibility gene targeted by TalC and genetically redundant with OsSWEET14 [35].

This study was primarily conducted to explore the possible existence of such a redundant TalC target. Revisiting available rice transcriptomics data with an original strategy helped identify a particularly interesting candidate TalC target locus. In addition, we provide evidence that talC triggers transcription of two divergent transcripts through the recognition of an EBE that is shown to be necessary for induction at this target locus. However, both gain- and loss-of-function approaches could not establish that this rice locus contributes to BLB susceptibility. Overall, this work identifies a second TalC target and provides the first unambiguous results indicating that a native TALE can trigger bidirectional transcription during infection.

2. Results

2.1. Reassessing TalC Target Predictions in the Context of a BAI3-Infected Rice Experimental Transcriptome Yields New Candidate TalC Target Loci

Prior searches for TalC targets in rice identified OsSWEET14 [31]. A subsequent Nipponbare-genome-wide inquiry of African TALE targets combined both in silico target predictions on upstream sequences of annotated genes and lists of differentially expressed genes upon African Xoo infection based on transcriptomic data [25]. In addition to OsSWEET14, the latter approach identified LOC_Os03g58790, a putative ATPase, as a candidate TalC target. However, the induction ratio and prediction rank values are moderately supportive for a TALE target (~4× and 264, respectively). Furthermore, the predicted TalC EBE for this locus starts with a T for position 0 of the RVD array, which contrasts with the functional TalC EBE on the OsSWEET14 promoter, which uncommonly starts with a C, as previously pointed out [31]. Considering both that TALEs may trigger the synthesis of transcripts overlooked in reference annotations and that previous TalC target searches did not incorporate the possibility of “antisense EBEs,” in order to maximize the likelihood of identifying candidate redundant susceptibility genes, we decided to reassess TalC target predictions.

For that, we analyzed a recently documented mRNA-seq dataset [36]. First, an experimental transcriptome (see our dataverse) was inferred with the StringTie software by combining mRNA-seq read alignments to the Nipponbare genome and the reference MSU7 annotation [37]. This experimental transcriptome served to compute read counts at the gene level, which permitted the assessment of differential expression (annotation and expression data are available on our dataverse). From this, we listed 261 loci deemed as differentially induced in comparisons between African-Xoo-strain-BAI3-infected leaves and a BAI3 T3SS defective mutant strain (BAI3H) or mock inoculated Nipponbare leaves. To predict potential TalC EBEs, we ran Talvez [38] using the sequences located 500 bp upstream and 200 bp downstream of the start of the experimental BAI3 T3SS upregulated loci. The resulting list of candidate loci with the 10 highest Talvez prediction scores is shown in Table 1. When available, TalC-dependent induction data (derived from 31) were included in this table. Annotated genome browser views of those candidate TalC target loci with read coverage profiles and predicted EBEs are provided in Figure S1.

As expected, the OsSWEET14 locus (LOC_Os11g31190, also previously designated as Os11N3) carries the top-scoring predicted EBE and exhibits elevated T3SS- and talC-dependent induction ratios (Table 1). Apart from this positive control locus, the INCREASED LEAF INCLINATION 2 (ILI2) locus (LOC_Os11g39000) is the only MSU7 reference annotation locus with a similar profile, being notably induced in a talC-dependent manner (adjusted p-value = 6.29 × 10⁻³; Table 1). In contrast, however, the predicted TalC EBE is positioned in an antisense orientation relative to the ILI2 coding strand.

Interestingly, our analytical approach also uncovered two novel loci, PIX.7679 and PIX.6359, without MSU7 annotation, which appeared to generate mRNA-seq reads only if the bacteria were able to deliver type III effectors (Figure S1 and Figure 1a,b). They were, therefore, noteworthy TalC target candidates even if they could not be assessed immediately with regard to talC-dependent induction because their expression is not measured in microarray data.

2.2. ATAC, A Candidate Alternative TalC-Target Locus Is Induced Bidirectionally and Comprises a Long ncRNA and the ILI2 Gene

Upon further scrutiny of the genomic positions of the predicted EBEs in Table 1, we unexpectedly realized that the ILI2 and PIX.6359 loci actually share the same candidate EBE and are arranged in opposite, head-to-head, orientations on Chromosome 11, with the TalC EBE positioned on the same strand as PIX.6359 (Figure 1a). The interval between the PIX.6359 or ILI2 start site and the proximal end of this EBE is, respectively, 52 bp or 115 bp. The distance from the PIX.6359 transcript start is consistent with distances usually observed between a TALE EBE and its target transcription start (40–60 bp). This configuration is reminiscent of the bidirectional transcription activation by TALEs previously observed in transient plant expression assays [16]. Because of this rather unique feature and because the observed expression profile of the ILI2 and PIX.6359 transcript was at this stage consistent with their being bidirectional targets of TalC (Figure 1b), we tentatively refer to this locus as ATAC (Alternative TalC Target locus), with PIX.6359 designated as ATAC1 and ILI2 as ATAC2.

In order to examine if the sequence of the predicted EBE in this region (EBE_TalC@ATAC) is likely to be targeted by TalC, we aligned this sequence with EBE_TalC on the OsSWEET14 promoter (EBE_TalC@OsSWEET14) and with an “optimal” EBE_TalC sequence based on the Talvez nucleotide–RVD association matrix (Figure 1c). The nine first 5′ nucleotides of EBE_TalC@ATAC and EBE_TalC@OsSWEET14 are strictly conserved. This includes the cytosine at position 0, an atypical feature of EBE_TalC@OsSWEET14. EBE_TalC@SWEET14 displays only one position conflicting with the RVD–nucleotide code. EBE_TalC@ATAC contains two such positions but possibly with compensatory better matches at other positions.

The StringTie-inferred ATAC1 transcript is 456 nt long and spans positions 23,222,261 to 23,222,716 on Chromosome 11. It is predicted to be composed of a single exon. The longest predicted ORF is 46aa long (Figure S2), and preliminary searches with blastx on NCBI for proteins similar to putative ATAC1 translation products yielded no meaningful hit outside of O. sativa hypothetical proteins. Because this suggested that the ATAC1 transcript does not code for a protein, we used PLncPRO [39] to examine if the ATAC1 transcript is a potential long noncoding RNA (lncRNA). Predicting the ATAC1 transcript status with PlncPRO returned a noncoding probability of 0.965, a value well above 0.8, which is considered the lower bound for high-confidence lncRNAs [39].

The ATAC2/ILI2 locus codes for a putative protein belonging to an atypical group of the bHLH transcription factor family that lacks the basic region of the domain and does not bind DNA. In rice, the ILI clade contains seven homologous genes, including ILI1, its founding member. These genes are phylogenetically and functionally related to members of the Arabidopsis PRE family [40]. Characterized proteins of this group positively regulate cell elongation in response to environmental and hormonal cues, but this signaling module has also been implicated in the regulation of Arabidopsis PAMP-triggered immunity (PTI) [41].

Under the LOC_Os11g39000 locus identifier, ATAC2/ILI2 has previously been shown by Yu et al. [31] to be talC-regulated based both on microarray and independent RT-PCR data. To ascertain that ATAC1 lncRNA induction is both T3SS and talC dependent, we profiled gene expression in response to wild-type BAI3, its talC-mutant derivative BAI3-1-1 [31], and complemented strains in leaves of wild-type Kitaake or the EBE_TalC@OsSWEET14-edited line sweet14–32 [35]. The qRT-PCR data in Figure 2 indicate that the OsSWEET14 control gene is induced on wild-type plants but not on the edited line, in response to the BAI3 strain and the BAI3-1-1 mutant complemented with talC. ATAC transcripts expression follow a similar pattern except that they also appear to be talC induced on the EBE_TalC@OsSWEET14-edited line. This demonstrates that both ATAC1 and ATAC2 are talC regulated and that this effect occurs independently of OsSWEET14 induction.

Altogether, these data support the hypothesis that TalC recognizes a second sequence in the rice genome at the ATAC locus and triggers bidirectional transcription of two opposite units.

2.3. Artificial TALEs Targeting the ATAC Locus Do Not Enhance the Virulence of a Xoo Strain Unable to Induce an OsSWEET Gene

To explain why, in contrast to Asian Xoo, African Xoo virulence is marginally affected on TALE-unresponsive OsSWEET14 promoter-edited lines, we postulated the existence of another TalC susceptibility gene target, redundant with OsSWEET14. Under this model, induction of the redundant factor is expected to compensate for the loss of the susceptibility functions provided by OsSWEET14 induction. If ATAC acted as such a redundant TalC target, we reasoned that the induction of the ATAC locus should rescue the virulence of an Asian strain that is unable to induce OsSWEET functions.

To address this possibility, we assembled a set of four “designer” or “artificial” TALEs (ArTALEs) as before [26]. Their respective EBEs were selected to maximize specificity toward the ATAC loci and as shown in Figure 3a and Figure S3, are located near EBE_TalC@ATAC on the same strand (ArTALE05 and 06) or the opposite, ATAC2, strand (ArTALE01 and 04). These constructs, together with the empty vector or plasmids expressing the control TalF_MAI1 or TalC_BAI3 native TALEs, were introduced into the ME2 strain, a non-virulent PXO99A mutant derivative deprived of the PthXo1 TALE [42] that is routinely employed for testing virulence complementation by OsSWEET targeting TALEs [8,16]. The resulting strains are described in Table 2. To functionally validate these ArTALEs, the transformed strains were infiltrated into Kitaake leaves and their ability to induce gene expression was monitored with qRT-PCR (Figure 3b). The positive control strain expressing TalC induced OsSWEET14, ATAC1, and ATAC2 relative to the empty vector strain, whereas the ME2 strain expressing TalF induced OsSWEET14 only. None of the ArTALEs were capable of dramatically affecting OsSWEET14 expression relative to water-infiltrated plants. In contrast, all four of them induced ATAC2 above empty control levels. The situation regarding ATAC1 is more contrasted. ArTALE05 and 06 induced this transcript at the same level as TalC, whereas ArTALE01 and, to a lesser extent, ArTALE04 were similar to the empty control strain. Overall, it appeared that all our ArTALEs are able to specifically activate ATAC transcription and that those designed to target an EBE with the same orientation as EBE_TalC@ATAC, similar to TalC, acted in a bidirectional manner.

Next, to test if ATAC transcript induction could confer enhanced virulence on the ME2 mutant, we screened those strains in leaf clipping virulence assays on Kitaake plants. As shown in Figure 3c, the TalC- and TalF-expressing strains caused disease lesions longer than 10 cm, well above the empty vector control strain and consistent with previous similar experiments [7,8]. However, with the exception of the ME2 derivative expressing ArTALE04, lesion lengths recorded for other ArTALE-expressing strains were not different from those produced with the empty vector control strain.

We, therefore, conclude that while our ArTALEs could simultaneously induce both ATAC1 and ATAC2 transcripts or ATAC2 alone, they did not consistently compensate for the loss of clade III OsSWEET induction and failed to confer increased virulence on the ME2 mutant strain in our conditions.

2.4. Multiplex Editing of the OsSWEET14 Promoter and the ATAC Locus Using CRISPR/Cas9

Negative results in these gain-of-function experiments could stem from a number of reasons. If ATAC is indeed redundant with OsSWEET14, another prediction of our working model is that TalC-unresponsive plants should display enhanced resistance to African Xoo strains. To investigate the significance of ATAC targeting by TalC with a complementary loss-of-function approach, we therefore applied CRISPR/Cas9 genome editing to inactivate TalC target loci in multiplex. To this end, we employed a previously described system [43] to generate a set of binary constructs with cas9 and an array of 4 gRNAs (Figure S4) for the genetic transformation of the Nipponbare variety. To abrogate TALE-mediated OsSWEET14 induction, all our three constructs bore a gRNA targeting EBE_TalF and another gRNA targeting EBE_TalC@OsSWEET14. In order to query different regions of the ATAC loci, they differed at the level of the two remaining gRNAs: a pair of gRNA targeting EBE_TalC@ATAC or a nearby sequence was included in one construct with the goal of abrogating TalC responsiveness. Additional pairs of gRNAs targeting either ATAC1 or ATAC2 transcript sequences for the insertion of deleterious mutations in the underlying locus, including the deletion of intervening sequences, were cloned in the second or third constructs (Figure S4 and Table S1). Finally, two complementary control constructs designed for single locus editing and comprising only the OsSWEET14 promoter gRNAs or the pair of gRNAs targeting ATAC2 were also included in transformation experiments.

Ultimately, after transformation, genotyping, and propagation of edited plants, we gathered a set of seven homozygous, non-transgenic T3 lines with a range of edition events as described in Figure 4a and Figure S5. The efficiency of EBE_TalF editing was too low to select consistently modified lines. However, all sweet14 alleles carry mutations that are likely to prevent TalC recognition. The atac-1 allele consists in a 100 bp deletion of the ATAC1 transcript, including a portion of the longest predicted ORF. Furthermore, we obtained several alleles with insertions or deletions upstream of position 8 of EBE_TalC@ATAC (attac2, atac4, and atac5) (see Figure 4a and Figure S5). We were, however, unable to recover edited lines combining an unresponsive allele at EBE_TalC@OsSWEET14 and a modified allele of ATAC2.

2.5. Bidirectional Induction of Transcription at ATAC Is Dependent on the Predicted TalC EBE, but TalC-Unresponsive Plants Remain Essentially Susceptible to an African Xoo Strain

Next, we examined TalC-responsiveness of edited lines using qRT-PCR following leaf infiltration of the BAI3 strain (Figure 4b). As opposed to wild-type Nipponbare or the atac-6 line, and consistent with our previous results [35], OsSWEET14 induction was abrogated in all lines carrying a mutation affecting the EBE_TalC on its promoter.

Figure 4

ATAC-locus-edited lines confirm that EBE_TalC@ATAC is functional but does not support a role as a BLB disease susceptibility locus. (a) A schematic diagram describing the nature of the editing event in the modified Nipponbare homozygous lines characterized below. A cross on top of a functional feature indicates that the sequence of the corresponding allele has been altered. Red arrows represent the OsSWEET14 or ATAC genes. Colored segments depict TALE target EBEs. Values in parentheses aligned with a line identifier inform about the extent and approximate location of DNA editions as in Figure S5. Features are not drawn to scale. (b) TalC target transcript expression ratios in multiplex edited lines 48 h post infiltration with BAI3 relative to water-infiltrated control leaves. The plotted data (gray dots) derive from three independent experiments (shades of gray; n = 6–9). The black dots and line ranges denote the mean and the standard error, respectively. Treatments with the same letter in an individual transcript panel are not significantly different (Dunn’s test for post hoc pairwise multiple comparisons; p-value < 0.05). For a couple of genotypes, sequence alterations precluded qRT-PCR on a specific target, and no data are reported (no mean difference letter; see the main text for details). (c) Leaf clipping susceptibility assays performed with the edited lines described in panel A challenged with the BAI3 strain. Lesion lengths were measured 14 days post inoculation. Additional water control leaves did not show BLB symptoms. The plotted data derive from three independent experiments (dots with different gray levels). The boxplots represent the median and the first and third quartiles, whereas the whisker extends from the hinge to the largest or the lowest value but no further than 1.5 times the distance between the first and third quartiles. Treatment comparison was performed separately depending on the wild-type rice background, and those with the same letter are not significantly different (Dunn’s test for post hoc pairwise multiple comparisons; p-value < 0.05). Numbers in parentheses represent the count of values for the corresponding treatment.

[Figure omitted. See PDF]

Plants containing the atac-1 allele were not measured for ATAC1 levels because of the large deletion in the genomic sequence of the transcript. Likewise, the atac-6 plants were not assayed for ATAC2 expression because of the deletion of most of the gene sequences, including EBE_TalC@ATAC (Figure S5). With these potential caveats in mind, plants with the atac-1 allele were not compromised for ATAC2 induction that was similar to the wild-type plants. In contrast, in atac-6 plants, ATAC1 was not responsive to BAI3 inoculations, possibly because EBE_TalC@ATAC is missing in those individuals. In agreement with this view, induction was essentially abrogated for both ATAC transcripts in other lines with more subtle modifications at EBE_TalC@ATAC (atac2, atac4, and atac5 alleles). This is indicative that an intact EBE_TalC@ATAC is necessary for TalC recognition and subsequent divergent transcription activation.

Having determined that our edition events panel encompasses a range of effects relative to the potential function and the regulation of the ATAC loci combined with TalC-unresponsive OsSWEET14 promoters, we asked if this resulted in altered susceptibility to the Xoo strain BAI3. As shown in Figure 4c, while in our leaf clipping assays with BAI3, and as reported previously [7], the double OsSWEET14/OsSWEET13 knockout plants (KO-S13-S14) displayed strongly reduced lesions relative to the wild-type Kitaake controls, the resistance level of our edited plants in the Nipponbare background was not commensurately affected. Even if the lesion lengths for two lines edited in multiplex at both ATAC and OsSWEET14 loci (sweet14–38 atac-1 and sweet14–42 atac-5) were significantly different than wild-type and single edited lines (sweet14–36 and sweet14–37), this effect was modest and not consistent across modified lines with similar anticipated genetic defects.

Overall, the production and phenotypical analysis of a series of ATAC-edited Nipponbare lines revealed that the putative TalC EBE is required for bidirectional induction of transcription at this locus by BAI3 but simultaneous abrogation of TalC responsiveness at both OsSWEET14 and ATAC loci did not confer dramatically enhanced resistance to disease caused by the African Xoo BAI3 strain.

3. Discussion

3.1. TalC Functions and the Dependence of Af Xoo on OsSWEET Induction for Infection

In this study, following up on our earlier results, we investigated potential alternative TalC targets in the rice genome and identified a second TalC-responsive locus coined ATAC. Recent data indicate that multi-EBE-edited Kitaake lines with an altered EBE_TalC@OsSWEET14 are resistant to the ME2 derivative complemented with talC, which is otherwise fully virulent on wild-type Kitaake rice [8]. Similarly, talC no longer provides increased lesions to ME2 on the OsSWEET14 knockout line versus wild-type Kitaake [7]. Thus, at least in the ME2 background, talC contribution to virulence is abrogated when either host OsSWEET14 gene activity or transcriptional induction is compromised. This is consistent with the view that talC virulence activity is exclusively channeled via an OsSWEET14 pathway and contradicts our hypothesis that alternative TalC rice targets, including the ATAC locus, could redundantly contribute to susceptibility.

Considering these observations and our inability to demonstrate a contribution of ATAC transcripts using either ArTALEs (Figure 3c) or edited lines (Figure 4c), we provisionally conclude that transcriptional targeting of ATAC is likely fortuitous and does not contribute to the virulence of African strains causing BLB. In this context, the mild ME2-virulence-enhancing activity of ArTALE04 only (Figure 3c) can only be interpreted as an off-target effect of this particular construct.

A growing list of TALEs has been shown to simultaneously target more than one host genomic locus for transcription activation. Most documented examples, such as AvrXa7 on OsSWEET14 and the transposon Os04g19960 [5], the X. oryzae pv. orzicola Tal2g on rice OsSULTR3;6 and Os06g46500 [44], or the X. citri subsp. citri TALE PthA4 on CsLOB1 and CsSWEET1 [45], involve a genuine susceptibility gene and a second host locus that does not measurably affect disease, as is likely the case with TalC. Nonetheless, dual upregulation of independent bona fide susceptibility genes by single TALEs has been demonstrated. For example, the African Xoo TALE TalB targets OsTFX1 and OsERF#123, both of which increase BLB symptoms upon induction [25]. The tomato (Solanum lycopersicum) transcription factor genes bHLH3 (Solyc03g097820) and bHLH6 (Solyc06g072520) are both regulated by AvrHah1, a TALE from X. gardneri and can individually promote water soaking, a typical disease readout [46].

Incidentally, regardless of whether this impinges on plant–pathogenic bacteria interactions, TALEs appear to repeatedly regulate genes of the large plant basic/helix–loop–helix (bHLH) transcription factor family [47]. In addition to AvrHah1 (mentioned above) and TalC, with the newly discovered ILI2/ATAC2 target reported in this study, an archetypal example of this is AvrBs3 from X. campestris pv. vesicatoria and its bHLH target UPA20, which controls pepper (Capsicum annuum) cell hypertrophy [13]. Many strains of the Xoo-related oryzicola pathovar of X. oryzae broadly upregulate the rice Os03g07540 locus, which also encodes a bHLH protein [48]. With the BLS256 strain, this effect is mediated by the Tal3c TALE [44]. Intriguingly, Os03g07540, which corresponds to ILI3 [40], is therefore a close homolog of ILI2/ATAC2, the alternative TalC target identified here.

3.2. Bidirectional Induction of Transcription by TALEs

Arguably, one of the most important outcomes of our work, with relevance to the molecular mechanisms of TALE-mediated transcription activation in plants, is the observation that the TalC bacterial effector triggers divergent bidirectional transcription at ATAC. Both ATAC1 and ATAC2 transcripts are coordinately induced in a TalC-dependent manner (Figure 2), and this requires an intact putative TalC EBE at the ATAC locus (Figure 4).

ATAC1 is a previously undocumented short (456 nt) rice transcript (Figure 1) whose putative translation products are less than 100 amino acids (Figure S2) and that is predicted to be an lncRNA by a classifier trained on a rice-specific dataset. Plant lncRNAs are comparatively less understood than animal lncRNAs, and they comprise a relatively diverse set of molecules having functions in plant development and adaptation to both abiotic and biotic stresses through a variety of mechanisms [49]. Whether ATAC1 has an intrinsic biological function or is just an accidental by-product of TalC activity is an interesting question that may be illuminated in future work. In line with another study [50], the discovery of the ATAC1 transcript in a reconstructed transcriptome is, nonetheless, a strong hint that the search for TALE target sequences should not solely rely on incomplete reference plant genome annotations and must incorporate experimental RNA-seq information to enable the identification of novel candidates.

Two previous studies have demonstrated that designer TALEs can activate unidirectional transcription of sequences upstream of the EBE on the opposite strand (antisense transcription). This was shown to occur both in transient assays in Nicotiana benthamiana and during leaf infection. ArTALEs designed to target an EBE in an antisense orientation relative to OsSWEET14 [15,16] or OsSULTR3;6 [16] susceptibility genes were able to complement the virulence of strains deficient for the cognate natural TALE. Wang et al. [16] provided additional insight and showed that both artificial and natural TALEs can activate bidirectional transcription in Agrobacterium-mediated transient expression assays in N. benthamiana with a T-DNA carrying a dual reporter construct. The same authors also described two functional EBEs in opposite orientation for the natural X. oryzae pv. orzicola TALE Tal3c upstream of the Os09g39810 locus, but they detected no sign of bidirectional transcription in mRNA-seq expression data during infection. Here, in contrast, our results represent the first tangible evidence that a native TALE can activate bidirectional transcription when delivered by the T3SS during infection.

ATAC1 is reminiscent of enhancer-associated lncRNAs (eRNAs) [51,52]. While in animals, enhancers and gene promoters frequently undergo divergent RNA-polymerase-II-mediated transcription, in plants, the evidence for bidirectional transcription from promoters of protein-coding genes is relatively scarce and previously defined bidirectional plant promoters [53,54,55,56] may entail a different phenomenon. However, hypomorphic mutations in the Arabidopsis exosome complex or its associated helicase system have recently enabled the detection of about a hundred loci with antisense RNAs upstream of mRNA transcripts and an equivalent number of noncoding regions with signs of bidirectional transcription reminiscent of eRNAs [57]. TALE-driven divergent transcription may be a consequence of an idiosyncratic mode of action of these bacterially encoded transcription factors. This property probably depends on the nature of the additional plant transcription initiation machinery components that are recruited and that are just beginning to be deciphered [58]. The current picture is, however, that TALEs predominantly activate unidirectional transcription downstream of their EBEs; for example, there is no signal of antisense transcription at the OsSWEET14 (LOC_Os11g31190) promoter (Figure S1). Thus, as suggested before [15,16], other promoter elements (e.g., the TATA-box) but also DNA conformational properties (dictated in part by chromatin modifications) conceivably determine the local mode of transcription activation (i.e., bidirectional, sense or antisense only) by stimulating the release of halted RNA polymerase II complexes or by inhibiting the exosomal processing of nascent transcripts. In this regard, it is probably not fortuitous that our ArTALE05 and 06 with an EBE on the same strand as EBE_TalC@ATAC trigger divergent transcription whereas ArTALE01 and 04 with an EBE oriented toward ATAC2 activate unidirectional sense transcription only (Figure 3).

In conclusion, pursuing the hypothesis that TalC targeted a second redundant susceptibility gene in rice, we did identify a second locus but failed to demonstrate its contribution to disease susceptibility. Hence, the basis of the unique ability of African Xoo to persistently colonize rice in the absence of clade III OsSWEET gene induction remains enigmatic and shall be investigated further. In this endeavor, we nonetheless documented the first example of bidirectional transcription activation by a native TALE during infection, which will be a valuable system for a better understanding of the molecular mechanisms underlying TALE activity.

4. Materials and Methods

4.1. Bacterial Strains, Plant Inoculations, and Growth of Bacteria in Planta

The wild-type rice genotypes Kitaake and Nipponbare varieties used in this study were originally obtained from the International Rice Research Institute and the IRGSP consortium [37], respectively. The OsSWEET14/OsSWEET13 double-knockout edited line corresponds to the sweet13;14 +1bp/+1bp line from Eom et al. [7]. All rice plants were grown under cycles of 12 h of light at 28 °C and 80% relative humidity (RH) and 12 h of darkness at 25 °C and 70% RH in the greenhouse.

The plasmids and Xoo bacterial strains used in this study are described in Tables S1 and S2, respectively. Bacterial strains stored in 15% glycerol at −80 °C were cultivated on PSA (10 g/L peptone, 10 g/L sucrose, 1 g/L glutamic acid, and 16 g/L Bacto Agar). The strains containing plasmids were grown on PSA supplemented as needed with kanamycin, gentamicin, or rifampicin at concentrations of 50 μg/mL, 20 μg/mL, and 50 μg/mL, respectively. The cultures were incubated at 28 °C for 48 h.

Inoculations were conducted as before [31]. Briefly, bacteria grown on plates for 48 h were resuspended in sterile water at an OD600 = 0.2 for leaf clipping assay inoculum preparation or OD600 = 0.5 for leaf infiltration. Leaves of 6-week-old plants were cut with scissors previously dipped in bacterial suspensions, whereas leaves of 4-week-old plants were pressure infiltrated with a 1 mL syringe. Lesion lengths were measured 14 days post inoculation.

4.2. RNA Isolation and qRT-PCR

The leaves of Kitaake, Nipponbare, or edited rice line plants were infiltrated with bacterial suspensions or water, and three independent replicate samples were processed at 48 h after inoculation. Total RNA was extracted using TRI reagent (EUROMEDEX, Souffelweyersheim, France) and further purified using the RNA Clean-Up & Concentration kit (Zymo Research, Irvine, CA, USA) according to the manufacturer’s instructions. A total of 1 μg of RNA per sample was treated with the TURBO DNA-free kit (Thermo Fisher Scientific, Waltham, MA, USA) and reverse-transcribed into cDNA using SuperScriptIII (Thermo Fisher Scientific, Waltham, MA, USA) with oligo-dT primers. PCR reactions were performed with SYBR-based Mesa Blue qPCR Mastermix (EUROGENTEC, Liège, Belgium).

For the OsSWEET14 gene, qRT-PCR oligos were reported before [35]. For the ATAC locus, a primer pair specific to the ATAC1 transcript (GAGGTGGAGCTGGAGCTATG and TACAACGGGCTACAACCACA) and a primer pair specific to the ATAC2 transcript (CGTCCGATCTCGCCATCAAC and GAGAACATTTTAGGCCGACAGCA) were used. The elongation gene factor-1 alpha (GenBank: GQ848072.1) was used for normalization of the expression values of complementary DNA in each sample as before [35]. Normalized expression values were calculated with the 2^−ΔCt methods, whereas expression ratios were calculated with the 2^−ΔΔCt method.

4.3. Plasmid Vector Construction

The design and construction of artificial TAL effector (ArTALEs) expression vectors involved the use of the method of Streubel et al. [26].

For rice genome editing, first, the pHD46 vector was generated by PCR amplifying the dual gRNA cassette from pENTR4:gRNA4 [43] with a pair of primers wherein the forward primer (ACTTAAGCTTGAATTCagaacgaactaagccggaca) comprises the HindIII and EcoRI restriction sites and the reverse primer (ACCTTTCTAGAgcccttcgaagggacaaaaa) comprises an XbaI site. The amplified sgRNA cassette was digested and gel-purified using the DNA Gel Recovery Kit (Zymo Research, Irvine, CA, USA). This new cassette was cloned into the original pENTR4:gRNA4 vector digested at the EcoRI and XbaI restriction sites that are located in between the attL sequences.

The E-CRISPR software [59] was used to design synthetic guide RNAs (sgRNAs). The sgRNA sequences are indicated in Figure S4. When the designed gRNA sequence 5′ nucleotide was not a guanine (G), a starting G was appended to this sequence to promote U6 promoter transcription initiation at this site. Complementary oligos spanning sgRNA sequences with appropriate 4 bp overhangs were thermally annealed to each other, diluted (1:200) with distilled water, and cloned into pENTR4:gRNA4 or pHD46 vectors using the procedure described by Zhang et al. [43].

For the assembly of an array of four gRNA expression units into a Gateway entry vector, pHD46-derived EcoRI-XbaI cassettes with ATAC locus targeting gRNAs were cloned into pHD23, a pENTR4:gRNA4-based vector bearing the OsSWEET14 promoter targeting gRNAs, digested with the same enzymes. This array was subsequently transferred into a modified pH-Ubi-cas9-7 [60] vector documented by Fayos et al. [61] in a Gateway LR reaction (Thermo Fisher Scientific, Waltham, MA, USA). The gRNA arrays of the resulting binary vectors (see Table S1) were confirmed again by sequencing before Agrobacterium transformation was continued.

4.4. Genetic Transformation of Rice and Modified Allele Genotyping

For rice genetic transformation, Agrobacterium tumefaciens EHA105 was transformed with the set of binary vectors by electroporation. The embryo calli derived from Nipponbare were used and co-cultured with bacteria as described [35,62]. T0 plants were screened to characterize multiplex editing events. The MATAB protocol [63] was used to extract DNA from leaf samples. The amplification products of the OsSWEET14 promoter were obtained as described before [35]. The ATAC locus was PCR amplified with specific forward (GACTGTAGCTTTGCCTAGCT) and reverse (CTACAGATGAAGAAGCTAGAG) primers. The resulting amplicons were Sanger sequenced and analyzed with DSDecode software [64] to identify mutations within the OsSWEET14 promoter. For simultaneous detection of mutation of the OsSWEET14 promoter and ATAC locus, the PCR products were cloned into the pGMT-T vector according to the instructions of the manufacturer (Promega, Madison, WI, USA). The clones were sequenced and analyzed with DSDecode software to decode the mutations.

4.5. Computational Methods for TalC Candidate Target Identification

The mRNA-seq data used in this study for experimental transcriptome inference and differential expression analysis are available in SRA (NCBI BioProject accession: PRJNA679478) and have been described before [36].

Paired-end reads were trimmed with AdapterRemoval version 2.2.2 [65] using the following relevant parameters: “--combined-output --gzip --trimns --trimqualities --mm 5 --minadapteroverlap 3 --minlength 25”. Trimmed reads were subsequently mapped onto the nuclear rice genome. Rice Nipponbare IRGSP-1.0 pseudomolecule sequences and MSU7 gene annotation data [37] were obtained from MSU (http://rice.uga.edu/ (accessed on 10 March 2022)). For mapping, HISAT2 (v.2.1.0) was run with the following parameters: “-q --new-summary –dta”. Experimental transcriptome reconstruction essentially followed the recommended procedure [66]. Transcripts were first inferred based on mRNA-seq mappings and the MSU7 reference annotation on an individual sample basis with StringTie (v.1.3.4d) and the following parameters: “-l GXY-1 --fr -m 100 -t -f 0.1 -c 5 -j 2”. Individual gtf files were subsequently merged with the reference annotation in a single step with the StringTie --merge mode and default parameter values.

Gene-level and strand-aware mated read pair counts were obtained with R (v.3.5.1) using the summarizeOverlaps function from the GenomicAlignments (v.1.18.1) Bioconductor package [67], discarding reads overlapping multiple features (“Union” mode and “inter.feature” true). StringTie gene loci that accumulated less than 5 total reads across all libraries in the dataset were excluded. Likewise, loci that overlapped rice loci with feature types tRNA, rRNA, or pseudogenic_tRNA in the gff annotation file for NCBI RefSeq assembly accession GCF_001433935.1 were excluded from differential expression (DE) analysis.

Differential expression was estimated with the DESeq2 (v.1.22.2) R Bioconductor package [68]. For normalization, size factors were estimated using the iterate method and the shorth function to compute a location for a sample. Dispersions were fitted using the parametric method. Two-tailed Wald tests were used for DE tests where the alternative hypothesis was that the absolute value of the log2 transform of the fold change ratio is greater than 0.5. Adjusted p-values were computed with the Benjamini & Hochberg method and considered significant when equal to or below 0.1. Outlier filtering based on Cook’s distance was implemented.

To generate a list of upstream sequences for EBE prediction, the DESeq2 gene list was arbitrarily further filtered to retain only the loci with a log2(FC) value ≥ 2 and an adjusted p-value ≤ 0.01 in both BAI3 vs. BAI3H and BAI3 vs. H₂O comparisons. EBE predictions for the TalC RVD sequence on this set of upstream sequences (−500 bp to +200 bp relative to annotated gene start) used Talvez (v.3.2) with parameters “-t 0 -l 19” and default RVD–nucleotide association matrixes [38].

The microarray expression data used to define TalC-induced genes were published earlier [31] and analyzed as previously described [38]. Only probesets with a log2 transform of the fold change ratio ≥ 1 in the XooBAI3-XooBAI3_Dtalc comparison were considered.

The PlncPRO [39] model for O. sativa was built using data and instructions from https://github.com/urmi-21/PlncPRO/blob/master/examples/Build%20and%20test%20monocot%20model.ipynb (accessed on 23 July 2021). The plncpro predict (“-v -r”) command was run using the StringTie-predicted ATAC1 transcript sequence.

Data plots were generated in R (v. 3.6.3) with the ggplot2 (v.3.3.2) or Gviz (v.1.30.3) packages, and statistical tests were conducted with functions from the rstatix (v.0.6.0) package.

Author Contributions

Conceptualization, H.D., F.A., S.B.-B. and S.C.; methodology, H.D., F.A., S.B.-B., S.F., C.S., B.S., V.V. and S.C.; software, M.G. and S.C.; validation, H.D., F.A. and S.C.; formal analysis, H.D., F.A., S.B.-B. and S.C; investigation, H.D., F.A., S.B.-B., S.F., E.T., F.D. and C.S.; resources, B.S.; data curation, H.D., F.A., S.B.-B. and S.C.; writing—original draft preparation, H.D., F.A. and S.C.; writing—review and editing, H.D., F.A., S.B.-B., S.F., M.G., E.T., F.D., C.S., B.S., O.K., V.V. and S.C.; visualization, H.D., F.A. and S.C.; supervision, O.K., V.V. and S.C.; project administration, O.K., V.V. and S.C.; funding acquisition, H.D., O.K., V.V. and S.C. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and supplementary material on our dataverse at https://doi.org/10.23708/XBTRUV, last accessed on 23 July 2021.

Acknowledgments

The authors acknowledge the ISO 9001-certified IRD i-Trop HPC (South Green Platform) at IRD Montpellier for providing high-performance computing resources.

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Transcripts at a second candidate target locus are potentially induced bidirectionally by TalC. (a) Genome browser snapshot of the ATAC locus with transcript annotation from MSU7, the experimental transcriptome annotation inferred by StringTie, the predicted TalC EBE and read coverage tracks (reads per million) for various experimental treatments in our mRNA-seq data. For the MSU7 annotation, shorter exon boxes correspond to untranslated regions. For coverage tracks, only the reads corresponding to the transcribed strand in paired-end data for a specific treatment were considered. Sequencing libraries were generated with RNA extracted 24 h after infiltration of Nipponbare leaves with the strains BAI3, BAI3H (T3SS-), and H2O. Coverage of reads mapping to the top genomic strand is represented with positive values and a golden-colored area, whereas coverage of reads mapping to the opposite strand is represented with negative values and a violet-colored area. (b) Expression of TalC target genes in our mRNA-seq data. Black dots correspond to individual replicate normalized count values, while the red dots and line ranges denote the mean and the standard error, respectively. (c) Alignment of the TalC RVD sequence and putative DNA target sequences. The numbers on top refer to the classical numbering of repeats in the TALE central repeat domain for TalC. The EBE Opti sequence was generated by concatenating the most probable nucleotide for each RVD in the RVD–nucleotide association matrix used by Talvez. In the alignment, nucleotides are colored according to their match to the corresponding TalC RVD in this matrix: red, best score; orange, intermediate; and blue, worst score.

Figure 2. ATAC transcript induction in response to BAI3 is talC dependent. Transcript expression of ATAC1, ATAC2, and SWEET14 on the wild-type Kitaake and the edited sweet14–32 lines measured 48 h after inoculation with derivatives of the BAI3 strain or water. The strains are represented on the x-axis: BAI3.1.1 is a BAI3 derivative defective for talC. “_TalC” or “_EV” suffixes indicate, respectively, that the talC strain carried a plasmid expressing TalC or the corresponding empty vector. Expression levels of the different transcripts measured by qRT-PCR are expressed in arbitrary units and calculated by the ∆Ct method. Individual biological replicate values (n = 3) are denoted with colored symbols as a function of the inoculated plant line. The black dots and line ranges denote the mean and the standard error, respectively.

Figure 3. Independent ArTALEs targeting the ATAC loci are functional but do not compensate the loss of clade III OsSWEET induction for the virulence of the Asian PXO99 strain pthXo1 mutant derivative. (a) Diagram depicting the orientations (color and arrowhead) and approximate locations of the different EBEs targeted by ArTALEs relative to the ATAC locus transcripts. Numbers on top of ArTALE EBEs identify individual ArTALEs in panels B and C below. (b) TalC target transcript expression ratios in response to ArTALEs designed to target the ATAC locus. qRT-PCR was performed with leaf samples collected 48 h post infiltration of the Kitaake wild-type rice line with derivatives of the ME2 (PXO99 mutated in the major TALE gene pthXo1) strain carrying an empty vector (EV) or related plasmids encoding the indicated TALEs. The plotted data (gray dots) derive from two independent experiments (n = 4–6). The black dots and line ranges denote the mean and the standard error, respectively. Treatments with the same letter in an individual transcript panel are not significantly different (Dunn’s test for post hoc pairwise multiple comparisons; p-value < 0.05). (c) Lesion length measured on wild-type rice line Kitaake 14 days post leaf clipping inoculation with the indicated strains. The plotted data derive from three independent experiments (dots with different gray levels). Additional water control leaves did not show BLB symptoms. The boxplots represent the median and the first and third quartiles, whereas the whisker extends from the hinge to the largest or the lowest value but no further than 1.5 times the distance between the first and third quartiles. Treatments with the same letter in an individual transcript panel are not significantly different (Dunn’s test for post hoc pairwise multiple comparisons; p-value < 0.05). Numbers in parentheses represent the count of values for the corresponding treatment.

Table 1

Candidate TalC targets in the Nipponbare genome.

EBE Position on IRGSP-1.0				EBE on Subject [−500, +200]			Talvez Score	Candidate Target Gene Annotation				mRNA-Seq Diff. Exp.BAI3 vs. BAI3H			Micro-Array Diff. Exp BAI3 vs. BAI3-1-1 *
Seq name	start	end	strand	start	end	pol ‡		id	ref_gene_id	name_Gramene	description_Gramene	log2FC	pvalue	padj	log2FC	padj
Chr11	18,174,547	18,174,569	-	410	432	S	14.45	PIX.6085	LOC_Os11g31190	11N3	Sugar transporter, TAL effector-mediated susceptibility to bacterial pathogen (Os11t0508600-01)	6.869	2.55 × 10−45	8.05 × 10−42	5.97	1.03 × 10−1
Chr12	17,463,497	17,463,519	+	577	599	S	13.86	PIX.7679	NA §	NA	NA	5.446	1.13 × 10−12	4.82 × 10−10	NA	NA
Chr5	29,688,055	29,688,077	-	106	128	AS	12.36	PIX.19089	LOC_Os05g51750	NA	Peptidase A1 domain containing protein (Os05t0596000-01)	7.077	5.47 × 10−37	9.09 × 10−34	NA	NA
Chr11	23,222,769	23,222,791	-	426	448	S	11.27	PIX.6359	NA	NA	NA	4.623	2.81 × 10−3	2.19 × 10−1	NA	NA
Chr11	23,222,769	23,222,791	-	363	385	AS	11.27	PIX.6360	LOC_Os11g39000	INCREASED LEAF INCLINATION 2	Basic helix–loop–helix dimerization region bHLH domain containing protein (Os11t0603000-01)	6.998	1.09 × 10−13	5.08 × 10−11	4.26	6.29 × 10−3
Chr7	21,798,885	21,798,907	-	496	518	S	10.83	PIX.22524	LOC_Os07g36430	NA	Conserved hypothetical protein (Os07t0549100-01)	4.66	3.82 × 10−3	2.83 × 10−1	NA	NA
Chr11	28,760,402	28,760,424	+	84	106	AS	10.74	PIX.6663	LOC_Os11g47600	NA	Similar to Class III chitinase homologue (OsChib3H-h) (Fragment) (Os11t0702100-01)	2.432	7.94 × 10−7	1.77 × 10−4	NA	NA
Chr10	18,168,970	18,168,992	+	369	391	S	10.46	LOC_Os10g34040	LOC_Os10g34040	FLOT1	Similar to Flotillin-like protein 1 (Os10t0481500-00)	4.599	1.68 × 10−3	1.39 × 10−1	NA	NA
Chr4	25,888,105	25,888,127	+	659	681	AS	10.39	PIX.16030	LOC_Os04g43730	WALL-ASSOCIATED KINASE GENE 51	Similar to OSIGBa0145M07.8 protein (Os04t0517700-01)	7.301	2.80 × 10−3	2.19 × 10−1	NA	NA
Chr4	15,500,717	15,500,739	-	360	382	S	10.22	PIX.15197	LOC_Os04g26550	NA	Conserved hypothetical protein (Os04t0332700-01)	3.877	5.94 × 10−13	2.61 × 10−10	NA	NA

* Dataset GSE19844; ^§ Not Available; ^‡ Polarity; S, sense; AS, antisense.

Table 2

List of X. oryzae strains used in this study.

Designation	Strain ID	Description	References
BAI3	CIX151	A rifampicin-resistant and fully virulent derivative of the BAI3 wild-type Xoo strain from Burkina Faso	[33]
BAI3H	CIX3878	CIX151 deleted for a portion of the hrcC gene	[31]
BAI3.1.1	CIX531	Also referred to as BAI3R ΔtalC; corresponds to CIX151 with a suicide plasmid insertion in talC	[31]
BAI3.1.1_pEV	CIX2533	CIX531 carrying the empty pSKX1 plasmid	[26]
BAI3.1.1_pTalC	CIX2529	CIX531 carrying the pSKX1plasmid containing the talC gene from strain BAI3	[26]
ME2	CIX4497	PXO99A derivative with a knockout mutation in pthXo1	[42]
ME2_EV	CIX4498	CIX4497 carrying the empty pSKX1 plasmid	This study
ME2_pTALC	CIX4499	CIX4497 carrying the pSKX1plasmid containing the talC gene from strain BAI3	This study
ME2_pTALF	CIX4500	CIX4497 carrying the pSKX1plasmid containing the talF gene from strain MAI1	This study
ME2_pArTAL01	CIX4501	CIX4497 carrying the pHD62 plasmid expressing ArtTale1	This study
ME2_pArTAL04	CIX4502	CIX4497 carrying the pHD65 plasmid expressing ArtTale4	This study
ME2_pArTAL05	CIX4503	CIX4497 carrying the pHD66 plasmid expressing ArtTale5	This study
ME2_pArTAL06	CIX4504	CIX4497 carrying the pHD67 plasmid expressing ArtTale6	This study

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms23105559/s1.

References

1. Niño-Liu, D.O.; Ronald, P.C.; Bogdanove, A.J. Xanthomonas oryzae pathovars: Model pathogens of a model crop. Mol. Plant Pathol.; 2006; 7, pp. 303-324. [DOI: https://dx.doi.org/10.1111/j.1364-3703.2006.00344.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20507449]

2. Savary, S.; Willocquet, L.; Pethybridge, S.J.; Esker, P.; McRoberts, N.; Nelson, A. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol.; 2019; 3, pp. 430-439. [DOI: https://dx.doi.org/10.1038/s41559-018-0793-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30718852]

3. Jiang, N.; Yan, J.; Liang, Y.; Shi, Y.; He, Z.; Wu, Y.; Zeng, Q.; Liu, X.; Peng, J. Resistance Genes and their Interactions with Bacterial Blight/Leaf Streak Pathogens (Xanthomonas oryzae) in Rice (Oryza sativa L.)—An Updated Review. Rice; 2020; 13, 3. [DOI: https://dx.doi.org/10.1186/s12284-019-0358-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31915945]

4. Kumar, A.; Kumar, R.; Sengupta, D.; Das, S.N.; Pandey, M.K.; Bohra, A.; Sharma, N.K.; Sinha, P.; Sk, H.; Ghazi, I.A. et al. Deployment of Genetic and Genomic Tools Toward Gaining a Better Understanding of Rice-Xanthomonasoryzae pv. oryzae Interactions for Development of Durable Bacterial Blight Resistant Rice. Front. Plant Sci.; 2020; 11, 1152. [DOI: https://dx.doi.org/10.3389/fpls.2020.01152] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32849710]

5. Li, T.; Liu, B.; Spalding, M.H.; Weeks, D.P.; Yang, B. High-efficiency TALEN-based gene editing produces disease-resistant rice. Nat. Biotechnol.; 2012; 30, pp. 390-392. [DOI: https://dx.doi.org/10.1038/nbt.2199] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22565958]

6. Duy, P.N.; Lan, D.T.; Thu, H.P.; Thu, H.P.T.; Thanh, H.N.; Pham, N.P.; Auguy, F.; Manh, T.B.; Cunnac, S.; Pham, X.H. Improved bacterial leaf blight disease resistance in the major elite Vietnamese rice cultivar TBR225 via editing of the OsSWEET14 promoter. PLoS ONE; 2021; 16, e0255470. [DOI: https://dx.doi.org/10.1371/journal.pone.0255470]

7. Eom, J.S.; Luo, D.; Atienza-Grande, G.; Yang, J.; Ji, C.; Thi Luu, V.; Huguet-Tapia, J.C.; Char, S.N.; Liu, B.; Nguyen, H. et al. Diagnostic kit for rice blight resistance. Nat. Biotechnol.; 2019; 37, pp. 1372-1379. [DOI: https://dx.doi.org/10.1038/s41587-019-0268-y]

8. Oliva, R.; Ji, C.; Atienza-Grande, G.; Huguet-Tapia, J.C.; Perez-Quintero, A.; Li, T.; Eom, J.S.; Li, C.; Nguyen, H.; Liu, B. et al. Broad-spectrum resistance to bacterial blight in rice using genome editing. Nat. Biotechnol.; 2019; 37, pp. 1344-1350. [DOI: https://dx.doi.org/10.1038/s41587-019-0267-z]

9. Xu, Z.; Xu, X.; Gong, Q.; Li, Z.; Li, Y.; Wang, S.; Yang, Y.; Ma, W.; Liu, L.; Zhu, B. et al. Engineering Broad-Spectrum Bacterial Blight Resistance by Simultaneously Disrupting Variable TALE-Binding Elements of Multiple Susceptibility Genes in Rice. Mol. Plant; 2019; 12, pp. 1434-1446. [DOI: https://dx.doi.org/10.1016/j.molp.2019.08.006]

10. Zafar, K.; Khan, M.Z.; Amin, I.; Mukhtar, Z.; Yasmin, S.; Arif, M.; Ejaz, K.; Mansoor, S. Precise CRISPR-Cas9 Mediated Genome Editing in Super Basmati Rice for Resistance Against Bacterial Blight by Targeting the Major Susceptibility Gene. Front. Plant Sci.; 2020; 11, 575. [DOI: https://dx.doi.org/10.3389/fpls.2020.00575]

11. Zeng, X.; Luo, Y.; Vu, N.T.Q.; Shen, S.; Xia, K.; Zhang, M. CRISPR/Cas9-mediated mutation of OsSWEET14 in rice cv. Zhonghua11 confers resistance to Xanthomonas oryzae pv. oryzae without yield penalty. BMC Plant Biol.; 2020; 20, 313. [DOI: https://dx.doi.org/10.1186/s12870-020-02524-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32620100]

12. Kay, S.; Hahn, S.; Marois, E.; Wieduwild, R.; Bonas, U. Detailed analysis of the DNA recognition motifs of the Xanthomonas type III effectors AvrBs3 and AvrBs3Deltarep16. Plant J.; 2009; 59, pp. 859-871. [DOI: https://dx.doi.org/10.1111/j.1365-313X.2009.03922.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19473322]

13. Kay, S.; Hahn, S.; Marois, E.; Hause, G.; Bonas, U. A Bacterial Effector Acts as a Plant Transcription Factor and Induces a Cell Size Regulator. Science; 2007; 318, pp. 648-651. [DOI: https://dx.doi.org/10.1126/science.1144956] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17962565]

14. Römer, P.; Strauss, T.; Hahn, S.; Scholze, H.; Morbitzer, R.; Grau, J.; Bonas, U.; Lahaye, T. Recognition of AvrBs3-Like Proteins Is Mediated by Specific Binding to Promoters of Matching Pepper Bs3 Alleles. Plant Physiol.; 2009; 150, pp. 1697-1712. [DOI: https://dx.doi.org/10.1104/pp.109.139931] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19448036]

15. Streubel, J.; Baum, H.; Grau, J.; Stuttmann, J.; Boch, J. Dissection of TALE-dependent gene activation reveals that they induce transcription cooperatively and in both orientations. PLoS ONE; 2017; 12, e0173580. [DOI: https://dx.doi.org/10.1371/journal.pone.0173580]

16. Wang, L.; Rinaldi, F.C.; Singh, P.; Doyle, E.L.; Dubrow, Z.E.; Tran, T.T.; Perez-Quintero, A.L.; Szurek, B.; Bogdanove, A.J. TAL Effectors Drive Transcription Bidirectionally in Plants. Mol. Plant; 2017; 10, pp. 285-296. [DOI: https://dx.doi.org/10.1016/j.molp.2016.12.002]

17. Boch, J.; Scholze, H.; Schornack, S.; Landgraf, A.; Hahn, S.; Kay, S.; Lahaye, T.; Nickstadt, A.; Bonas, U. Breaking the Code of DNA Binding Specificity of TAL-Type III Effectors. Science; 2009; 326, pp. 1509-1512. [DOI: https://dx.doi.org/10.1126/science.1178811]

18. Moscou, M.J.; Bogdanove, A.J. A Simple Cipher Governs DNA Recognition by TAL Effectors. Science; 2009; 326, 1501. [DOI: https://dx.doi.org/10.1126/science.1178817]

19. Streubel, J.; Blücher, C.; Landgraf, A.; Boch, J. TAL effector RVD specificities and efficiencies. Nat. Biotechnol.; 2012; 30, pp. 593-595. [DOI: https://dx.doi.org/10.1038/nbt.2304]

20. de Lange, O.; Schreiber, T.; Schandry, N.; Radeck, J.; Braun, K.H.; Koszinowski, J.; Heuer, H.; Strauß, A.; Lahaye, T. Breaking the DNA-binding code of Ralstonia solanacearum TAL effectors provides new possibilities to generate plant resistance genes against bacterial wilt disease. New Phytol.; 2013; 199, pp. 773-786. [DOI: https://dx.doi.org/10.1111/nph.12324]

21. Doyle, E.L.; Hummel, A.W.; Demorest, Z.L.; Starker, C.G.; Voytas, D.F.; Bradley, P.; Bogdanove, A.J. TAL effector specificity for base 0 of the DNA target is altered in a complex, effector- and assay-dependent manner by substitutions for the tryptophan in cryptic repeat -1. PLoS ONE; 2013; 8, e82120. [DOI: https://dx.doi.org/10.1371/journal.pone.0082120] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24312634]

22. Gao, H.; Wu, X.; Chai, J.; Han, Z. Crystal structure of a TALE protein reveals an extended N-terminal DNA binding region. Cell Res.; 2012; 22, pp. 1716-1720. [DOI: https://dx.doi.org/10.1038/cr.2012.156] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23147789]

23. Mak, A.N.-S.; Bradley, P.; Cernadas, R.A.; Bogdanove, A.J.; Stoddard, B.L. The Crystal Structure of TAL Effector PthXo1 Bound to Its DNA Target. Science; 2012; 335, pp. 716-719. [DOI: https://dx.doi.org/10.1126/science.1216211] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22223736]

24. White, F.F.; Yang, B. Host and Pathogen Factors Controlling the Rice-Xanthomonas oryzae Interaction. Plant Physiol.; 2009; 150, pp. 1677-1686. [DOI: https://dx.doi.org/10.1104/pp.109.139360] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19458115]

25. Tran, T.T.; Pérez-Quintero, A.L.; Wonni, I.; Carpenter, S.C.D.; Yu, Y.; Wang, L.; Leach, J.E.; Verdier, V.; Cunnac, S.; Bogdanove, A.J. et al. Functional analysis of African Xanthomonas oryzae pv. oryzae TALomes reveals a new susceptibility gene in bacterial leaf blight of rice. PLoS Pathog.; 2018; 14, e1007092. [DOI: https://dx.doi.org/10.1371/journal.ppat.1007092] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29864161]

26. Streubel, J.; Pesce, C.; Hutin, M.; Koebnik, R.; Boch, J.; Szurek, B. Five phylogenetically close rice SWEET genes confer TAL effector-mediated susceptibility to Xanthomonas oryzae pv. oryzae. New Phytol.; 2013; 200, pp. 808-819. [DOI: https://dx.doi.org/10.1111/nph.12411]

27. Bezrutczyk, M.; Yang, J.; Eom, J.S.; Prior, M.; Sosso, D.; Hartwig, T.; Szurek, B.; Oliva, R.; Vera-Cruz, C.; White, F.F. et al. Sugar flux and signaling in plant-microbe interactions. Plant J.; 2018; 93, pp. 675-685. [DOI: https://dx.doi.org/10.1111/tpj.13775]

28. Yang, B.; Sugio, A.; White, F.F. Os8N3 is a host disease-susceptibility gene for bacterial blight of rice. Proc. Natl. Acad. Sci. USA; 2006; 103, pp. 10503-10508. [DOI: https://dx.doi.org/10.1073/pnas.0604088103]

29. Zhou, J.; Peng, Z.; Long, J.; Sosso, D.; Liu, B.; Eom, J.-S.; Huang, S.; Liu, S.; Cruz, C.V.; Frommer, W. et al. Gene targeting by the TAL effector PthXo2 reveals cryptic resistance gene for bacterial blight of rice. Plant J.; 2015; 82, pp. 632-643. [DOI: https://dx.doi.org/10.1111/tpj.12838]

30. Antony, G.; Zhou, J.; Huang, S.; Li, T.; Liu, B.; White, F.; Yang, B. Rice xa13 Recessive Resistance to Bacterial Blight Is Defeated by Induction of the Disease Susceptibility Gene Os-11N3. Plant Cell; 2010; 22, pp. 3864-3876. [DOI: https://dx.doi.org/10.1105/tpc.110.078964]

31. Yu, Y.; Streubel, J.; Balzergue, S.; Champion, A.; Boch, J.; Koebnik, R.; Feng, J.; Verdier, V.; Szurek, B. Colonization of rice leaf blades by an African strain of Xanthomonas oryzae pv. oryzae depends on a new TAL effector that induces the rice nodulin-3 Os11N3 gene. Mol. Plant Microbe Interact.; 2011; 24, pp. 1102-1113. [DOI: https://dx.doi.org/10.1094/MPMI-11-10-0254] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21679014]

32. Doucouré, H.; Pérez-Quintero, A.L.; Reshetnyak, G.; Tekete, C.; Auguy, F.; Thomas, E.; Koebnik, R.; Szurek, B.; Koita, O.; Verdier, V. et al. 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. Front. Microbiol.; 2018; 9, 1657. [DOI: https://dx.doi.org/10.3389/fmicb.2018.01657] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30127769]

33. Gonzalez, C.; Szurek, B.; Manceau, C.; Mathieu, T.; Séré, Y.; Verdier, V. Molecular and Pathotypic Characterization of New Xanthomonas oryzae Strains from West Africa. Mol. Plant Microbe Interact.; 2007; 20, pp. 534-546. [DOI: https://dx.doi.org/10.1094/MPMI-20-5-0534] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17506331]

34. Poulin, L.; Grygiel, P.; Magne, M.; Gagnevin, L.; Rodriguez-R, L.M.; Serna, N.F.; Zhao, S.; El Rafii, M.; Dao, S.; Tekete, C. et al. New Multilocus Variable-Number Tandem-Repeat Analysis Tool for Surveillance and Local Epidemiology of Bacterial Leaf Blight and Bacterial Leaf Streak of Rice Caused by Xanthomonas oryzae. Appl. Environ. Microbiol.; 2015; 81, pp. 688-698. [DOI: https://dx.doi.org/10.1128/AEM.02768-14]

35. Blanvillain-Baufumé, S.; Reschke, M.; Solé, M.; Auguy, F.; Doucoure, H.; Szurek, B.; Meynard, D.; Portefaix, M.; Cunnac, S.; Guiderdoni, E. et al. Targeted promoter editing for rice resistance to Xanthomonas oryzae pv. oryzae reveals differential activities for SWEET14-inducing TAL effectors. Plant Biotechnol. J.; 2017; 15, pp. 306-317. [DOI: https://dx.doi.org/10.1111/pbi.12613]

36. Reshetnyak, G.; Jacobs, J.M.; Auguy, F.; Sciallano, C.; Claude, L.; Medina, C.; Perez-Quintero, A.L.; Comte, A.; Thomas, E.; Bogdanove, A. et al. An atypical class of non-coding small RNAs is produced in rice leaves upon bacterial infection. Sci Rep.; 2021; 11, 24141. [DOI: https://dx.doi.org/10.1038/s41598-021-03391-9]

37. Kawahara, Y.; de la Bastide, M.; Hamilton, J.P.; Kanamori, H.; McCombie, W.R.; Ouyang, S.; Schwartz, D.C.; Tanaka, T.; Wu, J.; Zhou, S. et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice; 2013; 6, 4. [DOI: https://dx.doi.org/10.1186/1939-8433-6-4]

38. Perez-Quintero, A.L.; Rodriguez-R, L.M.; Dereeper, A.; López, C.; Koebnik, R.; Szurek, B.; Cunnac, S. An Improved Method for TAL Effectors DNA-Binding Sites Prediction Reveals Functional Convergence in TAL Repertoires of Xanthomonas oryzae Strains. PLoS ONE; 2013; 8, e68464. [DOI: https://dx.doi.org/10.1371/journal.pone.0068464]

39. Singh, U.; Khemka, N.; Rajkumar, M.S.; Garg, R.; Jain, M. PLncPRO for prediction of long non-coding RNAs (lncRNAs) in plants and its application for discovery of abiotic stress-responsive lncRNAs in rice and chickpea. Nucleic Acids Res.; 2017; 45, e183. [DOI: https://dx.doi.org/10.1093/nar/gkx866]

40. Zhang, L.Y.; Bai, M.Y.; Wu, J.; Zhu, J.Y.; Wang, H.; Zhang, Z.; Wang, W.; Sun, Y.; Zhao, J.; Sun, X. et al. Antagonistic HLH/bHLH transcription factors mediate brassinosteroid regulation of cell elongation and plant development in rice and Arabidopsis. Plant Cell; 2009; 21, pp. 3767-3780. [DOI: https://dx.doi.org/10.1105/tpc.109.070441]

41. Fan, M.; Bai, M.Y.; Kim, J.G.; Wang, T.; Oh, E.; Chen, L.; Park, C.H.; Son, S.H.; Kim, S.K.; Mudgett, M.B. et al. The bHLH transcription factor HBI1 mediates the trade-off between growth and pathogen-associated molecular pattern-triggered immunity in Arabidopsis. Plant Cell; 2014; 26, pp. 828-841. [DOI: https://dx.doi.org/10.1105/tpc.113.121111] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24550223]

42. Yang, B.; White, F.F. Diverse members of the AvrBs3/PthA family of type III effectors are major virulence determinants in bacterial blight disease of rice. Mol. Plant Microbe Interact.; 2004; 17, pp. 1192-1200. [DOI: https://dx.doi.org/10.1094/MPMI.2004.17.11.1192] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15553245]

43. Zhou, H.; Liu, B.; Weeks, D.P.; Spalding, M.H.; Yang, B. Large chromosomal deletions and heritable small genetic changes induced by CRISPR/Cas9 in rice. Nucleic Acids Res.; 2014; 42, pp. 10903-10914. [DOI: https://dx.doi.org/10.1093/nar/gku806] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25200087]

44. Cernadas, R.A.; Doyle, E.L.; Niño-Liu, D.O.; Wilkins, K.E.; Bancroft, T.; Wang, L.; Schmidt, C.L.; Caldo, R.; Yang, B.; White, F.F. et al. Code-assisted discovery of TAL effector targets in bacterial leaf streak of rice reveals contrast with bacterial blight and a novel susceptibility gene. PLoS Pathog.; 2014; 10, e1003972. [DOI: https://dx.doi.org/10.1371/journal.ppat.1003972]

45. Hu, Y.; Zhang, J.; Jia, H.; Sosso, D.; Li, T.; Frommer, W.B.; Yang, B.; White, F.F.; Wang, N.; Jones, J.B. Lateral organ boundaries 1 is a disease susceptibility gene for citrus bacterial canker disease. Proc. Natl. Acad. Sci. USA; 2014; 111, pp. E521-E529. [DOI: https://dx.doi.org/10.1073/pnas.1313271111]

46. Schwartz, A.R.; Morbitzer, R.; Lahaye, T.; Staskawicz, B.J. TALE-induced bHLH transcription factors that activate a pectate lyase contribute to water soaking in bacterial spot of tomato. Proc. Natl. Acad. Sci. USA; 2017; 114, pp. E897-E903. [DOI: https://dx.doi.org/10.1073/pnas.1620407114]

47. Toledo-Ortiz, G.; Huq, E.; Quail, P.H. The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell; 2003; 15, pp. 1749-1770. [DOI: https://dx.doi.org/10.1105/tpc.013839]

48. Wilkins, K.E.; Booher, N.J.; Wang, L.; Bogdanove, A.J. TAL effectors and activation of predicted host targets distinguish Asian from African strains of the rice pathogen Xanthomonas oryzae pv. oryzicola while strict conservation suggests universal importance of five TAL effectors. Front. Plant Sci.; 2015; 6, 536. [DOI: https://dx.doi.org/10.3389/fpls.2015.00536]

49. Yu, Y.; Zhang, Y.; Chen, X.; Chen, Y. Plant Noncoding RNAs: Hidden Players in Development and Stress Responses. Annu. Rev. Cell Dev. Biol.; 2019; 35, pp. 407-431. [DOI: https://dx.doi.org/10.1146/annurev-cellbio-100818-125218]

50. Mücke, S.; Reschke, M.; Erkes, A.; Schwietzer, C.-A.; Becker, S.; Streubel, J.; Morgan, R.D.; Wilson, G.G.; Grau, J.; Boch, J. Transcriptional Reprogramming of Rice Cells by Xanthomonas oryzae TALEs. Front. Plant Sci.; 2019; 10, 162. [DOI: https://dx.doi.org/10.3389/fpls.2019.00162]

51. Arnold, P.R.; Wells, A.D.; Li, X.C. Diversity and Emerging Roles of Enhancer RNA in Regulation of Gene Expression and Cell Fate. Front. Cell Dev. Biol.; 2020; 7, 377. [DOI: https://dx.doi.org/10.3389/fcell.2019.00377] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31993419]

52. Chen, L.; Zhu, Q.H.; Kaufmann, K. Long non-coding RNAs in plants: Emerging modulators of gene activity in development and stress responses. Planta; 2020; 252, 92. [DOI: https://dx.doi.org/10.1007/s00425-020-03480-5]

53. Banerjee, J.; Sahoo, D.K.; Dey, N.; Houtz, R.L.; Maiti, I.B. An Intergenic Region Shared by At4g35985, and At4g35987, in Arabidopsis thaliana Is a Tissue Specific and Stress Inducible Bidirectional Promoter Analyzed in Transgenic Arabidopsis and Tobacco Plants. PLoS ONE; 2013; 8, e79622. [DOI: https://dx.doi.org/10.1371/journal.pone.0079622] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24260266]

54. Kourmpetli, S.; Lee, K.; Hemsley, R.; Rossignol, P.; Papageorgiou, T.; Drea, S. Bidirectional promoters in seed development and related hormone/stress responses. BMC Plant Biol.; 2013; 13, 187. [DOI: https://dx.doi.org/10.1186/1471-2229-13-187] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24261334]

55. Liu, X.; Yang, W.; Li, Y.; Li, S.; Zhou, X.; Zhao, Q.; Fan, Y.; Lin, M.; Chen, R. The intergenic region of the maize defensin-like protein genes Def1 and Def2 functions as an embryo-specific asymmetric bidirectional promoter. J. Exp. Bot.; 2016; 67, pp. 4403-4413. [DOI: https://dx.doi.org/10.1093/jxb/erw226] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27279278]

56. Wang, Q.; Wan, L.; Li, D.; Zhu, L.; Qian, M.; Deng, M. Searching for bidirectional promoters in Arabidopsis thaliana. BMC Bioinform.; 2009; 10, S29. [DOI: https://dx.doi.org/10.1186/1471-2105-10-S1-S29] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19208129]

57. Thieffry, A.; Vigh, M.L.; Bornholdt, J.; Ivanov, M.; Brodersen, P.; Sandelin, A. Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways. Plant Cell; 2020; 32, pp. 1845-1867. [DOI: https://dx.doi.org/10.1105/tpc.19.00815]

58. Yuan, M.; Ke, Y.; Huang, R.; Ma, L.; Yang, Z.; Chu, Z.; Xiao, J.; Li, X.; Wang, S. A host basal transcription factor is a key component for infection of rice by TALE-carrying bacteria. eLife; 2016; 5, e19605. [DOI: https://dx.doi.org/10.7554/eLife.19605]

59. Heigwer, F.; Kerr, G.; Boutros, M. E-CRISP: Fast CRISPR target site identification. Nat. Methods; 2014; 11, pp. 122-123. [DOI: https://dx.doi.org/10.1038/nmeth.2812]

60. Miao, J.; Guo, D.; Zhang, J.; Huang, Q.; Qin, G.; Zhang, X.; Wan, J.; Gu, H.; Qu, L.-J. Targeted mutagenesis in rice using CRISPR-Cas system. Cell Res.; 2013; 23, pp. 1233-1236. [DOI: https://dx.doi.org/10.1038/cr.2013.123]

61. Fayos, I.; Meunier, A.C.; Vernet, A.; Navarro-Sanz, S.; Portefaix, M.; Lartaud, M.; Bastianelli, G.; Périn, C.; Nicolas, A.; Guiderdoni, E. Assessment of the roles of SPO11-2 and SPO11-4 in meiosis in rice using CRISPR/Cas9 mutagenesis. J. Exp. Bot.; 2020; 71, pp. 7046-7058. [DOI: https://dx.doi.org/10.1093/jxb/eraa391] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32842152]

62. Sallaud, C.; Meynard, D.; van Boxtel, J.; Gay, C.; Bès, M.; Brizard, J.P.; Larmande, P.; Ortega, D.; Raynal, M.; Portefaix, M. et al. Highly efficient production and characterization of T-DNA plants for rice (Oryza sativa L.) functional genomics. Theor. Appl. Genet.; 2003; 106, pp. 1396-1408. [DOI: https://dx.doi.org/10.1007/s00122-002-1184-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12677401]

63. Romero, L.E.; Lozano, I.; Garavito, A.; Carabali, S.J.; Triana, M.; Villareal, N.; Reyes, L.; Duque, M.C.; Martinez, C.P.; Calvert, L. et al. Major QTLs control resistance to Rice Hoja Blanca Virus and its vector Tagosodes orizicolus. G3 Genes Genomes Genet.; 2014; 4, pp. 133-142.

64. Liu, W.; Xie, X.; Ma, X.; Li, J.; Chen, J.; Liu, Y.G. DSDecode: A Web-Based Tool for Decoding of Sequencing Chromatograms for Genotyping of Targeted Mutations. Mol. Plant; 2015; 8, pp. 1431-1433. [DOI: https://dx.doi.org/10.1016/j.molp.2015.05.009]

65. Schubert, M.; Lindgreen, S.; Orlando, L. AdapterRemoval v2: Rapid adapter trimming, identification, and read merging. BMC Res. Notes; 2016; 9, 88. [DOI: https://dx.doi.org/10.1186/s13104-016-1900-2]

66. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc.; 2016; 11, pp. 1650-1667. [DOI: https://dx.doi.org/10.1038/nprot.2016.095]

67. Lawrence, M.; Huber, W.; Pagès, H.; Aboyoun, P.; Carlson, M.; Gentleman, R.; Morgan, M.; Carey, V.J. Software for Computing and Annotating Genomic Ranges. PLoS Comput. Biol.; 2013; 9, e1003118. [DOI: https://dx.doi.org/10.1371/journal.pcbi.1003118]

68. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.; 2014; 15, 550. [DOI: https://dx.doi.org/10.1186/s13059-014-0550-8]