Conserved alternative splicing of Arabidopsis transthyretin-like determines protein localization and S-allantoin synthesis in peroxisomes

Non-canonical peroxisome targeting signals: identification of novel PTS1 tripeptides and characterization of enhancer elements by computational permutation analysis

BACKGROUND

High-accuracy prediction tools are essential in the post-genomic era to define organellar proteomes in their full complexity. We recently applied a discriminative machine learning approach to predict plant proteins carrying peroxisome targeting signals (PTS) type 1 from genome sequences. For Arabidopsis thaliana 392 gene models were predicted to be peroxisome-targeted. The predictions were extensively tested in vivo, resulting in a high experimental verification rate of Arabidopsis proteins previously not known to be peroxisomal.

RESULTS

In this study, we experimentally validated the predictions in greater depth by focusing on the most challenging Arabidopsis proteins with unknown non-canonical PTS1 tripeptides and prediction scores close to the threshold. By in vivo subcellular targeting analysis, three novel PTS1 tripeptides (QRL>, SQM>, and SDL>) and two novel tripeptide residues (Q at position -3 and D at pos. -2) were identified. To understand why, among many Arabidopsis proteins carrying the same C-terminal tripeptides, these proteins were specifically predicted as peroxisomal, the residues upstream of the PTS1 tripeptide were computationally permuted and the changes in prediction scores were analyzed. The newly identified Arabidopsis proteins were found to contain four to five amino acid residues of high predicted targeting enhancing properties at position -4 to -12 in front of the non-canonical PTS1 tripeptide. The identity of the predicted targeting enhancing residues was unexpectedly diverse, comprising besides basic residues also proline, hydroxylated (Ser, Thr), hydrophobic (Ala, Val), and even acidic residues.

CONCLUSIONS

Our computational and experimental analyses demonstrate that the plant PTS1 tripeptide motif is more diverse than previously thought, including an increasing number of non-canonical sequences and allowed residues. Specific targeting enhancing elements can be predicted for particular sequences of interest and are far more diverse in amino acid composition and positioning than previously assumed. Machine learning methods become indispensable to predict which specific proteins, among numerous candidate proteins carrying the same non-canonical PTS1 tripeptide, contain sufficient enhancer elements in terms of number, positioning and total strength to cause peroxisome targeting.

Plant peroxisomes: biogenesis and function

Peroxisomes are eukaryotic organelles that are highly dynamic both in morphology and metabolism. Plant peroxisomes are involved in numerous processes, including primary and secondary metabolism, development, and responses to abiotic and biotic stresses. Considerable progress has been made in the identification of factors involved in peroxisomal biogenesis, revealing mechanisms that are both shared with and diverged from non-plant systems. Furthermore, recent advances have begun to reveal an unexpectedly large plant peroxisomal proteome and have increased our understanding of metabolic pathways in peroxisomes. Coordination of the biosynthesis, import, biochemical activity, and degradation of peroxisomal proteins allows for highly dynamic responses of peroxisomal metabolism to meet the needs of a plant. Knowledge gained from plant peroxisomal research will be instrumental to fully understanding the organelle's dynamic behavior and defining peroxisomal metabolic networks, thus allowing the development of molecular strategies for rational engineering of plant metabolism, biomass production, stress tolerance, and pathogen defense.

Structural and functional insights into (S)-ureidoglycine aminohydrolase, key enzyme of purine catabolism in Arabidopsis thaliana

The ureide pathway has recently been identified as the metabolic route of purine catabolism in plants and some bacteria. In this pathway, uric acid, which is a major product of the early stage of purine catabolism, is degraded into glyoxylate and ammonia via stepwise reactions of seven different enzymes. Therefore, the pathway has a possible physiological role in mobilization of purine ring nitrogen for further assimilation. (S)-Ureidoglycine aminohydrolase enzyme converts (S)-ureidoglycine into (S)-ureidoglycolate and ammonia, providing the final substrate to the pathway. Here, we report a structural and functional analysis of this enzyme from Arabidopsis thaliana (AtUGlyAH). The crystal structure of AtUGlyAH in the ligand-free form shows a monomer structure in the bicupin fold of the β-barrel and an octameric functional unit as well as a Mn(2+) ion binding site. The structure of AtUGlyAH in complex with (S)-ureidoglycine revealed that the Mn(2+) ion acts as a molecular anchor to bind (S)-ureidoglycine, and its binding mode dictates the enantioselectivity of the reaction. Further kinetic analysis characterized the functional roles of the active site residues, including the Mn(2+) ion binding site and residues in the vicinity of (S)-ureidoglycine. These analyses provide molecular insights into the structure of the enzyme and its possible catalytic mechanism.

Mining Arabidopsis thaliana RNA-seq data with Integrated Genome Browser reveals stress-induced alternative splicing of the putative splicing regulator SR45a

PREMISE OF THE STUDY

High-throughput sequencing of cDNA libraries prepared from diverse samples (RNA-seq) can reveal genome-wide changes in alternative splicing. Using RNA-seq data to assess splicing at the level of individual genes requires the ability to visualize read alignments alongside genomic annotations. To meet this need, we added RNA-seq visualization capability to Integrated Genome Browser (IGB), a free desktop genome visualization tool. To illustrate this capability, we present an in-depth analysis of abiotic stresses and their effects on alternative splicing of SR45a (AT1G07350), a putative splicing regulator from Arabidopsis thaliana.

METHODS

cDNA libraries prepared from Arabidopsis plants that were subjected to heat and dehydration stresses were sequenced on an Illumina GAIIx sequencer, yielding more than 511 million high-quality 75-base, single-end sequence reads. Reads were aligned onto the reference genome and visualized in IGB.

KEY RESULTS

Using IGB, we confirmed exon-skipping alternative splicing in SR45a. Exon-skipped variant AT1G07350.1 encodes full-length SR45a protein with intact RS and RNA recognition motifs, while nonskipped variant AT1G07350.2 lacks the C-terminal RS region due to a frameshift in the alternative exon. Heat and drought stresses increased both transcript abundance and the proportion of exon-skipped transcripts encoding the full-length protein. We identified new splice sites and observed frequent intron retention flanking the alternative exon.

CONCLUSIONS

This study underlines the importance of visual inspection of RNA-seq alignments when investigating alternatively spliced genes. We showed that heat and dehydration stresses increase overall abundance of SR45a mRNA while also increasing production of transcripts encoding the full-length SR45a protein relative to other splice variants.

Peroxisomal localisation of the final steps of the mevalonic acid pathway in planta

In plants, the mevalonic acid (MVA) pathway provides precursors for the formation of triterpenes, sesquiterpenes, phytosterols and primary metabolites important for cell integrity. Here, we have cloned the cDNA encoding enzymes catalysing the final three steps of the MVA pathway from Madagascar periwinkle (Catharanthus roseus), mevalonate kinase (MVK), 5-phosphomevalonate kinase (PMK) and mevalonate 5-diphosphate decarboxylase (MVD). These cDNA were shown to functionally complement MVA pathway deletion mutants in the yeast Saccharomyces cerevisiae. Transient transformations of C. roseus cells with yellow fluorescent protein (YFP)-fused constructs reveal that PMK and MVD are localised to the peroxisomes, while MVK was cytosolic. These compartmentalisation results were confirmed using the Arabidopsis thaliana MVK, PMK and MVD sequences fused to YFP. Based on these observations and the arguments raised here we conclude that the final steps of the plant MVA pathway are localised to the peroxisome.

The biochemistry of nitrogen mobilization: purine ring catabolism

The enzymatic route of purine ring catabolism has recently been completed by the discovery of several novel enzymes identified through comparative genome analyses. Here, we review these recent discoveries and present an overview of purine ring catabolism in plants. Xanthine is oxidized to urate in the cytosol, followed by three enzymatic steps taking place in the peroxisome and four reactions in the endoplasmic reticulum releasing the four ring nitrogen as ammonia. Although the main physiological function of purine degradation might lie in the remobilization of nitrogen resources, it has also emerged that catabolic intermediates, the ureides allantoin and allantoate, are likely to be involved in protecting plants against abiotic stress. Conserved alternative splicing mediating the peroxisomal as well as cytosolic localization of allantoin synthase potentially links purine ring catabolism to brassinosteroid signaling.

Probing the evolution of hydroxyisourate hydrolase into transthyretin through active-site redesign

5-Hydroxyisourate hydrolase (HIUase) and transthyretin (TTR) are closely related phylogenetically and structurally, while performing quite different functions. The former catalyzes the hydrolysis of 5-hydroxyisourate within the urate degradation pathway, and the latter is a carrier protein involved in the extracellular transport of thyroid hormones and in the cotransport of retinol. The evolution of HIUase into TTR represents a remarkable example of adaptation of a new function by active-site modification of an enzyme. On the basis of phylogenetic reconstructions and structural comparison of HIUase and TTR, two mutations (Y116T and I16A) were likely to be crucial events in order to induce, after a gene duplication event, the conversion of the enzyme into a binding protein. By rational reshaping of the active sites of HIUase and functional analyses of its mutant forms, we have provided insights into how its neofunctionalization could be achieved. We show here that the two mutations at the active sites of HIUase open up the two ends of the channel that transverses the entire tetrameric protein, generating two cavities accessible to the thyroxine molecule and abrogating, at the same time, the enzymatic activity. Our data indicate that a small number of critical mutations affecting the active site of an enzyme may be sufficient to generate a drastically different function, while a large number of additional mutations may be required for the fine-tuning of the structural and functional features of new proteins.

Toward a definition of the complete proteome of plant peroxisomes: Where experimental proteomics must be complemented by bioinformatics

In the past few years, proteome analysis of Arabidopsis peroxisomes has been established by the complementary efforts of four research groups and has emerged as the major unbiased approach to identify new peroxisomal proteins on a large scale. Collectively, more than 100 new candidate proteins from plant peroxisomes have been identified, including long-awaited low-abundance proteins. More than 50 proteins have been validated as peroxisome targeted, nearly doubling the number of established plant peroxisomal proteins. Sequence homologies of the new proteins predict unexpected enzyme activities, novel metabolic pathways and unknown non-metabolic peroxisome functions. Despite this remarkable success, proteome analyses of plant peroxisomes remain highly material intensive and require major preparative efforts. Characterization of the membrane proteome or post-translational protein modifications poses major technical challenges. New strategies, including quantitative mass spectrometry methods, need to be applied to allow further identifications of plant peroxisomal proteins, such as of stress-inducible proteins. In the long process of defining the complete proteome of plant peroxisomes, the prediction of peroxisome-targeted proteins from plant genome sequences emerges as an essential complementary approach to identify additional peroxisomal proteins that are, for instance, specific to peroxisome variants from minor tissues and organs or to abiotically stressed model and crop plants.

Defining the plant peroxisomal proteome: from Arabidopsis to rice

Peroxisomes are small subcellular organelles mediating a multitude of processes in plants. Proteomics studies over the last several years have yielded much needed information on the composition of plant peroxisomes. In this review, the status of peroxisome proteomics studies in Arabidopsis and other plant species and the cumulative advances made through these studies are summarized. A reference Arabidopsis peroxisome proteome is generated, and some unique aspects of Arabidopsis peroxisomes that were uncovered through proteomics studies and hint at unanticipated peroxisomal functions are also highlighted. Knowledge gained from Arabidopsis was utilized to compile a tentative list of peroxisome proteins for the model monocot plant, rice. Differences in the peroxisomal proteome between these two model plants were drawn, and novel facets in rice were expounded upon. Finally, we discuss about the current limitations of experimental proteomics in decoding the complete and dynamic makeup of peroxisomes, and complementary and integrated approaches that would be beneficial to defining the peroxisomal metabolic and regulatory roadmaps. The synteny of genomes in the grass family makes rice an ideal model to study peroxisomes in cereal crops, in which these organelles have received much less attention, with the ultimate goal to improve crop yield.