Inferring Hypotheses on Functional Relationships of Genes: Analysis of the Arabidopsis thaliana Subtilase Gene Family

The gene family of subtilisin-like serine proteases (subtilases) in Arabidopsis thaliana comprises 56 members, divided into six distinct subfamilies. Whereas the members of five subfamilies are similar to pyrolysins, two genes share stronger similarity to animal kexins. Mutant screens confirmed 144 T-DNA insertion lines with knockouts for 55 out of the 56 subtilases. Apart from SDD1, none of the confirmed homozygous mutants revealed any obvious visible phenotypic alteration during growth under standard conditions. Apart from this specific case, forward genetics gave us no hints about the function of the individual 54 non-characterized subtilase genes. Therefore, the main objective of our work was to overcome the shortcomings of the forward genetic approach and to infer alternative experimental approaches by using an integrative bioinformatics and biological approach. Computational analyses based on transcriptional co-expression and co-response pattern revealed at least two expression networks, suggesting that functional redundancy may exist among subtilases with limited similarity. Furthermore, two hubs were identified, which may be involved in signalling or may represent higher-order regulatory factors involved in responses to environmental cues. A particular enrichment of co-regulated genes with metabolic functions was observed for four subtilases possibly representing late responsive elements of environmental stress. The kexin homologs show stronger associations with genes of transcriptional regulation context. Based on the analyses presented here and in accordance with previously characterized subtilases, we propose three main functions of subtilases: involvement in (i) control of development, (ii) protein turnover, and (iii) action as downstream components of signalling cascades. Supplemental material is available in the Plant Subtilase Database (PSDB) (http://csbdb.mpimp-golm.mpg.de/psdb.html), as well as from the CSB.DB (http://csbdb.mpimp-golm.mpg.de)

GO-PROMTO Illuminates Protein Membrane Topologies of Glycan Biosynthetic Enzymes in the Golgi Apparatus of Living Tissues

The Golgi apparatus is the main site of glycan biosynthesis in eukaryotes. Better understanding of the membrane topology of the proteins and enzymes involved can impart new mechanistic insights into these processes. Publically available bioinformatic tools provide highly variable predictions of membrane topologies for given proteins. Therefore we devised a non-invasive experimental method by which the membrane topologies of Golgi-resident proteins can be determined in the Golgi apparatus in living tissues. A Golgi marker was used to construct a series of reporters based on the principle of bimolecular fluorescence complementation. The reporters and proteins of interest were recombinantly fused to split halves of yellow fluorescent protein (YFP) and transiently co-expressed with the reporters in the Nicotiana benthamiana leaf tissue. Output signals were binary, showing either the presence or absence of fluorescence with signal morphologies characteristic of the Golgi apparatus and endoplasmic reticulum (ER). The method allows prompt and robust determinations of membrane topologies of Golgi-resident proteins and is termed GO-PROMTO (for GOlgi PROtein Membrane TOpology). We applied GO-PROMTO to examine the topologies of proteins involved in the biosynthesis of plant cell wall polysaccharides including xyloglucan and arabinan. The results suggest the existence of novel biosynthetic mechanisms involving transports of intermediates across Golgi membranes

A Genome Wide Association Study of arabinoxylan content in 2-row spring barley grain

In barley endosperm arabinoxylan (AX) is the second most abundant cell wall polysaccharide and in wheat it is the most abundant polysaccharide in the starchy endosperm walls of the grain. AX is one of the main contributors to grain dietary fibre content providing several health benefits including cholesterol and glucose lowering effects, and antioxidant activities. Due to its complex structural features, AX might also affect the downstream applications of barley grain in malting and brewing. Using a high pressure liquid chromatography (HPLC) method we quantified AX amounts in mature grain in 128 spring 2-row barley accessions. Amounts ranged from ~ 5.2 μg/g to ~ 9 μg/g. We used this data for a Genome Wide Association Study (GWAS) that revealed three significant quantitative trait loci (QTL) associated with grain AX levels which passed a false discovery threshold (FDR) and are located on two of the seven barley chromosomes. Regions underlying the QTLs were scanned for genes likely to be involved in AX biosynthesis or turnover, and strong candidates, including glycosyltransferases from the GT43 and GT61 families and glycoside hydrolases from the GH10 family, were identified. Phylogenetic trees of selected gene families were built based on protein translations and were used to examine the relationship of the barley candidate genes to those in other species. Our data reaffirms the roles of existing genes thought to contribute to AX content, and identifies novel QTL (and candidate genes associated with them) potentially influencing the AX content of barley grain. One potential outcome of this work is the deployment of highly associated single nucleotide polymorphisms markers in breeding programs to guide the modification of AX abundance in barley grain

Evolution of Plant Nucleotide-Sugar Interconversion Enzymes

Nucleotide-diphospho-sugars (NDP-sugars) are the building blocks of diverse polysaccharides and glycoconjugates in all organisms. In plants, 11 families of NDP-sugar interconversion enzymes (NSEs) have been identified, each of which interconverts one NDP-sugar to another. While the functions of these enzyme families have been characterized in various plants, very little is known about their evolution and origin. Our phylogenetic analyses indicate that all the 11 plant NSE families are distantly related and most of them originated from different progenitor genes, which have already diverged in ancient prokaryotes. For instance, all NSE families are found in the lower land plant mosses and most of them are also found in aquatic algae, implicating that they have already evolved to be capable of synthesizing all the 11 different NDP-sugars. Particularly interesting is that the evolution of RHM (UDP-L-rhamnose synthase) manifests the fusion of genes of three enzymatic activities in early eukaryotes in a rather intriguing manner. The plant NRS/ER (nucleotide-rhamnose synthase/epimerase-reductase), on the other hand, evolved much later from the ancient plant RHMs through losing the N-terminal domain. Based on these findings, an evolutionary model is proposed to explain the origin and evolution of different NSE families. For instance, the UGlcAE (UDP-D-glucuronic acid 4-epimerase) family is suggested to have evolved from some chlamydial bacteria. Our data also show considerably higher sequence diversity among NSE-like genes in modern prokaryotes, consistent with the higher sugar diversity found in prokaryotes. All the NSE families are widely found in plants and algae containing carbohydrate-rich cell walls, while sporadically found in animals, fungi and other eukaryotes, which do not have or have cell walls with distinct compositions. Results of this study were shown to be highly useful for identifying unknown genes for further experimental characterization to determine their functions in the synthesis of diverse glycosylated molecules

Response of cell wall composition and RNA-seq transcriptome to methyl-jasmonate in Brachypodium distachyon callus

Main conclusion: Methyl-jasmonate induces large increases in p-coumarate linked to arabinoxylan in Brachypodium and in abundance of GT61 and BAHD family transcripts consistent with a role in synthesis of this linkage. Jasmonic acid (JA) signalling is required for many stress responses in plants, inducing large changes in the transcriptome, including up-regulation of transcripts associated with lignification. However, less is known about the response to JA of grass cell walls and the monocot-specific features of arabinoxylan (AX) synthesis and acylation by ferulic acid (FA) and para-coumaric acid (pCA). Here, we show that methyl-jasmonate (MeJA) induces moderate increases in FA monomer, > 50% increases in FA dimers, and five–sixfold increases in pCA ester-linked to cell walls in Brachypodium callus. Direct measurement of arabinose acylated by pCA (Araf-pCA) indicated that most or all the increase in cell-wall pCA was due to pCA ester-linked to AX. Analysis of the RNA-seq transcriptome of the callus response showed that these cell-wall changes were accompanied by up-regulation of members of the GT61 and BAHD gene families implicated in AX decoration and acylation; two BAHD paralogues were among the most up-regulated cell-wall genes (seven and fivefold) after 24 h exposure to MeJA. Similar responses to JA of orthologous BAHD and GT61 transcripts are present in the RiceXPro public expression data set for rice seedlings, showing that they are not specific to Brachypodium or to callus. The large response of AX-pCA to MeJA may, therefore, indicate an important role for this linkage in response of primary cell walls of grasses to JA signalling

Absolute Quantitation of In Vitro Expressed Plant Membrane Proteins by Targeted Proteomics (MRM) for the Determination of Kinetic Parameters.

The purification of a functional soluble protein from biological or in vitro expression systems can be problematic and the enrichment of a functional membrane protein for biochemical analyses can be a serious technical challenge. Recently we have been characterizing plant endomembrane nucleotide sugar transporters using a yeast expression system. However, rather than enriching these in vitro expressed proteins to homogeneity, we have been conducting biochemical characterization of these transport proteins in yeast microsomal fractions. While this approach has enabled us to estimate a variety of kinetic parameters, the accurate determination of the turnover number of an enzyme-substrate complex (k cat) requires that the catalytic site concentration (amount of protein) in the total reaction volume is known. As a result, we have been employing targeted proteomics (multiple reaction monitoring) with peptide standards and a triple quadrupole mass spectrometer to estimate the absolute amount of protein in a mixed protein microsomal fraction. The following method details the steps required to define the absolute quantitation of an in vitro expressed membrane protein to define complete kinetic parameters