Root Hair Single Cell Type Specific Profiles of Gene Expression and Alternative Polyadenylation Under Cadmium Stress

Time-Course Transcriptome Analysis of Aquilegia vulgaris Root Reveals the Cell Wall's Roles in Salinity Tolerance

Salt stress has a considerable impact on the development and growth of plants. The soil is currently affected by salinisation, a problem that is becoming worse every year. This means that a significant amount of salt-tolerant plant material needs to be added. Aquilegia vulgaris has aesthetically pleasing leaves, unique flowers, and a remarkable tolerance to salt. In this study, RNA-seq technology was used to sequence and analyse the transcriptome of the root of Aquilegia vulgaris seedlings subjected to 200 mM NaCl treatment for 12, 24, and 48 h. In total, 12 Aquilegia vulgaris seedling root transcriptome libraries were constructed. At the three time points of salt treatment compared with the control, 3888, 1907, and 1479 differentially expressed genes (DEGs) were identified, respectively. Various families of transcription factors (TFs), mainly AP2, MYB, and bHLH, were identified and might be linked to salt tolerance. Gene Ontology (GO) analysis of DEGs revealed that the structure and composition of the cell wall and cytoskeleton may be crucial in the response to salt stress. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs showed a significant enrichment of the pentose and glucuronate interconversion pathway, which is associated with cell wall metabolism after 24 and 48 h of salt treatment. Based on GO and KEGG analyses of DEGs, the pentose and glucuronate interconversion pathway was selected for further investigation. AP2, MYB, and bHLH were found to be correlated with the functional genes in this pathway based on a correlation network. This study provides the groundwork for understanding the key pathways and gene networks in response to salt stress, thereby providing a theoretical basis for improving salt tolerance in Aquilegia vulgaris.

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

Inhibitor AN3661 reveals biological functions of Arabidopsis CLEAVAGE and POLYADENYLATION SPECIFICITY FACTOR 73

Cleavage and polyadenylation specificity factor (CPSF) is a protein complex that plays an essential biochemical role in mRNA 3'-end formation, including poly(A) signal recognition and cleavage at the poly(A) site. However, its biological functions at the organismal level are mostly unknown in multicellular eukaryotes. The study of plant CPSF73 has been hampered by the lethality of Arabidopsis (Arabidopsis thaliana) homozygous mutants of AtCPSF73-I and AtCPSF73-II. Here, we used poly(A) tag sequencing to investigate the roles of AtCPSF73-I and AtCPSF73-II in Arabidopsis treated with AN3661, an antimalarial drug with specificity for parasite CPSF73 that is homologous to plant CPSF73. Direct seed germination on an AN3661-containing medium was lethal; however, 7-d-old seedlings treated with AN3661 survived. AN3661 targeted AtCPSF73-I and AtCPSF73-II, inhibiting growth through coordinating gene expression and poly(A) site choice. Functional enrichment analysis revealed that the accumulation of ethylene and auxin jointly inhibited primary root growth. AN3661 affected poly(A) signal recognition, resulted in lower U-rich signal usage, caused transcriptional readthrough, and increased the distal poly(A) site usage. Many microRNA targets were found in the 3' untranslated region lengthened transcripts; these miRNAs may indirectly regulate the expression of these targets. Overall, this work demonstrates that AtCPSF73 plays important part in co-transcriptional regulation, affecting growth, and development in Arabidopsis.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

Alternative Polyadenylation Is a Novel Strategy for the Regulation of Gene Expression in Response to Stresses in Plants

Precursor message RNA requires processing to generate mature RNA. Cleavage and polyadenylation at the 3'-end in the maturation of mRNA is one of key processing steps in eukaryotes. The polyadenylation (poly(A)) tail of mRNA is an essential feature that is required to mediate its nuclear export, stability, translation efficiency, and subcellular localization. Most genes have at least two mRNA isoforms via alternative splicing (AS) or alternative polyadenylation (APA), which increases the diversity of transcriptome and proteome. However, most previous studies have focused on the role of alternative splicing on the regulation of gene expression. In this review, we summarize the recent advances concerning APA in the regulation of gene expression and in response to stresses in plants. We also discuss the mechanisms for the regulation of APA for plants in the adaptation to stress responses, and suggest that APA is a novel strategy for the adaptation to environmental changes and response to stresses in plants.

Integration of transcriptome and metabolome analyses reveals sorghum roots responding to cadmium stress through regulation of the flavonoid biosynthesis pathway

Cadmium (Cd) pollution is a serious threat to plant growth and human health. Although the mechanisms controlling the Cd response have been elucidated in other species, they remain unknown in Sorghum (Sorghum bicolor (L.) Moench), an important C₄ cereal crop. Here, one-week-old sorghum seedlings were exposed to different concentrations (0, 10, 20, 50, 100, and 150 μM) of CdCl₂ and the effects of these different concentrations on morphological responses were evaluated. Cd stress significantly decreased the activities of the enzymes peroxidase (POD), superoxide dismutase (SOD), glutathione S-transferase (GST) and catalase (CAT), and increased malondialdehyde (MDA) levels, leading to inhibition of plant height, decreases in lateral root density and plant biomass production. Based on these results, 10 μM Cd concentration was chosen for further transcription and metabolic analyses. A total of 2683 genes and 160 metabolites were found to have significant differential abundances between the control and Cd-treated groups. Multi-omics integrative analysis revealed that the flavonoid biosynthesis pathway plays a critical role in regulating Cd stress responses in sorghum. These results provide new insights into the mechanism underlying the response of sorghum to Cd.

Review: Mechanisms underlying alternative polyadenylation in plants - looking in the right places

Recent years have seen an explosion of interest in the subject of alternative polyadenylation in plants. Connections between the polyadenylation complex and numerous developmental and stress responses are well-established. However, those that link stimuli with the functioning of the polyadenylation complex are less well understood. To this end, it is imperative to clearly delineate the roles of the polyadenylation complex in both plant growth AND alternative polyadenylation. It is also necessary to understand the ways by which other molecular processes may contribute to alternative polyadenylation. This review discusses these issues, with a focus on instances that reveal mechanisms by which mRNA polyadenylation may be regulated. Insights from from characterizations of mutants affected in the polyadenylation complex are discussed, as are the limitations of such characterizations when it comes to teasing out cause and effect. These limitations encourage explorations to other processes that are beyond the core polyadenylation complex. Two such processes that sculpt the plant transcriptome - transcription termination and the epigenetic control of transposon activity - also contribute to regulated poly(A) site choice. These subjects define "the right places" - molecular mechanisms that contribute to the wide-ranging control of gene expression via mRNA polyadenylation.

Root hair-specific transcriptome reveals response to low phosphorus in Cicer arietinum

Root hairs (RH) are a single-cell extension of root epidermal cells. In low phosphorus (LP) availability, RH length and density increase thus expanding the total root surface area for phosphate (Pi) acquisition. However, details on genes involved in RH development and response to LP are missing in an agronomically important leguminous crop, chickpea. To elucidate this response in chickpea, we performed tissue-specific RNA-sequencing and analyzed the transcriptome modulation for RH and root without RH (Root-RH) under LP. Root hair initiation and cellular differentiation genes like RSL TFs and ROPGEFs are upregulated in Root-RH, explaining denser, and ectopic RH in LP. In RH, genes involved in tip growth processes and phytohormonal biosynthesis like cell wall synthesis and loosening (cellulose synthase A catalytic subunit, CaEXPA2, CaGRP2, and CaXTH2), cytoskeleton/vesicle transport, and ethylene biosynthesis are upregulated. Besides RH development, genes involved in LP responses like lipid and/or pectin P remobilization and acid phosphatases are induced in these tissues summarizing a complete molecular response to LP. Further, RH displayed preferential enrichment of processes involved in symbiotic interactions, which provide an additional benefit during LP. In conclusion, RH shows a multi-faceted response that starts with molecular changes for epidermal cell differentiation and RH initiation in Root-RH and later induction of tip growth and various LP responses in elongated RH.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

Divergence in the Regulation of the Salt Tolerant Response Between Arabidopsis thaliana and Its Halophytic Relative Eutrema salsugineum by mRNA Alternative Polyadenylation

Salt tolerance is an important mechanism by which plants can adapt to a saline environment. To understand the process of salt tolerance, we performed global analyses of mRNA alternative polyadenylation (APA), an important regulatory mechanism during eukaryotic gene expression, in Arabidopsis thaliana and its halophytic relative Eutrema salsugineum with regard to their responses to salt stress. Analyses showed that while APA occurs commonly in both Arabidopsis and Eutrema, Eutrema possesses fewer APA genes than Arabidopsis (47% vs. 54%). However, the proportion of APA genes was significantly increased in Arabidopsis under salt stress but not in Eutrema. This indicated that Arabidopsis is more sensitive to salt stress and that Eutrema exhibits an innate response to such conditions. Both species utilized distal poly(A) sites under salt stress; however, only eight genes were found to overlap when their 3' untranslated region (UTR) lengthen genes were compared, thus revealing their distinct responses to salt stress. In Arabidopsis, genes that use distal poly(A) sites were enriched in response to salt stress. However, in Eutrema, the use of poly(A) sites was less affected and fewer genes were enriched. The transcripts with upregulated poly(A) sites in Arabidopsis showed enriched pathways in plant hormone signal transduction, starch and sucrose metabolism, and fatty acid elongation; in Eutrema, biosynthetic pathways (stilbenoid, diarylheptanoid, and gingerol) and metabolic pathways (arginine and proline) showed enrichment. APA was associated with 42% and 29% of the differentially expressed genes (DE genes) in Arabidopsis and Eutrema experiencing salt stress, respectively. Salt specific poly(A) sites and salt-inducible APA events were identified in both species; notably, some salt tolerance-related genes and transcription factor genes exhibited differential APA patterns, such as CIPK21 and LEA4-5. Our results suggest that adapted species exhibit more orderly response at the RNA maturation step under salt stress, while more salt-specific poly(A) sites were activated in Arabidopsis to cope with salinity conditions. Collectively, our findings not only highlight the importance of APA in the regulation of gene expression in response to salt stress, but also provide a new perspective on how salt-sensitive and salt-tolerant species perform differently under stress conditions through transcriptome diversity.

QuantifyPoly(A): reshaping alternative polyadenylation landscapes of eukaryotes with weighted density peak clustering

The dynamic choice of different polyadenylation sites in a gene is referred to as alternative polyadenylation, which functions in many important biological processes. Large-scale messenger RNA 3' end sequencing has revealed that cleavage sites for polyadenylation are presented with microheterogeneity. To date, the conventional determination of polyadenylation site clusters is subjective and arbitrary, leading to inaccurate annotations. Here, we present a weighted density peak clustering method, QuantifyPoly(A), to accurately quantify genome-wide polyadenylation choices. Applying QuantifyPoly(A) on published 3' end sequencing datasets from both animals and plants, their polyadenylation profiles are reshaped into myriads of novel polyadenylation site clusters. Most of these novel polyadenylation site clusters show significantly dynamic usage across different biological samples or associate with binding sites of trans-acting factors. Upstream sequences of these clusters are enriched with polyadenylation signals UGUA, UAAA and/or AAUAAA in a species-dependent manner. Polyadenylation site clusters also exhibit species specificity, while plants ones generally show higher microheterogeneity than that of animals. QuantifyPoly(A) is broadly applicable to any types of 3' end sequencing data and species for accurate quantification and construction of the complex and dynamic polyadenylation landscape and enables us to decode alternative polyadenylation events invisible to conventional methods at a much higher resolution.

Alternative Polyadenylation in response to temperature stress contributes to gene regulation in Populus trichocarpa

BACKGROUND

Genome-wide change of polyadenylation (polyA) sites (also known as alternative polyadenylation, APA) is emerging as an important strategy of gene regulation in response to stress in plants. But little is known in woody perennials that are persistently dealing with multiple abiotic stresses.

RESULTS

Here, we performed a genome-wide profiling of polyadenylation sites under heat and cold treatments in Populus trichocarpa. Through a comprehensive analysis of polyA tail sequences, we identified 25,919 polyA-site clusters (PACs), and revealed 3429 and 3139 genes shifted polyA sites under heat and cold stresses respectively. We found that a small proportion of genes possessed APA that affected the open reading frames; and some shifts were commonly identified. Functional analysis of genes displaying shifted polyA tails suggested that pathways related to RNA metabolism were linked to regulate the APA events under both heat and cold stresses. Interestingly, we found that the heat stress induced a significantly more antisense PACs comparing to cold and control conditions. Furthermore, we showed that a unique cis-element (AAAAAA) was predominately enriched downstream of PACs in P. trichocarpa genes; and this sequence signal was only absent in shifted PACs under the heat condition, indicating a distinct APA mechanism responsive to heat tolerance.

CONCLUSIONS

This work provides a comprehensive picture of global polyadenylation patterns in response to temperatures stresses in trees. We show that the frequent change of polyA tail is a potential mechanism of gene regulation responsive to stress, which are associated with distinctive sequence signatures.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

Cadmium Immobilization in the Rhizosphere and Plant Cellular Detoxification: Role of Plant-Growth-Promoting Rhizobacteria as a Sustainable Solution

Food is the major cadmium (Cd)-exposure pathway from agricultural soils to humans and other living entities and must be reduced in an effective way. A plant can select beneficial microbes, like plant-growth-promoting rhizobacteria (PGPR), depending upon the nature of root exudates in the rhizosphere, for its own benefits, such as plant growth promotion as well as protection from metal toxicity. This review intends to seek out information on the rhizo-immobilization of Cd in polluted soils using the PGPR along with plant nutrient fertilizers. This review suggests that the rhizo-immobilization of Cd by a combination of PGPR and nanohybrid-based plant nutrient fertilizers would be a potential and sustainable technology for phytoavailable Cd immobilization in the rhizosphere and plant cellular detoxification, by keeping the plant nutrition flow and green dynamics of plant nutrition and boosting the plant growth and development under Cd stress.

scDAPA: detection and visualization of dynamic alternative polyadenylation from single cell RNA-seq data

MOTIVATION

Alternative polyadenylation (APA) plays a key post-transcriptional regulatory role in mRNA stability and functions in eukaryotes. Single cell RNA-seq (scRNA-seq) is a powerful tool to discover cellular heterogeneity at gene expression level. Given 3' enriched strategy in library construction, the most commonly used scRNA-seq protocol-10× Genomics enables us to improve the study resolution of APA to the single cell level. However, currently there is no computational tool available for investigating APA profiles from scRNA-seq data.

RESULTS

Here, we present a package scDAPA for detecting and visualizing dynamic APA from scRNA-seq data. Taking bam/sam files and cell cluster labels as inputs, scDAPA detects APA dynamics using a histogram-based method and the Wilcoxon rank-sum test, and visualizes candidate genes with dynamic APA. Benchmarking results demonstrated that scDAPA can effectively identify genes with dynamic APA among different cell groups from scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

The scDAPA package is implemented in Shell and R, and is freely available at https://scdapa.sourceforge.io.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Discovery of alternative polyadenylation dynamics from single cell types

Alternative polyadenylation (APA) occurs in the process of mRNA maturation by adding a poly(A) tail at different locations, resulting increased diversity of mRNA isoforms and contributing to the complexity of gene regulatory network. Benefit from the development of high-throughput sequencing technologies, we could now delineate APA profiles of transcriptomes at an unprecedented pace. Especially the single cell RNA sequencing (scRNA-seq) technologies provide us opportunities to interrogate biological details of diverse and rare cell types. Despite increasing evidence showing that APA is involved in the cell type-specific regulation and function, efficient and specific laboratory methods for capturing poly(A) sites at single cell resolution are underdeveloped to date. In this review, we summarize existing experimental and computational methods for the identification of APA dynamics from diverse single cell types. A future perspective is also provided.

The Position and Complex Genomic Architecture of Plant T-DNA Insertions Revealed by 4SEE

The integration of T-DNA in plant genomes is widely used for basic research and agriculture. The high heterogeneity in the number of integration events per genome, their configuration, and their impact on genome integrity highlight the critical need to detect the genomic locations of T-DNA insertions and their associated chromosomal rearrangements, and the great challenge in doing so. Here, we present 4SEE, a circular chromosome conformation capture (4C)-based method for robust, rapid, and cost-efficient detection of the entire scope of T-DNA locations. Moreover, by measuring the chromosomal architecture of the plant genome flanking the T-DNA insertions, 4SEE outlines their associated complex chromosomal aberrations. Applying 4SEE to a collection of confirmed T-DNA lines revealed previously unmapped T-DNA insertions and chromosomal rearrangements such as inversions and translocations. Uncovering such events in a feasible, robust, and cost-effective manner by 4SEE in any plant of interest has implications for accurate annotation and phenotypic characterization of T-DNA insertion mutants and transgene expression in basic science applications as well as for plant biotechnology.

Alternative Polyadenylation and Salicylic Acid Modulate Root Responses to Low Nitrogen Availability

Nitrogen (N) is probably the most important macronutrient and its scarcity limits plant growth, development and fitness. N starvation response has been largely studied by transcriptomic analyses, but little is known about the role of alternative polyadenylation (APA) in such response. In this work, we show that N starvation modifies poly(A) usage in a large number of transcripts, some of them mediated by FIP1, a component of the polyadenylation machinery. Interestingly, the number of mRNAs isoforms with poly(A) tags located in protein-coding regions or 5'-UTRs significantly increases in response to N starvation. The set of genes affected by APA in response to N deficiency is enriched in N-metabolism, oxidation-reduction processes, response to stresses, and hormone responses, among others. A hormone profile analysis shows that the levels of salicylic acid (SA), a phytohormone that reduces nitrate accumulation and root growth, increase significantly upon N starvation. Meta-analyses of APA-affected and fip1-2-deregulated genes indicate a connection between the nitrogen starvation response and salicylic acid (SA) signaling. Genetic analyses show that SA may be important for preventing the overgrowth of the root system in low N environments. This work provides new insights on how plants interconnect different pathways, such as defense-related hormonal signaling and the regulation of genomic information by APA, to fine-tune the response to low N availability.

Alternative polyadenylation of the stacyose synthase gene mediates source-sink regulation in cucumber

Alternative polyadenylation (APA) is a pervasive mechanism for gene regulation in eukaryotes. Stachyose is the main assimilate translocated in the cucumber phloem. Stachyose synthase (CsSTS) catalyzes the last step of stachyose biosynthesis in cucumber leaves and plays a key role in the regulation of assimilate partitioning between source and sink. In this study, three CsSTS mRNAs with the same open reading frame and the 5`untranslated region (UTR), but differing in their 3`UTRs, named CsSTS1 (short), CsSTS2 (medium), and CsSTS3 (long), were identified. Southern blot and sequence analysis of the cucumber genome confirmed that these transcripts are regulated through APA from a single gene. No significant difference of in vitro translation efficiency was found among three mRNAs. However, the relative stabilities of three transcripts varied among different tissues and different leaf development stages of cucumber. CsSTS1 expression in cucumber calli was up-regulated by the raffinose (substrate of CsSTS) and down-regulated by stachyose (product of CsSTS), respectively. In cucumber plants, all three isoforms have considerable expression in non-fruit node leaves. However, in fruit-carrying node leaves, the expression of CsSTS2 and CsSTS3 was severely inhibited and only CsSTS1 was highly expressed, indicating fruit setting has a remarkable effect on the relative expression level of three transcripts. This "fruit setting" effect could be observed until at least 36 h after the fruit was removed from the node. Our results suggest that abundant expression of CsSTS1 is beneficial for stachyose loading in source leaves, and APA is a delicate mechanism for CsSTS to regulate cucumber source-sink balance.