The Use of Nanopore Sequencing to Analyze the Chloroplast Transcriptome Part I: Library Preparation.

Global understanding of plastid gene expression has always been impaired by its complexity. RNA splicing, editing, and intercistronic processing create multiple transcripts isoforms that can hardly be resolved using traditional molecular biology techniques. During the last decade, the wide adoption of RNA-seq-based techniques has, however, allowed an unprecedented understanding of all the different steps of chloroplast gene expression, from transcription to translation. Current strategies are nonetheless unable to identify and quantify full length transcripts isoforms, a limitation that can now be overcome using Nanopore Sequencing. We here provide a complete protocol to produce, from total leaf RNA, cDNA libraries ready for Nanopore sequencing. While most Nanopore protocols take advantage of the mRNA polyA tail we here first ligate an RNA adapter to the 3' ends of the RNAs and use it to initiate the template switching reverse transcription. The cDNA is then prepared and indexed for use with the regular Oxford Nanopore v14 chemistry. This protocol is of particular interest to researchers willing to simultaneously study the multiple post-transcriptional processes prevalent in the chloroplast.

The Use of Nanopore Sequencing to Analyze the Chloroplast Transcriptome Part II: Bioinformatic Analyzes and Virtual RNA Blots.

Nanopore sequencing of full-length cDNAs offers unprecedented details of the plastid RNA metabolism. After the generation of the nanopore reads, several bioinformatic steps are required to analyze the data. In this chapter, we describe in a few simple command lines the processing and mapping of the reads as well as the generation of virtual Northern blots as a simple and familiar way to visualize Nanopore sequencing data.

DiffSegR: an RNA-seq data driven method for differential expression analysis using changepoint detection.

To fully understand gene regulation, it is necessary to have a thorough understanding of both the transcriptome and the enzymatic and RNA-binding activities that shape it. While many RNA-Seq-based tools have been developed to analyze the transcriptome, most only consider the abundance of sequencing reads along annotated patterns (such as genes). These annotations are typically incomplete, leading to errors in the differential expression analysis. To address this issue, we present DiffSegR - an R package that enables the discovery of transcriptome-wide expression differences between two biological conditions using RNA-Seq data. DiffSegR does not require prior annotation and uses a multiple changepoints detection algorithm to identify the boundaries of differentially expressed regions in the per-base log2 fold change. In a few minutes of computation, DiffSegR could rightfully predict the role of chloroplast ribonuclease Mini-III in rRNA maturation and chloroplast ribonuclease PNPase in (3'/5')-degradation of rRNA, mRNA and tRNA precursors as well as intron accumulation. We believe DiffSegR will benefit biologists working on transcriptomics as it allows access to information from a layer of the transcriptome overlooked by the classical differential expression analysis pipelines widely used today. DiffSegR is available at https://aliehrmann.github.io/DiffSegR/index.html.

An mTRAN-mRNA interaction mediates mitochondrial translation initiation in plants.

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that "mitochondrial translation factors" mTRAN1 and mTRAN2 are land plant-specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)-like RNA binding proteins of the mitoribosomal "small" subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5' regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

Cell specialization and coordination in Arabidopsis leaves upon pathogenic attack revealed by scRNA-seq.

Plant defense responses involve several biological processes that allow plants to fight against pathogenic attacks. How these different processes are orchestrated within organs and depend on specific cell types is poorly known. Here, using single-cell RNA sequencing (scRNA-seq) technology on three independent biological replicates, we identified several cell populations representing the core transcriptional responses of wild-type Arabidopsis leaves inoculated with the bacterial pathogen Pseudomonas syringae DC3000. Among these populations, we retrieved major cell types of the leaves (mesophyll, guard, epidermal, companion, and vascular S cells) with which we could associate characteristic transcriptional reprogramming and regulators, thereby specifying different cell-type responses to the pathogen. Further analyses of transcriptional dynamics, on the basis of inference of cell trajectories, indicated that the different cell types, in addition to their characteristic defense responses, can also share similar modules of gene reprogramming, uncovering a ubiquitous antagonism between immune and susceptible processes. Moreover, it appears that the defense responses of vascular S cells, epidermal cells, and mesophyll cells can evolve along two separate paths, one converging toward an identical cell fate, characterized mostly by lignification and detoxification functions. As this divergence does not correspond to the differentiation between immune and susceptible cells, we speculate that this might reflect the discrimination between cell-autonomous and non-cell-autonomous responses. Altogether our data provide an upgraded framework to describe, explore, and explain the specialization and the coordination of plant cell responses upon pathogenic challenge.

High-Throughput Protein-Protein Interactions Screening Using Pool-Based Liquid Yeast Two-Hybrid Pipeline.

Because of its adaptability to high-throughput approaches and a low operating cost, the yeast two-hybrid (Y2H) assay remains the most widely used one for high-throughput protein-protein interactions (PPI) mapping experiments. Here we provide a detailed protocol for a liquid culture-based high-throughput binary protein-protein Y2H screen pipeline of a pool of 50 proteins used as baits against a collection of ~12,000 Arabidopsis proteins encoded by sequence-verified open reading frames (ORFs).

Limiting etioplast gene expression induces apical hook twisting during skotomorphogenesis of Arabidopsis seedlings.

When covered by a layer of soil, seedling development follows a dark-specific program (skotomorphogenesis). In the dark, seedlings consist of small, non-green cotyledons, a long hypocotyl, and an apical hook to protect meristematic cells. We recently highlighted the role played by mitochondria in the high energy-consuming reprogramming of Arabidopsis skotomorphogenesis. Here, the role played by plastids, another energy-supplying organelle, in skotomorphogenesis is investigated. This study was conducted in dark conditions to exclude light signals so as to better focus on those produced by plastids. It was found that limitation of plastid gene expression (PGE) induced an exaggerated apical hook bending. Inhibition of PGE was obtained at the levels of transcription and translation using the antibiotics rifampicin (RIF) and spectinomycin, respectively, as well as plastid RPOTp RNA polymerase mutants. RIF-treated seedlings also showed expression induction of marker nuclear genes for mitochondrial stress, perturbation of mitochondrial metabolism, increased ROS levels, and an augmented capacity of oxygen consumption by mitochondrial alternative oxidases (AOXs). AOXs act to prevent overreduction of the mitochondrial electron transport chain. Previously, we reported that AOX1A, the main AOX isoform, is a key component in the developmental response to mitochondrial respiration deficiency. In this work, we suggest the involvement of AOX1A in the response to PGE dysfunction and propose the importance of signaling between plastids and mitochondria. Finally, it was found that seedling architecture reprogramming in response to RIF was independent of canonical organelle retrograde pathways and the ethylene signaling pathway.

Evolutionary analyses and expression patterns of TCP genes in Ranunculales.

TCP transcription factors play a role in a large number of developmental processes and are at the crossroads of numerous hormonal biosynthetic and signaling pathways. The complete repertoire of TCP genes has already been characterized in several plant species, but not in any species of early diverging eudicots. We focused on the order Ranunculales because of its phylogenetic position as sister group to all other eudicots and its important morphological diversity. Results show that all the TCP genes expressed in the floral transcriptome of Nigella damascena (Ranunculaceae) are the orthologs of the TCP genes previously identified from the fully sequenced genome of Aquilegia coerulea. Phylogenetic analyses combined with the identification of conserved amino acid motifs suggest that six paralogous genes of class I TCP transcription factors were present in the common ancestor of angiosperms. We highlight independent duplications in core eudicots and Ranunculales within the class I and class II subfamilies, resulting in different numbers of paralogs within the main subclasses of TCP genes. This has most probably major consequences on the functional diversification of these genes in different plant clades. The expression patterns of TCP genes in Nigella damascena were consistent with the general suggestion that CIN and class I TCP genes may have redundant roles or take part in same pathways, while CYC/TB1 genes have more specific actions. Our findings open the way for future studies at the tissue level, and for investigating redundancy and subfunctionalisation in TCP genes and their role in the evolution of morphological novelties.

Genome-Wide Transcriptomic Analysis of the Effects of Infection with the Hemibiotrophic Fungus Colletotrichum lindemuthianum on Common Bean.

Bean anthracnose caused by the hemibiotrophic fungus Colletotrichum lindemuthianum is one of the most important diseases of common bean (Phaseolus vulgaris) in the world. In the present study, the whole transcriptome of common bean infected with C. lindemuthianum during compatible and incompatible interactions was characterized at 48 and 72 hpi, corresponding to the biotrophy phase of the infection cycle. Our results highlight the prominent role of pathogenesis-related (PR) genes from the PR10/Bet vI family as well as a complex interplay of different plant hormone pathways including Ethylene, Salicylic acid (SA) and Jasmonic acid pathways. Gene Ontology enrichment analysis reveals that infected common bean seedlings responded by down-regulation of photosynthesis, ubiquitination-mediated proteolysis and cell wall modifications. In infected common bean, SA biosynthesis seems to be based on the PAL pathway instead of the ICS pathway, contrarily to what is described in Arabidopsis. Interestingly, ~30 NLR were up-regulated in both contexts. Overall, our results suggest that the difference between the compatible and incompatible reaction is more a question of timing and strength, than a massive difference in differentially expressed genes between these two contexts. Finally, we used RT-qPCR to validate the expression patterns of several genes, and the results showed an excellent agreement with deep sequencing.

Analysis of the Plant Mitochondrial Transcriptome.

Transcriptome analyses are widely performed, but the analysis of organelle genome expression is often overlooked. In this chapter, we describe three methods to analyse the accumulation, splicing, editing and processing of plant mitochondrial transcript expression: a classical RT-qPCR assay, an RNA-Seq approach with its bioinformatical and statistical analysis pipeline, as well as a useful complementary technique, the Northern-blot analysis, using short biotinylated oligonucleotides as probes.

Full Length Transcriptome Highlights the Coordination of Plastid Transcript Processing.

Plastid gene expression involves many post-transcriptional maturation steps resulting in a complex transcriptome composed of multiple isoforms. Although short-read RNA-Seq has considerably improved our understanding of the molecular mechanisms controlling these processes, it is unable to sequence full-length transcripts. This information is crucial, however, when it comes to understanding the interplay between the various steps of plastid gene expression. Here, we describe a protocol to study the plastid transcriptome using nanopore sequencing. In the leaf of Arabidopsis thaliana, with about 1.5 million strand-specific reads mapped to the chloroplast genome, we could recapitulate most of the complexity of the plastid transcriptome (polygenic transcripts, multiple isoforms associated with post-transcriptional processing) using virtual Northern blots. Even if the transcripts longer than about 2500 nucleotides were missing, the study of the co-occurrence of editing and splicing events identified 42 pairs of events that were not occurring independently. This study also highlighted a preferential chronology of maturation events with splicing happening after most sites were edited.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L.

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

The Genomic Impact of Mycoheterotrophy in Orchids.

Mycoheterotrophic plants have lost the ability to photosynthesize and obtain essential mineral and organic nutrients from associated soil fungi. Despite involving radical changes in life history traits and ecological requirements, the transition from autotrophy to mycoheterotrophy has occurred independently in many major lineages of land plants, most frequently in Orchidaceae. Yet the molecular mechanisms underlying this shift are still poorly understood. A comparison of the transcriptomes of Epipogium aphyllum and Neottia nidus-avis, two completely mycoheterotrophic orchids, to other autotrophic and mycoheterotrophic orchids showed the unexpected retention of several genes associated with photosynthetic activities. In addition to these selected retentions, the analysis of their expression profiles showed that many orthologs had inverted underground/aboveground expression ratios compared to autotrophic species. Fatty acid and amino acid biosynthesis as well as primary cell wall metabolism were among the pathways most impacted by this expression reprogramming. Our study suggests that the shift in nutritional mode from autotrophy to mycoheterotrophy remodeled the architecture of the plant metabolism but was associated primarily with function losses rather than metabolic innovations.

BdERECTA controls vasculature patterning and phloem-xylem organization in Brachypodium distachyon.

BACKGROUND: The vascular system of plants consists of two main tissue types, xylem and phloem. These tissues are organized into vascular bundles that are arranged into a complex network running through the plant that is essential for the viability of land plants. Despite their obvious importance, the genes involved in the organization of vascular tissues remain poorly understood in grasses. RESULTS: We studied in detail the vascular network in stems from the model grass Brachypodium distachyon (Brachypodium) and identified a large set of genes differentially expressed in vascular bundles versus parenchyma tissues. To decipher the underlying molecular mechanisms of vascularization in grasses, we conducted a forward genetic screen for abnormal vasculature. We identified a mutation that severely affected the organization of vascular tissues. This mutant displayed defects in anastomosis of the vascular network and uncommon amphivasal vascular bundles. The causal mutation is a premature stop codon in ERECTA, a LRR receptor-like serine/threonine-protein kinase. Mutations in this gene are pleiotropic indicating that it serves multiple roles during plant development. This mutant also displayed changes in cell wall composition, gene expression and hormone homeostasis. CONCLUSION: In summary, ERECTA has a pleiotropic role in Brachypodium. We propose a major role of ERECTA in vasculature anastomosis and vascular tissue organization in Brachypodium.

The Consequences of a Disruption in Cyto-Nuclear Coadaptation on the Molecular Response to a Nitrate Starvation in Arabidopsis.

Mitochondria and chloroplasts are important actors in the plant nutritional efficiency. So, it could be expected that a disruption of the coadaptation between nuclear and organellar genomes impact plant response to nutrient stresses. We addressed this issue using two Arabidopsis accessions, namely Ct1 and Jea, and their reciprocal cytolines possessing the nuclear genome from one parent and the organellar genomes of the other one. We measured gene expression, and quantified proteins and metabolites under N starvation and non-limiting conditions. We observed a typical response to N starvation at the phenotype and molecular levels. The phenotypical response to N starvation was similar in the cytolines compared to the parents. However, we observed an effect of the disruption of genomic coadaptation at the molecular levels, distinct from the previously described responses to organellar stresses. Strikingly, genes differentially expressed in cytolines compared to parents were mainly repressed in the cytolines. These genes encoded more mitochondrial and nuclear proteins than randomly expected, while N starvation responsive ones were enriched in genes for chloroplast and nuclear proteins. In cytolines, the non-coadapted cytonuclear genomic combination tends to modulate the response to N starvation observed in the parental lines on various biological processes.

Landscape of the Noncoding Transcriptome Response of Two Arabidopsis Ecotypes to Phosphate Starvation.

Root architecture varies widely between species; it even varies between ecotypes of the same species, despite strong conservation of the coding portion of their genomes. By contrast, noncoding RNAs evolve rapidly between ecotypes and may control their differential responses to the environment, since several long noncoding RNAs (lncRNAs) are known to quantitatively regulate gene expression. Roots from ecotypes Columbia and Landsberg erecta of Arabidopsis (Arabidopsis thaliana) respond differently to phosphate starvation. Here, we compared transcriptomes (mRNAs, lncRNAs, and small RNAs) of root tips from these two ecotypes during early phosphate starvation. We identified thousands of lncRNAs that were largely conserved at the DNA level in these ecotypes. In contrast to coding genes, many lncRNAs were specifically transcribed in one ecotype and/or differentially expressed between ecotypes independent of phosphate availability. We further characterized these ecotype-related lncRNAs and studied their link with small interfering RNAs. Our analysis identified 675 lncRNAs differentially expressed between the two ecotypes, including antisense RNAs targeting key regulators of root-growth responses. Misregulation of several lincRNAs showed that at least two ecotype-related lncRNAs regulate primary root growth in ecotype Columbia. RNA-sequencing analysis following deregulation of lncRNA NPC48 revealed a potential link with root growth and transport functions. This exploration of the noncoding transcriptome identified ecotype-specific lncRNA-mediated regulation in root apexes. The noncoding genome may harbor further mechanisms involved in ecotype adaptation of roots to different soil environments.

The Analysis of the Editing Defects in the dyw2 Mutant Provides New Clues for the Prediction of RNA Targets of Arabidopsis E+-Class PPR Proteins.

C to U editing is one of the post-transcriptional steps which are required for the proper expression of chloroplast and mitochondrial genes in plants. It depends on several proteins acting together which include the PLS-class pentatricopeptide repeat proteins (PPR). DYW2 was recently shown to be required for the editing of many sites in both organelles. In particular almost all the sites associated with the E+ subfamily of PPR proteins are depending on DYW2, suggesting that DYW2 is required for the function of E+-type PPR proteins. Here we strengthened this link by identifying 16 major editing sites controlled by 3 PPR proteins: OTP90, a DYW-type PPR and PGN and MEF37, 2 E+-type PPR proteins. A re-analysis of the DYW2 editotype showed that the 49 sites known to be associated with the 18 characterized E+-type PPR proteins all depend on DYW2. Considering only the 288 DYW2-dependent editing sites as potential E+-type PPR sites, instead of the 795 known editing sites, improves the performances of binding predictions systems based on the PPR code for E+-type PPR proteins. However, it does not compensate for poor binding predictions.

Thirteen New Plastid Genomes from Mixotrophic and Autotrophic Species Provide Insights into Heterotrophy Evolution in Neottieae Orchids.

Mixotrophic species use both organic and mineral carbon sources. Some mixotrophic plants combine photosynthesis and a nutrition called mycoheterotrophy, where carbon is obtained from fungi forming mycorrhizal symbiosis with their roots. These species can lose photosynthetic abilities and evolve full mycoheterotrophy. Besides morphological changes, the latter transition is associated with a deep alteration of the plastid genome. Photosynthesis-related genes are lost first, followed by housekeeping genes, eventually resulting in a highly reduced genome. Whether relaxation of selective constraints already occurs for the plastid genome of mixotrophic species, which remain photosynthetic, is unclear. This is partly due to the difficulty of comparing plastid genomes of autotrophic, mixotrophic, and mycoheterotrophic species in a narrow phylogenetic framework. We address this question in the orchid tribe Neottieae, where this large assortment of nutrition types occurs. We sequenced 13 new plastid genomes, including 9 mixotrophic species and covering all 6 Neottieae genera. We investigated selective pressure on plastid genes in each nutrition type and conducted a phylogenetic inference of the group. Surprisingly, photosynthesis-related genes did not experience selection relaxation in mixotrophic species compared with autotrophic relatives. Conversely, we observed evidence for selection intensification for some plastid genes. Photosynthesis is thus still under purifying selection, maybe because of its role in fruit formation and thus reproductive success. Phylogenetic analysis resolved most relationships, but short branches at the base of the tree suggest an evolutionary radiation at the beginning of Neottieae history, which, we hypothesize, may be linked to mixotrophy emergence.

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

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.