# Mitochondrial Transcriptome Control and Intercompartment Cross-Talk During Plant Development

Adnan Khan Niazi, E. Delannoy, Rana Khalid Iqbal, Daria Mileshina, R. Val, M. Gabryelska, E. Wyszko, L. Soubigou-Taconnat, M. Szymanski, J. Barciszewski, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Adnan Khan Niazi, E. Delannoy, Rana Khalid Iqbal, Daria Mileshina, R. Val, et al., Mitochondrial Transcriptome Control and Intercompartment Cross-Talk During Plant Development. Cells, 2019, 8 (6), pp.583. 10.3390/cells8060583. [hal-02184302](https://hal.science/hal-02184302)

HAL Id: hal-02184302

https://hal.science/hal-02184302

Submitted on 15 Jul 2019

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Mitochondrial Transcriptome Control and Intercompartment Cross-Talk During Plant Development

Adnan Khan Niazi <fnref n="1" />,<fnref n="2" />,† ID, Etienne Delannoy <fnref n="3" />,†, Rana Khalid Iqbal <fnref n="1" />,†, Daria Mileshina <fnref n="1" />, Romain Val <fnref n="1" />, Marta Gabryelska <fnref n="4" />, Eliza Wyszko <fnref n="4" /> ID, Ludivine Soubigou-Taconnat <fnref n="3" />, Maciej Szymanski <fnref n="5" />, Jan Barciszewski <fnref n="4" />,<fnref n="6" /> ID, Frédérique Weber-Lotfi <fnref n="1" />, José Manuel Gualberto <fnref n="1" /> and André Dietrich <fnref n="1" />,*

Received: 21 March 2019; Accepted: 11 June 2019; Published: 13 June 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** We address here organellar genetic regulation and intercompartment genome coordination. We developed earlier a strategy relying on a tRNA-like shuttle to mediate import of nuclear transgene-encoded custom RNAs into mitochondria in plants. In the present work, we used this strategy to drive *trans*-cleaving hammerhead ribozymes into the organelles, to knock down specific mitochondrial RNAs and analyze the regulatory impact. In a similar approach, the tRNA mimic was used to import into mitochondria in *Arabidopsis thaliana* the *orf77*, an RNA associated with cytoplasmic male sterility in maize and possessing sequence identities with the *atp9* mitochondrial RNA. In both cases, inducible expression of the transgenes allowed to characterise early regulation and signaling responses triggered by these respective manipulations of the organellar transcriptome. The results imply that the mitochondrial transcriptome is tightly controlled by a "buffering" mechanism at the early and intermediate stages of plant development, a control that is released at later stages. On the other hand, high throughput analyses showed that knocking down a specific mitochondrial mRNA triggered a retrograde signaling and an anterograde nuclear transcriptome response involving a series of transcription factor genes and small RNAs. Our results strongly support transcriptome coordination mechanisms within the organelles and between the organelles and the nucleus.

**Keywords:** anterograde regulation; cytoplasmic male sterility (CMS); plant mitochondria; retrograde regulation; ribozyme; RNA trafficking; signaling

<footnote n="1">¹ Institute of Plant Molecular Biology (IBMP), CNRS and University of Strasbourg, 12 rue du Général Zimmer, 67084 Strasbourg, France; adnan1753@yahoo.com (A.K.N.); rkiqbal@etu.unistra.fr (R.K.I.); d.mileshina@gmail.com (D.M.); r.val@arvalis.fr (R.V.); lotfif@unistra.fr (F.W.-L.); jose.gualberto@ibmp-cnrs.unistra.fr (J.M.G.)</footnote>

<footnote n="2">² Centre of Agricultural Biochemistry and Biotechnology (CABB), University of Agriculture, Faisalabad 38000, Pakistan</footnote>

<footnote n="3">³ Institute of Plant Sciences Paris-Saclay IPS2, CNRS, INRA, Université Paris-Sud, Université Evry, Université Paris-Saclay, Paris Diderot, Sorbonne Paris-Cité, 91405 Orsay, France; etienne.delannoy@inra.fr (E.D.); ludivine.soubigou-taconnat@inra.fr (L.S.-T.)</footnote>

<footnote n="4">⁴ Institute of Bioorganic Chemistry, Polish Academy of Sciences, Ul. Z. Noskowskiego 12/14, 61-704 Poznan, Poland; marta.gabryelska@gmail.com (M.G.); eliza.wyszko@ibch.poznan.pl (E.W.); Jan.Barciszewski@ibch.poznan.pl (J.B.)</footnote>

<footnote n="5">⁵ Department of Computational Biology, Institute of Molecular Biology and Biotechnology, A. Mickiewicz University Poznan, Ul. Umultowska 89, 61-614 Poznan, Poland; mszyman@amu.edu.pl</footnote>

<footnote n="6">⁶ NanoBioMedical Centre of the Adam Mickiewicz University, Umultowska 85, 61614 Poznan, Poland</footnote>

<footnote>* Correspondence: diandid.5926@orange.fr</footnote>

<footnote>† These authors contributed equally to this project.</footnote>
# 1. Introduction

Mitochondria massively import proteins encoded by the nuclear genome, but still contain their own genetic system, distinct from that of the nucleo-cytosolic compartment. Mammalian mitochondria possess an extremely compact 16.5 kb circular genome with a single non-coding region that carries the promoters for transcription in both orientations [1]. In plants, mitochondrial genomes are large (usually 200–700 kb) and have a low gene density. Nevertheless, due to the occurrence of multiple promoters spread all along, they are still almost entirely transcribed, including the large non-coding regions [2]. The mitochondrial DNA (mtDNA) encodes essential subunits of the different oxidative phosphorylation (OXPHOS) complexes and one can thus presume that the mitochondrial transcriptome needs to be tightly kept under control. Studies on the regulation of mammalian mitochondrial gene expression so far highlighted a major role of post-transcriptional processes, including coordination of mitochondrial and cytosolic translation [3–6]. Translation regulation and efficiency, rather than mtDNA copy number or transcription rate, were proposed to be crucial for the synthesis of the OXPHOS subunits [7]. Methylation of the mtDNA and epigenetic-like modifications of the mitochondrial transcription factor A (TFAM) were suggested to participate in establishing the mitochondrial gene expression profiles [8]. Known nuclear gene expression regulators turned out to be also involved in mtDNA transcription [9]. Conversely, models for mtDNA transcriptional regulation by regular nuclear-encoded mitochondrial factors were proposed [1]. In plant mitochondria as well, the major genetic control points are generally considered to be post-transcriptional, although there is some transcriptional control of gene expression [4,10]. A number of processes driven by nuclear-encoded protein factors, such as 5'- and 3'-end maturation of transcripts, intron splicing and RNA editing are all potential sites of regulation [11].

As a further level of complexity, mitochondrial biogenesis and functioning involve both nuclear gene expression and mitochondrial gene expression, especially to build the protein complexes of the electron transport chain (ETC). This is likely to require a coordinated communication between the organelles and the nucleus [12]. Intercompartment coordination mechanisms include both anterograde (nucleus to organelle) and retrograde (organelle to nucleus) signals [13]. Anterograde mechanisms control organelle gene expression in response to endogenous or exogenous signals that are perceived by the nucleus. Indeed, expression of the mitochondrial genome, from DNA replication and transcription to transcript processing, editing and translation, is entirely carried out by proteins encoded by nuclear-located genes [3,14]. Retrograde regulation refers to signals sent by the organelles to communicate their functional and developmental state to the nucleus, which can then modulate anterograde control and cellular metabolism accordingly [11,15]. Mitochondrial retrograde regulation (MRR) is an important mechanism of communication between mitochondria and the nucleus and is conserved among yeasts, mammals and plants. However, the signaling molecules and the signal transduction mechanisms are diverse [15].

The present work aims to further explore these pathways with a special interest in mitochondrial genetic expression control and signaling during development. Plants offer in this respect a relevant and easily accessible working model. Establishment of the plant mitochondrial transcriptome has been documented in detail along seed germination in *A. thaliana*, but little is known about the spatial and temporal control of the organellar transcripts during post-germination development [16–18]. Our experiments were based on the analysis and manipulation of the organellar transcriptome in *Arabidopsis thaliana* through a novel approach set up in our laboratory and based on the existence of a natural process of transfer RNA (tRNA) import from the cytosol into mitochondria in plant cells (reviewed in [19]). Circumventing the lack of methodologies for mitochondrial transformation, we showed previously that a tRNA mimic can be used in vivo as a shuttle for importing into plant mitochondria cargo RNAs expressed from nuclear transgenes [20]. Taking *trans*-cleaving hammerhead ribozymes as cargo sequences allowed to specifically knock down mitochondrial RNAs and demonstrate that the mt-DNA-encoded MATR protein is a mitochondrial maturase involved in the splicing of multiple organellar introns [20,21]. We now used tRNA mimic-mediated
shuttling of cargo RNAs to affect the mitochondrial transcriptome in *A. thaliana* plants at different developmental stages, so as to investigate the resulting regulation processes. Two strategies were applied. First, individual mitochondrial mRNAs encoding subunits of OXPHOS complexes were chosen as targets for specific *trans*-cleaving hammerhead ribozymes driven into the organelles by the tRNA mimic. In a second strategy, the tRNA mimic was used to drive into mitochondria a sequence associated with cytoplasmic male sterility (CMS). The results altogether imply that the mitochondrial transcriptome is tightly controlled by a buffering mechanism at early and intermediate stages of development. On the other hand, knocking down a major mitochondrial mRNA triggered a retrograde signaling and an anterograde response from the nucleus. Whereas it is currently considered that mitochondrial regulation processes occur essentially at the post-transcriptional stage (see above), our results support mRNA coordination mechanisms within the organelles and between the organelles and the nucleus.

## 2. Materials and Methods

### 2.1. Preparation of Gene Constructs for In Vivo Expression

For the catalytic RNA-mediated knockdown strategy, the complete sequences encoding the *trans*-ribozymes, the linker L, the *Turnip yellow mosaic virus* (TYMV) PKTLS import shuttle and the *Hepatitis delta virus* (HDV) antigenomic *cis*-cleaving ribozyme (cHDV) were assembled by polymerase chain reaction (PCR) amplification from the pCK-PSTYPKTLScHDV plasmid [20]. Direct megaprimers contained a 5' *Xho*I site followed by the desired ribozyme and linker sequences upstream of a sequence corresponding to the 5' region of the PKTLS. The reverse primer contained a 5' *Spe*I site followed by a sequence annealing to the 3' region of the cHDV [22]. The obtained PCR products were cloned into the pGEM-T vector (Promega, Madison, WI, USA) for sequencing, re-excized with the *Xho*I and *Spe*I restriction enzymes and cloned into the *Xho*I and *Spe*I sites of the inducible transcription unit (O~LexA~-46 estradiol-inducible promoter) of the pER8 vector [23], raising a series of pER8-Rzxxxx-L-PKTLS-cHDV plasmids. These were used to transform *A. thaliana* or *Nicotiana tabacum*. Sequence details can be found in Figure 1. As previously [20,21], the linker (L) to provide spacing between the *trans*-ribozyme and PKTLS motifs was selected from a pool of random sequences, based on MFOLD [24] predictions of weak RNA secondary structures unlikely to interfere either with PKTLS formation or *trans*-ribozyme binding to the target sequence.

For the CMS RNA strategy, a PCR product comprising the complete *orf*77 sequence with four nucleotides of 5'-UTR and 92 nucleotides of 3'-UTR was amplified using total DNA from CMS-S maize as a template. The forward primer comprised a *Hind*III site and a sequence encoding a 6xHis tag. The reverse primer had a *Bam*HI site. The resulting product was cloned into the *Hind*III and *Bam*HI sites upstream of the sequences encoding the TYMV PKTLS and the HDV *cis*-ribozyme in an available derivative of the pUCAP plasmid [25] carrying these sequences. After sequencing of the recombinant plasmids, the assembled *orf*77-PKTLS-HDV construct was re-excized upon digestion with *Asc*I and *Pac*I and cloned into the *Asc*I and *Pac*I restriction sites of the estradiol-inducible transcription unit of the pER8 vector [23]. The resulting plasmid called pER8-*orf*77-PKTLS-cHDV was used for transformation of *A. thaliana*. For control, plants were transformed with empty pER8 plasmid. Sequence details can be found in Figure 2.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Structure of the ribozyme-PKTLS chimeric RNAs targeted to mitochondria in the present study. (a) Structure of the *Rznad9*-L-PKTLS RNA. The *trans*-cleaving ribozyme *Rznad9* directed against the *A. thaliana* mitochondrial *nad9* mRNA (*Atnad9*) is attached to the 5'-end of the PKTLS shuttle via a 40 nucleotide linker (L) selected from a pool of random sequences. The *Rznad9* hammerhead sequence is annealed to its target sequence motif in the *Atnad9* mRNA. The *Rzsdh3* (b), *Rzcox3* (c) and *Rzatp9* (d) ribozymes, directed against the *N. tabacum* *sdh3* (*Ntsdh3*) mRNA, the *A. thaliana* *cox3* (*Atcox3*) mRNA and the *A. thaliana* *atp9* (*Atatp9*) mRNA are attached to the same L-PKTLS moiety to generate the *Rzsdh3*-L-PKTLS, *Rzcox3*-L-PKTLS and *Rzatp9*-L-PKTLS, respectively. Ribozyme cleavage sites are indicated in the target sequence motifs by thick arrows. Their precise location is after position 421 in the *A. thaliana* *nad9* coding sequence (24662-25234 in accession JF729201), after position 216 in the *N. tabacum* *sdh3* coding sequence (77198-77524 in accession BA000042, complementary strand), after position 685 in the *A. thaliana* *cox3* coding sequence (328926-329723 in accession JF729201, complementary strand), and after position 99 in the *A. thaliana* *atp9* coding sequence (269920-270177 in accession JF729201, complementary strand). The *Rzsdh3* and *Rzatp9* ribozymes are as described earlier in Sultan et al. [21] and Val et al. [20], respectively.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Structure of the *orf77-PKTLS* RNA targeted to mitochondria in the present study. The *Zea mays* CMS-S-specific *orf77* coding sequence (initiation codon and termination codon in bold) with its four nucleotide upstream and 92 nucleotide downstream sequences (italics) is directly attached to the 5'-end of the PKTLS shuttle. Short additional sequences deriving from the *Hind*III and *Bam*HI restriction sites introduced for cloning purposes are underlined. The chimeric RNA was expressed from a nuclear transgene and driven into the organelles by the PKTLS moiety.</figcaption></figure>

## 2.2. Nuclear Transformations

*A. thaliana* plants (ecotype Col-0) were transformed via *Agrobacterium tumefaciens* through floral dip [26]. Seeds from transformed plants were recovered and germinated on hygromycin medium under long day conditions (16 h light/8 h dark) for transformant selection. The presence and integrity of the constructs was confirmed by PCR and homozygous lines of the relevant transformants were selected.

To generate tobacco transformants, leaf discs from *N. tabacum* plants grown in vitro on germination medium (4.3 g/L DUCHEFA micro- and macroelements M0238, 10 g/L agar, 1% *w/v* sucrose, pH 5.7) were exposed to the relevant recombinant *A. tumefaciens* suspension and subsequently cultured in vitro on coculture medium (4.3 g/L DUCHEFA micro- and macroelements M0238, 825 mg/L NH₄NO₃, 2 mg/mL glycine, 100 mg/L myo-inositol, 0.5 mg/L nicotinic acid, 0.5 mg/L putrescine, 0.1 mg/L thiamine, 3% *w/v* sucrose, 0.8% *w/v* bacto-agar, 2 mg/mL benzylaminopurine (BAP), 0.05 mg/L naphthaleneacetic acid (NAA), pH 5.8) for 3-4 days at 28 °C (12 h light/12 h dark). The leaf discs were then transferred to selection medium (4.3 g/L DUCHEFA micro- and macroelements M0238, 825 mg/L NH₄NO₃, 2 mg/mL glycine, 100 mg/L myo-inositol, 0.5 mg/L nicotinic acid, 0.5 mg/L putrescine, 0.1 mg/L thiamine, 3% *w/v* sucrose, 0.8% *w/v* bacto-agar, 350 mg/L cefotaxime, 20 mg/L hygromycin, pH 5.8) at 25 °C (12 h light/12 h dark). The medium was changed every two to three weeks. Calli appearing on the edges were transferred to individual pots with selection medium when young leaves were formed, and planted into mineral-rich soil when roots had emerged. After acclimatization in growth chambers (25 °C, 12 h light/12 h dark) the plantlets were grown in the greenhouse. Recovered seeds were further selected through germination on hygromycin medium.

## 2.3. Seedling Growth and Transgene Induction

*A. thaliana* or *N. tabacum* transgenic seedlings were grown under long day light conditions (16h light/8h dark) or in the dark on solid agar medium. For transgene induction, the agar layer was overlaid with liquid medium containing 10 µM β-estradiol [23] Alternatively, plants grown on solid
agar in the light were transferred to culture plates containing liquid medium supplemented with 10 µM β-estradiol. Only the roots were dipping into the medium. Developmental stages considered are defined in the “Results” section. At all stages and in all growth conditions, samples of 3 to 4 entire plants were harvested each day till day 4 after induction. To avoid variations that would be due to the circadian clock, plants were always induced and collected at the same time of the day. Whole plant samples were immediately frozen in liquid nitrogen. Control plants (see the “Results” section) were grown and treated in the same way to take into account any side effects.

## 2.4. RNA Extraction, Northern Bloting and RT-qPCR Analyses

Total RNA was extracted from *A. thaliana* or *N. tabacum* seedlings following standard TRI Reagent protocols (Molecular Research Center). Northern blot analyses were carried out through standard protocols using non-radioactive digoxygenin-labeled probes prepared with the PCR DIG Probe Synthesis kit (Roche, Basel, Switzerland). Prior to reverse transcription, RNA samples were treated up to 3 times with RNase-free DNase 1 (Thermo Scientific, Waltham, MA, USA). Reverse transcription was carried out with SuperScript III (Invitrogen, Carlsbad, CA, USA), RevertAid (Fermentas, Waltham, MA, USA) or GoScript (Promega, Madison, WI, USA) reverse transcriptase according to the corresponding recommended protocols, using 2 pmol of a specific primer or 250 ng of a random hexanucleotide mixture. Reactions were stopped by 15 min incubation at 70 °C and the RT samples served directly for standard or real-time PCR. The absence of residual DNA contamination in the original RNA samples was assessed by the absence of PCR product amplification from cDNAs resulting from reverse transcription assays run without the enzyme.

Real time PCR (qPCR) was performed in 384-well optical plates on an iCycler (BioRad, Hercules, CA, USA) or a LightCycler 480 II (Roche) using 5 µL of PCR master mix containing 480 SYBER Green I fluorescent reporter (Roche) with 2.5 µM forward and reverse specific primers. DNA or cDNA template and water were added to a total volume of 10 µL. Each sample was performed in triplicate. Reactions were run in the following conditions: pre-heating at 95 °C for 10 min, followed by 40 cycles of 15 sec at 95 °C, 30 sec at 60 °C and 15 sec at 72 °C. Melting curves were run to assess the specificity of the products.

The nuclear genes encoding actin 2 (accession AT3G18780; AB158612), glyceraldehyde 3-phosphate dehydrogenase c (GAPDH, accession AT1G13440; AY049259) and RGS1-HXK1 interacting protein 1 (RHIP1, accession AT4G26410; BT002964), as well as the mitochondrial gene encoding mitochondrial ribosomal protein L2 (rpl2, accession BA000042, region 361051-363948) served as reference genes. To design qPCR primers, we used the Universal ProbeLibrary and the ProbeFinder software (Roche) (https://lifescience.roche.com/en_fr/articles/Universal-ProbeLibrary-System-Assay-Design.html#ProbeFinder). For short sequences that were unsuccessful with such an approach, we used the Primer 3 software [27].

RT-qPCR data obtained with RNAs from transformants expressing ribozyme-PKTLS transcripts were analyzed with the Student's t-test. Statistical analyses of the mitochondrial gene expression RT-qPCR data obtained with RNAs from the *orf77-PKTLS* expressing transformants were carried out using the GraphPad Prism version 7.01 software (www.graphpad.com/). Correlation analyses, one way ANOVA, multiple comparison tests and Tukey tests were used to calculate the *p* values. *p* values < 0.05 were considered statistically significant.

## 2.5. Isolation of Mitochondria and Import Assessment

To confirm that the *orf77-PKTLS* RNA expressed in the transformants was targeted to the organelles, mitochondria were isolated according to established protocols [28] from light-grown plants at day 2 after induction of transgene expression with β-estradiol. Isolated organelles were extensively treated with RNase in buffer containing 100 µg/mL RNase A and 750 U/mL RNase T1. Mitochondrial RNA was extracted following standard TRI Reagent protocols (Molecular Research Center), with 150 µL of Trizol reagent (Invitrogen) added to an amount of mitochondria equivalent to 200 µg of proteins.
Total RNA was prepared in parallel from the same plants. Mitochondrial and total RNAs were probed by RT-qPCR for *orf77-PKTLS*, for specific nuclear transcripts and for selected mitochondrial transcripts, so as to characterize mitochondrial enrichment and confirm the absence of significant contamination.

## 2.6. Microarray Assays

Microarray analyses were carried out on the INRA transcriptomic platform, first at the Unité de Recherche en Génomique Végétale (URGV, Evry, France) and subsequently at the Institute of Plant Sciences Paris-Saclay (IPS2, Orsay, France), using CATMAv6.2 (Roche-NimbleGen technology, Madison, WI, USA) and CATMAv7 (AGILENT technology, Santa Clara, CA, USA) arrays. High density CATMAv6.2 microarray slides contained per chamber 270,000 primers representing all *A. thaliana* genes, i.e., 30,834 probes referring to the TAIRv8 annotation (including 476 probes of mitochondrial and chloroplast genes) (www.arabidopsis.org/), plus 1289 probes corresponding to EUGENE software predictions, 5352 probes corresponding to repeated elements, 658 probes for miRNAs, 342 probes for other non-coding RNAs (rRNAs, tRNAs, snRNAs, soRNAs) and finally 36 control probes. CATMAv7 slides contained per chamber 180,000 primers corresponding to 38,360 probes referring to the TAIRv8 annotation (including primers for mitochondrial and chloroplast genes), plus primers corresponding to EUGENE software specific predictions, primers corresponding to repeats, primers for miRNAs and other ncRNAs (rRNAs, tRNAs, snRNAs, snoRNAs) and the 36 controls. Each of these primers was designed in both orientations (forward strand in triplicate and reverse strand without duplication).

RNA samples from plants collected each day from Day 0 to Day 4 post-induction (i.e., 5 test samples and 5 control samples) were submitted to RT-qPCR analysis to assess *nad9* knockdown and individually used for microarray assays. Two to three independent biological replicates were analyzed and fluorochrome reversal (dye swap) was systematically applied. Total RNA extracted from *A. thaliana* plants following TRI-reagent protocols (Molecular Research Center) and treated twice with RNase-free DNase I (Fermentas) was finally purified using NucleoSpin RNA kits (Macherey-Nagel). Further steps were mostly adapted from Lurin et al. [29]. RNA samples were checked for quality on a bioanalyser (Agilent), quantified with a Quant-iT RiboGreen RNA Assay Kit (Thermo Fisher Scientific) and amplified with the Complete Whole Transcriptome Amplification kit (WTA2, Sigma-Aldrich) using random hexamer primers. After purification, 30 pmoles of cDNA per sample and per slide were hybridized overnight at 42 °C in the presence of formamide. Two micron scanning was performed with an InnoScan900 scanner (InnopsysR, Carbonne, France) and raw data were extracted using the MapixR software (InnopsysR, Carbonne, France).

## 2.7. Microarray Data Analyses

Microarray raw data comprised the logarithm of median feature pixel intensity at wavelengths 635 nm (red) and 532 nm (green). To correct the dye bias, global intensity-dependent normalization was run using the LOESS procedure [30]. Differential analysis was based on log-ratio averaging over the duplicate probes and over the technical replicates. The data were submitted to *limma* moderated t-statistics [31,32] using the R software package (www.r-project.org). No evidence that the specific variances would vary between probes was raised by *limma* and consequently the moderated t-statistics was assumed to follow a standard normal distribution. The squeezeVar function of the *limma* library was used to smooth the specific variances by computing empirical Bayes posterior means. To control the false discovery rate, adjusted *p*-values were calculated using the optimized FDR (false discovery rate) approach and the kerfdr library [33,34]. Probes with an adjusted *p*-value ≤ 0.05 were considered as differentially expressed. Complementary analyses were developed with the MapMan software (www.mapman.gabipd.org). To link the levels of the *nad9* transcript to nuclear transcripts, a sparse PLS (Partial Least Squares) regression approach was performed using the spls (v2.2) R package with the default parameters, except *eta* = 0.5, *K* = 1, *kappa* = 0.5 [35]. The 99% confidence interval of the coefficient of each gene was calculated by bootstrapping (N = 1000). Only genes for
which the confidence interval did not include 0 were considered as associated with *nad9* expression. As we have directly manipulated the *nad9* transcript levels, this association can be interpreted as a regulation of these genes by the *nad9* transcript levels. The sparse PLS analyses were performed separately for the etiolated and the light-grown seedlings using the log2 fold-changes. Microarray data from this manuscript were deposited into the CATdb database [36,37] (tools.ips2.u-psud.fr/CATdb/, Projects: RS13-02_Mitomanip and RS14-01_Mitomanip2) and into the Gene Expression Omnibus (GEO) repository [38] at the National Center for Biotechnology Information (NCBI) (accession numbers GSE127756 and GSE93122), according to the "Minimum Information About a Microarray Experiment" standards [39].

### 3. Results

#### 3.1. Trans-Ribozyme-Mediated Modulation of Mitochondrial RNA Steady State Levels Depends on the Plant Developmental Stage

Four mitochondrial mRNAs, *nad9*, *sdh3*, *cox3* and *atp9*, were initially chosen as targets for specific *trans*-cleaving hammerhead ribozymes. The catalytic RNAs, *Rznad9*, *Rzsdh3*, *Rzcox3* and *Rzatp9* were designed on the same basis as in previous studies [20,21], i.e., with only 2 base pairs in helix II and a UUUU tetraloop. The *Rzatp9* ribozyme was already described before [20]. The ribozyme sequences were attached as 5'-trailors to the previously set up PKTLS mitochondrial RNA shuttle [20] through a weakly structured linker (L). The PKTLS sequence corresponded to the last 120 nucleotides of the *Turnip yellow mosaic virus* (TYMV) genomic RNA and included the tRNA-like structure and an upstream pseudoknot considered to improve the interaction with the aminoacyl-tRNA synthetase [40], which is a prerequisite for mitochondrial import [41,42]. Figure 1 shows the resulting design. The specifically designed *trans*-ribozymes associated with the linker and the PKTLS shuttle were expressed from nuclear transgenes in stably transformed *A. thaliana* plants (for cleavage of the *nad9*, *cox3* and *atp9* mitochondrial targets) and *N. tabacum* plants (for cleavage of the *sdh3* mitochondrial target). To master the expression, the sequences were placed under the control of the commonly used estradiol-inducible promoter system of the pER8 vector [23]. In this system, a chimeric transcription activator (XVE), taking together the DNA-binding domain of the bacterial repressor LexA (X), the acidic transactivating domain of the *Herpes simplex virus* VP16 factor (V) and the regulatory region of the human estrogen receptor (E), is expressed under the G10-90 strong constitutive promoter. The *trans*-activating activity of the chimeric XVE factor is triggered upon estradiol binding. The XVE-estradiol complex activates the target promoter, which consists of eight copies of the LexA operator fused upstream of a minimal 35S promoter. The activated target promoter drives the expression of the transgene, which is thus strictly dependent on estradiol delivery to the cells [23]. The XVE factor is constitutively expressed, with possible side effects. In our hands, estradiol delivery to the plants showed by itself a slight activating effect on mitochondrial genes. Also, our ribozyme-PKTLS constructs include a tRNA-like moiety that is recognized by the cognate aminoacyl-tRNA synthetase and enters tRNA pathways. Untransformed plants not exposed to estradiol were thus not appropriate controls to take into account all possible side effects and we took advantage of the fact that the SDH3 subunit is not encoded by a mitochondrial gene in *A. thaliana*, so that the *Rzsdh3* chimeric ribozyme directed against the *N. tabacum* mitochondrial *sdh3* mRNA has no RNA target in *A. thaliana* mitochondria. Transformant *A. thaliana* expressing the *Rzsdh3*-L-PKTLS chimeric ribozyme upon estradiol treatment was thus used as a control, as in previous studies [21]. Expression of the ribozyme transgenes was induced with estradiol at different developmental stages of light-grown plants and in dark-grown seedlings. Test lines and the *Rzsdh3* "no target" control line were grown, treated and harvested in parallel in strictly identical conditions. Expected knockdown of the corresponding target RNAs in test lines was analyzed every day for up to four days after onset of estradiol induction. For both the test lines and the control line, the level of the target RNA of interest at Day 1, Day 2, Day 3 and Day 4 post-induction was compared to the level at Day 0. At each growth stage, the kinetics obtained with the test lines were subsequently compared day by day to the kinetics of the *Rzsdh3* control line. Initial comparison to
Day 0 and day-by-day normalization against the control at the same growth stage always gave the same profiles. *Rznad9* and *Rzatp9* were the most active in vivo and were exploited at all growth stages. In all our assays, expression of the ribozymes had no visually detectable phenotypic effects within the four-day time frame of the experiments (See for example, Figure S1).

Induced expression of the previously described *Rzatp9*-L-PKTLS RNA in *A. thaliana* seedlings at an early stage of growth (maximum 4 true leaves) was efficient (Figure 3a), but failed to trigger a significant decrease in the steady-state level of the *atp9* target mRNA within four days in a day by day analysis (Figure 3b). A more detailed kinetic analysis along 16 h of initial phase of ribozyme expression led to the same absence of significant decrease of the target RNA (Figure S3). Similarly, expression of the *Rznad9*-L-PKTLS RNA at the early stage of growth was unable to knock down the level of *nad9* mRNA in a 4-day (Figure 3c) or 16 h (Figure S3) analysis. Thus, at an early stage of development, *A. thaliana* seedlings seemed to be able to maintain the steady-state level of mRNAs like *atp9* and *nad9* in the presence of the corresponding *trans*-cleaving ribozymes. This suggests that mitochondrial function is of primary importance at that stage of growth and that a still active transcription might efficiently "buffer" transcriptome changes through replacement of cleaved RNAs.

At an intermediate stage of *A. thaliana* development (up to 10 true leaves), expression (Figure 3d) and mitochondrial import of the *Rznad9*-L-PKTLS RNA resulted in a striking "bounce-back" profile of the *nad9* transcript level. Figure 3e shows that, following initial efficient knockdown of the *nad9* mRNA at Day 1 after onset of *Rznad9*-L-PKTLS expression, the steady-state level of the target RNA went back to normal at Day 2, before dropping again at Day 3 and 4. The profile established by RT-qPCR was confirmed by northern blot (Figure S2). The *atp9* transcript level also showed a marked bounce-back profile upon expression and mitochondrial import of the *Rzatp9*-L-PKTLS RNA at the intermediate stage of *A. thaliana* development (Figure 3f). Further analyses showed a similar bounce-back kinetics for the steady-state level of the *cox3* mitochondrial mRNA (Figure 3h) upon induced expression of the corresponding *Rzcox3*-L-PKTLS ribozyme in light-grown *A. thaliana* plants at the intermediate stage of development (Figure 3g). Thus, the bounce-back profile of the RNA steady-state level was not restricted to a given target, but appeared to be potentially a general response of mitochondrial RNAs to the expression and organellar import of specific ribozymes at intermediate stages of growth. Moreover, the phenomenon was not either restricted to *A. thaliana*, as the mitochondrial *sdh3* mRNA level responded in the same way to the expression of the *Rzsdh3*-L-PKTLS RNA in *N. tabacum* plants at the equivalent developmental stage (Figure 3i,j). These observations altogether imply that the buffering mechanism that might operate at the early stage of development (Figure 3b,c and Figure S3) becomes looser at a later stage, but disturbance of the mitochondrial transcriptome through the knockdown of an individual mitochondrial RNA is still sensed and transiently complemented. Notably, ribozyme-mediated knockdown of the mitochondrial *matR* mRNA encoding an organellar splicing factor was previously shown to be efficient and continuous in *A. thaliana* plants at an intermediate developmental stage [21]. Transcriptomic buffering might thus mainly apply to mRNAs encoding subunits of the OXPHOS complexes.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Chimeric ribozyme expression and knockdown of steady-state levels of mitochondrial target RNAs in transformed seedlings at different stages of growth. (a–h) *A. thaliana* control seeds and seeds carrying the *Rzatp9*-L-PKTLS, *Rznad9*-L-PKTLS or *Rzcox3*-L-PKTLS transgene were germinated in the light on solid MS-agar medium. Plants at early stage (upper panel) or intermediate stage (lower panel) of development were transferred at Day 0 to wells in culture plates containing liquid medium supplemented with estradiol for transgene induction. Kinetics of induced expression of the *Rzatp9*-L-PKTLS (*a*), *Rznad9*-L-PKTLS (*d*) or *Rzcox3*-L-PKTLS (*g*) RNA and of the steady-state level of the mitochondrial *atp9* (*b*,*f*), *nad9* (*c*,*e*) or *cox3* (*h*) target RNA were analyzed by RT-qPCR with total RNA from plant samples collected each day from Day 0 to Day 4 post-induction. (i,j) Transformant *N. tabacum* carrying the *Rzsdh3*-L-PKTLS transgene was germinated in the light on solid MS-agar medium and transferred at intermediate stage of development to liquid medium supplemented with estradiol for transgene induction. Kinetics of induced expression of the *Rzsdh3*-L-PKTLS RNA (*i*) and of the steady-state level of the mitochondrial *sdh3* target RNA (*j*) were analyzed by RT-qPCR with total RNA from transformed plant samples collected each day from Day 0 to Day 4 post-induction. Data from three independent biological replicates were analyzed with the Student's *t*-test; NS = not significant; \* = *p* ≤ 0.05; \*\* = *p* ≤ 0.01; \*\*\* = *p* ≤ 0.001; \*\*\*\* = *p* ≤ 0.0001.</figcaption></figure>
Continuous knockdown of the *nad9* target RNA (Figure 4b), or of the *atp9* target RNA (Figure 4c) was recorded upon *Rznad9*-L-PKTLS (Figure 4a) or *Rzatp9*-L-PKTLS expression in *A. thaliana* plants at the bolting stage of development (up to 18 true leaves). Mitochondrial transcriptome surveillance and buffering thus seems to be progressively released at later developmental stages. Notably, continuous and strong knockdown of the *atp9* mRNA was also obtained earlier in flowering plants expressing the *Rzatp9*-L-PKTLS RNA [20].

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Chimeric ribozyme expression and knockdown of steady-state levels of mitochondrial target RNAs in transformed seedlings in different conditions. (a–e) *A. thaliana* control seedlings and seedlings carrying the *Rznad9*-L-PKTLS, or the *Rzatp9*-L-PKTLS transgene were grown in the light (upper panel) or in the dark (lower panel) on solid MS-agar medium. Plants at bolting stage of development (upper panel) were transferred at Day 0 to wells in culture plates containing liquid medium supplemented with estradiol for transgene induction. Kinetics of induced expression of the *Rznad9*-L-PKTLS RNA (a), of the steady-state level of the mitochondrial *nad9* target RNA (b) and of the steady-state level of the mitochondrial *atp9* target RNA (c) were analyzed by RT-qPCR with total RNA from plant samples collected each day from Day 0 to Day 4 post-induction. Plates with ten-day-old seedlings grown in the dark (lower panel) were overlayed with liquid medium supplemented with estradiol for transgene induction. Kinetics of induced expression of the *Rznad9*-L-PKTLS RNA (d) and of the steady-state level of the mitochondrial *nad9* target RNA (e) were analyzed by RT-qPCR with total RNA from plant samples collected each day from Day 0 to Day 4 post-induction. Data from three independent biological replicates were analyzed with the Student's *t*-test; NS = not significant; \* = *p* ≤ 0.05; \*\* = *p* ≤ 0.01; \*\*\* = *p* ≤ 0.001; \*\*\*\* = *p* ≤ 0.0001.</figcaption></figure>

Remarkably, a bounce-back profile of the *nad9* mRNA level (Figure 4e) was also observed when expressing the *Rznad9*-L-PKTLS transcript in 10-day-old *A. thaliana* seedlings grown in the dark
(Figure 4d), indicating that the phenomenon was related to the age and importance of the mitochondrial function rather than to the growth conditions *per se*.

## 3.2. CMS RNA-Triggered Modulation of Mitochondrial RNA Steady State Levels Depends on the Plant Developmental Stage

The possibility that physiological conditions might impair ribozyme cleavage activity at some developmental stages cannot be completely ruled out. In particular, even when optimized, ribozymes require a minimal free Mg⁺⁺ concentration. We thus developed a parallel strategy to affect the mitochondrial transcriptome and test its reactivity. The concept was to use the PKTLS shuttle to drive a CMS RNA into the mitochondria of *A. thaliana* plants. Mitochondria of maize lines with the CMS-S cytoplasm express characteristic transcripts that carry a specific open reading frame called *orf77* ([43] and references therein). The *orf77* is a chimera combining sequences of unknown origin with sequences derived from the mitochondrial *atp9* gene (Figure S4). Due to these similarities, it is hypothesized that the *orf77* transcript can impair *atp9* expression at the RNA or protein level. ATP9 is a core component of the F₀ moiety of the mitochondrial ATP synthase located at the inner membrane. Ten or more ATP9 subunits assemble to form an inner membrane-embedded, hydrophobic ring that functions as a rotor pumping protons from the intermembrane space to the matrix side. ATP9 is thus a major protein for mitochondrial function.

In the present work, we imported the maize *orf77* CMS-S RNA into mitochondria in *A. thaliana* plants, so as to analyze the resulting impact on the whole mitochondrial transcriptome. Similarly, to the above *trans*-cleaving ribozymes, the *orf77* sequence was attached to the PKTLS shuttle as a 5'-trailor (Figure 2) and placed under the control of the estradiol-inducible promoter system in the pER8 expression vector [23]. The *orf77*-PKTLS RNA was expressed from a nuclear transgene in stably transformed *A. thaliana* plants (Figure S6). Expression was induced with estradiol at different developmental stages of light-grown plants and in dark-grown plants. In this strategy, plants transformed with a pER8 plasmid deprived of transgene served as a control. The resulting impact on the whole mitochondrial transcriptome was analyzed day by day for up to four days after onset of estradiol induction. At all stages, *orf77*-PKTLS expression had no visually detectable phenotypic effects on the plants within this time frame (Figure 5 and Figure S5). In this case, the cargo sequence attached to the PKTLS was longer and we could not rely on a well-defined activity to confirm organellar targeting. Purified, RNase-treated mitochondria were thus prepared from *A. thaliana* plants expressing the transgene and organellar import of the *orf77*-PKTLS transcript was confirmed by RT-qPCR probing, which showed about tenfold enrichment of the *orf77*-PKTLS RNA in the mitochondrial fraction *versus* the actin and GAPDH nuclear transcripts (Figure S7).
<table><caption>**Table 1.** Mitochondrial RNAs whose level increased (in black) or decreased (in purple) significantly in relation with the expression and mitochondrial import of the *orf77-PKTLS* RNA.</caption><thead><tr><th>Conditions and<br>Developmental<br>Stage</th><th>Day 1</th><th>Day 2</th><th>Day 3</th><th>Day 4</th></tr></thead><tbody><tr><td rowspan="4">Light Early stage</td><td></td><td>↓nad5a*<br>↓atp4**<br>↑rps7*<br>↑AOX1d*</td><td>↓AOX1d***</td><td>↓AOX1d***</td></tr><tr><td rowspan="4">Light Intermediate<br>stage</td><td>↑nad5b**<br>↓ccmFN2**</td><td>↑nad5a**<br>↑atp4*<br>↑rpl2*<br>↑rpl5*</td><td>↑cob**<br>↓ccmFN2***</td><td>↑rps7*<br>↑rpl2**<br>↑rpl5*</td></tr><tr><td></td><td></td><td>↑AOX1d***</td><td>↑AOX1d***</td></tr><tr><td>↓nad2b**<br>↑nad5a**<br>↓nad5b**<br>↓nad9**<br>↓cox1**<br>↓ccmB**<br>↓mttB*<br>↓rpl5*</td><td>↓atp4*<br>↓atp9****<br>↓rps4**<br>↓rps7****<br>↓rpl2****<br>↓rpl5**<br>↓rrn18**</td><td>↑nad1b*<br>↓nad2b*<br>↓nad4**<br>↑nad5a***<br>↓nad5b*<br>↓nad6*<br>↓nad9**<br>↓cob***<br>↓cox1*<br>↓cox2**<br>↓cox3****<br>↑atp9*<br>↓ccmB***<br>↓ccmFC*<br>↓ccmFN1***<br>↓ccmFN2**<br>↓mttB**<br>↓rrn18***<br>↓rrn26**</td><td></td></tr><tr><td rowspan="4">Light Bolting stage</td><td></td><td></td><td>↓AOX1a*<br>↓AOX1d*</td><td>↓ccmFC***<br>↑matR*<br>↓matR*<br>↑rps7**<br>↓rpl5*<br>↓rrn26*</td></tr><tr><td></td><td></td><td></td><td>↓AOX1d**</td></tr><tr><td>↑nad2b**<br>↑nad5a**<br>↑nad5b****<br>↑nad9*<br>↓cox2*<br>↓matR*<br>↑rps7**<br>↓rpl5*<br>↓rrn26*</td><td>↑nad4*<br>↑nad5b*<br>↑cox1*<br>↓ccmFN2**</td><td>↓nad1b***<br>↓nad2b**<br>↑nad5b**<br>↓ccmFN2**<br>↓mttB**<br>↓rrn26**</td><td>↓ccmFC***<br>↓matR*<br>↓mttB**<br>↓rrn18****<br>↑rrn26**</td></tr><tr><td>↑AOX1d**</td><td></td><td></td><td></td></tr></tbody></table>

Significant variations of nuclear RNAs coding for the mitochondrially targeted enzyme alternative oxidase (AOX1, underlined) are also included. Expression of the *orf77-PKTLS* transgene was induced at Day 0 and sampling was run every day until day 4 post-induction. Analyses were done by RT-qPCR with total RNA. Mitochondrial RNAs showing level variations versus control at Day 0 were not further considered. Data from three independent biological replicates were analyzed with the Graph Pad Prism software; * = *p* ≤ 0.05; ** = *p* ≤ 0.01; *** = *p* ≤ 0.001; **** = *p* ≤ 0.0001.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Synthetic scheme representing the impact of *orf77-PKTLS* expression on the mitochondrial transcriptome in transformed *A. thaliana* plants at different stages of growth in the light and at 10 days of growth in the dark. Plants were grown to the appropriate stage and the transgene was induced with estradiol as in Figures 3 and 4. Following transgene induction, samples were collected each day from Day 0 to Day 4 post-induction and RNAs were analyzed by RT-qPCR. The total numbers of positively or negatively affected transcripts combining all daily samples for a given growth stage are indicated and represented by proportional circular areas. Detailed day-by-day and gene-by-gene results are given in Table 1. Data from three independent biological replicates were analyzed with the Graph Pad Prism version 7.01 software.</figcaption></figure>

As for ribozymes, the effect of *orf77-PKTLS* expression and mitochondrial import on the levels of essentially all mitochondrial RNAs was characterized in transformed *A. thaliana* plants at different growth stages in the light and in dark-grown seedlings. A synthetic scheme of the results is given in Figure 5, while detailed data are compiled in Table 1. The analyses confirmed the transcriptome buffering process highlighted above in the case of *trans*-ribozyme-mediated RNA knockdown. Induced expression of the *orf77-PKTLS* RNA in transgenic *A. thaliana* seedlings at an early stage of development (maximum 4 true leaves) had almost no effect on the mitochondrial transcriptome, with only 3 RNAs positively or negatively affected and only at Day 2. At an intermediate stage of development, i.e., the stage where ribozyme-mediated target knockdown showed a bounce-back profile (Figure 3e,f,h,j), transcriptome buffering became somehow looser and up to eight mitochondrial RNAs were affected as a whole. Finally, at the bolting stage buffering was essentially released, as up to 27 mitochondrial RNAs were positively or negatively affected in response to *orf77-PKTLS* expression and organellar import. Expressing the *orf77-PKTLS* RNA in dark-grown seedlings affected up to 17 mitochondrial RNAs.

As mentioned, *orf77* contains sequences from the *atp9* gene (Figure S4), so that the *orf77-PKTLS* transcript was expected to interfere with the *atp9* mRNA. Remarkably, such an effect was indeed observed, but only at the bolting stage, i.e., at a stage where the transcriptome buffering was released. The level of *atp9* decreased at Day 2 and Day 4, while increasing at Day 3 (Table 1), further suggesting a dynamic process. Like *atp9*, *mttB* and *rpl5* showed some modulation at the bolting stage, with an initial decrease and subsequent increase at Day 4.
Except in early stage light-grown seedlings, all genetic functions encoded by the mtDNA, i.e., OXPHOS subunits, biogenesis co-factors (especially cytochrome c maturation factors), as well as ribosomal proteins and rRNAs, were represented in the organellar transcriptome response to *orf77-PKTLS* expression and mitochondrial import. Strikingly, the mitochondrial genes affected at the bolting stage were essentially downregulated, while the limited set of reactive genes at the intermediate stage was upregulated, with the exception of the cytochrome c maturation gene *ccmFN2* (Table 1). It is tempting to speculate that upregulation at the intermediate stage might be representative for transcriptome buffering, as proposed above, while general downregulation at the bolting stage would reflect the release of such a compensation process. Due to these differences, *nad5b*, *cob*, *atp4*, *rpl2*, *rpl5* and *rps7* showed an opposite behaviour between the intermediate and bolting stages. For instance, *atp4* was upregulated by about 50% at the intermediate stage and downregulated by 50% at the bolting stage.

The response to *orf77-PKTLS* expression and mitochondrial import in 10-day-old dark-grown seedlings was more balanced between up- and downregulation. As for light-grown intermediate stage samples, this led to an opposite behaviour of a number of transcripts, like for example the different *nad* mRNAs or *cox1*, *versus* the bolting stage of light-grown plants (Table 1). Notably, as intermediate stage light-grown samples, dark-grown seedlings also showed a bounce-back profile in ribozyme-mediated knockdown assays (see above, § 3.1. and Figure 4e). When affected, *ccmFN2* was always downregulated, whichever the growth conditions or growth stage (Table 1).

### 3.3. Ribozyme-Mediated Knockdown of a Mitochondrial mRNA Triggers a Nuclear Transcriptome Response

Microarray analyses revealed that continuous ribozyme-directed knockdown of the *nad9* mitochondrial mRNA in light-grown *A. thaliana* plants at the bolting stage of development affected the level of multiple nuclear-encoded RNAs. However, the extent of *nad9* knockdown differed from one kinetics to the other for the samples from biological replicates dedicated to the microarray assays, so that we decided to perform a sparse PLS analysis of the data (Supplementary File 2). The latter highlighted the significant regulation of 208 genes by the ribozyme-directed knockdown of *nad9* (Supplementary File 3). Notably, these included as many as 12 upregulated (Table 2) and 7 downregulated (Table 3) transcription factor genes, as well as 12 upregulated (Table 2) and 8 downregulated (Table 3) small nucleolar RNA (snoRNA) genes. Downregulation of the questioned miR414 microRNA was also noticed (Table 3), as well as up- or downregulation of genes involved in signaling, RNA metabolism or translation (Tables 2 and 3). Successful knockdown of *nad9*, i.e., a major mitochondrial mRNA encoding an OXPHOS subunit, in light-grown mature plants thus appeared to trigger a complex nuclear response driven by retrograde/anterograde regulation processes.
<table><caption>**Table 2.** Selected *A. thaliana* nuclear-encoded RNAs upregulated upon knockdown of the mitochondrial *nad9* mRNA.</caption><thead><tr><th>Modulation</th><th>Annotation / Function / Organellar Localization</th><th>Gene</th></tr></thead><tbody><tr><td rowspan="15">Upregulated (opposite to nad9)</td><td>Basic helix-loop-helix (bHLH) DNA-binding superfamily protein; transcription factor; response to ethylene</td><td>AT1G05710</td></tr><tr><td>ABI3-interacting protein 3, AIP3, PFD4, PREFOLDIN 4; protein chaperone; ABI3 is an auxin-inducible transcription factor</td><td>AT1G08780</td></tr><tr><td>ATMYB60, Myb domain protein 60, MYB60; transcription factor; response to abscisic acid, jasmonic acid, salicylic acid</td><td>AT1G08810</td></tr><tr><td>Agamous-like 87, AGL87, MADS-box family protein; transcription factor</td><td>AT1G22590</td></tr><tr><td>GL2, GLABRA 2, HD-ZIP IV family of homeobox-leucine zipper protein with lipid-binding START domain; transcription factor</td><td>AT1G79840</td></tr><tr><td>DREB subfamily A-6 of ERF/AP2 transcription factor family; one AP2 domain; ethylene-activated</td><td>AT2G22200</td></tr><tr><td>AP2/B3-like transcriptional factor family protein; transcription factor</td><td>AT2G33720</td></tr><tr><td>FMA (FAMA), basic helix-loop-helix (bHLH) DNA-binding superfamily protein; transcription factor/ transcriptional activator</td><td>AT3G24140</td></tr><tr><td>EDF3, ethylene response DNA-binding factor 3; transcription factor; AP2 domain; ethylene responding</td><td>AT3G25730</td></tr><tr><td>IAA30, indole-3-acetic acid inducible protein 30; transcription factor; response to auxin</td><td>AT3G62100</td></tr><tr><td>B-BOX domain protein 23, BBX23; transcription factor</td><td>AT4G10240</td></tr><tr><td>PUCHI, ethylene response factor (ERF) subfamily B-1 of ERF/AP2 transcription factor family; one AP2 domain; ethylene response</td><td>AT5G18560</td></tr><tr><td>ARGOS, Auxin-regulated gene involved in organ size; response to ethylene, auxin; membrane, cytoplasm, mitochondrion</td><td>AT3G59900</td></tr><tr><td>EMB3103, Embryo-defective 3103, PDM2, Pigment-Defective Mutant2; pentatricopeptide repeat (PPR) superfamily protein; endonuclease; chloroplast/mitochondrion</td><td>AT1G10910</td></tr><tr><td>RNH1C, RNase H family protein, RNase H domain-containing protein; chloroplast</td><td>AT1G24090</td></tr><tr><td>RING-finger, DEAD-like helicase, PHD and SNF2 domain-containing protein</td><td>AT2G40770</td></tr><tr><td>AT-SR34B, Serine/arginine-rich protein splicing factor 34B, SR34B</td><td>AT4G02430</td></tr><tr><td>Mitochondrial nuclease 1, MNU1; putative endonuclease or glycosyl hydrolase; mitochondrial RNA 5'-end processing; chloroplast/mitochondrion</td><td>AT5G64710</td></tr><tr><td>EMB2394, Embryo-defective 2394; structural constituent of chloroplast ribosome; response to cytokinin</td><td>AT1G05190</td></tr><tr><td>rRNA, cytosolic small ribosomal subunit</td><td>AT2G01010</td></tr><tr><td>5.8S rRNA, cytosolic large ribosomal subunit</td><td>AT2G01020</td></tr><tr><td>40S ribosomal protein S23 (RPS23A), ribosomal protein S12/S23 family protein; cytosolic small ribosomal subunit</td><td>AT3G09680</td></tr></tbody></table>
<table><caption>**Table 2.** *Cont.*</caption><thead><tr><th>Modulation</th><th>Annotation / Function / Organellar Localization</th><th>Gene</th></tr></thead><tbody><tr><td></td><td>Ribosomal protein L15, RPL15; large subunit of the chloroplast<br>ribosome; response to cytokinin</td><td>AT3G25920</td></tr><tr><td></td><td>EMB3126, Embryo-defective 3126, plastid ribosomal protein<br>L1, PRPL1; ribosomal protein L1p/L10e family; chloroplast<br>large ribosomal subunit</td><td>AT3G63490</td></tr><tr><td></td><td>60S ribosomal protein L31, RPL31B; cytosolic large<br>ribosomal subunit</td><td>AT4G26230</td></tr><tr><td></td><td>U3 ribonucleoprotein, Utp family protein; rRNA processing</td><td>AT5G08600</td></tr><tr><td></td><td>MA3 domain-containing translation regulatory factor 1, MRF1;<br>colocalizes with cytosolic large ribosomal subunit;<br>isomerase activity</td><td>AT5G63190</td></tr><tr><td>snoRNA</td><td></td><td>AT1G03743; AT1G19373;<br>AT1G19376; AT1G75166;<br>AT2G35387; AT3G27865;<br>AT3G47342; AT3G47347;<br>AT3G58193; AT3G58196;<br>AT4G39366; AT5G44286</td></tr></tbody></table>

Data taken from sparse PLS analysis of microarray results obtained with RNAs from light-grown *A. thaliana* plants at the bolting stage. Individual Day 0 to Day 4 RNA samples from three independent biological replicates were analyzed (i.e., 15 test samples and 15 control samples). Organellar localization is mentioned where annotated. Selection was for genes of potential relevance in regulation mechanisms, including transcription factors, hormone signaling or RNA metabolism.

To cast off the effects of photosynthesis processes and chloroplast pathways, and try to simplify the picture, microarray assays were continued with RNAs from 10-day-old dark-grown *A. thaliana* seedlings expressing the *Rznad9*-L-PKTLS chimeric ribozyme. In this case, the mitochondrial transcriptome was constrained by buffering, so that the response was also expected to be representative for the moderate bounce-back profile of *nad9* knockdown in dark-grown plants (see Figure 4e). Microarray analyses of the RNAs from dark-grown *A. thaliana* seedlings expressing the *Rznad9*-L-PKTLS ribozyme (Supplementary File 4) revealed that the number of differentially expressed genes reached 3046 at Day 1 and then dropped to 1622 and 416 at Day 2 and Day 3, before re-increasing to 959 at Day 4. Sparse PLS sorting of these data (Supplementary File 5) showed a regulation of 49 genes by the ribozyme-directed knockdown of *nad9* (Supplementary File 3), including 3 downregulated transcription factor genes, as well as genes involved in signaling or RNA metabolism (Table 4). In parallel, we mined the original microarray data for modulation of transcripts encoding proteins shown or predicted to localize to mitochondria (Supplementary File 6). Interestingly, the analysis pointed to upregulation of a series of members of the large “small auxin upregulated RNAs” (SAUR) family, as well as to upregulation of transcripts for cytokinin and auxin response factors (Table 5), suggesting a hormone signaling-related process. MapMan analysis of the data also pointed to hormone metabolism pathways and especially to auxin-regulated genes. Conversely, downregulation of a jasmonate-responding factor, of a set of pentatricopeptide repeat (PPR) proteins and of further RNA metabolism factors was noticed (Table 5). Also, to be noted, the knockdown of *nad9* in dark-grown *A. thaliana* seedlings negatively affected the level of a series of mitochondrial mRNAs encoding subunits of OXPHOS complexes or ribosomal proteins, while mRNAs coding for biogenesis co-factors like CCM factors or MATR were little affected (Figure S8).
<table><caption>**Table 3.** Selected *A. thaliana* nuclear-encoded RNAs downregulated upon knockdown of the mitochondrial *nad9* mRNA.</caption><thead><tr><th>Modulation</th><th>Annotation/Function/Organellar Localization</th><th>Gene</th></tr></thead><tbody><tr><td rowspan="15">Downregulated<br>(as nad9)</td><td>MicroRNA414, mir414, short open reading frame 16, SORF16; miRNA; identified as a translated small open reading frame by ribosome profiling</td><td>AT1G67195</td></tr><tr><td>B-BOX domain protein 25, BBX25, Salt tolerance homologue, STH; transcription factor; zinc ion binding</td><td>AT2G31380</td></tr><tr><td>ATIBH1, IBH1, ILI1 binding BHLH 1 ILI1 binding bHLH 1; transcription factor; brassinosteroid signaling; gibberellic acid signaling</td><td>AT2G43060</td></tr><tr><td>Arabidopsis thaliana response regulator 2, ARR5, ATRR2, IBC6, induced by cytokinin 6, response regulator 5, RR5; transcription repressor; cytokinin signaling</td><td>AT3G48100</td></tr><tr><td>Basic helix-loop-helix (bHLH) DNA-binding superfamily protein; transcription factor</td><td>AT4G01460</td></tr><tr><td>KELP; homodimers or heterodimers with the kiwi protein; transcriptional co-activator</td><td>AT4G10920</td></tr><tr><td>WRKY DNA-binding protein 24, ATWRKY24, WRKY24; WRKY transcription factor group II-c</td><td>AT5G41570</td></tr><tr><td>BOA, Brother of lux ARRHYTHMO; transcription factor; circadian clock; mRNA cell-to-cell mobile</td><td>AT5G59570</td></tr><tr><td>SAUR53, small auxin-upregulated RNA 53, SAUR-like auxin-responsive protein family; <u>mitochondrion</u></td><td>AT1G19840</td></tr><tr><td>SAUR65, small auxin-upregulated RNA 65, SAUR-like auxin-responsive protein family; membrane; <u>mitochondrion</u></td><td>AT1G29460</td></tr><tr><td>ATCLE19, CLAVATA3/ESR-related 19, CLE19, embryo surrounding region 19, ESR19; receptor binding; signal transduction; <u>mitochondrion</u></td><td>AT3G24225</td></tr><tr><td>VQ motif-containing protein 29, VQ29; response to hypoxia</td><td>AT4G37710</td></tr><tr><td>RNA-binding (RRM/RBD/RNP motifs) family protein</td><td>AT1G33470</td></tr><tr><td>ATRNS1, Ribonuclease 1, RNS1; endoribonuclease</td><td>AT2G02990</td></tr><tr><td>Eukaryotic translation initiation factor 2 (eIF-2) family protein; cytosol</td><td>AT1G76820</td></tr><tr><td>EMB3113, Embryo-defective 3113, ribosomal protein S5, RPS5, SCA1, SCABRA 1; structural component of the 70S chloroplast ribosome; <u>mitochondrion</u> small ribosomal subunit</td><td>AT2G33800</td></tr><tr><td>U3 containing 90S pre-ribosomal complex subunit</td><td>AT2G43110</td></tr><tr><td>ATRAB8D, ATRABE1B, RAB GTPase homolog E1B, RABE1B; translation elongation factor; membrane; <u>chloroplast</u></td><td>AT4G20360</td></tr><tr><td>snoRNA</td><td>AT2G35382; AT2G43137; AT2G43138; AT2G43139; AT2G43141; AT4G02550; AT4G02555; AT4G13245</td></tr></tbody></table>

Data taken from sparse PLS analysis of microarray results obtained with RNAs from light-grown *A. thaliana* plants at the bolting stage. Individual Day 0 to Day 4 RNA samples from three independent biological replicates were analyzed (i.e., 15 test samples and 15 control samples). Organellar localization is mentioned where annotated. Selection was for genes of potential relevance in regulation mechanisms, including transcription factors, hormone signaling or RNA metabolism.
<table><caption>**Table 4.** Selected *A. thaliana* nuclear-encoded RNAs up- or downregulated upon knockdown of the mitochondrial *nad9* mRNA.</caption><thead><tr><th>Modulation</th><th>Annotation/Function/Organellar Localization</th><th>Gene</th></tr></thead><tbody><tr><td rowspan="3">Upregulated (opposite to nad9)</td><td>Small nuclear RNA U6acat, mRNA splicing</td><td>AT5G40395</td></tr><tr><td>ANAC028, NAC domain-containing protein 28, NAC028; transcription factor</td><td>AT1G65910</td></tr><tr><td>Transcription elongation factor Spt5; KOW domain</td><td>AT2G34210</td></tr><tr><td rowspan="3">Downregulated (as nad9)</td><td>MYR2, homeodomain-like superfamily protein; transcription factor</td><td>AT3G04030</td></tr><tr><td>HVA22-like protein F, HVA22F; membrane protein; response to abscisic acid</td><td>AT2G42820</td></tr><tr><td>ALY2 RNA-binding (RRM/RBD/RNP motifs) family protein; mRNA transport</td><td>AT5G02530</td></tr></tbody></table>

Data taken from sparse PLS analysis of microarray results obtained with RNAs from 10-day-old dark-grown *A. thaliana* seedlings. Individual Day 0 to Day 4 RNA samples from two independent biological replicates were analyzed (i.e., 10 test samples and 10 control samples). Selection was for genes of potential relevance in regulation mechanisms, including transcription factors, hormone signaling or RNA metabolism.

<table><caption>**Table 5.** Selected *A. thaliana* nuclear RNAs up- or downregulated upon knockdown of the mitochondrial *nad9* mRNA and coding for proteins predicted to be mitochondrion-targeted.</caption><thead><tr><th>Modulation</th><th>Annotation/Function/Organellar Localization</th><th>Gene</th></tr></thead><tbody><tr><td rowspan="15">Upregulated (opposite to nad9)</td><td>auxin-responsive family protein (SAUR72)</td><td>AT3G12830</td></tr><tr><td>auxin-responsive protein, putative (SAUR9)</td><td>AT4G36110</td></tr><tr><td>auxin-responsive family protein (SAUR41)</td><td>AT1G16510</td></tr><tr><td>auxin-responsive protein-related (SAUR77)</td><td>AT1G17345</td></tr><tr><td>auxin-responsive family protein (SAUR53); <u>mitochondrion</u></td><td>AT1G19840</td></tr><tr><td>auxin-responsive family protein (SAUR71)</td><td>AT1G56150</td></tr><tr><td>auxin-responsive family protein (SAUR52); <u>mitochondrion</u></td><td>AT1G75590</td></tr><tr><td>auxin-responsive protein-related (SAUR36); <u>mitochondrion/nucleus</u> (nucleus confirmed by GFP targeting experiments)</td><td>AT2G45210</td></tr><tr><td>auxin-responsive family protein (SAUR59); <u>mitochondrion</u></td><td>AT3G60690</td></tr><tr><td>auxin-responsive protein, putative (SAUR25); <u>mitochondrion</u></td><td>AT4G13790</td></tr><tr><td>auxin-responsive family protein (SAUR1); <u>chloroplast</u></td><td>AT4G34770</td></tr><tr><td>auxin-responsive protein, putative (SAUR23)</td><td>AT5G18060</td></tr><tr><td>MIF1 (MINI ZINC FINGER 1); transcription factor; response to abscisic acid, auxin, brassinosteroid, cytokinin, gibberellin; <u>nucleus</u> (cytosol reported from GFP data)</td><td>AT1G74660</td></tr><tr><td>PLS (POLARIS), cytokinin and auxin responses; <u>mitochondrion</u></td><td>AT4G39403</td></tr><tr><td>Tetratricopeptide repeat (TPR)-like superfamily protein; <u>mitochondrion</u></td><td>AT1G28690</td></tr></tbody></table>
<table><caption>**Table 5. Cont.**</caption><thead><tr><th>Modulation</th><th>Annotation/Function/Organellar Localization</th><th>Gene</th></tr></thead><tbody><tr><td rowspan="10">Downregulated (as nad9)</td><td>S-RBP11, SMALL RNA-BINDING PROTEIN 11; salt stress response; chloroplast</td><td>AT5G06210</td></tr><tr><td>ARGAH2, Arginine amidohydrolase 2, response to jasmonate; chloroplast/mitochondrion (confirmed by MS data)</td><td>AT4G08870</td></tr><tr><td>Tetratricopeptide repeat (TPR)-like superfamily protein; chloroplast/mitochondrion (chloroplast confirmed by MS data)</td><td>AT2G37230</td></tr><tr><td>pentatricopeptide (PPR) repeat-containing protein</td><td>AT3G62470</td></tr><tr><td>pentatricopeptide (PPR) repeat-containing protein</td><td>AT4G01030</td></tr><tr><td>EMB1417 (embryo-defective 1417), PPR protein; RNA binding; endonuclease activity</td><td>AT4G21190</td></tr><tr><td>pentatricopeptide (PPR) repeat-containing protein</td><td>AT4G21880</td></tr><tr><td>pentatricopeptide (PPR) repeat-containing protein</td><td>AT5G65560</td></tr><tr><td>EMB1586, Embryo-defective 1586, increased size exclusion limit 1, ISE1; DEAD-box RNA helicase; chloroplast/mitochondrion (mitochondrion confirmed by GFP data)</td><td>AT1G12770</td></tr><tr><td>AGS1, AHG2-1 suppressor 1, bacterial-type poly(A) polymerase; mRNA polyadenylation; chloroplast/mitochondrion</td><td>AT2G17580</td></tr><tr><td>ATTRM2A, TRM2A, tRNA methyltransferase 2A, RNA methyltransferase family protein; cytoplasm</td><td>AT3G21300</td></tr><tr><td>ATP-dependent RNA helicase</td><td>AT5G39840</td></tr><tr><td>ARFB1A, ATARFB1A (ADP-ribosylation factor B1A); GTP binding; protein transport; Golgi apparatus</td><td>AT2G15310</td></tr></tbody></table>

Data taken from microarray analysis of RNAs from 10-day-old dark-grown *A. thaliana* seedlings. Individual Day 0 to Day 4 RNA samples from two independent biological replicates were analyzed (i.e., 10 test samples and 10 control samples). All proteins included have a prediction consensus for mitochondrial targeting in the SUBA database [44] (http://suba.plantenergy.uwa.edu.au/). Experimentally reported localization compiled in the SUBA database is also indicated. Selection was for genes of potential relevance in regulation mechanisms, including hormone signaling or RNA metabolism.

## 4. Discussion

All complexes of the mitochondrial OXPHOS chain consist of subunits encoded by the nuclear genome and subunits expressed from mitochondrial genes, hence the idea of a coordinated expression of the two genomes. Nevertheless, according to the available data mtDNA transcription and mitochondrial transcript levels did not seem to be major regulation issues by themselves, so that mitochondrial genetic control is considered to occur mainly at the post-transcriptional level [3,4]. However, over the past decades, there have been a number of reports describing changes in the abundance of specific mitochondrial transcripts during plant growth [45], such as during the development of wheat leaves [46] or seed germination in maize and rice [47–49]. Variations in the steady state levels of a number of mitochondrial mRNAs were also observed with *A. thaliana* nuclear mutants deficient for the editing or splicing of a given organellar RNA [50–53]. The question of a control of mitochondrial gene expression at the transcriptome level, connected with mitochondrial retrograde pathways, can thus still be raised. So-called intergenomic signaling has been characterized in yeast [54]. Unaffected by respiration and metabolic signals, the process requires mtDNA and appears as a direct communication between the mitochondrial and nuclear genetic systems. Similarly, the idea of an involvement of chloroplast transcription in signaling to the nucleus has been put forward [55,56].

The present work follows this line. We disturbed the mitochondrial transcriptome in plants at different stages of development through specific in vivo knockdown of different organellar target mRNAs or introduction of a CMS RNA into mitochondria. The ribozymes and the CMS RNA were expressed from an inducible promoter, so that at each growth stage the plants were in a regular physiological configuration before induction. This enabled to reveal early and specific events,
a notable advantage over the use of mutants that would have to cope with the mutation and adapt
permanently their pathways. Also, to be noted, transgene expression upon single estradiol induction
is transient and lasts only a few days. The question of the existence of transcriptome control within
mitochondria and of intergenic coordination could then be addressed by analyzing the mitochondrial
and nuclear responses. Remarkably, the observations highlighted RNA control mechanisms in
mitochondria with relation to the plant developmental stage. The mitochondrial transcriptome
appeared to be tightly constrained in young, actively developing seedlings, as ribozyme-directed
target RNA knockdown was not observed at that stage and the *orf77* CMS RNA triggered almost no
effect. This suggests in particular that in young seedlings high mitochondrial transcription rates [57]
might compensate negative variations of RNA levels and "protect" the transcriptome at a growth
stage where efficient mitochondrial activity is essential. A corollary of such a hypothesis would be
that changes affecting the organellar transcriptome can be "sensed". Compensation became only
transient at an intermediate stage of development, perhaps due to limited RNA synthesis, hence the
bounce-back effect. Finally, at the bolting stage, ribozyme-mediated knockdown of target RNAs became
efficient. Both *nad9* knockdown and organellar targeting of the *orf77* CMS RNA (presumed competitor
of *atp9*) were then associated with a decrease of many mitochondrial RNAs, especially of mRNAs
encoding subunits of OXPHOS complexes, suggesting a coordination and cohesion of the mitochondrial
transcriptome. According to complementary analyses, the observed decrease of mitochondrial RNAs
was not associated with variations in mtDNA copy number or in the number of mitochondria per
cell. Notably, directed knockdown of *nad9*, i.e., a mitochondrial mRNA encoding a subunit of an
OXPHOS complex, negatively affected the level of many other mitochondrial mRNAs encoding
OXPHOS subunits or ribosome subunits. By contrast, in earlier experiments knockdown of the *matR*
mitochondrial mRNA, encoding a maturase, affected intron splicing, but had limited effect on the
level of intronless mitochondrial RNAs [21], suggesting that the mitochondrial transcriptome might be
somehow coordinated to optimize expression of OXPHOS complex subunits.

As revealed by microarray analyses, successful ribozyme-mediated knockdown of *nad9*, i.e., a mitochondrial mRNA encoding a subunit of an OXPHOS complex, positively or negatively affected the level of multiple nuclear-encoded RNAs. This implies that directed alteration of the mitochondrial transcriptome resulted in a retrograde/anterograde response. The effect was fast, starting right away from the first day of ribozyme expression and knockdown of the target RNA, which implies a very fast retrograde signaling and raises the question of the nature of such a signaling. So far described retrograde signals are essentially metabolic signals sent to inform the nucleo-cytosolic compartment about the status of the organelles or resulting from impaired organellar function or biogenesis [54,58–60]. In particular, different signaling pathways respond to OXPHOS chain dysfunction. Our earlier studies showed that the profile of the mitochondrial OXPHOS complexes was still unaffected in plant cells four days after onset of ribozyme-mediated knockdown of the *atp9* mRNA [20]. Also, the level of MATR protein started to decrease significantly only from day 4 after onset of ribozyme-mediated knockdown of the *matR* mRNA [21]. With such a slow turnover of mitochondrial proteins, it seems unlikely that mitochondrial respiratory and metabolic dysfunction would occur fast enough to generate the signals underlying the observed retrograde response to the knockdown of the *nad9* mRNA. The data presented here thus support the idea of a fast communication between the mitochondrial and nuclear genomes through yet to be identified signals. Mitochondrial-nuclear communication through intergenomic signaling described in yeast was suggested to involve the Abf1p transcription factor and seemed indeed to be distinct from regulation activated by respiration dysfunction [54]. On the other hand, it was shown that sigma factors utilized by chloroplast-encoded RNA polymerase to transcribe specific sets of plastid genes are a source of retrograde signals to the nucleus [56]. However, no significant up- or downregulation of the nuclear genes for mitochondrial RNA polymerases ever appeared in our RT-qPCR and microarray analyses in response to *nad9* knockdown. It can still be hypothesized that other regular nuclear-encoded factors involved in plant mitochondrial transcription play a regulatory role.
Whichever the mitochondrial signal, ribozyme-mediated *nad9* knockdown triggered an anterograde response involving a series of nuclear transcription factors, either in light-grown or dark-grown seedlings. A number of these belonged to families highlighted in previous studies as taking part in mitochondrion-nucleus communication during plant growth, development, defense or stress (MYB, WRKY, ARR, AP2/EREBP, IAA, bHLH, ANAC/NAM) [14,61,62]. Many of them are responsive to hormone signaling. Notably, a meta-analysis of available datasets identified the plant hormones abscissic acid, auxin, cytokinin, jasmonate, and salicylic acid as main regulators of mitochondrial functions [61]. Modulation of the miR414 micro RNA (miRNA) level in response to *nad9* knockdown was also in line. miR414 was proposed earlier to be involved in the regulation of the genes participating in auxin response pathways and energy metabolism during plant embryonic development [63] and a number of its target genes were reported to play an important role in hormone signal transduction, especially in the auxin and ethylene pathways [64]. miR164 is a further miRNA highlighted in previous mitochondrion-nucleus communication analyses [61]. Also striking is the up- or downregulation of a number of small nucleolar RNAs (snoRNAs) in response to *nad9* knockdown. Beyond their fundamental role in the modification and processing of ribosomal RNAs, snoRNAs recently appeared to be a source of short regulatory RNAs [65,66].

Mining the data to identify proteins potentially sent to the mitochondria upon *nad9* knockdown in dark-grown seedlings in turn highlighted some interesting features. Firstly, the results suggest an upregulation of putative mitochondrial auxin-responsive proteins encoded by small auxin upregulated RNAs (SAUR), as well as further hormone-responsive factors (Table 5). The family of SAUR auxin-responsive genes has about 60–140 members in most higher plant species [67] and plays a central role in fine-tuning of growth in response to auxin or other hormones. Whether the SAUR proteins pointed out in our analyses are indeed targeted into mitochondria and have a specific functional or regulatory activity in there remains to be investigated in the future. The second remarkable feature highlighted by the analyses was a downregulation of a number of putative mitochondrial RNA metabolism factors, especially pentatricopeptide repeat (PPR) proteins (Table 5). PPR proteins in plants are encoded by a large family of nuclear genes with up to 800 members and most of them are targeted to the organelles [68–70]. Aiding in RNA editing, maturation, stabilisation, intron splicing, transcription and translation of organellar genes, they are major RNA metabolism factors in both mitochondria and chloroplasts. The AT4G01030 PPR gene highlighted in our plant study was annotated as coding for an editing factor, while endonuclease activity was attributed to the protein encoded by the AT4G21190 gene. The function of the remaining ones (Table 5) remains to be investigated and experimental evidence is needed to confirm mitochondrial targeting.

**5. Conclusions**

In conclusion, establishing directed manipulation of the mitochondrial transcriptome brought evidence that RNA levels are controlled and coordinated in plant mitochondria, and that interfering with specific mitochondrial transcripts causes a retrograde response that affects the nuclear transcriptome and triggers an anterograde response. The exact mechanisms underlying the observed transcriptome control and intercompartment cross-talk are still to be deciphered, but they do not seem to rely on signals that would be generated by mitochondrial metabolic dysfunction. While the control of mitochondrial gene expression had so far mostly been documented at the post-transcriptional level, our observations introduce a broader and more integrated view of the genetic regulation in the organelles.

**6. Patents**

Commercial use of the ribozyme/PKTLS strategy to manipulate mitochondrial transcripts is covered by patent Number WO 2010/031918.

**Supplementary Materials:** The following materials are available online at http://www.mdpi.com/2073-4409/8/6/583/s1, Supplementary File 1 containing: Figure S1. Phenotype of *A. thaliana* transgenic seedlings expressing the target-deprived *Rzsdh3* ribozyme or the *Rznad9* ribozyme; Figure S2. *Rznad9*-L-PKTLS-mediated knockdown of
*nad9* steady-state levels in transformed *A. thaliana* seedlings at intermediate stage of growth; Figure S3. Short time kinetics of ribozyme expression and target RNA steady-state level in *A. thaliana* plants at the early stage of development; Figure S4. Common regions between the *orf77* and *atp9* sequences; Figure S5. Phenotype of *A. thaliana* transgenic seedlings expressing the *orf77-PKTLS* RNA *versus* control seedlings; Figure S6. Kinetics of *orf77-PKTLS* RNA expression at the early stage of development in light-grown transformed *A. thaliana* plants; Figure S7. The *orf77-PKTLS* RNA is recovered in the mitochondrial fraction; Figure S8. Knockdown of *nad9* negatively impacts the level of further mitochondrial mRNAs coding for subunits of oxidative phosphorylation complexes; Supplementary File 2: sparse PLS analysis of microarray data from light-grown seedlings; Supplementary File 3: significant results from sparse PLS analyses of microarray data from light-grown and dark-grown seedlings; Supplementary File 4: microarray analysis of samples from dark-grown seedlings; Supplementary File 5: sparse PLS analysis of microarray data from dark-grown seedlings; Supplementary File 6: mining microarray data from dark-grown seedlings for transcripts encoding proteins shown or predicted to localize to mitochondria

**Author Contributions:** Conceptualization, A.D.; Formal analysis, E.D.; Investigation, A.K.N., R.K.I., D.M., R.V., M.G., E.W. and L.S.-T.; Methodology, M.S.; Project administration, A.D.; Supervision, J.B., F.W.-L., J.M.G. and A.D.; Writing—original draft, A.D.

**Funding:** This research was funded by grants from the French State program "Investments for the Future" (LABEX ANR-11-LABX-0057_MITOCROSS and ANR-10-LABX-0040-SPS), from the French National Research Agency (ANR-06-MRAR-037-02, ANR-09-BLAN-0240-01), from the Polish National Science Centre (UMO-2013/09/B/NZ1/03359), as well as by regular funding from the French National Center for Scientific Research (CNRS-UPR2357) and the University of Strasbourg. R.V. was supported by a fellowship shared by the CNRS and the Région Alsace. Cooperation between the IBMP and the CABB is supported by the PERIDOT Program developed in France by the Ministère de l'Europe et des Affaires Etrangères (MEAE) and the Ministère de l'Enseignement Supérieur, de la Recherche et de l'Innovation (MESRI) and in Pakistan by the Higher Education Commission (HEC).

**Acknowledgments:** We thank Anne Cosset, Anne-Marie Duchêne and Clara Chicois for their help at some stages of this work.

**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. Shokolenko, I.N.; Alexeyev, M.F. Mitochondrial transcription in mammalian cells. *Front. Biosci.* **2017**, *22*, 835-853.

2. Holec, S.; Lange, H.; Kühn, K.; Alioua, M.; Börner, T.; Gagliardi, D. Relaxed transcription in *Arabidopsis* mitochondria is counterbalanced by RNA stability control mediated by polyadenylation and polynucleotide phosphorylase. *Mol. Cell. Biol.* **2006**, *26*, 2869-2876. [CrossRef] [PubMed]

3. Rackham, O.; Mercer, T.R.; Filipovska, A. The human mitochondrial transcriptome and the RNA-binding proteins that regulate its expression. *Wiley Interdiscip. Rev. RNA* **2012**, *3*, 675-695. [CrossRef] [PubMed]

4. Small, I.D.; Rackham, O.; Filipovska, A. Organelle transcriptomes: Products of a deconstructed genome. *Curr. Opin. Microbiol.* **2013**, *16*, 652-658. [CrossRef] [PubMed]

5. Pearce, S.F.; Rebelo-Guiomar, P.; D'Souza, A.R.; Powell, C.A.; Van Haute, L.; Minczuk, M. Regulation of mammalian mitochondrial gene expression: Recent advances. *Trends Biochem. Sci.* **2017**, *42*, 625-639. [CrossRef] [PubMed]

6. Rebelo-Guiomar, P.; Powell, C.A.; Van Haute, L.; Minczuk, M. The mammalian mitochondrial epitranscriptome. *Biochim. Biophys. Acta Gene Regul. Mech.* **2019**, *1862*, 429-446. [CrossRef] [PubMed]

7. Pohjoismaki, J.L.; Kruger, M.; Al-Furoukh, N.; Lagerstedt, A.; Karhunen, P.J.; Braun, T. Postnatal cardiomyocyte growth and mitochondrial reorganization cause multiple changes in the proteome of human cardiomyocytes. *Mol. BioSyst.* **2013**, *9*, 1210-1219. [CrossRef] [PubMed]

8. Mposhi, A.; Van der Wijst, M.G.; Faber, K.N.; Rots, M.G. Regulation of mitochondrial gene expression, the epigenetic enigma. *Front. Biosci. (Landmark Ed)* **2017**, *22*, 1099-1113.

9. Barshad, G.; Marom, S.; Cohen, T.; Mishmar, D. Mitochondrial DNA transcription and its regulation: An evolutionary perspective. *Trends Genet.* **2018**, *34*, 682-692. [CrossRef]

10. Robles, P.; Quesada, V. Transcriptional and post-transcriptional regulation of organellar gene expression (OGE) and its roles in plant salt tolerance. *Int. J. Mol. Sci.* **2019**, *20*, 1056. [CrossRef]

11. Woodson, J.D.; Chory, J. Coordination of gene expression between organellar and nuclear genomes. *Nat. Rev. Genet.* **2008**, *9*, 383-395. [CrossRef] [PubMed]

12. Cogliati, S.; Lorenzi, I.; Rigoni, G.; Caicci, F.; Soriano, M.E. Regulation of mitochondrial electron transport chain assembly. *J. Mol. Biol.* **2018**, *430*, 4849-4873. [CrossRef] [PubMed]
13. Weinhouse, C. Mitochondrial-epigenetic crosstalk in environmental toxicology. *Toxicology* **2017**, *391*, 5–17.
[CrossRef] [PubMed]

14. Ng, S.; De Clercq, I.; Van Aken, O.; Law, S.R.; Ivanova, A.; Willems, P.; Giraud, E.; Van Breusegem, F.; Whelan, J. Anterograde and retrograde regulation of nuclear genes encoding mitochondrial proteins during growth, development, and stress. *Mol. Plant* **2014**, *7*, 1075–1093.
[CrossRef] [PubMed]

15. Liu, Z.; Butow, R.A. Mitochondrial retrograde signaling. *Annu. Rev. Genet.* **2006**, *40*, 159–185.
[CrossRef]

16. Dekkers, B.J.; Pearce, S.; van Bolderen-Veldkamp, R.P.; Marshall, A.; Widera, P.; Gilbert, J.; Drost, H.G.; Bassel, G.W.; Müller, K.; King, J.R.; et al. Transcriptional dynamics of two seed compartments with opposing roles in *Arabidopsis* seed germination. *Plant Physiol.* **2013**, *163*, 205–215.
[CrossRef]

17. Law, S.R.; Narsai, R.; Taylor, N.L.; Delannoy, E.; Carrie, C.; Giraud, E.; Millar, A.H.; Small, I.; Whelan, J. Nucleotide and RNA metabolism prime translational initiation in the earliest events of mitochondrial biogenesis during *Arabidopsis* germination. *Plant Physiol.* **2012**, *158*, 1610–1627.
[CrossRef]

18. Law, S.R.; Narsai, R.; Whelan, J. Mitochondrial biogenesis in plants during seed germination. *Mitochondrion* **2014**, *19*, 214–221.
[CrossRef]

19. Salinas, T.; Duchêne, A.M.; Maréchal-Drouard, L. Recent advances in tRNA mitochondrial import. *Trends Biochem. Sci.* **2008**, *33*, 320–329.
[CrossRef] [PubMed]

20. Val, R.; Wyszko, E.; Valentin, C.; Szymanski, M.; Cosset, A.; Alioua, M.; Dreher, T.W.; Barciszewski, J.; Dietrich, A. Organelle trafficking of chimeric ribozymes and genetic manipulation of mitochondria. *Nucleic Acids Res.* **2011**, *39*, 9262–9274.
[CrossRef] [PubMed]

21. Sultan, L.D.; Mileshina, D.; Grewe, F.; Rolle, K.; Abudraham, S.; Glodowicz, P.; Niazi, A.K.; Keren, I.; Shevtsov, S.; Klipcan, L.; et al. The reverse transcriptase/RNA maturase protein MatR is required for the splicing of various group II introns in *Brassicaceae* mitochondria. *Plant Cell* **2016**, *28*, 2805–2829.
[CrossRef] [PubMed]

22. Perrotta, A.T.; Been, M.D. A pseudoknot-like structure required for efficient self-cleavage of hepatitis delta virus RNA. *Nature* **1991**, *350*, 434–436.
[CrossRef] [PubMed]

23. Zuo, J.; Niu, Q.W.; Chua, N.H. Technical advance: An estrogen receptor-based transactivator XVE mediates highly inducible gene expression in transgenic plants. *Plant J.* **2000**, *24*, 265–273.
[CrossRef] [PubMed]

24. Mathews, D.H.; Sabina, J.; Zuker, M.; Turner, D.H. Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. *J. Mol. Biol.* **1999**, *288*, 911–940.
[CrossRef] [PubMed]

25. Van Engelen, F.A.; Molthoff, J.W.; Conner, A.J.; Nap, J.P.; Pereira, A.; Stiekema, W.J. pBINPLUS: An improved plant transformation vector based on pBIN19. *Transgenic Res.* **1995**, *4*, 288–290.
[CrossRef] [PubMed]

26. Clough, S.J.; Bent, A.F. Floral dip: A simplified method for *Agrobacterium*-mediated transformation of *Arabidopsis thaliana*. *Plant J.* **1998**, *16*, 735–743.
[CrossRef] [PubMed]

27. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3–new capabilities and interfaces. *Nucleic Acids Res.* **2012**, *40*, e115.
[CrossRef]

28. Sweetlove, L.J.; Taylor, N.L.; Leaver, C.J. Isolation of intact, functional mitochondria from the model plant *Arabidopsis thaliana*. *Methods Mol. Biol.* **2007**, *372*, 125–136.

29. Lurin, C.; Andrés, C.; Aubourg, S.; Bellaoui, M.; Bitton, F.; Bruyère, C.; Caboche, M.; Debast, C.; Gualberto, J.; Hoffmann, B.; et al. Genome-wide analysis of *Arabidopsis* pentatricopeptide repeat proteins reveals their essential role in organelle biogenesis. *Plant Cell* **2004**, *16*, 2089–2103.
[CrossRef]

30. Yang, Y.H.; Dudoit, S.; Luu, P.; Lin, D.M.; Peng, V.; Ngai, J.; Speed, T.P. Normalization for cDNA microarray data: A robust composite method addressing single and multiple slide systematic variation. *Nucleic Acids Res.* **2002**, *30*, e15.
[CrossRef]

31. Matthew, E.R.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. *limma* powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* **2015**, *43*, e47.

32. Smyth, G.K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat. Appl. Genet. Mol. Biol.* **2004**, *3*, 3.
[CrossRef] [PubMed]

33. Guedj, M.; Robin, S.; Celisse, A.; Nuel, G. Kerfdr: A semi-parametric kernel-based approach to local false discovery rate estimation. *BMC Bioinform.* **2009**, *10*, 84.
[CrossRef] [PubMed]

34. Storey, J.D.; Tibshirani, R. Statistical significance for genomewide studies. *Proc. Natl. Acad. Sci. USA* **2003**, *100*, 9440–9445.
[CrossRef] [PubMed]
35. Chun, H.; Keleş, S. Sparse partial least squares regression for simultaneous dimension reduction and variable selection. *J. R. Stat. Soc. Series B Stat. Methodol.* **2010**, *72*, 3–25. [CrossRef] [PubMed]

36. Gagnot, S.; Tamby, J.P.; Martin-Magniette, M.L.; Bitton, F.; Taconnat, L.; Balzergue, S.; Aubourg, S.; Renou, J.P.; Lecharny, A.; Brunaud, V. CATdb: A public access to *Arabidopsis* transcriptome data from the URGV-CATMA platform. *Nucleic Acids Res.* **2008**, *36*, D986–D990. [CrossRef] [PubMed]

37. Zaag, R.; Tamby, J.P.; Guichard, C.; Tariq, Z.; Rigaill, G.; Delannoy, E.; Renou, J.P.; Balzergue, S.; Mary-Huard, T.; Aubourg, S.; et al. GEM2Net: From gene expression modeling to–Omics networks, a new CATdb module to investigate *Arabidopsis thaliana* genes involved in stress response. *Nucleic Acids Res.* **2015**, *43*, D1010–D1017. [CrossRef]

38. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for functional genomics data sets–update. *Nucleic Acids Res.* **2013**, *41*, D991–D995. [CrossRef]

39. Brazma, A.; Hingamp, P.; Quackenbush, J.; Sherlock, G.; Spellman, P.; Stoeckert, C.; Aach, J.; Ansorge, W.; Ball, C.A.; Causton, H.C.; et al. Minimum information about a microarray experiment (MIAME)–toward standards for microarray data. *Nat. Genet.* **2001**, *29*, 365–371. [CrossRef]

40. Matsuda, D.; Dreher, T.W. The tRNA-like structure of Turnip yellow mosaic virus RNA is a 3'-translational enhancer. *Virology* **2004**, *321*, 36–46. [CrossRef]

41. Delage, L.; Duchêne, A.M.; Zaepfel, M.; Maréchal-Drouard, L. The anticodon and the D-domain sequences are essential determinants for plant cytosolic tRNA(Val) import into mitochondria. *Plant J.* **2003**, *34*, 623–633. [CrossRef] [PubMed]

42. Dietrich, A.; Maréchal-Drouard, L.; Carneiro, V.; Cosset, A.; Small, I. A single base change prevents import of cytosolic tRNA(Ala) into mitochondria in transgenic plants. *Plant J.* **1996**, *10*, 913–918. [CrossRef]

43. Matera, J.T.; Monroe, J.; Smelser, W.; Gabay-Laughnan, S.; Newton, K.J. Unique changes in mitochondrial genomes associated with reversions of S-type cytoplasmic male sterility in maizemar. *PLoS ONE* **2011**, *6*, e23405. [CrossRef] [PubMed]

44. Heazlewood, J.L.; Verboom, R.E.; Tonti-Filippini, J.; Small, I.; Millar, A.H. SUBA: The *Arabidopsis* subcellular database. *Nucleic Acids Res.* **2007**, *35*, D213–D218. [CrossRef]

45. Millar, A.H.; Whelan, J.; Soole, K.L.; Day, D.A. Organization and regulation of mitochondrial respiration in plants. *Annu. Rev. Plant Biol.* **2011**, *62*, 79–104. [CrossRef]

46. Topping, J.F.; Leaver, C.J. Mitochondrial gene expression during wheat leaf development. *Planta* **1990**, *182*, 399–407. [CrossRef]

47. Howell, K.A.; Millar, A.H.; Whelan, J. Ordered assembly of mitochondria during rice germination begins with pro-mitochondrial structures rich in components of the protein import apparatus. *Plant Mol. Biol.* **2006**, *60*, 201–223. [CrossRef] [PubMed]

48. Howell, K.A.; Cheng, K.; Murcha, M.W.; Jenkin, L.E.; Millar, A.H.; Whelan, J. Oxygen initiation of respiration and mitochondrial biogenesis in rice. *J. Biol. Chem.* **2007**, *282*, 15619–15631. [CrossRef] [PubMed]

49. Ishizaki, K.; Larson, T.R.; Schauer, N.; Fernie, A.R.; Graham, I.A.; Leaver, C.J. The critical role of *Arabidopsis* electron-transfer flavoprotein:ubiquinone oxidoreductase during dark-induced starvation. *Plant Cell* **2005**, *17*, 2587–2600. [CrossRef]

50. Sung, T.Y.; Tseng, C.C.; Hsieh, M.H. The SLO1 PPR protein is required for RNA editing at multiple sites with similar upstream sequences in *Arabidopsis* mitochondria. *Plant J.* **2010**, *63*, 499–511. [CrossRef] [PubMed]

51. Hsieh, W.Y.; Liao, J.C.; Chang, C.Y.; Harrison, T.; Boucher, C.; Hsieh, M.H. The SLOW GROWTH3 pentatricopeptide repeat protein is required for the splicing of mitochondrial NADH dehydrogenase subunit 7 intron 2 in *Arabidopsis*. *Plant Physiol.* **2015**, *168*, 490–501. [CrossRef] [PubMed]

52. Koprivova, A.; Des Francs-Small, C.C.; Calder, G.; Mugford, S.T.; Tanz, S.; Lee, B.R.; Zechmann, B.; Small, I.; Kopriva, S. Identification of a pentatricopeptide repeat protein implicated in splicing of intron 1 of mitochondrial *nad7* transcripts. *J. Biol. Chem.* **2010**, *285*, 32192–32199. [CrossRef] [PubMed]

53. Nakagawa, N.; Sakurai, N. A mutation in At-nMat1a, which encodes a nuclear gene having high similarity to group II intron maturase, causes impaired splicing of mitochondrial NAD4 transcript and altered carbon metabolism in *Arabidopsis thaliana*. *Plant Cell Physiol.* **2006**, *47*, 772–783. [CrossRef] [PubMed]

54. Woo, D.K.; Phang, T.L.; Trawick, J.D.; Poyton, R.O. Multiple pathways of mitochondrial-nuclear communication in yeast: Intergenomic signaling involves ABF1 and affects a different set of genes than retrograde regulation. *Biochim. Biophys. Acta* **2009**, *1789*, 135–145. [CrossRef] [PubMed]
55. Börner, T. The discovery of plastid-to-nucleus retrograde signaling–A personal perspective. *Protoplasma* **2017**, *254*, 1845–1855. [CrossRef] [PubMed]

56. Woodson, J.D.; Perez-Ruiz, J.M.; Schmit, R.J.; Ecker, J.R.; Chory, J. Sigma factor-mediated plastid retrograde signals control nuclear gene expression. *Plant J.* **2013**, *73*, 1–13. [CrossRef] [PubMed]

57. Holec, S.; Lange, H.; Canaday, J.; Gagliardi, D. Coping with cryptic and defective transcripts in plant mitochondria. *Biochim. Biophys. Acta* **2008**, *1779*, 566–573. [CrossRef]

58. Nam, M.; Akie, T.E.; Sanosaka, M.; Craige, S.M.; Kant, S.; Keaney, J.F., Jr.; Cooper, M.P. Mitochondrial retrograde signaling connects respiratory capacity to thermogenic gene expression. *Sci. Rep.* **2017**, *7*, 2013. [CrossRef]

59. da Cunha, F.M.; Torelli, N.Q.; Kowaltowski, A.J. Mitochondrial retrograde signaling: Triggers, pathways, and outcomes. *Oxid. Med. Cell. Longev.* **2015**, *2015*, 482582. [CrossRef]

60. Jazwinski, S.M. The retrograde response: A conserved compensatory reaction to damage from within and from without. *Prog. Mol. Biol. Transl. Sci.* **2014**, *127*, 133–154.

61. Berkowitz, O.; De Clercq, I.; Van Breusegem, F.; Whelan, J. Interaction between hormonal and mitochondrial signalling during growth, development and in plant defence responses. *Plant Cell Environ.* **2016**, *39*, 1127–1139. [CrossRef] [PubMed]

62. Kühn, K.; Yin, G.; Duncan, O.; Law, S.R.; Kubiszewski-Jakubiak, S.; Kaur, P.; Meyer, E.; Wang, Y.; Colas des Francs-Small, C.; Giraud, E.; et al. Decreasing electron flux through the cytochrome and/or alternative respiratory pathways triggers common and distinct cellular responses dependent on growth conditions. *Plant Physiol.* **2015**, *167*, 228–250. [CrossRef]

63. Zhang, F.; Zhao, J.; Xu, S.; Fang, W.; Chen, F.; Teng, N. MicroRNA and putative target discoveries in *Chrysanthemum* polyploidy breeding. *Int. J. Genomics* **2017**, *2017*, 6790478. [CrossRef] [PubMed]

64. Ma, C.; Yang, J.; Cheng, Q.; Mao, A.; Zhang, J.; Wang, S.; Weng, Y.; Wen, C. Comparative analysis of miRNA and mRNA abundance in determinate cucumber by high-throughput sequencing. *PLoS ONE* **2018**, *13*, e0190691.

65. Scott, M.S.; Ono, M. From snoRNA to miRNA: Dual function regulatory non-coding RNAs. *Biochimie* **2011**, *93*, 1987–1992. [CrossRef] [PubMed]

66. Stepanov, G.A.; Filippova, J.A.; Komissarov, A.B.; Kuligina, E.V.; Richter, V.A.; Semenov, D.V. Regulatory role of small nucleolar RNAs in human diseases. *Biomed. Res. Int.* **2015**, *2015*, 206849. [CrossRef] [PubMed]

67. Stortenbeker, N.; Bemer, M. The SAUR gene family: The plant's toolbox for adaptation of growth and development. *J. Exp. Bot.* **2019**, *70*, 17–27. [CrossRef] [PubMed]

68. Chen, L.; Li, Y.X.; Li, C.; Shi, Y.; Song, Y.; Zhang, D.; Li, Y.; Wang, T. Genome-wide analysis of the pentatricopeptide repeat gene family in different maize genomes and its important role in kernel development. *BMC Plant Biol.* **2018**, *18*, 366. [CrossRef] [PubMed]

69. Gorchs Rovira, A.; Smith, A.G. PPR proteins–Orchestrators of organelle RNA metabolism. *Physiol. Plant* **2019**. [CrossRef]

70. Xing, H.; Fu, X.; Yang, C.; Tang, X.; Guo, L.; Li, C.; Xu, C.; Luo, K. Genome-wide investigation of pentatricopeptide repeat gene family in poplar and their expression analysis in response to biotic and abiotic stresses. *Sci. Rep.* **2018**, *8*, 2817. [CrossRef]

<figure><img src="image_9.png" /></figure>