Phenotypic effects of changes in the FTVTxK region of an Arabidopsis secondary wall cellulose synthase compared with results from analogous mutations in other isoforms.

Understanding protein structure and function relationships in cellulose synthase (CesA), including divergent isomers, is an important goal. Here, we report results from mutant complementation assays that tested the ability of sequence variants of AtCesA7, a secondary wall CesA of Arabidopsis thaliana, to rescue the collapsed vessels, short stems, and low cellulose content of the irx3-1 AtCesA7 null mutant. We tested a catalytic null mutation and seven missense or small domain changes in and near the AtCesA7 FTVTSK motif, which lies near the catalytic domain and may, analogously to bacterial CesA, exist within a substrate "gating loop." A low-to-high gradient of rescue occurred, and even inactive AtCesA7 had a small positive effect on stem cellulose content but not stem elongation. Overall, secondary wall cellulose content and stem length were moderately correlated, but the results were consistent with threshold amounts of cellulose supporting particular developmental processes. Vibrational sum frequency generation microscopy allowed tissue-specific analysis of cellulose content in stem xylem and interfascicular fibers, revealing subtle differences between selected genotypes that correlated with the extent of rescue of the collapsing xylem phenotype. Similar tests on PpCesA5 from the moss Physcomitrium (formerly Physcomitrella) patens helped us to synergize the AtCesA7 results with prior results on AtCesA1 and PpCesA5. The cumulative results show that the FTVTxK region is important for the function of an angiosperm secondary wall CesA as well as widely divergent primary wall CesAs, while differences in complementation results between isomers may reflect functional differences that can be explored in further work.

Novel transcriptional responses to heat revealed by turning up the heat at night.

KEY MESSAGE: The circadian clock controls many molecular activities, impacting experimental interpretation. We quantify the genome-wide effects of time-of-day on the heat-shock response and the effects of "diurnal bias" in stress experiments. Heat stress has significant adverse effects on plant productivity worldwide. Most experiments examining heat stress are performed during daytime hours, generating a 'diurnal bias' in the pathways and regulatory mechanisms identified. Such bias may confound downstream interpretations and limit our understanding of the full response to heat stress. Here we show that the transcriptional and physiological responses to a sudden heat shock in Arabidopsis are profoundly sensitive to the time of day. We observe that plant tolerance and acclimation to heat shock vary throughout the day and are maximal at dusk. Consistently, over 75% of heat-responsive transcripts show a time of day-dependent response, including many previously characterized heat-response genes. This temporal sensitivity implies a complex interaction between time and temperature where daily variations in basal transcription influence thermotolerance. When we examined these transcriptional responses, we uncovered novel night-response genes and cis-regulatory elements, underpinning new aspects of heat stress responses not previously appreciated. Exploiting this temporal variation can be applied to most environmental responses to understand the underlying network wiring. Therefore, we propose that using time as a perturbagen is an approach that will enhance our understanding of plant regulatory networks and responses to environmental stresses.

Analysis of differential gene expression and alternative splicing is significantly influenced by choice of reference genome.

RNA-seq analysis has enabled the evaluation of transcriptional changes in many species including nonmodel organisms. However, in most species only a single reference genome is available and RNA-seq reads from highly divergent varieties are typically aligned to this reference. Here, we quantify the impacts of the choice of mapping genome in rice where three high-quality reference genomes are available. We aligned RNA-seq data from a popular productive rice variety to three different reference genomes and found that the identification of differentially expressed genes differed depending on which reference genome was used for mapping. Furthermore, the ability to detect differentially used transcript isoforms was profoundly affected by the choice of reference genome: Only 30% of the differentially used splicing features were detected when reads were mapped to the more commonly used, but more distantly related reference genome. This demonstrated that gene expression and splicing analysis varies considerably depending on the mapping reference genome, and that analysis of individuals that are distantly related to an available reference genome may be improved by acquisition of new genomic reference material. We observed that these differences in transcriptome analysis are, in part, due to the presence of single nucleotide polymorphisms between the sequenced individual and each respective reference genome, as well as annotation differences between the reference genomes that exist even between syntenic orthologs. We conclude that even between two closely related genomes of similar quality, using the reference genome that is most closely related to the species being sampled significantly improves transcriptome analysis.

TGNap1 is required for microtubule-dependent homeostasis of a subpopulation of the plant trans-Golgi network.

Defining convergent and divergent mechanisms underlying the biogenesis and function of endomembrane organelles is fundamentally important in cell biology. In all eukaryotes, the Trans-Golgi Network (TGN) is the hub where the exocytic and endocytic pathways converge. To gain knowledge in the mechanisms underlying TGN biogenesis and function, we characterized TGNap1, a protein encoded by a plant gene of unknown function conserved with metazoans. We demonstrate that TGNap1 is a TGN protein required for the homeostasis of biosynthetic and endocytic traffic pathways. We also show that TGNap1 binds Rab6, YIP4 and microtubules. Finally, we establish that TGNap1 contributes to microtubule-dependent biogenesis, tracking and function of a TGN subset, likely through interaction with Rab6 and YIP4. Our results identify an important trafficking determinant at the plant TGN and reveal an unexpected reliance of post-Golgi traffic homeostasis and organelle biogenesis on microtubules in plants.

Neural Net Classification Combined With Movement Analysis to Evaluate Setaria viridis as a Model System for Time of Day of Anther Appearance.

In many plant species, the time of day at which flowers open to permit pollination is tightly regulated. Proper time of flower opening, or Time of Day of Anther Appearance (TAA), may coordinate flowering opening with pollinator activity or may shift temperature sensitive developmental processes to cooler times of the day. The genetic mechanisms that regulate the timing of this process in cereal crops are unknown. To address this knowledge gap, it is necessary to establish a monocot model system that exhibits variation in TAA. Here, we examine the suitability of Setaria viridis, the model for C4 photosynthesis, for such a role. We developed an imaging system to monitor the temporal regulation of growth, flower opening time, and other physiological characteristics in Setaria. This system enabled us to compare Setaria varieties Ames 32254, Ames 32276, and PI 669942 variation in growth and daily flower opening time. We observed that TAA occurs primarily at night in these three Setaria accessions. However, significant variation between the accessions was observed for both the ratio of flowers that open in the day vs. night and the specific time of day where the rate is maximal. Characterizing this physiological variation is a requisite step toward uncovering the molecular mechanisms regulating TAA. Leveraging the regulation of TAA could provide researchers with a genetic tool to improve crop productivity in new environments.

Cellulose synthase 'class specific regions' are intrinsically disordered and functionally undifferentiated.

Cellulose synthases (CESAs) are glycosyltransferases that catalyze formation of cellulose microfibrils in plant cell walls. Seed plant CESA isoforms cluster in six phylogenetic clades, whose non-interchangeable members play distinct roles within cellulose synthesis complexes (CSCs). A 'class specific region' (CSR), with higher sequence similarity within versus between functional CESA classes, has been suggested to contribute to specific activities or interactions of different isoforms. We investigated CESA isoform specificity in the moss, Physcomitrella patens (Hedw.) B. S. G. to gain evolutionary insights into CESA structure/function relationships. Like seed plants, P. patens has oligomeric rosette-type CSCs, but the PpCESAs diverged independently and form a separate CESA clade. We showed that P. patens has two functionally distinct CESAs classes, based on the ability to complement the gametophore-negative phenotype of a ppcesa5 knockout line. Thus, non-interchangeable CESA classes evolved separately in mosses and seed plants. However, testing of chimeric moss CESA genes for complementation demonstrated that functional class-specificity is not determined by the CSR. Sequence analysis and computational modeling showed that the CSR is intrinsically disordered and contains predicted molecular recognition features, consistent with a possible role in CESA oligomerization and explaining the evolution of class-specific sequences without selection for class-specific function.

Comparative Structural and Computational Analysis Supports Eighteen Cellulose Synthases in the Plant Cellulose Synthesis Complex.

A six-lobed membrane spanning cellulose synthesis complex (CSC) containing multiple cellulose synthase (CESA) glycosyltransferases mediates cellulose microfibril formation. The number of CESAs in the CSC has been debated for decades in light of changing estimates of the diameter of the smallest microfibril formed from the ß-1,4 glucan chains synthesized by one CSC. We obtained more direct evidence through generating improved transmission electron microscopy (TEM) images and image averages of the rosette-type CSC, revealing the frequent triangularity and average cross-sectional area in the plasma membrane of its individual lobes. Trimeric oligomers of two alternative CESA computational models corresponded well with individual lobe geometry. A six-fold assembly of the trimeric computational oligomer had the lowest potential energy per monomer and was consistent with rosette CSC morphology. Negative stain TEM and image averaging showed the triangularity of a recombinant CESA cytosolic domain, consistent with previous modeling of its trimeric nature from small angle scattering (SAXS) data. Six trimeric SAXS models nearly filled the space below an average FF-TEM image of the rosette CSC. In summary, the multifaceted data support a rosette CSC with 18 CESAs that mediates the synthesis of a fundamental microfibril composed of 18 glucan chains.

The valine and lysine residues in the conserved FxVTxK motif are important for the function of phylogenetically distant plant cellulose synthases.

Cellulose synthases (CESAs) synthesize the ß-1,4-glucan chains that coalesce to form cellulose microfibrils in plant cell walls. In addition to a large cytosolic (catalytic) domain, CESAs have eight predicted transmembrane helices (TMHs). However, analogous to the structure of BcsA, a bacterial CESA, predicted TMH5 in CESA may instead be an interfacial helix. This would place the conserved FxVTxK motif in the plant cell cytosol where it could function as a substrate-gating loop as occurs in BcsA. To define the functional importance of the CESA region containing FxVTxK, we tested five parallel mutations in Arabidopsis thaliana CESA1 and Physcomitrella patens CESA5 in complementation assays of the relevant cesa mutants. In both organisms, the substitution of the valine or lysine residues in FxVTxK severely affected CESA function. In Arabidopsis roots, both changes were correlated with lower cellulose anisotropy, as revealed by Pontamine Fast Scarlet. Analysis of hypocotyl inner cell wall layers by atomic force microscopy showed that two altered versions of Atcesa1 could rescue cell wall phenotypes observed in the mutant background line. Overall, the data show that the FxVTxK motif is functionally important in two phylogenetically distant plant CESAs. The results show that Physcomitrella provides an efficient model for assessing the effects of engineered CESA mutations affecting primary cell wall synthesis and that diverse testing systems can lead to nuanced insights into CESA structure-function relationships. Although CESA membrane topology needs to be experimentally determined, the results support the possibility that the FxVTxK region functions similarly in CESA and BcsA.

ER network homeostasis is critical for plant endosome streaming and endocytosis.

Eukaryotic cells internalize cargo at the plasma membrane via endocytosis, a vital process that is accomplished through a complex network of endosomal organelles. In mammalian cells, the ER is in close association with endosomes and regulates their fission. Nonetheless, the physiological role of such interaction on endocytosis is yet unexplored. Here, we probed the existence of ER-endosome association in plant cells and assayed its physiological role in endocytosis. Through live-cell imaging and electron microscopy studies, we established that endosomes are extensively associated with the plant ER, supporting conservation of interaction between heterotypic organelles in evolutionarily distant kingdoms. Furthermore, by analyzing ER-endosome dynamics in genetic backgrounds with defects in ER structure and movement, we also established that the ER network integrity is necessary for homeostasis of the distribution and streaming of various endosome populations as well as for efficient endocytosis. These results support a novel model that endocytosis homeostasis depends on a spatiotemporal control of the endosome dynamics dictated by the ER membrane network.

Computational and genetic evidence that different structural conformations of a non-catalytic region affect the function of plant cellulose synthase.

The ß-1,4-glucan chains comprising cellulose are synthesized by cellulose synthases in the plasma membranes of diverse organisms including bacteria and plants. Understanding structure-function relationships in the plant enzymes involved in cellulose synthesis (CESAs) is important because cellulose is the most abundant component in the plant cell wall, a key renewable biomaterial. Here, we explored the structure and function of the region encompassing transmembrane helices (TMHs) 5 and 6 in CESA using computational and genetic tools. Ab initio computational structure prediction revealed novel bi-modal structural conformations of the region between TMH5 and 6 that may affect CESA function. Here we present our computational findings on this region in three CESAs of Arabidopsis thaliana (AtCESA1, 3, and 6), the Atcesa3(ixr1-2) mutant, and a novel missense mutation in AtCESA1. A newly engineered point mutation in AtCESA1 (Atcesa1(F954L) ) that altered the structural conformation in silico resulted in a protein that was not fully functional in the temperature-sensitive Atcesa1(rsw1-1) mutant at the restrictive temperature. The combination of computational and genetic results provides evidence that the ability of the TMH5-6 region to adopt specific structural conformations is important for CESA function.

Cellulose synthases: new insights from crystallography and modeling.

Detailed information about the structure and biochemical mechanisms of cellulose synthase (CelS) proteins remained elusive until a complex containing the catalytic subunit (BcsA) of CelS from Rhodobacter sphaeroides was crystalized. Additionally, a 3D structure of most of the cytosolic domain of a plant CelS (GhCESA1 from cotton, Gossypium hirsutum) was produced by computational modeling. This predicted structure contributes to our understanding of how plant CelS proteins may be similar and different as compared with BcsA. In this review, we highlight how these structures impact our understanding of the synthesis of cellulose and other extracellular polysaccharides. We show how the structures can be used to generate hypotheses for experiments testing mechanisms of glucan synthesis and translocation in plant CelS.

Membrane-tethered transcription factors provide a connection between stress response and developmental pathways.

Membrane-tethered transcription factors (MTTFs) are proteins that are targeted to membranes and are capable of regulating gene expression. In this way, they are physically restrained from entering the nucleus and are innately dormant. Upon specific signal recognition cues, MTTFs are activated through cleavage by a protease that releases the transcription factor domain into the cytosol thus allowing it to translocate to the nucleus where it can regulate gene expression. MTTFs are classically thought to provide an advantage to an organism by allowing for rapid signal transduction in response to cellular and environmental stresses. However, recent findings suggest that MTTFs may not only act as a means to respond quickly to stress but also are able to regulate developmental pathways, illustrating a point of interaction between stress and development. 

Control of root hair development in Arabidopsis thaliana by an endoplasmic reticulum anchored member of the R2R3-MYB transcription factor family.

The evolution of roots and root hairs was a crucial innovation that contributed to the adaptation of plants to a terrestrial environment. Initiation of root hairs involves transcriptional cues that in part determine cell patterning of the root epidermis. Once root hair initiation has occurred, elongation of the root hair takes place. Although many genes have been identified as being involved in root hair development, many contributors remain uncharacterized. In this study we report on the involvement of a member (here dubbed maMYB) of the plant-specific R2R3-MYB family of transcription factors in root hair elongation in Arabidopsis. We show that maMYB is associated with the endoplasmic reticulum membrane with the transcription factor domain exposed to the cytosol, suggesting that it may function as a membrane-tethered transcription factor. We demonstrate that a truncated form of maMYB (maMYB84⁻³°9), which contains the R2R3-MYB transcription factor domain, is localized and retained in the nucleus, where it regulates gene expression. Silencing of maMyb resulted in plants with significantly shorter root hairs but similar root hair density compared with wild type, implying a role of the protein in root hair elongation. 2,4-D (2,4-dichlorophenoxyacetic acid), an exogenous auxin analog that promotes root hair elongation, rescued the short root hair phenotype and maMyb mRNA was induced in the presence of 2,4-D and IAA (indole-3-acetic acid). These results indicate a functional role of maMYB, which is integrated with auxin, in root hair elongation in Arabidopsis.

The targeting of the oxysterol-binding protein ORP3a to the endoplasmic reticulum relies on the plant VAP33 homolog PVA12.

In plants, sterols play fundamental roles as membrane constituents in the biosynthesis of steroid hormones, and act as precursors for cell wall deposition. Sterols are synthesized in the endoplasmic reticulum (ER), but mainly accumulate in the plasma membrane. How sterols are trafficked in plant cells is largely unknown. In non-plant systems, oxysterol-binding proteins have been involved in sterol trafficking and homeostasis. There are at least twelve homologs of oxysterol-binding proteins in the Arabidopsis genome, but the biology of these proteins remains for the most part obscure. Here, we report our analysis of the targeting requirements and the sterol-binding properties of a small Arabidopsis oxysterol-binding protein, ORP3a. We have determined that ORP3a is a bona fide sterol-binding protein with sitosterol-binding properties. Live-cell imaging analyses revealed that ORP3a is localized at the ER, and that binding to this organelle depends on a direct interaction with PVA12, a member of the largely uncharacterized VAP33 family of plant proteins. Molecular modeling analyses and site-directed mutagenesis led to the identification of a novel protein domain that is responsible for the PVA12-ORP3a interaction. Disruption of the integrity of this domain caused redistribution of ORP3a to the Golgi apparatus, suggesting that ORP3a may cycle between the ER and the Golgi. These results represent new insights into the biology of sterol-binding proteins in plant cells, and elucidate a hitherto unknown relationship between members of oxysterol-binding protein and VAP33 families of plant proteins in the early plant secretory pathway.

Membrane-tethered transcription factors in Arabidopsis thaliana: novel regulators in stress response and development.

Membrane-tethered transcription factors (MTTFs) differ from cytosolic transcription factors (TF) in that they are innately membrane-bound. To attain TF activity, MTTFs are released from the membrane anchor as a result of proteolytic cleavage. This enables MTTFs to travel to the nucleus and modulate gene expression. Arabidopsis MTTFs characterized to date belong to either the bZIP or the NAC family. In this review, we highlight the most recent findings on Arabidopsis MTTFs that ascribe different yet important roles to these proteins: the MTTFs in the bZIP family appear to regulate stress signaling pathways, whereas members of the NAC family are involved in both development and stress response.