Phylogenetic analysis and stress response of the plant U2 small nuclear ribonucleoprotein B″ gene family
BMC Genomics volume 23, Article number: 744 (2022)
Alternative splicing (AS) is an important channel for gene expression regulation and protein diversification, in addition to a major reason for the considerable differences in the number of genes and proteins in eukaryotes. In plants, U2 small nuclear ribonucleoprotein B″ (U2B″), a component of splicing complex U2 snRNP, plays an important role in AS. Currently, few studies have investigated plant U2B″, and its mechanism remains unclear.
Phylogenetic analysis, including gene and protein structures, revealed that U2B″ is highly conserved in plants and typically contains two RNA recognition motifs. Subcellular localisation showed that OsU2B″ is located in the nucleus and cytoplasm, indicating that it has broad functions throughout the cell. Elemental analysis of the promoter region showed that it responded to numerous external stimuli, including hormones, stress, and light. Subsequent qPCR experiments examining response to stress (cold, salt, drought, and heavy metal cadmium) corroborated the findings. The prediction results of protein–protein interactions showed that its function is largely through a single pathway, mainly through interaction with snRNP proteins.
U2B″ is highly conserved in the plant kingdom, functions in the nucleus and cytoplasm, and participates in a wide range of processes in plant growth and development.
In contrast to the model of the 1 gene-1 polypeptide in prokaryotes, eukaryotic genes contain exons and introns. They generate mature mRNAs via splicing, which is an important post-transcriptional regulatory link during pre-mRNA splicing . In the splicing process, the U1 snRNP identifies and binds to the 5'-terminal splicing site of the intron . Subsequently, U2 cofactor identifies the 3' terminal splicing site to make the U2 snRNP combine with the branching point and form a pre-spliceosome (complex A). Afterwards, U4/U6 and U5 snRNPs bind to complex A and form complex B. Subsequently, the U4/U6 snRNP dissociates, and the U2 snRNP combines with the U6 snRNP to form a catalytic centre [3, 4]. Shortly afterwards, U1 and U4 snRNPs are removed from the complex, and U6 snRNPs bind to the 5' terminal splicing site to form a complex C-spliceosome . Finally, a splicing reaction occurs. Different transcripts are produced by alternative splicing (AS) and are eventually translated into different proteins . In eukaryotes, AS is an important channel of protein diversification and gene expression regulation and includes five modes. In addition, polyadenylation sites and variable selection of promoters can increase mRNA diversity . Splicing is regulated by trans-acting factors, splicing regulatory elements (SREs) and RNA secondary structures . SREs can be classified into silencers and enhancers. Silencers include intron splicing silencers (ISSs) and exon splicing silencers (ESSs), and enhancers include intron splicing enhancers (ISEs) and exon splicing enhancers (ESEs) . Most activating proteins that bind to ISEs and ESEs are splicing regulatory (SR) proteins, including one or two RNA recognition motifs (RRMs) and one SR domain containing abundant arginine and serine. RRM recognised by two amino acid sequences (RNP1 and RNP2) is the most commonly used RNA-binding domain in eukaryotes . A typical view of RNA–RRM interactions is that single-stranded RNA binds to a β-folded surface. Good electrostatic interactions, hydrogen bonding, and stacking between RNA bases and aromatic residues located in the RNP motif are considered major factors involved in RNA binding . Splicing is also regulated by the trans-acting factor-spliceosome, which includes five small nuclear RNAs and more than 100 core proteins .
Spliceosomes can be divided into two categories: major and minor spliceosomes . Major spliceosomes, including U1, U2, U4, U5, and U6 snRNPs , remove introns whose 5' termini have a “GU” splice site and whose 3' termini have an “AG” splicing site. Minor spliceosomes , including U5, U11, U12, U4, and U6 , are splice introns whose 5' termini have an “AU” splice site and whose 3' termini have an “AC” splicing site. The spliceosome is a multi-subunit RNA–protein complex. Currently, the identified splicing-related proteins are divided into five types: major snRNP proteins, splicing factors, splicing regulation factors, novel spliceosome proteins, and possible splicing-related proteins. Major snRNP proteins are divided into seven protein families: Sm core proteins, U1 snRNP-specific proteins, 17S U2 snRNP-specific proteins, U5 snRNP-specific proteins, U4/U6 snRNP-specific proteins, Tri-snRNP-specific proteins, and 18S U11/U2 snRNP-specific proteins. The major types of spliceosome-U2-dependent proteins contain several guanosine-rich U-snRNPs (U1, U2, U5, and U4/U6) and an array of non-snRNP proteins. U4 and U6 typically form complexes. Each snRNP type , including a conserved U-snRNA, binds to the ring containing seven Sm core proteins (B/B″, D3, D2, D1, E, F, and G) and a specific protein of each snRNP. Sm proteins are highly conserved RNA-binding proteins classified into seven Sm core proteins and like-Sm (LSm) proteins. The heptamer ring formed by the seven core Sm proteins is the basic structure that binds most snRNPs. Only U6 snRNP binds to LSm2-8 to form heptamers .
Recent research has shown that the structure of human 17S U2 snRNP is consistent with the electron microscope structure of previously isolated U2 snRNP [19, 20] and its entire structure in the human spliceosome [21,22,23,24]. Human 17S U2 snRNP presents a bipartite three-dimensional (3D) structure, as observed through high-resolution cryo-electron microscopy. The 17S U2 snRNP is located in the nucleus and is the active form of the U2 snRNP. It binds to the branching point of the pre-mRNA during the splicing assembly. The U2 snRNP plays an important role in selecting the mRNA precursor branching site adenosine, a nucleophile, in the first step of splicing . The mature 17S U2 snRNP includes 12S U2 snRNP formed by U2 snRNA, Sm protein, U2A', and U2B″ , and two splicing factors (SF3a and SF3b) . The first step in 17S U2 snRNP formation involves the Sm protein complex and U2A'-U2B″ dimer binding to the Sm site and stem-loop IV of U2 snRNA, respectively, forming 12S U2 (Fig. S1). In addition, the binding of U2B″ to the corresponding sites requires the participation of U2A', and the specific binding of U2B″ may be attributed to this common binding. In vitro, U2B″ is able to bind human stem-loop IV, Drosophila U2 snRNA stem-loop IV, and stem-loop II of human U1 snRNA. Yeast U2 snRNA nucleotides that form base pairs with branching sites were initially isolated in a fulcrum-interacting stem-loop ; however, it is not clear whether human U2 snRNA folds in a similar manner. The second step is to combine SF3b with 12S U2 to form 15S U2 snRNP and then combine with SF3a to form functional 17S U2 snRNP . As a spliceosomal protein, U2B″, which is an important tool for studying the evolution of RNA-binding patterns and RNA-binding specificity, belongs to the U1A/U2B″/SNF family. Although a high degree of sequence and structural conservation has been found in the family , proteins from modern members of this family have unique RNA-binding properties . Three proteins recognise their RNA targets using RNA recognition motifs (RRMs); however, their N-terminal RRMs differ from most RRMs in their high affinity and specificity for RNA. Notably, these homologous proteins contain two RRMs, whereas yeast contains only one . These phenomena were probably caused by evolution. The RRMs of the three proteins showed a high degree of conservation, similar to the two RNA stem loops that they recognised. When comparing modern RRMs, it has become apparent that each has unique RNA-binding properties, while U1A binds exclusively to U1 stem-loop II and U2B″ binds exclusively to U2 stem-loop IV in humans . Although they are components of snRNPs, how U1A and U2B″ participate in splicing remains unclear; in fact, snRNP recombination had no effect on splicing in the absence of U1A in vitro . Other experiments used fly homologous SNF mutations to exclude the protein from either U1 or U2 snRNP, resulting in relatively mild phenotypic consequences . In contrast, mutants that simultaneously knock out U1A and U2B″ in C. elegans are embryonically lethal, similar to knocking out SNF in Drosophila . It is possible that these proteins have other functions. In worms, U1A and U2B″, which are functionally redundant, interact with various snRNPs. The knockout of these two genes is necessary for a lethal phenotype. However, the absence of one does not lead to a lethal phenotype, and the other replaces the missing phenotype and binds to its corresponding snRNP. U2B″ homologous genes have been identified in vertebrates, yeast, and plants. Although the genes play an important role in the formation of the spliceosome in yeast, their mechanisms of action remain unclear. Analysis of the C. elegans genome has revealed two members of the U1A/U2B″ family. Notably, the two genes are present in a single operon. None of them has unique sequence features . The phenomenon is difficult to explain from a molecular perspective, especially based on protein structure. Therefore, analysing the function of proteins from an evolutionary perspective could provide a novel avenue for enhancing our understanding protein function.
Especially in plants, little is known about the specific functions of U2B″. Consequently, there is a growing interest in the elucidation of plant U2B″ proteins and how they mediate U2 snRNP, which could also involve numerous U2B″-interacting proteins. The present study provides a phylogenetic description of plant U2B″ and provides a clue for future functional studies in plant AS regulation. In this work, we studied the U2B″ gene family across different plants and analysed their phylogenetics, gene and protein structures, spatiotemporal gene expression profiles under different stimulation conditions, and subcellular localisation. The results of the present study will provide basic information on the phylogeny, structure, and expression of the gene family and lay a foundation for further functional identification of plants. The design process for the entire study is illustrated in Fig. S2.
Phylogenetic analysis of plant U2B″ gene family
To identify U2B″ genes in different plant species, we conducted a BLAST search using the Arabidopsis U2B″ protein (AT2G30260) against the Phytozome database. After filtering the sequences, 117 putative U2B″ sequences from 80 different plant species were obtained. The sequences were divided into five subfamilies (Fig. S3), including 45 dicotyledons and two monocotyledons in blue, 17 monocotyledons in pink, one fern in white, four bryophytes in green, and nine algae in yellow. Regarding the distribution of the number of genes in the 80 species, 50 species contained one gene, 24 species contained two genes, five species contained three genes (Gossypium hirsutum, Malus domestica, Medicago truncatula, Populus trichocarpa, Salix purpurea), and one species (Triticum aestivum) contained four genes (Table S1).
Gene copy numbers and species were mapped to a phylogenetic tree (Fig. S3). There were 45 dicotyledons, including 23 species with a single copy of the gene, 17 species with a double copy, and five species with three copies, and two monocotyledons, including one species with a single copy and one species with a double copy in the blue subfamily. In addition, the pink subfamily contained only monocotyledons, with a total of 17 species, including 12 species with a single copy, four species with a double copy, and one species with four copies. The white subfamily contained only one fern with a single copy. The green subfamily includes three bryophytes and three monocotyledons. Monocotyledons were all single copies, whereas bryophytes included two species with single copies and one species with a double copy. The yellow subfamily consisted of nine algae and contained eight species with single copies and one species with a double copy.
The 117 sequences from 80 plant species could provide us with a more comprehensive view of the evolutionary and developmental relationships among U2B″gene family. The clear topological structure and high bootstrap values (red) supporting each branch indicate the validity of this phylogenetic reconstruction of the U2B″ gene family. In addition, the yellow subfamily constitutes the basic part of the phylogenetic tree, which is far from other subfamilies, indicating that it is significantly associated with U2B″ in other plants. The green subfamily represents bryophytes, and the white subfamily of ferns is closely related to the yellow subfamily, which represents algae. However, the pink subfamily, representing monocotyledons, is closely related to the blue subfamily, representing dicotyledons, and is located far from the three subfamilies mentioned above. Multiple sequences from the same species were closely clustered in the same subfamily. Notably, each subfamily contained not only the plants it represented but also plants represented by other subfamilies. For example, the green subfamily contained not only bryophytes but also monocotyledons (Ananas comosus, Spirodela polyrhiza, Zostera marina). The blue subfamily, which represents dicotyledons, also contained two monocotyledons (Ananas comosus, Musa acuminata). The phenomena may be attributed to the degrees of plant development and evolution. Excluding the example above, the subfamilies strictly corresponded to the plants they represented.
Analysis of gene structure
To understand the general function of plant U2B″, it is necessary to study its genomic structure and its conserved gene motifs. For visualisation, the gene structure and corresponding motifs were linked in a phylogenetic tree (Fig. 1). Surprisingly, more than 80 sequences showed a 5 exon-2 UTR structure (Fig. 1, right panel). The proportion was close to 70%, indicating that they are highly conserved in the plant genomic structure across plant genomes, which suggests strong functional conservation of U2B″. Most of the remaining sequences had 3–6 exons, with a few exceptions. For example, Chlamydomonas contains seven exons. It should be noted that there were a few sequences with a single intron, or even no introns. For example, one of the two sequences of Arabidopsis thaliana had two introns, and the other had only one intron, whereas one of the two sequences of Capsella rubella had two introns and the other had no introns. In addition, the gene structures of the yellow, green, and white subfamilies were smaller, the pink subfamily was larger, and the blue subfamily was different in size, and there was no obvious trend.
Although gene structure variation was minimal among U2B″, we determined whether motif composition in their cDNA sequences reflected any differences. Further investigation of conserved DNA motifs showed that 89% of the 104 sequences had similar sequence characteristics and contained 8–10 of the top 10 identified motifs (Fig. S4A) and were located at similar sequence positions (Fig. 1, middle panel). The remaining 13 sequences, whose motifs were less than eight in number, were evenly distributed in all subfamilies except the green subfamily. Therefore, the gene structure may have a subtle relationship with conserved motifs.
Analysis of protein structure
Conserved domains and motifs were analysed for phylogenetic tree construction (Fig. 2). Different subfamilies of proteins were highly conserved (Fig. S5). A total of 106 peptides were predicted to contain an N-terminal domain named RRM1 and C-terminal domain named RRM2 (Fig. 2, right panel). In addition, there were eight peptides with only one RRM2, and one peptide with only one RRM1. This may lead to a decline in RNA recognition capacity. It is worth noting that two peptides (Prunus persica, Salix purpurea) did not have RRMs, possibly because they had more than one copy of the gene. Although some domains are lost during development and evolution, other peptides can perform similar functions.
Since the differences in protein domains were very minor, we attempted to explore differences in the motifs (Fig. 2, middle panel) of the protein sequences, with obvious differences identified. The numbers of motifs in the yellow, green, and white subfamilies were lower than those in the pink and blue subfamilies. The differences could be related to the different degrees of development and evolution among different subfamilies.
Homology modelling and amino acid conservation estimation
Since the protein domains were highly conserved, we estimated amino acid conservation. First, amino acid multiple sequence alignment was performed (Fig. S5). As expected, there was a high degree of conservation. The yellow, green, and white subfamilies differed from the other two subfamilies. The number of red amino acids was obviously greater, and the degree of red was greater, indicating that the yellow, green, and white subfamilies were closer together, whereas the pink and blue subfamilies were closer to one another, consistent with the results of the phylogenetic analysis. However, the two RRM domains of U2B″ also showed differences. The amino acid sequences of the RRM1 domain were highly conserved, whereas those of RRM2 were not as highly conserved. This may be related to the functions of the proteins, even though they all recognise the RNA. In addition, to demonstrate this point in more detail, we selected some plants (A. thaliana, Zea mays, Dunaliella salina, Sphagnum fallax) and plotted a two-dimensional map of their U2B″ proteins, which also revealed a high degree of conservation (Fig. 3A). Generally, the RRM1 domain is located between the 12th and 83rd amino acids, while the amino acid range of the RRM2 domain varies slightly between the 154th and 230th amino acids. Both RRM domains belong to the cl17169 superfamily (SF) and all proteins have disordered structures between the two domains. Furthermore, two had a coiled-coil structure.
A 3D model of plant U2B″ was constructed based on this template (Fig. 3B). The results of the conservative estimation of the 3D model were consistent with those of the amino acid multiple sequence alignment. RRM1 showed high quality, whereas RRM2 showed low quality.
To further study the potential expression profiles of plant U2B'', the 1500-bp sequence upstream of the genes was analysed using PlantCARE. From 117 sequences, 12,885 elements were identified, including 1162 blank elements and 1339 unnamed elements, accounting for 9% and 10% of the total, respectively. Subsequently, the remaining elements were screened, and 1731 elements, accounting for 13% of the total, were removed. Finally, 8653 elements, accounting for 68% of the initial total, were included in the analysis. First, we classified the elements into four categories according to their function (Fig. S6A, Table S2), hormone response (879 elements), light response (541 elements), regulation of basic transcription (5932 elements), and stress response (1186 elements). An additional element was not clearly defined (115 AAGAA-motif). Hormone-responsive elements included 252 ABA-responsive elements (ABREs), 177 as-1 motifs, 177 CGTCA motifs, 96 estrogen response elements (EREs), and 177 TGACG motifs. Notably, three of five motifs had the same number, although whether they corresponded to each other remains to be further studied. The light-responsive elements comprised 162 BOX4, 288 G-box, and 91 TCT motifs. The most abundant regulatory elements of basic transcription included three common elements:2848 TATA boxes, 2768 CAAT boxes, and 316 AT-TATA boxes. The remaining stress-responsive elements were divided into four categories:143, 523, 316, and 204.
We then combined the distribution map of elements (Fig. 4, right panel) on the sequence with the phylogenetic tree (Fig. 4, left panel). Excluding the white subfamily, the distribution of various elements in different subfamilies was relatively uniform. In general, the blue and pink subfamilies were distributed more than the other two subfamilies (Table S3). Owing to the considerable differences among family members, the number of promoters in the blue and pink subfamilies is also discrete. This may be due to differences in the degrees of evolution and development of different subfamilies. Notably, the numbers of elements in some sequences (Oropetium thomaeum, Ostreococcus lucimarinus, Micromonas sp. RCC299) were significantly lower than those in other sequences, indicating that they can only respond to a small number of stimuli. By contrast, they would have a single means of regulating gene expression. Subsequently, we performed an enrichment analysis of the stimuli of the sequence response (Fig. 5 and Fig. S6B). Among the 117 sequences, 112 responded to hormones, 113 to light, and 115 to stress, accounting for 96%, 97%, and 98%, respectively. Among them, 109 sequences (Fig. S6A) responded to the three stimuli simultaneously, exhibiting very high correlation. Notably, each group had two sequences that could respond to two of the three stimuli. Surprisingly, three sequences from three different plants (Amaranthus hypochondriacus, Arabidopsis lyrata, Populus trichocarpa) could only respond to stress. A specific analysis was then performed (Fig. 5). Surprisingly, only four sequences contained all elements at the same time. However, there were 49 sequences with 9–11 elements and 15 sequences with 5–8 elements simultaneously. Overall, the findings suggest that the U2B″ can respond to diverse stimuli and has a complex response network.
Knowledge of the subcellular localisation of proteins is also essential for the understanding of their functions. To further understand the function of U2B″, we performed subcellular localisation analysis to identify where U2B″ functions inside the cell (Fig. 6). First, we used web tools to predict the subcellular localisation of U2B″ in several species. The results were intriguing. In lower plants, U2B″ tends to be located in the cytoplasm, whereas in higher plants, U2B″ tends to be located in the nucleus (Fig. 6A). This means that over the course of evolution, the functions of U2B″ may have diverged based on its distinct positions. Considering that it is a component of the spliceosome and that splicing is essentially completed in the nucleus, we speculated that it may be located in the nucleus. To test this, we selected rice, which is one of the world’s most important crops. NLS-mCherry (a nuclear localisation marker) was selected as the control, to observe whether U2B″ was present in the nucleus. As expected, both were located in the nucleus. U2B″ and mCherry exhibited strong fluorescence, and a yellow signal was observed when the signals were merged (Fig. 6B). Notably, mCherry also had a strong localisation signal in the nucleolus, but U2B″ did not. This indicates that U2B″ plays a role in the nucleoplasm or nuclear membrane. Sequence analysis showed that the subsequences "KRKK" and "KKRR" in the junction structure between the two RRMs mainly promoted U2B″ to locate in the nucleus (Fig. 6C). The two RRMs are both disordered structures, suggesting that the structure could influence the positioning behaviour of U2B″. In addition, U2B″ has cytoplasmic localisation signals. Analysis of the transmembrane domain showed that U2B″ was within the cell, indicating that it had no localisation signal on the cell membrane (Fig. 6D).
Expression of plant U2B″ genes
qPCR was performed to further study the levels of expression of U2B″ in plants. Based on the results of promoter analysis, we selected four stress treatments: cold, drought, salt, and heavy metal cadmium (in the text referred to as Cd or Cd2+) (Fig. 7, Table S4). In general, excluding Cd stress, the levels of expression in the aboveground parts were higher than those in the belowground parts. The results suggest that U2B″ has contrasting regulatory mechanisms in the roots and shoots. Notably, there were not only differences in the levels of expression in the two parts but also contrasting expression patterns. U2B″ was not sensitive to Cd2+ and drought in roots, but showed significant differences in the shoot. U2B″ expression began to increase significantly after 3 h of drought treatment, which was sustained for 6 h. At 6 h, the relative expression was twice that of the control. However, it began to decline after 12 h of drought stress treatment. In contrast, the level of expression was significantly downregulated after 3 h of Cd treatment, and the relative expression was only 59% that of the control. However, it began to increase after 6 h and reached a level of significant difference, which was 77% of that of the control. After 12 h of Cd treatment, it began to decrease and finally reached the same level as 3 h, which was 60% of the control. When rice was treated with cold and salt, both roots and shoots exhibited sensitivity. Although they all exhibited sensitivity to stress, their response trends were distinct. Root and shoot exhibited similar trends under the salt treatment; both decreased first (40% and 81% of the control, respectively), increased (85% and 111% of the control, respectively), and then declined (56% and 95% of the control, respectively). However, the difference was that U2B″ was more responsive in the root. The difference between the salt treatment and the control was at most 60% in the root, whereas in the shoot, it was only 19%. In addition, U2B″ showed contrasting response patterns in the roots and shoots when subjected to cold stress. The levels of expressions were the lowest in the roots after 6 h of chilling treatment, and the highest in shoots after 3 h of chilling treatment.
Alternative splicing profile analysis and splicing isoforms
Four representative monocotyledons, four dicotyledons, four algae, and three bryophytes were selected for this study. By comparing their gene structures, we observed that AS was not common in the U2B″ families (Fig. S7). AS occurred three times in the selected monocotyledons, twice in dicotyledons, and once in bryophytes; however, it did not occur in algae. In addition, more than two occurrences of AS were non-existent in the plants. Among the 15 selected plants, five plants had AS, including twice in Arabidopsis, while the other four had single occurrences. In Z. mays transcripts, only one of the five exons was preserved after AS. In the other four splicing isoforms, there was no significant difference between the gene structure and representative transcript forms, with some sequence differences observed in the UTRs. Overall, U2B″ may not undergo large splicing changes and may act uniformly on substrates with conserved RRM motifs across isoforms.
U2B″ needs to be translated into a protein to perform its function. We selected three common plants (Arabidopis thaliana, Oryza sativa, Z. mays) and constructed an interaction network (Fig. 8, Table S5) based on data obtained from STRING. We obtained 9, 10, and 10 interacting proteins in A. thaliana, O. sativa, and Z. mays. The numbers of interacting proteins showed clear consistency, suggesting that the U2B″ protein has a fairly complex interaction network in different plants. Regarding the genes encoding the corresponding proteins, there was an interesting phenomenon in which genes in Z. mays are evenly distributed on chromosomes, whereas those in A. thaliana and O. sativa were concentrated on a few chromosomes. Furthermore, among the nine proteins interacting with AT2G30260, there were nine snRNP proteins, one RNA-binding protein, and one putative splicing factor. Excluding a protein in Z. mays that has an unknown function, there were eight and seven snRNP proteins in O. sativa and Z. mays, respectively. In addition, both O. Sativa and Z. mays had an Sm-like protein, and it should be noted that the interaction protein in Z. mays also contains a zinc ion binding protein. In terms of the types of interacting proteins, the U2B″ protein showed a high degree of consistency in different plants. Notably, each plant has a unique interaction protein that differs from that of other plants. This may be caused by different degrees of evolution and development.
Structural analysis of the spliceosome has long been considered one of the most promising research fields in structural biology. Several human diseases can be attributed to incorrect splicing of genes or regulation of spliceosomes. In addition, 35% of genetic disorders in humans are caused by the variable splicing of single genes caused by gene mutations . Some diseases are caused by mutations in splice proteins that affect the splicing of many transcripts. Some cancers are also involved in incorrect regulation of splicing factors [38, 39]. In a previous study, a small percentage of medulloblastoma samples had the same non-coding mutations. Following a preliminary study, the authors were surprised to find that this mutation (A > G) affected more than a small number of people, as nearly 100% of the adult patients with a certain disease subtype had the mutation. In addition, another mutation (A > C) in the same genome position is mainly present in liver cancer and chronic lymphoid leukaemia. Surprisingly, the mutation is located in the RNA-binding site of an important splicing RNA (U1 snRNA). Although at present, research on splicing mainly focuses on humans, there are still many reports related to plants. Some studies have found that the deletion or mutation of the mu1a gene of Magnaporthe (M. oryzae U1A, MU1A) may lead to abnormal splicing of the precursor mRNA, affecting the normal expression of some proteins, normal growth and metabolism processes, and ultimately, the pathogenicity of Magnaporthe. In the defence response to fungal infection, AS may be as important as traditional transcriptional control in C. sublineola-inoculated sorghum seedlings, as reported in a study based on next-generation sequencing technologies . To screen for pathogenic effectors that regulate plant AS, a fluorescence reporting system was established to screen nine splicing regulation (SREs) effectors from 87 effectors of Phytophthora infestans . Further studies have shown that SRE3 in combination with U1-70 K physically manipulates the plant AS mechanism, thereby regulating AS-mediated plant immunity. The biological clock of Arabidopsis not only controls gene transcription but also influences its post-transcriptional regulation by influencing AS . Six differential AS events occurred in the defence response of germs in mock-inoculated and S. sclerotiorum-inoculated susceptible and tolerant B. napus plants, as determined by an analysis of 18 RNA-seq libraries . However, the potential mechanism via which splice proteins affect splicing regulation remains unclear. The work described in this study provides an introductory layout for such works in future.
In the present study, we identified 117 U2B″ genes in 80 plant species. In the phylogenetic tree, plants with higher degrees of evolution and development were closely clustered together, and plants with lower degrees of evolution and development were similarly closely clustered together. Surprisingly, U2B″ in different plants showed consistently high conservation in multiple analyses. Nearly 70% of the sequences showed a 5 exon-2 UTR structure, and 90% of the sequences were predicted to contain an N-terminal domain named RRM1 and a C-terminal domain annotated as RRM2. Generally, if a gene is highly conserved in many species, it is believed that these similarities across species indicate that the gene performs some basic functions necessary for many life forms, and therefore retains these sequences during evolution. These sequences are generally necessary for life activities, and mutations in these sequences often lead to death. Therefore, few mutations have been preserved during evolution, so that the genes are highly conserved. The same is true for U2B″, which may perform one or more functions essential for life. We prefer U2B″ to perform a single or a few functions essential for life. U2B'' knockout did not lead to a lethal phenotype in worms . U1A replaces snRNP and binds to its corresponding snRNP. However, we found that the N- and C-terminal domains of the protein were different in the multisequence alignment. Although both domains can recognise RNA, the N-terminus is more conserved than the C-terminus. We speculate that RRM1 performs the most basic functions required by all species, whereas RRM2 performs distinct functions in evolution and development.
The U2 snRNP is a subcomplex in early spliceosomal assembly. U2B″ plays an important role in the snRNP complex. However, very few studies have elucidated its biological function in plants. Promoter analysis revealed that according to the complex response elements, U2B″ could respond to various stimuli. First, many stress-responsive elements were identified, including AREs, MYBs, and STREs. As a necessary cis-regulatory element for anaerobic induction, AREs play an important role in plants. MYB-MYC is a plant-specific element widely involved in plant growth, development, and response to abiotic stress. STRE is reportedly involved in a variety of stress responses in Neurospora, including heat, osmotic, and oxidative stress, and the molecular mechanism mediated by STRE may not be conserved . In addition to many stress response components, there are a large number of hormone-responsive elements upstream of U2B″, including ABRE, activating sequence-1, CGTCA motif, ERE, and the TGACG motif. ABREs are response elements of ABA  that play a variety of important roles in plants. In higher plants, AS-1 of the cauliflower mosaic virus 35S promoter mediates SA- and IAA-inducible transcriptional activation . As MeJA-responsive elements, the CGTCA and TGACG motifs are found in a large number of plants. The estrogen response element (ERE) is a conserved DNA sequence in the promoters of the estrogen target genes. It can bind to estrogen receptors, is transcriptionally regulated, and is usually located in the promoters of its target genes . These results indicate that U2B″ is involved extensively in plant growth, development, and stress responses. We also confirmed this using qPCR. However, the interaction protein of U2B″ is relatively simple. These include the major snRNPs common to the three plants and other minor proteins specific to each plant. This may be caused by different degrees of evolution and development. Overall, the above studies show that U2B″ can respond to a variety of external stimuli but can only interact with a limited number of proteins.
In the present study, 117 U2B″ genes were identified in 80 plant types. Comprehensive bioinformatic analysis and partial experiments of the gene family were performed. Considering that U2B″ is highly conserved in different plants and is located in the nucleus, it can respond to a variety of external stimuli, indicating that it is extensively involved in plant growth and development. However, the mechanism by which U2B″ regulates plant growth and development requires further study.
Identification of plant U2B″
The A. thaliana U2B″ protein (AT2G30260) was used as the original sequence for carrying out a protein BLAST search with an e-value cutoff of 1e-10  against 80 plant genome sequences from Phytozome v12.1. Consequently, 117 putative U2B″/U2B″-like sequences were identified for subsequent research.
Protein sequences of the plant U2B″ identified above were obtained for phylogenetic tree construction to infer clustering patterns and evolutionary relationships. The transcript with the longest coding sequence was selected for loci with multiple splicing subtypes. Subsequently, multiple sequence alignment of proteins were performed using Muscle V3.8 with the default settings . A phylogenetic tree was constructed using MrBayes v3.2.2 (Jones model) .
Analysis of gene and protein structure
The structural information of genes and proteins was obtained from the Phytozome v12.1 and Pfam protein family databases. The top 10 detected motifs were obtained using default settings on the MEME server from the cDNA and protein sequences of plant U2B″ .
Promoter motif prediction and the protein–protein interaction network
The 1.5-kb promoter sequences of plant U2B″ were selected from the Phytozome database and used for the prediction of cis-elements using PlantCARE . Arabidopsis U2B″ protein was input to STRING and PlantSPEAD  to find the most interacting partners.
The CDS of OsU2B″ (LOC_Os03g18720.1) was amplified using gene-specific primers (5'-CCTGTTGTTTGGTGTTACTTAAGCTTATGTTGTCCGGCGACATACC-3' and 5'- TCCTCGCCCTTGCTCACCATGGATCCTCACTTCTTTGCGTAGGATATAGCC-3') and inserted into the GFP (pGreen II-UBI-GFP) vector to generate a fusion construct according to Lu et al. . NLS − mCherry (SV40 large − T antigen NLS, a signal peptide guiding GFP into the nucleus) was selected as a marker for nuclear localisation . Nipponbare rice seedlings grown at approximately 30 °C for one week were selected, and their stems were collected to extract protoplasts. Equal volumes of the constructed plasmids and extracted protoplasts were mixed with 40% PEG4000 solution at a ratio of 1:5 (v/v), cultured overnight at 30 °C, and observed under a laser confocal microscope (Leica TCS SP8) more than three times. The subcellular localisation of U2B″ in some species was predicted using the online PSORT database. NucPred was used to predict the nuclear localisation of protein sequences . DeepTMHMM was used to predict the protein transmembrane domains.
Plant growth and RT-qPCR analysis
Nipponbare seeds were first soaked in carbendazim for one day and then in water for another day. The seeds were then transferred to the climate chamber and subjected to stress when they grew two true leaves (22 d). Materials were grown at 22 °C/20 °C under a 16 h/ 8 h light/dark photoperiod and a light intensity of 120 μmol m −2 s −1 in an incubator (QY-14; Nanjing Quanyou Electronic Technology Co., Ltd., China). Hydroponics was used in all four stress treatments, including cold treatment at 8 °C and drought treatment, with 20% PEG6000, salt treatment with 100 mmol L−1 sodium chloride, and heavy metal cadmium treatment with 100 μmol L−1 cadmium sulfate. Total RNA was isolated using TRIzol reagent (Invitrogen) and converted to complementary DNA (cDNA) using PrimeScript™ RT Master Mix (TAKARA) according to the manufacturer’s protocol. RT-qPCR was performed using a 7500 Real-time PCR Detection System (Bio-Rad) in conjunction with SYBR Fast qPCR Mix (TAKARA) and was repeated at least three times. Forward and reverse primers of OsU2B″ were used to produce a single amplification (5'- GCAACCGAAGATGGTTCTACTG -3' and 5'- GCTTTGGGCAGCAGCATTAG -3'). Gene expression was presented as a heat map using TBtools v1.048 . OsACTIN (LOC_Os05g36290) was used as a reference gene.
Model construction and amino acid conservation assessment
Homology modelling was carried out using the A. thaliana U2B″ protein sequence (UniProtKB AC: O22922) based on the template 4pkd.1. B using SWISS-MODEL [57, 58]. Amino acid multiple sequence alignment and conservation scores based on frequency-based differences were obtained using the default values from NCBI.
Availability of data and materials
All the data generated during this study are included in this published article and its supplementary files. The data used and/or analyzed during the current study are available from the Phytozome (https://phytozome-next.jgi.doe.gov/pz/), PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/), STRING (https://string-db.org/), PlantSPEAD (http://chemyang.ccnu.edu.cn/ccb/database/PlantSPEAD/index.php/home/index), PSORT (https://www.genscript.com/psort.html), DeepTMHMM (https://dtu.biolib.com/DeepTMHMM), UniProtKB (https://www.uniprot.org/), SWISS-MODEL (https://swissmodel.expasy.org/) and NCBI (https://www.ncbi.nlm.nih.gov/) website.
Estrogen response element
Exon splicing enhancer
Intron splicing enhancer
RNA recognition motifs
Splicing regulatory elements
U2 small nuclear ribonucleoprotein B
Lu S, Gao C, Wang YZ, He YY, Du JR, Chen MX, Zhao H, Fang H, Wang BH, Cao YY. Phylogenetic analysis of the plant U2 snRNP auxiliary factor large subunit a gene family in response to developmental cues and environmental stimuli. Front Plant Sci. 2021;12: 739671.
Hilleren PJ, Kao HY, Siliciano PG. The amino-terminal domain of yeast U1–70K is necessary and sufficient for function. Mol Cell Biol. 1995;15:6341–50.
Urlaub H, Hartmuth K, Kostka S, Grelle G, Luhrmann R. A general approach for identification of RNA-protein cross-linking sites within native human spliceosomal small nuclear ribonucleoproteins (snRNPs) analysis Of RNA-protein contacts in native U1 and U4/U6. U5 snRNPs. J Biol Chem. 2000;275(52):41458–68.
Gornemann J, Kotovic KM, Hujer K, Neugebauer KM. Cotranscriptional spliceosome assembly occurs in a stepwise fashion and requires the cap binding complex. Mol Cell. 2005;19:53–63.
Lacadie SA, Rosbash M. Cotranscriptional spliceosome assembly dynamics and the role of U1 snRNA: 50 ss base pairing in yeast. Mol Cell. 2005;19:65–75.
Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M. Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 2010;152(4):1787–95.
Tazi J, Bakkour N, Stamm S. Alternative splicing and disease. Mol Basis Dis. 2008;1792:14–26.
Uhart M, Sirand-Pugnet P, Labarère J. Evolution of mitochondrial SSU-rDNA variable domain sequences and rRNA secondary structures, and phylogeny of the Agrocybe aegerita multispecies complex. Res Microbiol. 2007;158(3):203–12.
Vuong CK, Black DL, Zheng S. The neurogenetics of alternative splicing. Nat Rev Neurosci. 2016;17(5):265–81.
Maris C, Dominguez C, Allain FHT. The RNA recognition motif, a plastic RNA-binding platform to regulate post-transcriptional gene expression. FEBS J. 2005;272(9):2118–31.
Cléry A, Blatter M, Allain FH-T. RNA recognition motifs: boring? Not quite. Curr Opin Struct Biol. 2008;18(3):290–8.
Galej WP. Structural studies of the spliceosome: past, present and future perspectives. Biochem Soc T. 2018;46(6):1407–22.
Wahl MC, Will CL, Lührmann R. The spliceosome: Design principles of a dynamic RNP machine. Cell. 2009;136(4):701–18.
Kaida D, Berg MG, Younis I, Kasim M, Singh LN, Wan L, Dreyfuss G. U1 snRNP protects pre-mRNAs from premature cleavage and polyadenylation. Nature. 2010;468(7324):664–8.
König H, Matter N, Bader R, Thiele W, Müller F. Splicing segregation: the minor spliceosome acts outside the nucleus and controls cell proliferation. Cell. 2007;131(4):718–29.
Tidow H, Andreeva A, Rutherford TJ, Fersht AR. Solution structure of the U11–48K CHHC zinc-finger domain that specifically binds the 5’ splice site of U12-Type introns. Structure. 2009;17(2):294–302.
Strunov A, Boldyreva LV, Pavlova GA, Pindyurin AV, Gatti M, Kiseleva E. A simple and effective method for ultrastructural analysis of mitosis in Drosophila S2 cells. MethodsX. 2016;3:551–9.
Wahl MC, Lührmann R. SnapShot: spliceosome dynamics II. Cell. 2015;162(2):456–456.
Zhang Z, Will CL, Bertram K, Dybkov O, Hartmuth K, Agafonov DE, Hofele R, Hofele R, Urlaub H, Kastner B, Lührmann R, Stark H. Molecular architecture of the human 17S U2 snRNP. Nature. 2020;583(7815):310–3.
Krämer A, Grüter P, Gröning K, Kastner B. Combined biochemical and electron microscopic analyses reveal the architecture of the mammalian U2 snRNP. J Cell Biol. 1999;145(7):1355–68.
Bertram K, Agafonov DE, Dybkov O, Haselbach D, Leelaram MN, Will CL, Urlaub H, Kastner B, Lührmann R, Stark H. Cryo-EM structure of a pre-catalytic human spliceosome primed for activation. Cell. 2017;170(4):701–13.
Haselbach D, Komarov I, Agafonov DE, Hartmuth K, Graf B, Dybkov O, Urlaub H, Kastner B, Lührmann R, Stark H. Structure and conformational dynamics of the human spliceosomal Bact complex. Cell. 2018; https://doi.org/10.1016/j.cell.2018.01.010.
Zhan X, Yan C, Zhang X, Lei J, Shi Y. Structures of the human pre-catalytic spliceosome and its precursor spliceosome. Cell Res. 2018;28(12):1129–40.
Zhang X, Yan C, Zhan X, Li L, Lei J, Shi Y. Structure of the human activated spliceosome in three conformational states. Cell Res. 2018;28(3):307–22.
Will CL, Lührmann R. Spliceosome structure and function. Cold Spring Harb Perspect Biol. 2011;3(7):a003707.
Yun HJ, Kwon J, Seol W. Specific inhibition of transcriptional activity of the constitutive androstane receptor (CAR) by the splicing factor SF3a3. Biol Chem. 2008;389(10):1313–8.
Bosetti F, Seemann R, Rapoport SI. Chronic lithium chloride administration to rats decreases brain protein level of epsilon (ϵ) subunit of eukaryotic initiation factor-2B. Neurosci Lett. 2002;327(1):71–3.
Perriman R Jr, MA. Invariant U2 snRNA nucleotides form a stem loop to recognize the intron early in splicing. Mol Cell. 2010;38(3):416–27.
Krämer A, Ferfoglia F, Huang CJ, Mulhaupt F, Nesic D, Tanackovic G. Structure-function analysis of the U2 snRNP-associated splicing factor SF3a. Biochem Soc T. 2005;33(3):439–42.
Williams SG, Harms MJ, Hall KB. Resurrection of an urbilaterian U1A/U2B″/SNF protein. J Mol Biol. 2013;425(20):3846–62.
Delaney KJ, Williams SG, Lawler M, Hall KB. Climbing the vertebrate branch of U1A/U2B″ protein evolution. RNA. 2014;20(7):1035–45.
Saldi T, Wilusz C, MacMorris M, Blumenthal T. Functional redundancy of worm spliceosomal proteins U1A and U2B”. Proc Natl Acad Sci. 2007;104(23):9753–7.
Scherly D, Boelens W, Dathan NA, van Venrooij WJ, Mattaj IW. Major determinants of the specificity of interaction between small nuclear ribonucleoproteins U1A and U2B” and their cognate RNAs. Nature. 1990;345(6275):502–6.
Will CL, Rümpler S, Klein Gunnewiek J, van Venrooij WJ, Lührmann R. In vitro reconstitution of mammalian U1 snRNPs active in splicing: the U1-C protein enhances the formation of early (E) spliceosomal complexes. Nucleic Acids Res. 1996;24(23):4614–23.
Stitzinger SM, Conrad TR, Zachlin AM, Salz HK. Functional analysis of SNF, the Drosophila U1A/U2B" homolog: identification of dispensable and indispensable motifs for both snRNP assembly and function in vivo. RNA. 1999;5:1440–50.
Flickinger TW, Salz HK. The Drosophila sex determination gene snf encodes a nuclear protein with sequence and functional similarity to the mammalian U1A snRNP protein. Gene Dev. 1994;8(8):914–25.
Shi Y. Mechanistic insights into precursor messenger RNA splicing by the spliceosome. Nat Rev Mol Cell Bio. 2017;18(11):655–70.
Shuai S, Suzuki H, Diaz-Navarro A, Nadeu F, Kumar SA, Gutierrez-Fernandez A, Delgado J, Pinyol M, López-Otín C, Puente XS, Taylor MD, Campo E, Stein LD. The U1 spliceosomal RNA is recurrently mutated in multiple cancers. Nature. 2019;574(7780):712–6.
Suzuki H, Kumar SA, Shuai S, et al. (2019) Recurrent noncoding U1 snRNA mutations drive cryptic splicing in SHH medulloblastoma. Nature. 2019;574(7780):707–15.
Wang L, Chen M, Zhu F, Fan T, Zhang J, Lo C. Alternative splicing is a Sorghum bicolor defense response to fungal infection. Planta. 2020;251(1):14.
Huang J, Lu X, Wu H, Xie Y, Peng Q, Gu L, Wu J, Wang Y, Reddy ASN, Dong S. Phytophthora effectors modulate genome-wide alternative splicing of host mRNAs to reprogram plant immunity. Mol Plant. 2020;13(10):1470–84.
Yang Q, Yang B, Li J, Wang Y, Tao R, Yang F, Wu X, Yan X, Ahmad M, Shen J, Bai S, Teng Y. ABA-responsive ABRE-BINDING FACTOR3 activates DAM3 expression to promote bud dormancy in Asian pear. Plant Cell Environ. 2020;43(6):1360–75.
Ma JQ, Xu W, Xu F, Lin A, Sun W, Jiang HH, Lu K, Li JN, Wei LJ. Differential alternative splicing genes and isoform regulation networks of rapeseed (Brassica napus L.) infected with Sclerotinia sclerotiorum. Genes. 2020;11(7):784.
Freitas FZ, Virgilio S, Cupertino FB, Kowbel DJ, Fioramonte M, Gozzo FC, Glass NL, Bertoloni MC. The SEB-1 transcription factor binds to the STRE motif in Neurospora crassa and regulates a variety of cellular processes including the stress response and reserve carbohydrate metabolism. G3-Genes Genom Genet. 2016;6(5):1327–43.
Yang Y, Li Y, Sancar A, Oztas O. The circadian clock shapes the Arabidopsis transcriptome by regulating alternative splicing and alternative polyadenylation. J Biol Chem. 2020;295(22):7608–19.
Niggeweg R, Thurow C, Kegler C, Gatz C. Tobacco transcription factor TGA2.2 is the main component of as-1-binding factor ASF-1 and is involved in salicylic acid- and auxin-inducible expression of as-1-containing target promoters. J Biol Chem. 2000;275(26):19897–905.
Govind AP, Thampan RV. Proteins interacting with the mammalian estrogen receptor: proposal for an integrated model for estrogen receptor mediated regulation of transcription. J Cell Biochem. 2001;80(4):571–9.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;3(5):1792–7.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME Suite: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:W202–8.
Lescot M, Patrice Déhais P, Thijs G, Marchal K, Moreau Y, de Peer YV, Rouzé P, Rombauts S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.
Chen MX, Mei LC, Wang F, Boyagane Dewayalage IKW, Yang JF, Dai L, Yang GF, Gao B, Chen CL, Liu YG, Zhang JH, Hao GF. PlantSPEAD: a web resource towards comparatively analysing stress-responsive expression of splicing-related proteins in plant. Plant Biotechol J. 2020;9(2):227–9.
Kalderon D, Richardson WD, Markham AF, Smith AE. Sequence requirements for nuclear location of simian virus 40 large-T antigen. Nature. 1984;311:33–8.
Brameier M, Krings A, Maccallum RM. NucPred-Predicting Nuclear Localization of Proteins. Bioinformatics. 2007;23(9):1159–60.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.
Kondo Y, Oubridge C, Roon AMV, Nagai K. Crystal structure of human U1 snRNP, a small nuclear ribonucleoprotein particle, reveals the mechanism of 5’ splice site recognition. eLife. 2015;4:e04986.
Andrew W, Martino B, Stefan B, Gabriel S, Gerardo T, Rafal G, Florian HT, Tjaart APDB, Christine R, Lorenza B, Rosalba L, Torsten S. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(W1):W296–303.
We would like to thank Professor Jianhua Zhang of Hong Kong Baptist University for the tobacco and Nipponbare seeds and Editage (http://www.editage.com) for editing and reviewing this manuscript for English language.
This work was supported by the National Natural Science Foundation of China (32172104) and Large Instruments Open Foundation of Nantong University (KFJN2223). The funding agency did not participate in any aspects of research design, data analysis or manuscript writing.
Ethics approval and consent to participate
All methods using plant materials were conducted in accordance with the relevant guidelines and regulations.
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fig. S1. Formation of mature 17S U2 snRNP. The upper one is12S core particle, the middle one is 15S pre-mature particle, and the lower oneis 17S functional maturity particle. Continuous thin black lines represent U2 snRNA. Fig. S2. The design process of the whole article. Fig. S3. Circle phylogenetic tree representation of the available plant U2B″ gene family. Phylogenetic analysisof plant U2B″ gene family was carried out by using software MrBayes v3.2.2. The posterior probability values are labeled at each major branch. Blue for dicotyledons, pink for monocotyledons,white for ferns, green for bryophytes, and yellow for algae. Fig. S4. Motifs of genomic structure and protein structure analysis. (A) Consensus sequence of top ten identified DNA motifs are listed in ascending order. (B) Consensus sequence of top ten identified amino-acid motifs are listed in ascending order. Fig. S5. The multiple sequence alignment of RRM domains for the conservative analysis. The sequences are arranged from top to bottom in phylogenetic tree. Fig. S6. Promoter classification and enrichment analysis. (A) Statistics of motifs function and number. The x-axis represents the number of elements. (B)Overall enrichment statistics of motifs in response to stress, hormones and light.Fig. S7. AS profile analysis. Summary of annotated alternatively spliced transcript isoforms for identified U2B″ genes. Pink represents monocotyledons, blue represents dicotyledons, green represents bryophytes, and yellow represents algae.
Table S1. Sequence summary of plant U2B″ gene phylogenetic analysis.
Table S2. Specific data of enrichment analysis.
Table S3. Statistical analysis of promoter distribution.
Table S4. Expression of U2B″ in rice under different stresses.
Table S5. Sequence summary of plant U2B″ protein-protein interaction network.
About this article
Cite this article
Gao, C., Lu, S., Zhou, R. et al. Phylogenetic analysis and stress response of the plant U2 small nuclear ribonucleoprotein B″ gene family. BMC Genomics 23, 744 (2022). https://doi.org/10.1186/s12864-022-08956-0
- Splicing factor
- Subcellular localisation