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.

Asc-Seurat: analytical single-cell Seurat-based web application

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of transcriptomes, arising as a powerful tool for discovering and characterizing cell types and their developmental trajectories. However, scRNA-seq analysis is complex, requiring a continuous, iterative process to refine the data and uncover relevant biological information. A diversity of tools has been developed to address the multiple aspects of scRNA-seq data analysis. However, an easy-to-use web application capable of conducting all critical steps of scRNA-seq data analysis is still lacking. We present Asc-Seurat, a feature-rich workbench, providing an user-friendly and easy-to-install web application encapsulating tools for an all-encompassing and fluid scRNA-seq data analysis. Asc-Seurat implements functions from the Seurat package for quality control, clustering, and genes differential expression. In addition, Asc-Seurat provides a pseudotime module containing dozens of models for the trajectory inference and a functional annotation module that allows recovering gene annotation and detecting gene ontology enriched terms. We showcase Asc-Seurat's capabilities by analyzing a peripheral blood mononuclear cell dataset.

CONCLUSIONS

Asc-Seurat is a comprehensive workbench providing an accessible graphical interface for scRNA-seq analysis by biologists. Asc-Seurat significantly reduces the time and effort required to analyze and interpret the information in scRNA-seq datasets.

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.

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.

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.