The single-cell transcriptome program of nodule development cellular lineages in Medicago truncatula.

Legumes establish a symbiotic relationship with nitrogen-fixing rhizobia by developing nodules. Nodules are modified lateral roots that undergo changes in their cellular development in response to bacteria, but the transcriptional reprogramming that occurs in these root cells remains largely uncharacterized. Here, we describe the cell-type-specific transcriptome response of Medicago truncatula roots to rhizobia during early nodule development in the wild-type genotype Jemalong A17, complemented with a hypernodulating mutant (sunn-4) to expand the cell population responding to infection and subsequent biological inferences. The analysis identifies epidermal root hair and stele sub-cell types associated with a symbiotic response to infection and regulation of nodule proliferation. Trajectory inference shows cortex-derived cell lineages differentiating to form the nodule primordia and, posteriorly, its meristem, while modulating the regulation of phytohormone-related genes. Gene regulatory analysis of the cell transcriptomes identifies new regulators of nodulation, including STYLISH 4, for which the function is validated.

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula.

BACKGROUND: Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. RESULTS: We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. CONCLUSIONS: Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus.

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

Functional and comparative genomics reveals conserved noncoding sequences in the nitrogen-fixing clade.

Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen-fixing bacteria in the root nodule. This root-nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade. To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen-fixing clade (NFC). We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1. CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Spatiotemporal cytokinin response imaging and ISOPENTENYLTRANSFERASE 3 function in Medicago nodule development.

Most legumes can establish a symbiotic association with soil rhizobia that trigger the development of root nodules. These nodules host the rhizobia and allow them to fix nitrogen efficiently. The perception of bacterial lipo-chitooligosaccharides (LCOs) in the epidermis initiates a signaling cascade that allows rhizobial intracellular infection in the root and de-differentiation and activation of cell division that gives rise to the nodule. Thus, nodule organogenesis and rhizobial infection need to be coupled in space and time for successful nodulation. The plant hormone cytokinin (CK) contributes to the coordination of this process, acting as an essential positive regulator of nodule organogenesis. However, the temporal regulation of tissue-specific CK signaling and biosynthesis in response to LCOs or Sinorhizobium meliloti inoculation in Medicago truncatula remains poorly understood. In this study, using a fluorescence-based CK sensor (pTCSn::nls:tGFP), we performed a high-resolution tissue-specific temporal characterization of the sequential activation of CK response during root infection and nodule development in M. truncatula after inoculation with S. meliloti. Loss-of-function mutants of the CK-biosynthetic gene ISOPENTENYLTRANSFERASE 3 (IPT3) showed impairment of nodulation, suggesting that IPT3 is required for nodule development in M. truncatula. Simultaneous live imaging of pIPT3::nls:tdTOMATO and the CK sensor showed that IPT3 induction in the pericycle at the base of nodule primordium contributes to CK biosynthesis, which in turn promotes expression of positive regulators of nodule organogenesis in M. truncatula.

A robust method of nuclei isolation for single-cell RNA sequencing of solid tissues from the plant genus Populus.

Single-cell transcriptome analysis has been extensively applied in humans and animal models to uncover gene expression heterogeneity between the different cell types of a tissue or an organ. It demonstrated its capability to discover key regulatory elements that determine cell fate during developmental programs. Single-cell analysis requires the isolation and labeling of the messenger RNA (mRNA) derived from each cell. These challenges were primarily addressed in mammals by developing microfluidic-based approaches. For plant species whose cells contain cell walls, these approaches have generally required the generation of isolated protoplasts. Many plant tissues' secondary cell wall hinders enzymatic digestion required for individual protoplast isolation, resulting in an unequal representation of cell types in a protoplast population. This limitation is especially critical for cell types located in the inner layers of a tissue or the inner tissues of an organ. Consequently, single-cell RNA sequencing (scRNA-seq) studies using microfluidic approaches in plants have mainly been restricted to Arabidopsis roots, for which well-established procedures of protoplast isolation are available. Here we present a simple alternative approach to generating high-quality protoplasts from plant tissue by characterizing the mRNA extracted from individual nuclei instead of whole cells. We developed the protocol using two different plant materials with varying cellular complexity levels and cell wall structure, Populus shoot apices, and more lignified stems. Using the 10× Genomics Chromium technology, we show that this procedure results in intact mRNA isolation and limited leakage, with a broad representation of individual cell transcriptomes.

Guard cell redox proteomics reveals a role of lipid transfer protein in plant defense.

Redox-based post-translational modifications (PTMs) involving protein cysteine residues as redox sensors are important to various physiological processes. However, little is known about redox-sensitive proteins in guard cells and their functions in stomatal immunity. In this study, we applied an integrative protein labeling method cysTMTRAQ, and identified guard cell proteins that were altered by thiol redox PTMs in response to a bacterial flagellin peptide flg22. In total, eight, seven and 20 potential redox-responsive proteins were identified in guard cells treated with flg22 for 15, 30 and 60 min, respectively. The proteins fall into several functional groups including photosynthesis, lipid binding, oxidation-reduction, and defense. Among the proteins, a lipid transfer protein (LTP)-II was confirmed to be redox-responsive and involved in plant resistance to Pseudomonas syringe pv. tomato DC3000. This study not only creates an inventory of potential redox-sensitive proteins in flg22 signal transduction in guard cells, but also highlights the biological relevance of the lipid transfer protein in plant defense against bacterial pathogens. SIGNIFICANCE: Protein redox modifications play important roles in many physiological processes. However, redox proteomics has rarely been studied in plant single cell-types. In this study, isobaric tandem mass tag-based redox proteomics technology was applied to discover redox-sensitive proteins and corresponding cysteine residues in guard cell response to a bacterial flagellin peptide flg22. Many redox-responsive proteins related to photosynthesis, lipid binding, oxidation-reduction, and defense were identified. Using reverse genetics and biochemical analyses, a lipid transfer protein was functionally characterized to be involved in plant defense against pathogens. The study highlights the utility of redox proteomics in discovering new proteins and redox modifications in important stomatal guard cell functions. Furthermore, detailed functional characterization demonstrates the biological relevance of the redox-responsive lipid transfer protein in plant pathogen defense.

Accelerating forest tree breeding by integrating genomic selection and greenhouse phenotyping.

Breeding forest species can be a costly and slow process because of the extensive areas needed for field trials and the long periods (e.g., five years) that are required to measure economically and environmentally relevant phenotypes (e.g., adult plant biomass or plant height). Genomic selection (GS) and indirect selection using early phenotypes (e.g., phenotypes collected in greenhouse conditions) are two ways by which tree breeding can be accelerated. These approaches can both reduce the costs of field-testing and the time required to make selection decisions. Moreover, these approaches can be highly synergistic. Therefore, in this study, we used a data set comprising DNA genotypes and longitudinal measurements of growth collected from a population of Populus deltoides W. Bartram ex Marshall (eastern cottonwood) in the greenhouse and the field, to evaluate the potential impact of integrating large-scale greenhouse phenotyping with conventional GS. We found that the integration of greenhouse phenotyping and GS can deliver very early selection decisions that are moderately accurate. Therefore, we conclude that the adoption of these approaches, in conjunction with reproductive techniques that shorten the generation interval, can lead to an unprecedented acceleration of selection gains in P. deltoides and, potentially, other commercially planted tree species.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides.

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

The Genetic Regulation of Alternative Splicing in Populus deltoides.

Alternative splicing (AS) is a mechanism of regulation of the proteome via enabling the production of multiple mRNAs from a single gene. To date, the dynamics of AS and its effects on the protein sequences of individuals in a large and genetically unrelated population of trees have not been investigated. Here we describe the diversity of AS events within a previously genotyped population of 268 individuals of Populus deltoides and their putative downstream functional effects. Using a robust bioinformatics pipeline, the AS events and resulting transcript isoforms were discovered and quantified for each individual in the population. Analysis of the AS revealed that, as expected, most AS isoforms are conserved. However, we also identified a substantial collection of new, unannotated splice junctions and transcript isoforms. Heritability estimates for the expression of transcript isoforms showed that approximately half of the isoforms are heritable. The genetic regulators of these AS isoforms and splice junction usage were then identified using a genome-wide association analysis. The expression of AS isoforms was predominately cis regulated while splice junction usage was generally regulated in trans. Additionally, we identified 696 genes encoding alternatively spliced isoforms that changed putative protein domains relative to the longest protein coding isoform of the gene, and 859 genes exhibiting this same phenomenon relative to the most highly expressed isoform. Finally, we found that 748 genes gained or lost micro-RNA binding sites relative to the longest protein coding isoform of a given gene, while 940 gained or lost micro-RNA binding sites relative to the most highly expressed isoform. These results indicate that a significant fraction of AS events are genetically regulated and that this isoform usage can result in protein domain architecture changes.

The uncharacterized gene EVE contributes to vessel element dimensions in Populus.

The radiation of angiosperms led to the emergence of the vast majority of today's plant species and all our major food crops. Their extraordinary diversification occurred in conjunction with the evolution of a more efficient vascular system for the transport of water, composed of vessel elements. The physical dimensions of these water-conducting specialized cells have played a critical role in angiosperm evolution; they determine resistance to water flow, influence photosynthesis rate, and contribute to plant stature. However, the genetic factors that determine their dimensions are unclear. Here we show that a previously uncharacterized gene, ENLARGED VESSEL ELEMENT (EVE), contributes to the dimensions of vessel elements in Populus, impacting hydraulic conductivity. Our data suggest that EVE is localized in the plasma membrane and is involved in potassium uptake of differentiating xylem cells during vessel development. In plants, EVE first emerged in streptophyte algae, but expanded dramatically among vessel-containing angiosperms. The phylogeny, structure and composition of EVE indicates that it may have been involved in an ancient horizontal gene-transfer event.

Bicarbonate Induced Redox Proteome Changes in Arabidopsis Suspension Cells.

Climate change as a result of increasing atmospheric CO2 affects plant growth and productivity. CO2 is not only a carbon donor for photosynthesis but also an environmental signal that can perturb cellular redox homeostasis and lead to modifications of redox-sensitive proteins. Although redox regulation of protein functions has emerged as an important mechanism in several biological processes, protein redox modifications and how they function in plant CO2 response remain unclear. Here a new iodoTMTRAQ proteomics technology was employed to analyze changes in protein redox modifications in Arabidopsis thaliana suspension cells in response to bicarbonate (mimic of elevated CO2) in a time-course study. A total of 47 potential redox-regulated proteins were identified with functions in carbohydrate and energy metabolism, transport, ROS scavenging, cell structure modulation and protein turnover. This inventory of previously unknown redox responsive proteins in Arabidopsis bicarbonate responses lays a foundation for future research toward understanding the molecular mechanisms underlying plant CO2 responses.

A family of metal-dependent phosphatases implicated in metabolite damage-control.

DUF89 family proteins occur widely in both prokaryotes and eukaryotes, but their functions are unknown. Here we define three DUF89 subfamilies (I, II, and III), with subfamily II being split into stand-alone proteins and proteins fused to pantothenate kinase (PanK). We demonstrated that DUF89 proteins have metal-dependent phosphatase activity against reactive phosphoesters or their damaged forms, notably sugar phosphates (subfamilies II and III), phosphopantetheine and its S-sulfonate or sulfonate (subfamily II-PanK fusions), and nucleotides (subfamily I). Genetic and comparative genomic data strongly associated DUF89 genes with phosphoester metabolism. The crystal structure of the yeast (Saccharomyces cerevisiae) subfamily III protein YMR027W revealed a novel phosphatase active site with fructose 6-phosphate and Mg(2+) bound near conserved signature residues Asp254 and Asn255 that are critical for activity. These findings indicate that DUF89 proteins are previously unrecognized hydrolases whose characteristic in vivo function is to limit potentially harmful buildups of normal or damaged phosphometabolites.

Protein Phosphorylation and Redox Modification in Stomatal Guard Cells.

Post-translational modification (PTM) is recognized as a major process accounting for protein structural variation, functional diversity, and the dynamics and complexity of the proteome. Since PTMs can change the structure and function of proteins, they are essential to coordinate signaling networks and to regulate important physiological processes in eukaryotes. Plants are constantly challenged by both biotic and abiotic stresses that reduce productivity, causing economic losses in crops. The plant responses involve complex physiological, cellular, and molecular processes, with stomatal movement as one of the earliest responses. In order to activate such a rapid response, stomatal guard cells employ cellular PTMs of key protein players in the signaling pathways to regulate the opening and closure of the stomatal pores. Here we discuss two major types of PTMs, protein phosphorylation and redox modification that play essential roles in stomatal movement under stress conditions. We present an overview of PTMs that occur in stomatal guard cells, especially the methods and technologies, and their applications in PTM identification and quantification. Our focus is on PTMs that modify molecular components in guard cell signaling at the stages of signal perception, second messenger production, as well as downstream signaling events and output. Improved understanding of guard cell signaling will enable generation of crops with enhanced stress tolerance, and increased yield and bioenergy through biotechnology and molecular breeding.

New nodes and edges in the glucosinolate molecular network revealed by proteomics and metabolomics of Arabidopsis myb28/29 and cyp79B2/B3 glucosinolate mutants.

Glucosinolates present in Brassicales are important for human health and plant defense against insects and pathogens. Here we investigate the proteomes and metabolomes of Arabidopsis myb28/29 and cyp79B2/B3 mutants deficient in aliphatic glucosinolates and indolic glucosinolates, respectively. Quantitative proteomics of the myb28/29 and cyp79B2/B3 mutants led to the identification of 2785 proteins, of which 142 proteins showed significant changes in the two mutants compared to wild type (WT). By mapping the differential proteins using STRING, we detected 59 new edges in the glucosinolate metabolic network. These connections can be classified as primary with direct roles in glucosinolate metabolism, secondary related to plant stress responses, and tertiary involved in other biological processes. Gene Ontology analysis of the differential proteins showed high level of enrichment in the nodes belonging to metabolic process including glucosinolate biosynthesis and response to stimulus. Using metabolomics, we quantified 292 metabolites covering a broad spectrum of metabolic pathways, and 89 exhibited differential accumulation patterns between the mutants and WT. The changing metabolites (e.g., γ-glutamyl amino acids, auxins and glucosinolate hydrolysis products) complement our proteomics findings. This study contributes toward engineering and breeding of glucosinolate profiles in plants in efforts to improve human health, crop quality and productivity. BIOLOGICAL SIGNIFICANCE: Glucosinolates in Brassicales constitute an important group of natural metabolites important for plant defense and human health. Its biosynthetic pathways and transcriptional regulation have been well-studied. Using Arabidopsis mutants of important genes in glucosinolate biosynthesis, quantitative proteomics and metabolomics led to identification of many proteins and metabolites that are potentially related to glucosinolate metabolism. This study provides a comprehensive insight into the molecular networks of glucosinolate metabolism, and will facilitate efforts toward engineering and breeding of glucosinolate profiles for enhanced crop defense, and nutritional value.

Preparation of Epidermal Peels and Guard Cell Protoplasts for Cellular, Electrophysiological, and -Omics Assays of Guard Cell Function.

Bioassays are commonly used to study stomatal phenotypes. There are multiple options in the choice of plant materials and species used for observation of stomatal and guard cell responses in vivo. Here, detailed procedures for bioassays of stomatal responses to abscisic acid (ABA) in Arabidopsis thaliana are described, including ABA promotion of stomatal closure, ABA inhibition of stomatal opening, and ABA promotion of reaction oxygen species (ROS) production in guard cells. We also include an example of a stomatal bioassay for the guard cell CO2 response using guard cell-enriched epidermal peels from Brassica napus. Highly pure preparations of guard cell protoplasts can be produced, which are also suitable for studies on guard cell signaling, as well as for studies on guard cell ion transport. Small-scale and large-scale guard cell protoplast preparations are commonly used for electrophysiological and -omics studies, respectively. We provide a procedure for small-scale guard cell protoplasting from A. thaliana. Additionally, a general protocol for large-scale preparation of guard cell protoplasts, with specifications for three different species, A. thaliana, B. napus, and Vicia faba is also provided.

Redox proteomics of tomato in response to Pseudomonas syringae infection.

Unlike mammals with adaptive immunity, plants rely on their innate immunity based on pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) for pathogen defense. Reactive oxygen species, known to play crucial roles in PTI and ETI, can perturb cellular redox homeostasis and lead to changes of redox-sensitive proteins through modification of cysteine sulfhydryl groups. Although redox regulation of protein functions has emerged as an important mechanism in several biological processes, little is known about redox proteins and how they function in PTI and ETI. In this study, cysTMT proteomics technology was used to identify similarities and differences of protein redox modifications in tomato resistant (PtoR) and susceptible (prf3) genotypes in response to Pseudomonas syringae pv tomato (Pst) infection. In addition, the results of the redox changes were compared and corrected with the protein level changes. A total of 90 potential redox-regulated proteins were identified with functions in carbohydrate and energy metabolism, biosynthesis of cysteine, sucrose and brassinosteroid, cell wall biogenesis, polysaccharide/starch biosynthesis, cuticle development, lipid metabolism, proteolysis, tricarboxylic acid cycle, protein targeting to vacuole, and oxidation-reduction. This inventory of previously unknown protein redox switches in tomato pathogen defense lays a foundation for future research toward understanding the biological significance of protein redox modifications in plant defense responses.

cysTMTRAQ-An integrative method for unbiased thiol-based redox proteomics.

Protein redox regulation plays important roles in many biological processes. Protein cysteine thiols are sensitive to redox changes and may function as redox switches, which turn signaling and metabolic pathways on or off to ensure speedy responses to environmental stimuli or stresses. Here we report a novel integrative proteomics method called cysTMTRAQ that combines two types of isobaric tags, cysteine tandem mass tags and isobaric tag for relative and absolute quantification, in one experiment. The method not only enables simultaneous analysis of cysteine redox changes and total protein level changes, but also allows the determination of bona fide redox modified cysteines in proteins through the correction of protein turnover. Overlooking the factor of protein-level changes in the course of protein posttranslational modification experiments could lead to misleading results. The capability to analyze protein posttranslational modification dynamics and protein-level changes in one experiment will advance proteomic studies in many areas of biology and medicine.