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.

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.

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.

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.

In situ transcriptomic and metabolomic study of the loss of photosynthesis in the leaves of mixotrophic plants exploiting fungi.

Mycoheterotrophic plants have lost photosynthesis and obtain carbon through mycorrhizal fungi colonizing their roots. They are likely to have evolved from mixotrophic ancestors, which rely on both photosynthesis and fungal carbon for their development. Whereas our understanding of the ecological and genomic changes associated with the evolutionary shift to mycoheterotrophy is deepening, little information is known about the specific metabolic and physiological features driving this evolution. We investigated this issue in naturally occurring achlorophyllous variants of temperate mixotrophic orchids. We carried out an integrated transcriptomic and metabolomic analysis of the response to achlorophylly in the leaves of three mixotrophic species sampled in natura. Achlorophyllous leaves showed major impairment of their photosynthetic and mineral nutrition functions, strong accumulation of free amino acids, overexpression of enzymes and transporters related to sugars, amino acids and fatty acid catabolism, as well as induction of some autophagy-related and biotic stress genes. Such changes were reminiscent of these reported for variegated leaves and appeared to be symptomatic of a carbon starvation response. Rather than decisive metabolic innovations, we suggest that the evolution towards mycoheterotrophy in orchids is more likely to be reliant on the versatility of plant metabolism and an ability to exploit fungal organic resources, especially amino acids, to replace missing photosynthates.

Synthetic data sets for the identification of key ingredients for RNA-seq differential analysis.

Numerous statistical pipelines are now available for the differential analysis of gene expression measured with RNA-sequencing technology. Most of them are based on similar statistical frameworks after normalization, differing primarily in the choice of data distribution, mean and variance estimation strategy and data filtering. We propose an evaluation of the impact of these choices when few biological replicates are available through the use of synthetic data sets. This framework is based on real data sets and allows the exploration of various scenarios differing in the proportion of non-differentially expressed genes. Hence, it provides an evaluation of the key ingredients of the differential analysis, free of the biases associated with the simulation of data using parametric models. Our results show the relevance of a proper modeling of the mean by using linear or generalized linear modeling. Once the mean is properly modeled, the impact of the other parameters on the performance of the test is much less important. Finally, we propose to use the simple visualization of the raw P-value histogram as a practical evaluation criterion of the performance of differential analysis methods on real data sets.

Two interacting PPR proteins are major Arabidopsis editing factors in plastid and mitochondria.

RNA editing is converting hundreds of cytosines into uridines during organelle gene expression of land plants. The pentatricopeptide repeat (PPR) proteins are at the core of this posttranscriptional RNA modification. Even if a PPR protein defines the editing site, a DYW domain of the same or another PPR protein is believed to catalyze the deamination. To give insight into the organelle RNA editosome, we performed tandem affinity purification of the plastidial CHLOROPLAST BIOGENESIS 19 (CLB19) PPR editing factor. Two PPR proteins, dually targeted to mitochondria and chloroplasts, were identified as potential partners of CLB19. These two proteins, a P-type PPR and a member of a small PPR-DYW subfamily, were shown to interact in yeast. Insertional mutations resulted in embryo lethality that could be rescued by embryo-specific complementation. A transcriptome analysis of these complemented plants showed major editing defects in both organelles with a very high PPR type specificity, indicating that the two proteins are core members of E+-type PPR editosomes.

Homoeologous exchanges cause extensive dosage-dependent gene expression changes in an allopolyploid crop.

Structural variation is a major source of genetic diversity and an important substrate for selection. In allopolyploids, homoeologous exchanges (i.e. between the constituent subgenomes) are a very frequent type of structural variant. However, their direct impact on gene content and gene expression had not been determined. Here, we used a tissue-specific mRNA-Seq dataset to measure the consequences of homoeologous exchanges (HE) on gene expression in Brassica napus, a representative allotetraploid crop. We demonstrate that expression changes are proportional to the change in gene copy number triggered by the HEs. Thus, when homoeologous gene pairs have unbalanced transcriptional contributions before the HE, duplication of one copy does not accurately compensate for loss of the other and combined homoeologue expression also changes. These effects are, however, mitigated over time. This study sheds light on the origins, timing and functional consequences of homeologous exchanges in allopolyploids. It demonstrates that the interplay between new structural variation and the resulting impacts on gene expression, influences allopolyploid genome evolution.

Function of the Plant DNA Polymerase Epsilon in Replicative Stress Sensing, a Genetic Analysis.

Faithful transmission of the genetic information is essential in all living organisms. DNA replication is therefore a critical step of cell proliferation, because of the potential occurrence of replication errors or DNA damage when progression of a replication fork is hampered causing replicative stress. Like other types of DNA damage, replicative stress activates the DNA damage response, a signaling cascade allowing cell cycle arrest and repair of lesions. The replicative DNA polymerase ε (Pol ε) was shown to activate the S-phase checkpoint in yeast in response to replicative stress, but whether this mechanism functions in multicellular eukaryotes remains unclear. Here, we explored the genetic interaction between Pol ε and the main elements of the DNA damage response in Arabidopsis (Arabidopsis thaliana). We found that mutations affecting the polymerase domain of Pol ε trigger ATR-dependent signaling leading to SOG1 activation, WEE1-dependent cell cycle inhibition, and tolerance to replicative stress induced by hydroxyurea, but result in enhanced sensitivity to a wide range of DNA damaging agents. Using knock-down lines, we also provide evidence for the direct role of Pol ε in replicative stress sensing. Together, our results demonstrate that the role of Pol ε in replicative stress sensing is conserved in plants, and provide, to our knowledge, the first genetic dissection of the downstream signaling events in a multicellular eukaryote.

Role of the Polymerase ϵ sub-unit DPB2 in DNA replication, cell cycle regulation and DNA damage response in Arabidopsis.

Faithful DNA replication maintains genome stability in dividing cells and from one generation to the next. This is particularly important in plants because the whole plant body and reproductive cells originate from meristematic cells that retain their proliferative capacity throughout the life cycle of the organism. DNA replication involves large sets of proteins whose activity is strictly regulated, and is tightly linked to the DNA damage response to detect and respond to replication errors or defects. Central to this interconnection is the replicative polymerase DNA Polymerase ϵ (Pol ϵ) which participates in DNA replication per se, as well as replication stress response in animals and in yeast. Surprisingly, its function has to date been little explored in plants, and notably its relationship with DNA Damage Response (DDR) has not been investigated. Here, we have studied the role of the largest regulatory sub-unit of Arabidopsis DNA Pol ϵ: DPB2, using an over-expression strategy. We demonstrate that excess accumulation of the protein impairs DNA replication and causes endogenous DNA stress. Furthermore, we show that Pol ϵ dysfunction has contrasting outcomes in vegetative and reproductive cells and leads to the activation of distinct DDR pathways in the two cell types.

Characterization of CYCLOIDEA-like genes in Proteaceae, a basal eudicot family with multiple shifts in floral symmetry.

BACKGROUND AND AIMS: The basal eudicot family Proteaceae (approx. 1700 species) shows considerable variation in floral symmetry but has received little attention in studies of evolutionary development at the genetic level. A framework for understanding the shifts in floral symmetry in Proteaceae is provided by reconstructing ancestral states on an upated phylogeny of the family, and homologues of CYCLOIDEA (CYC), a key gene for the control of floral symmetry in both monocots and eudicots, are characterized. METHODS: Perianth symmetry transitions were reconstructed on a new species-level tree using parsimony and maximum likelihood. CYC-like genes in 35 species (31 genera) of Proteaceae were sequenced and their phylogeny was reconstructed. Shifts in selection pressure following gene duplication were investigated using nested branch-site models of sequence evolution. Expression patterns of CYC homologues were characterized in three species of Grevillea with different types of floral symmetry. KEY RESULTS: Zygomorphy has evolved 10-18 times independently in Proteaceae from actinomorphic ancestors, with at least four reversals to actinomorphy. A single duplication of CYC-like genes occurred prior to the diversification of Proteaceae, with putative loss or divergence of the ProtCYC1 paralogue in more than half of the species sampled. No shifts in selection pressure were detected in the branches subtending the two ProtCYC paralogues. However, the amino acid sequence preceding the TCP domain is strongly divergent in Grevillea ProtCYC1 compared with other species. ProtCYC genes were expressed in developing flowers of both actinomorphic and zygomorphic Grevillea species, with late asymmetric expression in the perianth of the latter. CONCLUSION: Proteaceae is a remarkable family in terms of the number of transitions in floral symmetry. Furthermore, although CYC-like genes in Grevillea have unusual sequence characteristics, they display patterns of expression that make them good candidates for playing a role in the establishment of floral symmetry.

GEM2Net: from gene expression modeling to -omics networks, a new CATdb module to investigate Arabidopsis thaliana genes involved in stress response.

CATdb (http://urgv.evry.inra.fr/CATdb) is a database providing a public access to a large collection of transcriptomic data, mainly for Arabidopsis but also for other plants. This resource has the rare advantage to contain several thousands of microarray experiments obtained with the same technical protocol and analyzed by the same statistical pipelines. In this paper, we present GEM2Net, a new module of CATdb that takes advantage of this homogeneous dataset to mine co-expression units and decipher Arabidopsis gene functions. GEM2Net explores 387 stress conditions organized into 18 biotic and abiotic stress categories. For each one, a model-based clustering is applied on expression differences to identify clusters of co-expressed genes. To characterize functions associated with these clusters, various resources are analyzed and integrated: Gene Ontology, subcellular localization of proteins, Hormone Families, Transcription Factor Families and a refined stress-related gene list associated to publications. Exploiting protein-protein interactions and transcription factors-targets interactions enables to display gene networks. GEM2Net presents the analysis of the 18 stress categories, in which 17,264 genes are involved and organized within 681 co-expression clusters. The meta-data analyses were stored and organized to compose a dynamic Web resource.

A chemical genetic strategy identify the PHOSTIN, a synthetic molecule that triggers phosphate starvation responses in Arabidopsis thaliana.

Plants display numerous strategies to cope with phosphate (Pi)-deficiency. Despite multiple genetic studies, the molecular mechanisms of low-Pi-signalling remain unknown. To validate the interest of chemical genetics to investigate this pathway we discovered and analysed the effects of PHOSTIN (PSN), a drug mimicking Pi-starvation in Arabidopsis. We assessed the effects of PSN and structural analogues on the induction of Pi-deficiency responses in mutants and wild-type and followed their accumulation in plants organs by high pressure liquid chromotography (HPLC) or mass-spectrophotometry. We show that PSN is cleaved in the growth medium, releasing its active motif (PSN11), which accumulates in plants roots. Despite the overaccumulation of Pi in the roots of treated plants, PSN11 elicits both local and systemic Pi-starvation effects. Nevertheless, albeit that the transcriptional activation of low-Pi genes by PSN11 is lost in the phr1;phl1 double mutant, neither PHO1 nor PHO2 are required for PSN11 effects. The range of local and systemic responses to Pi-starvation elicited, and their dependence on the PHR1/PHL1 function suggests that PSN11 affects an important and early step of Pi-starvation signalling. Its independence from PHO1 and PHO2 suggest the existence of unknown pathway(s), showing the usefulness of PSN and chemical genetics to bring new elements to this field.

Identification of phosphatin, a drug alleviating phosphate starvation responses in Arabidopsis.

Inorganic phosphate (Pi) is present in most soils at suboptimal concentrations, strongly limiting plant development. Plants have the ability to sense and adapt to the surrounding ionic environment, and several genes involved in the response to Pi starvation have been identified. However, a global understanding of the regulatory mechanisms involved in this process is still elusive. Here, we have initiated a chemical genetics approach and isolated compounds that inhibit the response to Pi starvation in Arabidopsis (Arabidopsis thaliana). Molecules were screened for their ability to inhibit the expression of a Pi starvation marker gene (the high-affinity Pi transporter PHT1;4). A drug family named Phosphatin (PTN; Pi starvation inhibitor), whose members act as partial suppressors of Pi starvation responses, was thus identified. PTN addition also reduced various traits of Pi starvation, such as phospholipid/glycolipid conversion, and the accumulation of starch and anthocyanins. A transcriptomic assay revealed a broad impact of PTN on the expression of many genes regulated by low Pi availability. Despite the reduced amount of Pi transporters and resulting reduced Pi uptake capacity, no reduction of Pi content was observed. In addition, PTN improved plant growth; this reveals that the developmental restrictions induced by Pi starvation are not a consequence of metabolic limitation but a result of genetic regulation. This highlights the existence of signal transduction pathway(s) that limit plant development under the Pi starvation condition.

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.

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.

Nucleotide and RNA metabolism prime translational initiation in the earliest events of mitochondrial biogenesis during Arabidopsis germination.

Mitochondria play a crucial role in germination and early seedling growth in Arabidopsis (Arabidopsis thaliana). Morphological observations of mitochondria revealed that mitochondrial numbers, typical size, and oval morphology were evident after 12 h of imbibition in continuous light (following 48 h of stratification). The transition from a dormant to an active metabolic state was punctuated by an early molecular switch, characterized by a transient burst in the expression of genes encoding mitochondrial proteins. Factors involved in mitochondrial transcription and RNA processing were overrepresented among these early-expressed genes. This was closely followed by an increase in the transcript abundance of genes encoding proteins involved in mitochondrial DNA replication and translation. This burst in the expression of factors implicated in mitochondrial RNA and DNA metabolism was accompanied by an increase in transcripts encoding components required for nucleotide biosynthesis in the cytosol and increases in transcript abundance of specific members of the mitochondrial carrier protein family that have previously been associated with nucleotide transport into mitochondria. Only after these genes peaked in expression and largely declined were typical mitochondrial numbers and morphology observed. Subsequently, there was an increase in transcript abundance for various bioenergetic and metabolic functions of mitochondria. The coordination of nucleus- and organelle-encoded gene expression was also examined by quantitative reverse transcription-polymerase chain reaction, specifically for components of the mitochondrial electron transport chain and the chloroplastic photosynthetic machinery. Analysis of protein abundance using western-blot analysis and mass spectrometry revealed that for many proteins, patterns of protein and transcript abundance changes displayed significant positive correlations. A model for mitochondrial biogenesis during germination is proposed, in which an early increase in the abundance of transcripts encoding biogenesis functions (RNA metabolism and import components) precedes a later cascade of gene expression encoding the bioenergetic and metabolic functions of mitochondria.

The cytoskeleton and the peroxisomal-targeted snowy cotyledon3 protein are required for chloroplast development in Arabidopsis.

Here, we describe the snowy cotyledon3 (sco3-1) mutation, which impairs chloroplast and etioplast development in Arabidopsis thaliana seedlings. SCO3 is a member of a largely uncharacterized protein family unique to the plant kingdom. The sco3-1 mutation alters chloroplast morphology and development, reduces chlorophyll accumulation, impairs thylakoid formation and photosynthesis in seedlings, and results in photoinhibition under extreme CO(2) concentrations in mature leaves. There are no readily apparent changes to chloroplast biology, such as transcription or assembly that explain the disruption to chloroplast biogenesis. Indeed, SCO3 is actually targeted to another organelle, specifically to the periphery of peroxisomes. However, impaired chloroplast development cannot be attributed to perturbed peroxisomal metabolic processes involving germination, fatty acid ß-oxidation or photorespiration, though there are so far undescribed changes in low and high CO(2) sensitivity in seedlings and young true leaves. Many of the chloroplasts are bilobed, and some have persistent membranous extensions that encircle other cellular components. Significantly, there are changes to the cytoskeleton in sco3-1, and microtubule inhibitors have similar effects on chloroplast biogenesis as sco3-1 does. The localization of SCO3 to the periphery of the peroxisomes was shown to be dependent on a functional microtubule cytoskeleton. Therefore, the microtubule and peroxisome-associated SCO3 protein is required for chloroplast development, and sco3-1, along with microtubule inhibitors, demonstrates an unexpected role for the cytoskeleton and peroxisomes in chloroplast biogenesis.

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).