Intron splicing and intron-mediated enhanced expression in monocots

Integration analysis of miRNA-mRNA pairs between two contrasting genotypes reveals the molecular mechanism of jujube (Ziziphus jujuba Mill.) response to high-temperature stress

With global warming, high temperature (HT) has become one of the most common abiotic stresses resulting in significant crop yield losses, especially for jujube (Ziziphus jujuba Mill.), an important temperate economic crop cultivated worldwide. This study aims to explore the coping mechanism of jujube to HT stress at the transcriptional and post-transcriptional levels, including identifying differentially expressed miRNAs and mRNAs as well as elucidating the critical pathways involved. High-throughput sequencing analyses of miRNA and mRNA were performed on jujube leaves, which were collected from "Fucumi" (heat-tolerant) and "Junzao" (heat-sensitive) cultivars subjected to HT stress (42 °C) for 0, 1, 3, 5, and 7 days, respectively. The results showed that 45 known miRNAs, 482 novel miRNAs, and 13,884 differentially expressed mRNAs (DEMs) were identified. Among them, integrated analysis of miRNA target genes prediction and mRNA-seq obtained 1306 differentially expressed miRNAs-mRNAs pairs, including 484, 769, and 865 DEMIs-DEMs pairs discovered in "Fucuimi", "Junzao" and two genotypes comparative groups, respectively. Furthermore, functional enrichment analysis of 1306 DEMs revealed that plant-pathogen interaction, starch and sucrose metabolism, spliceosome, and plant hormone signal transduction were crucial pathways in jujube leaves response to HT stress. The constructed miRNA-mRNA network, composed of 20 DEMIs and 33 DEMs, displayed significant differently expressions between these two genotypes. This study further proved the regulatory role of miRNAs in the response to HT stress in plants and will provide a theoretical foundation for the innovation and cultivation of heat-tolerant varieties.

Combination of the endogenous promoter-intron significantly improves salt and drought tolerance conferred by TdSHN1 transcription factor in transgenic tobacco

Recent years have witnessed a renewed interest in introns as a tool to increase gene expression. We previously isolated TdSHN1 gene encoding a transcription factor in durum wheat. Here we show that TdSHN1 intron contains many CT-stretches and the motif CGATT known to be important for IME. When subjected to bioinformatics analysis using IMEter software, TdSHN1 intron obtained a score of 17.04 which indicates that it can moderately enhance gene expression. TdSHN1 gene including its intron was placed under the control of TdSHN1 endogenous salt and drought-inducible promoter or the constitutive 35S promoter and transferred into tobacco. Transgenic lines were obtained and designated gD (with 35S promoter) and PI (with native promoter). A third construct was also used in which intron-less cDNA was driven by the 35S promoter (cD lines). Results showed that, gD lines exhibited lower stomatal density than cD lines. When subjected to drought and salt stresses, gD lines outperformed intron-less cD lines and WT. Indeed, gD lines exhibited longer roots, higher biomass production, retained more chlorophyll, produced less ROS and MDA and had higher antioxidant activity. qRT-PCR analysis revealed that gD lines had higher TdSHN1 expression levels than cD lines. In addition, expression of ROS-scavengering, stress-related and wax biosynthesis tobacco genes was higher in gD lines compared to cD lines and WT. Interestingly, under stress conditions, PI transgenic lines showed higher TdSHN1 expression levels and outperformed gD lines. These results suggest that TdSHN1 intron enhances gene expression when used alone or in combination with TdSHN1 endogenous promoter.

Intron-Mediated Enhancement: A Tool for Heterologous Gene Expression in Plants?

Many plant promoters were characterized and used for transgene expression in plants. Even though these promoters drive high levels of transgene expression in plants, the expression patterns are rarely constitutive but restricted to some tissues and developmental stages. In terms of crop improvement not only the enhancement of expression per se but, in particular, tissue-specific and spatial expression of genes plays an important role. Introns were used to boost expression in transgenic plants in the field of crop improvement for a long time. However, the mechanism behind this so called intron-mediated enhancement (IME) is still largely unknown. This review highlights the complexity of IME on the levels of its regulation and modes of action and gives an overview on IME methodology, examples in fundamental research and models of proposed mechanisms. In addition, the application of IME in heterologous gene expression is discussed.

A Leader Intron of a Soybean Elongation Factor 1A (eEF1A) Gene Interacts with Proximal Promoter Elements to Regulate Gene Expression in Synthetic Promoters

Introns, especially the first intron in the 5' untranslated region (5'UTR), can significantly impact gene expression via intron-mediated enhancement (IME). In this study, we demonstrate the leader intron of a soybean elongation factor 1A (eEF1A) gene (GmScreamM8) was essential for the high activity of the native promoter. Furthermore, the interaction of the GmScreamM8 leader intron with regulatory element sequences from several soybean eEF1A promoters was studied using synthetic promoters, which consisted of element tetramers upstream of a core promoter used to regulate a green fluorescent protein (gfp) reporter gene. Element tetramers, placed upstream of a GmScreamM8 core promoter, showed very high activity using both transient expression in lima bean cotyledons and stable expression in soybean hairy roots, only if the native leader intron was included, suggesting an interaction between intronic sequences and promoter elements. Partial deletions of the leader intron showed that a 222 bp intronic sequence significantly contributed to very high levels of GFP expression. Generation of synthetic intron variants with a monomeric or trimeric repeat of the 222 bp intronic sequence, yielded almost two-fold higher expression compared to the original intron, while partial deletion of the 222 bp intronic repeated sequence significantly decreased gene expression, indicating that this intronic sequence was essential for the intron-element interaction enhancement.

Isolation and characterization of "GmScream" promoters that regulate highly expressing soybean (Glycine max Merr.) genes

To increase our understanding of the regulatory components that control gene expression, it is important to identify, isolate and characterize new promoters. In this study, a group of highly expressed soybean (Glycine max Merr.) genes, which we have named "GmScream", were first identified from RNA-Seq data. The promoter regions were then identified, cloned and fused with the coding region of the green fluorescent protein (gfp) gene, for introduction and analysis in different tissues using 3 tools for validation. Approximately half of the GmScream promoters identified showed levels of GFP expression comparable to or higher than the Cauliflower Mosaic Virus 35S (35S) promoter. Using transient expression in lima bean cotyledonary tissues, the strongest GmScream promoters gave over 6-fold higher expression than the 35S promoter while several other GmScream promoters showed 2- to 3-fold higher expression. The two highest expressing promoters, GmScreamM4 and GmScreamM8, regulated two different elongation factor 1A genes in soybean. In stably transformed soybean tissues, GFP driven by the GmScreamM4 or GmScreamM8 promoter exhibited constitutive high expression in most tissues with preferentially higher expression in proliferative embryogenic tissues, procambium, vascular tissues, root tips and young embryos. Using deletion analysis of the promoter, two proximal regions of the GmScreamM8 promoter were identified as contributing significantly to high levels of gene expression.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1914-5) contains supplementary material, which is available to authorized users.

Testing the IMEter on rice introns and other aspects of intron-mediated enhancement of gene expression

In many eukaryotes, spliceosomal introns are able to influence the level and site of gene expression. The mechanism of this Intron Mediated Enhancement (IME) has not yet been elucidated, but regulation of gene expression is likely to occur at several steps during and after transcription. Different introns have different intrinsic enhancing properties, but the determinants of these differences remain unknown. Recently, an algorithm called IMEter, which is able to predict the IME potential of introns without direct testing, has been proposed. A computer program was developed for Arabidopsis thaliana and rice (Oryza sativa L.), but was only tested experimentally in Arabidopsis by measuring the enhancement effect on GUS expression of different introns inserted within otherwise identical plasmids. To test the IMEter potential in rice, a vector bearing the upstream regulatory sequence of a rice β-tubulin gene (OsTub6) fused to the GUS reporter gene was used. The enhancing intron interrupting the OsTub6 5'-UTR was precisely replaced by seven other introns carrying different features. GUS expression level in transiently transformed rice calli does not significantly correlate with the calculated IMEter score. It was also found that enhanced GUS expression was mainly due to a strong increase in the mRNA steady-state level and that mutations at the splice recognition sites almost completely abolished the enhancing effect. Splicing also appeared to be required for IME in Arabidopsis cell cultures, where failure of the OsTub6 5' region to drive high level gene expression could be rescued by replacing the poorly spliced rice intron with one from Arabidopsis.

Transgene expression systems in the Triticeae cereals

The control of transgene expression is vital both for the elucidation of gene function and for the engineering of transgenic crops. Given the dominance of the Triticeae cereals in the agricultural economy of the temperate world, the development of well-performing transgene expression systems of known functionality is of primary importance. Transgenes can be expressed either transiently or stably. Transient expression systems based on direct or virus-mediated gene transfer are particularly useful in situations where the need is to rapidly screen large numbers of genes. However, an unequivocal understanding of gene function generally requires that a transgene functions throughout the plant's life and is transmitted through the sexual cycle, since this alone allows its effect to be decoupled from the plant's response to the generally stressful gene transfer event. Temporal, spatial and quantitative control of a transgene's expression depends on its regulatory environment, which includes both its promoter and certain associated untranslated region sequences. While many transgenic approaches aim to manipulate plant phenotype via ectopic gene expression, a transgene sequence can be also configured to down-regulate the expression of its endogenous counterpart, a strategy which exploits the natural gene silencing machinery of plants. In this review, current technical opportunities for controlling transgene expression in the Triticeae species are described. Apart from protocols for transient and stable gene transfer, the choice of promoters and other untranslated regulatory elements, we also consider signal peptides, as they too govern the abundance and particularly the sub-cellular localization of transgene products.

The leader intron of AtMHX can elicit, in the absence of splicing, low-level intron-mediated enhancement that depends on the internal intron sequence

BACKGROUND

Introns stimulate gene expression in a wide range of organisms by increasing the levels of mature mRNA, without affecting mRNA stability. Although introns sometimes function as transcriptional enhancers, they usually stimulate expression by a process termed intron-mediated enhancement (IME). The mechanism of IME is largely unknown. While splicing per se is not sufficient for IME, as evident from the fact that not all introns increase expression, it is not clear yet whether splicing of the enhancing introns is essential for enhancement. The leader intron (LI) of the Arabidopsis AtMHX gene was previously shown to substantially increase the expression of the AtMHX promoter. Here we investigated whether this LI acts as a transcriptional enhancer and whether its splicing is essential for IME.

RESULTS

Expression in transformed Arabidopsis plants of an AtMHX::GUS construct from which the LI was eliminated was similar to a construct that included only the minimal promoter fused to GUS. Yet, almost no expression was seen in constructs that included the LI in addition to the minimal promoter or the LI inserted in various locations in the promoter. While the LI enhanced 272-fold the expression of the weak AtMHX promoter, only a 3-fold enhancement was observed for the strong CaMV 35S promoter. In the context of the AtMHX promoter, an unspliceable version of the LI that had mutated 5' and 3' splice sites mediated a low-level (5-fold) enhancement. Eliminating the internal 320 nt of the 416 nt unspliceable intron resulted in loss of ability to mediate low-level enhancement.

CONCLUSIONS

Although AtMHX promoter shows almost no expression in the absence of its LI, this intron does not act as a transcriptional enhancer and is unable to support expression in the absence of the enhancer elements of the promoter. It is also shown that the same intron can have very different contributions to expression of different promoters. Our results also demonstrate that while splicing is essential for substantial IME, in the absence of splicing low-level enhancement can be obtained. Notably, it is shown that the internal intron sequence plays a significant role in mediating the low-level enhancement of unspliced introns.

Two-dimensional liquid chromatography technique coupled with mass spectrometry analysis to compare the proteomic response to cadmium stress in plants

Plants are useful in studies of metal toxicity, because their physiological responses to different metals are correlated with the metal exposure dose and chemical state. Moreover a network of proteins and biochemical cascades that may lead to a controlled homeostasis of metals has been identified in many plant species. This paper focuses on the global protein variations that occur in a Populus nigra spp. clone (Poli) that has an exceptional tolerance to the presence of cadmium. Protein separation was based on a two-dimensional liquid chromatography technique. A subset of 20 out of 126 peaks were identified as being regulated differently under cadmium stress and were fingerprinted by MALDI-TOF. Proteins that were more abundant in the treated samples were located in the chloroplast and in the mitochondrion, suggesting the importance of these organelles in the response and adaptation to metal stress.

In trangenic rice, alpha- and beta-tubulin regulatory sequences control GUS amount and distribution through intron mediated enhancement and intron dependent spatial expression

The genomic upstream sequence of the rice tubulin gene OsTub6 has been cloned, sequenced and characterized. The 5'UTR sequence is interrupted by a 446 bp long leader intron. This feature is shared with two other rice beta-tubulin genes (OsTub4 and OsTub1) that, together with OsTub6, group in the same clade in the evolutionary phylogenetic tree of plant beta-tubulins. Similarly to OsTub4, the leader intron of OsTub6 is capable of sustaining intron mediated enhancement (IME) of gene expression, in transient expression assays. A general picture is drawn for three rice alpha-tubulin and two rice beta-tubulin genes in which the first intron of the coding sequence for the formers and the intron present in the 5'UTR for the latters, are important elements for controlling gene expression. We used OsTua2:GUS, OsTua3:GUS, OsTub4:GUS and OsTub6:GUS chimeric constructs to investigate the in vivo pattern of beta-glucuronidase (GUS) expression in transgenic rice plants. The influence of the regulatory introns on expression patterns was evaluated for two of them, OsTua2 and OsTub4. We have thus characterized distinct patterns of expression attributable to each tubulin isotype and we have shown that the presence of the regulatory intron can greatly influence both the amount and the actual site of expression. We propose the term Intron Dependent Spatial Expression (IDSE) to highlight this latter effect.

Plant spliceosomal introns: not only cut and paste

Spliceosomal introns in higher eukaryotes are present in a high percentage of protein coding genes and represent a high proportion of transcribed nuclear DNA. In the last fifteen years, a growing mass of data concerning functional roles carried out by such intervening sequences elevated them from a selfish burden carried over by the nucleus to important active regulatory elements. Introns mediate complex gene regulation via alternative splicing; they may act in cis as expression enhancers through IME (intron-mediated enhancement of gene expression) and in trans as negative regulators through the generation of intronic microRNA. Furthermore, some introns also contain promoter sequences for alternative transcripts. Nevertheless, such regulatory roles do not require long conserved sequences, so that introns are relatively free to evolve faster than exons: this feature makes them important tools for evolutionary studies and provides the basis for the development of DNA molecular markers for polymorphisms detection. A survey of introns functions in the plant kingdom is presented.

Analysis of regulatory elements of the promoter and the 3' untranslated region of the maize Hrgp gene coding for a cell wall protein

Hydroxyproline-rich glycoproteins (HRGP) are structural components of the plant cell wall. Hrgp genes from maize and related species have a conserved 500 bp sequence in the 5'-flanking region, and all Hrgp genes from monocots have an intron located in the 3' untranslated region. To study the role of these conserved regions, several deletions of the Hrgp gene were fused to the beta-glucuronidase ( GUS) gene and used to transform maize tissues by particle bombardment. The overall pattern of GUS activity directed by sequential deletions of the Hrgp promoter was different in embryos and young shoots. In embryos, the activity of the full-length Hrgp promoter was in the same range as that of the p35SI promoter construct, based on the strong 35S promoter, whereas in the fast-growing young shoots it was 20 times higher. A putative silencer element specific for young shoots was found in the -1,076/-700 promoter region. Other major cis elements for Hrgp expression are probably located in the regions spanning -699/-510 and -297/-160. Sequences close to the initial ATG and mRNA leader were also important since deletion of the region -52/+16 caused a 75% reduction in promoter activity. The presence of the Hrgp intron in the 3' untranslated region changed the levels of GUS activity directed by the Hrgp and the 35S promoters. This pattern of activity was complex, and was dependent on the promoter and cell type analysed.

Splicing of the maize Sh1 first intron is essential for enhancement of gene expression, and a T-rich motif increases expression without affecting splicing

Certain plant and animal introns increase expression of protein-coding sequences when placed in the 5' region of the transcription unit. The mechanisms of intron-mediated enhancement have not been defined, but are generally accepted to be post- or cotranscriptional in character. One of the most effective plant introns in stimulating gene expression is the 1,028-bp first intron of the Sh1 gene that encodes maize (Zea mays) sucrose synthase. To address the mechanisms of intron-mediated enhancement, we used reporter gene fusions to identify features of the Sh1 first intron required for enhancement in cultured maize cells. A 145-bp derivative conferred approximately the same 20- to 50-fold stimulation typical for the full-length intron in this transient expression system. A 35-bp motif contained within the intron is required for maximum levels of enhancement but not for efficient transcript splicing. The important feature of this redundant 35-bp motif is T-richness rather than the specific sequence. When transcript splicing was abolished by mutations at the intron borders, enhancement was reduced to about 2-fold. The requirement of splicing for enhancement was not because of upstream translation initiation codons contained in unspliced transcripts. On the basis of our current findings, we conclude that splicing of the Sh1 intron is integral to enhancement, and we hypothesize that transcript modifications triggered by the T-rich motif and splicing may link the mRNA with the trafficking system of the cell.

Intron-dependent stimulation of marker gene expression in cultured insect cells

We tested in a systematic fashion the effect of an intron on the level of luciferase expression in cultured C6/36 Aedes albopictus cells. The intron was inserted in both orientations, upstream and downstream of the luciferase coding region in two different luciferase expression vectors. The two parental luciferase expression vectors differed only in their promoters, one containing the Drosophila melanogaster actin5C promoter and the other the Autographa californica nuclear polyhedrosis virus hr5/ie1 enhancer/promoter. All resulting intron-containing constructs were tested for their ability to express luciferase in transient assays following electroporation into C6/36 cells. We found that the introns stimulate luciferase expression between twelve and sixtyfold, depending on the promoter. Enhanced expression was only seen when the intron was present in the correct orientation upstream of the luciferase ORF. When the 3' splice sites of the enhanced intron-containing constructs were mutated, the expression level dropped back to below the level of the intronless parental constructs, suggesting that the intron-dependent stimulation of luciferase expression is depending on splicing and is not due to other effects the intron may have on transcription or translation. The luciferase transcripts of all constructs were analysed by reverse transcription, PCR amplification and sequencing, and the results show a perfect correlation between efficient splicing of the intron and elevated levels of luciferase expression. Our findings have the potential to be very useful for boosting expression of foreign proteins in the widely used baculoviral or non-viral systems in insect cells.

A long leader intron of the Ostub16 rice beta-tubulin gene is required for high-level gene expression and can autonomously promote transcription both in vivo and in vitro

A 2 kb DNA fragment, upstream of the rice beta-tubulin isotype 16 (Ostub16) coding sequence, was isolated using inverse PCR and screening of a tubulin-enriched lambda library. An intron (863 bp) present in the 5' untranslated region (5' UTR) is spliced out to produce the most abundant mRNA species which corresponds to the previously cloned Ostub16 cDNA. Transient expression assays performed on rice embryogenic calluses with chimeric Ostub16::GUS constructs demonstrated that the entire 2 kb upstream sequence has a strong promoter activity, and that the 863 bp intron is required for high-level GUS expression. In addition, the intron sequence is capable per se of sustaining a weak but consistent GUS expression. Two rare Ostub16 transcripts, with a start site mapping within this intron sequence, were detected in rice coleoptile cells. The transcription start site mapped at position -290 with respect to the ATG codon, and the shorter molecule originated from splicing of the same precursor mRNA. Therefore transcriptional expression of rice beta-tubulin isotype 16 results in the synthesis of two premRNA molecules (I and II) encoding for three different mRNA species. We discuss these findings in terms of function and molecular evolution of the mechanisms that control plant beta-tubulin gene expression.

Intron-mediated enhancement of gene expression independent of unique intron sequences and splicing

Either of the first two introns of the Arabidopsis tryptophan pathway gene PAT1 elevates mRNA accumulation from a PAT1:beta-glucuronidase (GUS) fusion roughly 5-fold without affecting the rate of PAT1:GUS transcription. To further explore the mechanism of this intron-mediated enhancement of gene expression, we wanted to determine whether splicing or specific intron sequences were necessary. In-frame derivatives of PAT1 intron 1, whose splicing was prevented by a point mutation or large deletions, were able to increase mRNA accumulation from a PAT1:GUS fusion, demonstrating that splicing per se is not required. Furthermore, each of a series of introns containing overlapping deletions that together span PAT1 intron 1 increased PAT1:GUS mRNA accumulation as much as the full-length intron did, indicating that all intron sequences are individually dispensable for this phenomenon. These results eliminate the simple idea that this intron stimulates mRNA accumulation via a unique RNA-stabilizing sequence or through the completed act of splicing. However, they are consistent with a possible role for redundant intron sequence elements or an association of the pre-mRNA with the spliceosome.

Genetic variation within and among populations of the threatened lichen Lobaria pulmonaria in Switzerland and implications for its conservation

The foliose epiphytic lichen Lobaria pulmonaria has suffered a significant decline in European lowlands during the last decades and therefore is considered as endangered throughout Europe. An assessment of the genetic variability is necessary to formulate biologically sound conservation recommendations for this species. We investigated the genetic diversity of the fungal symbiont of L. pulmonaria using 143 specimens sampled from six populations (two small, one medium, three large) in the lowland, the Jura Mountains, the pre-Alps and the Alps of Switzerland. Among all nuclear and mitochondrial regions sequenced for this study, variability was found only in the internal transcribed spacer (ITS I), with three polymorphic sites, and in the nuclear ribosomal large subunit (nrLSU), with four polymorphic sites. The variable sites in the nrLSU are all located within a putative spliceosomal intron. We sequenced these two regions for 81 specimens and detected six genotypes. Two genotypes were common, two were found only in the more diverse populations and two were found only in one population each. There was no correlation between population size and genetic diversity. The highest genetic diversity was found in populations where the fungal symbiont is reproducing sexually. Populations with low genetic diversity included only the two same common genotypes. Our study provides evidence suggesting that L. pulmonaria is self-incompatible and heterothallic. Based on our results we give populations with sexually reproducing individuals a higher rank in terms of conservation priority than strictly asexual populations. The remaining lowland populations are so small, that one single catastrophic event such as a windthrow might destroy the entire population. Hence we suggest augmenting such populations in size and genetic diversity using small thallus fragments or vegetative diaspores collected in other populations. As we did not detect any locally adapted genotypes, these transplants can be taken from any other genetically diverse population in Switzerland.

Prediction of locally optimal splice sites in plant pre-mRNA with applications to gene identification in Arabidopsis thaliana genomic DNA

Prediction of splice site selection and efficiency from sequence inspection is of fundamental interest (testing the current knowledge of requisite sequence features) and practical importance (genome annotation, design of mutant or transgenic organisms). In plants, the dominant variables affecting splice site selection and efficiency include the degree of matching to the extended splice site consensus and the local gradient of U- and G+C-composition (introns being U-rich and exons G+C-rich). We present a novel method for splice site prediction, which was particularly trained for maize and Arabidopsis thaliana. The method extends our previous algorithm based on logitlinear models by considering three variables simultaneously: intrinsic splice site strength, local optimality and fit with respect to the overall splice pattern prediction. We show that the method considerably improves prediction specificity without compromising the high degree of sensitivity required in gene prediction algorithms. Applications to gene identification are illustrated for Arabidopsis and suggest that successful methods must combine scoring for splice sites, coding potential and similarity with potential homologs in non-trivial ways. A WWW version of the SplicePredictor program is available at http:/gnomic.stanford.edu/volker/SplicePredi ctor.html/

Prediction of splice sites in plant pre-mRNA from sequence properties

Heterologous introns are often inaccurately or inefficiently processed in higher plants. The precise features that distinguish the process of pre-mRNA splicing in plants from splicing in yeast and mammals are unclear. One contributing factor is the prominent base compositional contrast between U-rich plant introns and flanking G + C-rich exons. Inclusion of this contrast factor in recently developed statistical methods for splice site prediction from sequence inspection significantly improved prediction accuracy. We applied the prediction tools to re-analyze experimental data on splice site selection and splicing efficiency for native and more than 170 mutated plant introns. In almost all cases, the experimentally determined preferred sites correspond to the highest scoring sites predicted by the model. In native genes, about 90% of splice sites are the locally highest scoring sites within the bounds of the flanking exon and intron. We propose that, in most cases, local context (about 50 bases upstream and downstream from a potential intron end) is sufficient to account for intrinsic splice site strength, and that competition for transacting factors determines splice site selection in vivo. We suggest that computer-aided splice site prediction can be a powerful tool for experimental design and interpretation.

Logitlinear models for the prediction of splice sites in plant pre-mRNA sequences

Pre-mRNA splicing in plants, while generally similar to the processes in vertebrates and yeast, is thought to involve plant specific cis-acting elements. Both monocot and dicot introns are typically strongly enriched in U nucleotides, and AU- or U-rich segments are thought to be involved in intron recognition, splice site selection, and splicing efficiency. We have applied logitlinear models to find optimal combinations of splice site variables for the purpose of separating true splice sites from a large excess of potential sites. It is shown that plant splice site prediction from sequence inspection is greatly improved when compositional contrast between exons and introns is considered in addition to degree of matching to the splice site consensus (signal quality). The best model involves subclassification of splice sites according to the identity of the base immediately upstream of the GU and AG signals and gives substantial performance gains compared with conventional profile methods.

Splicing of precursors to mRNA in higher plants: mechanism, regulation and sub-nuclear organisation of the spliceosomal machinery

The removal of introns from pre-mRNA transcripts and the concomitant ligation of exons is known as pre-mRNA splicing. It is a fundamental aspect of constitutive eukaryotic gene expression and an important level at which gene expression is regulated. The process is governed by multiple cis-acting elements of limited sequence content and particular spatial constraints, and is executed by a dynamic ribonucleoprotein complex termed the spliceosome. The mechanism and regulation of pre-mRNA splicing, and the sub-nuclear organisation of the spliceosomal machinery in higher plants is reviewed here. Heterologous introns are often not processed in higher plants indicating that, although highly conserved, the process of pre-mRNA splicing in plants exhibits significant differences that distinguish it from splicing in yeast and mammals. A fundamental distinguishing feature is the presence of and requirement for AU or U-rich intron sequence in higher-plant pre-mRNA splicing. In this review we document the properties of higher-plant introns and trans-acting spliceosomal components and discuss the means by which these elements combine to determine the accuracy and efficiency of pre-mRNA processing. We also detail examples of how introns can effect regulated gene expression by affecting the nature and abundance of mRNA in plants and list the effects of environmental stresses on splicing. Spliceosomal components exhibit a distinct pattern of organisation in higher-plant nuclei. Effective probes that reveal this pattern have only recently become available, but the domains in which spliceosomal components concentrate were identified in plant nuclei as enigmatic structures some sixty years ago. The organisation of spliceosomal components in plant nuclei is reviewed and these recent observations are unified with previous cytochemical and ultrastructural studies of plant ribonuleoprotein domains.

Intron position affects expression from the tpi promoter in rice

A series of promoter-GUS fusion constructs containing a portion of the rice triosephosphate isomerase (tpi) promoter, the first tpi intron, and the gene encoding bacterial beta-glucuronidase (GUS) were made. These constructs were electroporated into rice protoplasts and transient expression was monitored. Inclusion of the first intron from the rice tpi gene enhanced expression of the GUS gene from the tpi promoter when it was placed 5' of the GUS gene. When the tpi intron was placed in the 3'-untranslated region no enhancement of GUS gene expression was observed, indicating the importance of position in intron-mediated enhancement of gene expression.

Intron-specific stimulation of anaerobic gene expression and splicing efficiency in maize cells

Most of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) genes characterized in plants and algae to date have one intron very close to the 5' end of the gene. To study the functional relevance of some of these introns for gene expression we have analysed the influence of three 5' introns on transient gene expression of the anaerobically inducible maize GapC4 promoter in maize cells. Under aerobic conditions, reporter gene expression is increased in the presence of the first introns of the GapC4 and GapC1 genes, and the first intron of the nuclear encoded chloroplast-specific GapA1 gene. In contrast, the GapC4 intron increases anaerobic gene expression above the level obtained for the intronless construct, while anaerobic expression of constructs harboring the GapA1 and GapC1 introns was similar to the anaerobic expression level of the intronless construct. Splicing analysis revealed that the GapC4 intron is processed more efficiently under anaerobic conditions, while no change in splicing efficiency is observed for the GapC1 and the GapA1 introns when subjected to anaerobic conditions. These results suggest that an increase in splicing efficiency contributes to the anaerobic induction of the maize GapC4 gene.

Characterization of intronic uridine-rich sequence elements acting as possible targets for nuclear proteins during pre-mRNA splicing in Nicotiana plumbaginifolia

Introns of nuclear pre-mRNAs in dicotyledonous plants, unlike introns in vertebrates or yeast, are distinctly rich in A+U nucleotides and this feature is essential for their processing. In order to define more precisely sequence elements important for intron recognition in plants, we investigated the effects of short insertions, either U-rich or A-rich, on splicing of synthetic introns in transfected protoplast of Nicotiana plumbaginifolia. It was found that insertions of U-rich (sequence UUUUUAU) but not A-rich (AUAAAAA) segments can activate splicing of a GC-rich synthetic infron, and that U-rich segments, or multimers thereof, can function irrespective of the site of insertion within the intron. Insertions of multiple U-rich segments, either at the same or different locations, generally had an additive, stimulatory effect on splicing. Mutational analysis showed that replacement of one or two U residues in the UUUUUAU sequence with A or C residues had only a small effect on splicing, but replacement with G residues was strongly inhibitory. Proteins that interact with fragments of natural and synthetic pre-mRNAs in vitro were identified in nuclear extracts of N.plumbaginifolia by UV cross- linking. The profile of cross-linked plant proteins was considerably less complex than that obtained with a HeLa cell nuclear extract. Two major cross-linkable plant proteins had apparent molecular mass of 50 and 54 kDa and showed affinity for oligouridilates present in synGC introns or for poly(U).

Expression of heat shock factor and heat shock protein 70 genes during maize pollen development

We have analysed the expression of heat shock protein 70 (HSP70) and heat shock factor (HSF) gene during maize pollen development, HSFs being the transcriptional activators of hsp genes. In order to eliminate the sporophytic tissues of anthers, we have isolated homogeneous cell populations corresponding to five stages of maize pollen development from microspores to mature pollen. We show that in the absence of heat stress, hsp70 genes are highly expressed late-bicellular pollen as compared to other stages. HSP70 transcripts are significantly accumulated in response to a heat shock at the late microspore stage but to a much lower extent than in vegetative tissues. The latest stages of pollen development, i.e. mid-tricellular and mature pollen, do not exhibit heat-induced accumulation of HSP70 transcripts. Therefore, we analysed the expression of hsf genes throughout pollen development. We demonstrate that at least three hsf genes are expressed in maize and that transcripts corresponding to one hsf gene, whose expression is independent of temperature in somatic as well as in microgametophytic tissues, are present at similar levels throughout pollen development. In addition, we show that the expression of the two other hsf genes is heat-inducible in maize vegetative tissues and is not significantly increased after heat shock at any stage of pollen development. These results indicate that the loss of hsp gene expression at late stages of pollen development is not due to a modification of hsf gene expression at the mRNA level and that hsf gene expression is differentially regulated in vegetative and microgametophytic tissues.

The plastid aldolase gene from Chlamydomonas reinhardtii: intron/exon organization, evolution, and promoter structure

Genomic clones encoding the plastidic fructose-1,6-bisphosphate aldolase of Chlamydomonas reinhardtii were isolated and sequenced. The gene contains three introns which are located within the coding sequence for the mature protein. No introns are located within or near the sequence encoding the transit-peptide, in contrast to the genes for plastidic aldolases of higher plants. Neither the number nor the positions of the three introns of the C. reinhardtii aldolase gene are conserved in the plastidic or cytosolic aldolase genes of higher plants and animals. The 5' border sequences of introns in the aldolase gene of C. reinhardtii exhibit the conserved plant consensus sequence. The 3' acceptor splice sites for introns 1 and 3 show much less similarity to the eukaryotic consensus sequences than do those of intron 2. The plastidic aldolase gene has two tandemly repeated CAAT box motifs in the promoter region. Genomic Southern blots indicate that the gene is encoded by a single locus in the C. reinhardtii genome.

Intron-dependent transient expression of the maize GapA1 gene

Transient expression experiments show that the maize GapA1 promoter exhibits a requirement for sequences contained within intron 1 and surrounding exon border regions for expression in maize Black Mexican Sweet cells. Maize GapA1-promoter constructs lacking intron 1 are inactive. Intron 1 and its exon border sequences, when reintroduced into constructs lacking introns, restore gene activity whereas intron 2 and its exon borders to not. The minimal promoter so defined encompasses roughly 250 bp upstream of the in vivo transcription start and appears also to include intron 1. An octameric sequence was identified in intron 1 of maize GapA1 which is similar to sequence motifs found in other maize introns known to increase transient expression. Partial restoration of gene expression in GapA1 constructs lacking intron 1 was achieved through insertion of the identified octameric sequence.

The marine red alga Chondrus crispus has a highly divergent beta-tubulin gene with a characteristic 5' intron: functional and evolutionary implications

We characterized a nuclear gene and its corresponding cDNA encoding beta-tubulin (gene TubB1) of the marine red alga Chondrus crispus. The deduced TubB1 protein is the most divergent beta-tubulin so far reported with only 64 to 69% amino acid identity relative to other beta-tubulins from higher and lower eukaryotes. Our analysis reveals that TubB1 has an accelerated evolutionary rate probably due to a release of functional constraints in connexion with a specialization of microtubular structures in rhodophytes. It further indicates that isoform diversity and functional differentiation of tubulins in eukaryotic cells may be controlled by independent selective constraints. TubB1 has a short spliceosomal intron at its 5' end which seems to be a characteristic feature of nuclear protein-coding genes from rhodophytes. The splice junctions of the four known rhodophyte introns comply well with the corresponding consensus sequences of higher plants in agreement with previous suggestions from phylogenetic inference that red algae and green plants may be sister groups. The paucity and asymmetrical location of introns in rhodophyte genes can be explained by differential intron loss due to conversion of genes by homologous recombination with cDNAs corresponding to reverse transcribed mRNAs or partially spliced pre-mRNAs, respectively. The identification of an intron containing TubB1 cDNA in C. crispus confirms that pre-mRNAs can escape both splicing and degradation in the nucleus prior to transport into the cytoplasm. Differential Southern hybridizations under non-stringent conditions with homologous and heterologous probes suggest that C. crispus contains a second degenerate beta-tubulin gene (or pseudogene?) which, however, is only distantly related to TubB1 as it is to the more conserved homologues of other organisms.

Arabidopsis heat shock factor: isolation and characterization of the gene and the recombinant protein

We have isolated the Arabidopsis heat shock factor gene Athsf1 as genomic and corresponding cDNA sequences via cross-hybridization with tomato clones. Sequence analysis indicates only a partial homology with the HSFs from tomato and other organisms which is confined to the DNA-binding and the oligomerization domains. The gene is constitutively expressed but the level of mRNA for Athsf1 increases two-fold upon heat shock. However, the putative promoter region lacks the canonical heat shock elements. After expression in Escherichia coli the recombinant Athsf1 protein binds specifically to a synthetic oligonucleotide containing five heat shock elements. The native size of recombinant ATHSF1 in vitro is consistent with a trimer as demonstrated by chemical cross-linking and pore exclusion limit analysis.

The GAPDH gene system of the red alga Chondrus crispus: promoter structures, intron/exon organization, genomic complexity and differential expression of genes

Our previous phylogenetic analysis based on cDNA sequences of chloroplast and cytosolic glyceraldehyde-3-phosphate dehydrogenases (GAPDH; genes GapA and GapC, respectively) of the red alga Chondrus crispus suggested that rhodophytes and green plants are sister groups with respect to plastids and mitochondria and diverged at about the same time or somewhat later than animals and fungi. Here we characterize the genomic sequences of genes GapC and GapA of C. crispus with respect to promotor structures, intron/exon organization, genomic complexity, G + C content, CpG suppression and their transcript levels in gametophytes and protoplasts, respectively. To our knowledge this is the first report on nuclear protein genes of red algae. The GapC gene is G + C-rich, contains no introns and displays a number of classic sequence motifs within its promotor region, such as TATA, CAAT, GC boxes and several elements resembling the plant-specific G-box palindrome. The GapA gene has a moderate G+C content, a single CAAT box motif in its promotor region and a single intron of 115 bp near its 5' end. This intron occupies a conserved position corresponding to that of intron 1 in the transit peptide region of chloroplast GAPDH genes (GapA and GapB) of higher plants. It has consensus sequences similar to those of yeast introns and folds into a conspicuous secondary structure of -61.3 kJ. CpG profiles of genes GapC and GapA and their flanking sequences show no significant CpG depletion suggesting that these genomic sequences are not methylated. Genomic Southern blots hybridized with generic and gene specific probes indicate that both genes are encoded by single loci composed of multiple polymorphic alleles. Northern hybridizations demonstrate that both genes are expressed in gametophytes but not in protoplasts where appreciable amounts of transcripts can only be detected for GapC.

Plant pre-mRNA splicing and splicing components

Pre-mRNA splicing or the removal of introns from precursor messenger RNAs depends on the accurate recognition of intron sequences by the plant splicing machinery. The major components of this machinery are small nuclear ribonucleoprotein protein particles (snRNPs) which consist of snRNAs and snRNP proteins. We have analysed various aspects of intron sequence and structure in relation to splice site selection and splicing efficiency and we have cloned snRNA genes and a gene encoding the snRNP protein, U2B". In the absence of an in vitro splicing system for plants, transient expression in protoplasts and stable plant transformations have been used to analyse splicing of intron constructs. We aim to address the function of the UsnRNP-specific protein, U2B", via the production of transgenic plants expressing antisense U2B" transcripts and epitope-tagged U2B" protein. In addition, we have cloned genes encoding other proteins which potentially interact with RNA, such as RNA helicases, and strategies involving transgenic plants are being developed to analyse their function.

Poly(A) tail length of a heat shock protein RNA is increased by severe heat stress, but intron splicing is unaffected

The small heat shock proteins (sHSPs) are induced in all eukaryotes in response to high temperature stress, but are most abundant among members of the plant kingdom where they accumulate in multiple subcellular compartments. We have analyzed the expression of the chloroplast-localized sHSP from Arabidopsis thaliana, HSP21, and characterized the structure of the gene encoding this protein to facilitate future genetic studies on the function of HSP21 in the heat shock response. HSP21 is encoded in Arabidopsis by a single gene whose coding region is interrupted by a single intron. Previous studies have shown that intron processing is disrupted by severe, abrupt heat stress but is protected by pretreatments that induce thermotolerance. The processing of the HSP21 transcript was investigated in response to an abrupt heat stress regime and a gradual heat stress regime, the latter of which is known to confer thermotolerance in plants. Under abrupt stress conditions the HSP21 transcript is somewhat longer than under gradual heat stress conditions. However, the molecular basis for the size difference is not impaired intron splicing, but rather a difference in the length of the poly(A) tail depending on the heat stress regime. The results suggest that an increase in poly(A) tail length may be a generalized response to severe, abrupt heat stress and that poly(A) tail metabolism may be one of numerous cellular processes normally protected in thermotolerant cells from the otherwise damaging effects of high temperature stress.

Characterization of two maize HSP90 heat shock protein genes: expression during heat shock, embryogenesis, and pollen development

We have isolated two genes from Zea mays encoding proteins of 82 and 81 kD that are highly homologous to the Drosophila 83-kD heat shock protein gene and have analyzed the structure and pattern of expression of these two genes during heat shock and development. Southern blot analysis and hybrid select translations indicate that the highly homologous hsp82 and hsp81 genes are members of a small multigene family composed of at least two and perhaps three or more gene family members. The deduced amino acid sequence of these proteins based on the nucleotide sequence of the coding regions shows 64-88% amino acid homology to other hsp90 family genes from human, yeast, Drosophila, and Arabidopsis. The promoter regions of both the hsp82 and hsp81 genes contain several heat shock elements (HSEs), which are putative binding sites for heat shock transcription factor (HSF) commonly found in the promoters of other heat shock genes. Gene-specific oligonucleotide probes were synthesized and used to examine the mRNA expression patterns of the hsp81 and hsp82 genes during heat shock, embryogenesis, and pollen development. The hsp81 gene is only mildly heat inducible in leaf tissue, but is strongly expressed in the absence of heat shock during the pre-meiotic and meiotic prophase stages of pollen development and in embryos, as well as in heat-shocked embryos and tassels. The hsp82 gene shows strong heat inducibility at heat-shock temperatures (37-42 degrees C) and in heat shocked embryos and tassels but is only weakly expressed in the absence of heat shock. Promoter-GUS reporter gene fusions made and analyzed by transient expression assays in Black Mexican Sweet (BMS) Maize protoplasts also indicate that the hsp82 and hsp81 are regulated differentially. The hsp82 promoter confers strong heat-inducible expression of the GUS reporter gene in heat-treated cells (60- to 80-fold over control levels), whereas the hsp81 promoter is only weakly heat inducible (5- to 10-fold over control levels).