Altered interactions between cis-regulatory elements partially resolve BLADE-ON-PETIOLE genetic redundancy in Capsella rubella.

Duplicated genes are thought to follow one of three evolutionary trajectories that resolve their redundancy: neofunctionalization, subfunctionalization, or pseudogenization. Differences in expression patterns have been documented for many duplicated gene pairs and interpreted as evidence of subfunctionalization and a loss of redundancy. However, little is known about the functional impact of such differences and about their molecular basis. Here, we investigate the genetic and molecular basis for the partial loss of redundancy between the two BLADE-ON-PETIOLE genes BOP1 and BOP2 in red shepherd's purse (Capsella rubella) compared to Arabidopsis (Arabidopsis thaliana). While both genes remain almost fully redundant in A. thaliana, BOP1 in C. rubella can no longer ensure wild-type floral organ numbers and suppress bract formation, due to an altered expression pattern in the region of the cryptic bract primordium. We use two complementary approaches, transgenic rescue of A. thaliana atbop1 atbop2 double mutants and deletions in the endogenous AtBOP1 promoter, to demonstrate that several BOP1 promoter regions containing conserved noncoding sequences interact in a nonadditive manner to control BOP1 expression in the bract primordium and that changes in these interactions underlie the evolutionary divergence between C. rubella and A. thaliana BOP1 expression and activity. Similarly, altered interactions between cis-regulatory regions underlie the divergence in functional promoter architecture related to the control of floral organ abscission by BOP1. These findings highlight the complexity of promoter architecture in plants and suggest that changes in the interactions between cis-regulatory elements are key drivers for evolutionary divergence in gene expression and the loss of redundancy.

Correction: Suppression of class I compensated cell enlargement by xs2 mutation is mediated by salicylic acid signaling.

[This corrects the article DOI: 10.1371/journal.pgen.1008873.].

RING/U-box E3 protein BIR1 interacts with and ubiquitinates barley growth repressor BROAD LEAF1.

Establishment of final leaf size in plants relies on the precise regulation of 2 interconnected processes, cell division and cell expansion. The barley (Hordeum vulgare) protein BROAD LEAF1 (BLF1) limits cell proliferation and leaf growth in the width direction. However, how the levels of this potent repressor of leaf growth are controlled remains unclear. Here, we used a yeast 2-hybrid screen to identify the BLF1-INTERACTING RING/U-BOX 1 (BIR1) E3 ubiquitin ligase that interacts with BLF1 and confirmed the interaction of the 2 proteins in planta. Inhibiting the proteasome caused overaccumulation of a BLF1-eGFP fusion protein when co-expressed with BIR1, and an in vivo ubiquitination assay in bacteria confirmed that BIR1 can mediate ubiquitination of BLF1 protein. Consistent with regulation of endogenous BLF1 in barley by proteasomal degradation, inhibition of the proteasome in BLF1-vYFP-expressing barley plants caused an accumulation of the BLF1 protein. The BIR1 protein co-localized with BLF1 in nuclei and appeared to reduce BLF1 protein levels. Analysis of bir1-1 knockout mutants suggested the involvement of BIR1 in leaf growth control, although mainly on leaf length. Together, our results suggest that proteasomal degradation, in part mediated by BIR1, helps fine-tune BLF1 protein levels in barley.

Poly(A) polymerase 1 contributes to competence acquisition of pollen tubes growing through the style in Arabidopsis thaliana.

Polyadenylation of mRNAs is critical for their export from the nucleus, stability, and efficient translation. The Arabidopsis thaliana genome encodes three isoforms of canonical nuclear poly(A) polymerase (PAPS) that redundantly polyadenylate the bulk of pre-mRNAs. However, previous studies have indicated that subsets of pre-mRNAs are preferentially polyadenylated by either PAPS1 or the other two isoforms. Such functional specialization raises the possibility of an additional level of gene-expression control in plants. Here we test this notion by studying the function of PAPS1 in pollen-tube growth and guidance. Pollen tubes growing through female tissue acquire the competence to find ovules efficiently and upregulate PAPS1 expression at the transcriptional, but not detectably at the protein level compared with in vitro grown pollen tubes. Using the temperature-sensitive paps1-1 allele we show that PAPS1 activity during pollen-tube growth is required for full acquisition of competence, resulting in inefficient fertilization by paps1-1 mutant pollen tubes. While these mutant pollen tubes grow almost at the wild-type rate, they are compromised in locating the micropyles of ovules. Previously identified competence-associated genes are less expressed in paps1-1 mutant than in wild-type pollen tubes. Estimating the poly(A) tail lengths of transcripts suggests that polyadenylation by PAPS1 is associated with reduced transcript abundance. Our results therefore suggest that PAPS1 plays a key role in the acquisition of competence and underline the importance of functional specialization between PAPS isoforms throughout different developmental stages.

The gain-of-function mutation blf13 in the barley orthologue of the rice growth regulator NARROW LEAF1 is associated with increased leaf width.

Canopy architecture in cereals plays an important role in determining yield. Leaf width represents one key aspect of this canopy architecture. However, our understanding of leaf width control in cereals remains incomplete. Classical mutagenesis studies in barely identified multiple morphological mutants, including those with differing leaf widths. Of these, we characterized the broad leaf13 (blf13) mutant in detail. Mutant plants form wider leaves due to increased post-initiation growth and cell proliferation. The mutant phenotype perfectly co-segregated with a missense mutation in the HvHNT1 gene which affected a highly conserved region of the encoded protein, orthologous to the rice NARROW LEAF1 (NAL1) protein. Causality of this mutation for the blf13 phenotype is further supported by correlative transcriptomic analyses and protein-protein interaction studies showing that the mutant HvNHT1 protein interacts more strongly with a known interactor than wild-type HvHNT1. The mutant HvHNT1 protein also showed stronger homodimerization compared with wild-type HvHNT1, and homology modelling suggested an additional interaction site between HvHNT1 monomers due to the blf13 mutation. Thus, the blf13 mutation parallels known gain-of-function NAL1 alleles in rice that increase leaf width and grain yield, suggesting that the blf13 mutation may have a similar agronomic potential in barley.

Molecular breeding of flower load related traits in dioecious autotetraploid Actinidia arguta.

Flowering plants exhibit a wide range of sexual reproduction systems, with the majority being hermaphroditic. However, some plants, such as Actinidia arguta (kiwiberry), have evolved into dioecious species with distinct female and male vines. In this study, we investigated the flower load and growth habits of female kiwiberry genotypes to identify the genetic basis of high yield with low maintenance requirements. Owing to the different selection approaches between female and male genotypes, we further extended our study to male kiwiberry genotypes. By combining both investigations, we present a novel breeding tool for dioecious crops. A population of A. arguta seedlings was phenotyped for flower load traits, in particular the proportion of non-floral shoots, proportion of floral shoots, and average number of flowers per floral shoot. Quantitative trait locus (QTL) mapping was used to analyse the genetic basis of these traits. We identified putative QTLs on chromosome 3 associated with flower-load traits. A pleiotropic effect of the male-specific region of the Y chromosome (MSY) on chromosome 3 affecting flower load-related traits between female and male vines was observed in an A. arguta breeding population. Furthermore, we utilized Genomic Best Linear Unbiased Prediction (GBLUP) to predict breeding values for the quantitative traits by leveraging genomic data. This approach allowed us to identify and select superior genotypes. Our findings contribute to the understanding of flowering and fruiting dynamics in Actinidia species, providing insights for kiwiberry breeding programs aiming to improve yield through the utilization of genomic methods and trait mapping. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01476-7.

Comparative Genomics Elucidates the Origin of a Supergene Controlling Floral Heteromorphism.

Supergenes are nonrecombining genomic regions ensuring the coinheritance of multiple, coadapted genes. Despite the importance of supergenes in adaptation, little is known on how they originate. A classic example of supergene is the S locus controlling heterostyly, a floral heteromorphism occurring in 28 angiosperm families. In Primula, heterostyly is characterized by the cooccurrence of two complementary, self-incompatible floral morphs and is controlled by five genes clustered in the hemizygous, ca. 300-kb S locus. Here, we present the first chromosome-scale genome assembly of any heterostylous species, that of Primula veris (cowslip). By leveraging the high contiguity of the P. veris assembly and comparative genomic analyses, we demonstrated that the S-locus evolved via multiple, asynchronous gene duplications and independent gene translocations. Furthermore, we discovered a new whole-genome duplication in Ericales that is specific to the Primula lineage. We also propose a mechanism for the origin of S-locus hemizygosity via nonhomologous recombination involving the newly discovered two pairs of CFB genes flanking the S locus. Finally, we detected only weak signatures of degeneration in the S locus, as predicted for hemizygous supergenes. The present study provides a useful resource for future research addressing key questions on the evolution of supergenes in general and the S locus in particular: How do supergenes arise? What is the role of genome architecture in the evolution of complex adaptations? Is the molecular architecture of heterostyly supergenes across angiosperms similar to that of Primula?

Testing the effect of individual scent compounds on pollinator attraction in nature using quasi-isogenic Capsella lines.

PREMISE: Floral scent, usually consisting of multiple compounds, is a complex trait, and its role in pollinator attraction has received increasing attention. However, disentangling the effect of individual floral scent compounds is difficult due to the complexity of isolating the effect of single compounds by traditional methods. METHODS: Using available quasi-isogenic lines (qILs) that were generated as part of the original mapping of the floral scent volatile-related loci CNL1 (benzaldehyde) and TPS2 (ß-ocimene) in Capsella, we generated four genotypes that should only differ in these two compounds. Plants of the four genotypes were introduced into a common garden outside the natural range of C. rubella or C. grandiflora, with individuals of a self-compatible C. grandiflora line as pollen donors, whose different genetic background facilitates the detection of outcrossing events. Visitors to flowers of all five genotypes were compared, and the seeds set during the common-garden period were collected for high-throughput amplicon-based sequencing to estimate their outcrossing rates. RESULTS: Benzaldehyde and ß-ocimene emissions were detected in the floral scent of corresponding genotypes. While some pollinator groups showed specific visitation preferences depending on scent compounds, the outcrossing rates in seeds did not vary among the four scent-manipulated genotypes. CONCLUSIONS: The scent-manipulated Capsella materials constructed using qILs provide a powerful system to study the ecological effects of individual floral scent compounds under largely natural environments. In Capsella, individual benzaldehyde and ß-ocimene emission may act as attractants for different types of pollinators.

Implementation of different relationship estimate methodologies in breeding value prediction in kiwiberry (Actinidia arguta).

In dioecious crops such as Actinidia arguta (kiwiberries), some of the main challenges when breeding for fruit characteristics are the selection of potential male parents and the long juvenile period. Currently, breeding values of male parents are estimated through progeny tests, which makes the breeding of new kiwiberry cultivars time-consuming and costly. The application of best linear unbiased prediction (BLUP) would allow direct estimation of sex-related traits and speed up kiwiberry breeding. In this study, we used a linear mixed model approach to estimate narrow sense heritability for one vine-related trait and five fruit-related traits for two incomplete factorial crossing designs. We obtained BLUPs for all genotypes, taking into consideration whether the relationship was pedigree-based or marker-based. Owing to the high cost of genome sequencing, it is important to understand the effects of different sources of relationship matrices on estimating breeding values across a breeding population. Because of the increasing implementation of genomic selection in crop breeding, we compared the effects of incorporating different sources of information in building relationship matrices and ploidy levels on the accuracy of BLUPs' heritability and predictive ability. As kiwiberries are autotetraploids, multivalent chromosome formation and occasionally double reduction can occur during meiosis, and this can affect the accuracy of prediction. This study innovates the breeding programme of autotetraploid kiwiberries. We demonstrate that the accuracy of BLUPs of male siblings, without phenotypic observations, strongly improved when a tetraploid marker-based relationship matrix was used rather than parental BLUPs and female siblings with phenotypic observations. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01419-8.

Supergene evolution via stepwise duplications and neofunctionalization of a floral-organ identity gene.

Heterostyly represents a fascinating adaptation to promote outbreeding in plants that evolved multiple times independently. While l-morph individuals form flowers with long styles, short anthers, and small pollen grains, S-morph individuals have flowers with short styles, long anthers, and large pollen grains. The difference between the morphs is controlled by an S-locus "supergene" consisting of several distinct genes that determine different traits of the syndrome and are held together, because recombination between them is suppressed. In Primula, the S locus is a roughly 300-kb hemizygous region containing five predicted genes. However, with one exception, their roles remain unclear, as does the evolutionary buildup of the S locus. Here we demonstrate that the MADS-box GLOBOSA2 (GLO2) gene at the S locus determines anther position. In Primula forbesii S-morph plants, GLO2 promotes growth by cell expansion in the fused tube of petals and stamen filaments beneath the anther insertion point; by contrast, neither pollen size nor male incompatibility is affected by GLO2 activity. The paralogue GLO1, from which GLO2 arose by duplication, has maintained the ancestral B-class function in specifying petal and stamen identity, indicating that GLO2 underwent neofunctionalization, likely at the level of the encoded protein. Genetic mapping and phylogenetic analysis indicate that the duplications giving rise to the style-length-determining gene CYP734A50 and to GLO2 occurred sequentially, with the CYP734A50 duplication likely the first. Together these results provide the most detailed insight into the assembly of a plant supergene yet and have important implications for the evolution of heterostyly.

Suppression of class I compensated cell enlargement by xs2 mutation is mediated by salicylic acid signaling.

The regulation of leaf size has been studied for decades. Enhancement of post-mitotic cell expansion triggered by impaired cell proliferation in Arabidopsis is an important process for leaf size regulation, and is known as compensation. This suggests a key interaction between cell proliferation and cell expansion during leaf development. Several studies have highlighted the impact of this integration mechanism on leaf size determination; however, the molecular basis of compensation remains largely unknown. Previously, we identified extra-small sisters (xs) mutants which can suppress compensated cell enlargement (CCE) via a specific defect in cell expansion within the compensation-exhibiting mutant, angustifolia3 (an3). Here we revealed that one of the xs mutants, namely xs2, can suppress CCE not only in an3 but also in other compensation-exhibiting mutants erecta (er) and fugu2. Molecular cloning of XS2 identified a deleterious mutation in CATION CALCIUM EXCHANGER 4 (CCX4). Phytohormone measurement and expression analysis revealed that xs2 shows hyper activation of the salicylic acid (SA) response pathway, where activation of SA response can suppress CCE in compensation mutants. All together, these results highlight the regulatory connection which coordinates compensation and SA response.

Multi-omics approach reveals the contribution of KLU to leaf longevity and drought tolerance.

KLU, encoded by a cytochrome P450 CYP78A family gene, generates an important-albeit unknown-mobile signal that is distinct from the classical phytohormones. Multiple lines of evidence suggest that KLU/KLU-dependent signaling functions in several vital developmental programs, including leaf initiation, leaf/floral organ growth, and megasporocyte cell fate. However, the interactions between KLU/KLU-dependent signaling and the other classical phytohormones, as well as how KLU influences plant physiological responses, remain poorly understood. Here, we applied in-depth, multi-omics analysis to monitor transcriptome and metabolome dynamics in klu-mutant and KLU-overexpressing Arabidopsis plants. By integrating transcriptome sequencing data and primary metabolite profiling alongside phytohormone measurements, our results showed that cytokinin signaling, with its well-established function in delaying leaf senescence, was activated in KLU-overexpressing plants. Consistently, KLU-overexpressing plants exhibited significantly delayed leaf senescence and increased leaf longevity, whereas the klu-mutant plants showed early leaf senescence. In addition, proline biosynthesis and catabolism were enhanced following KLU overexpression owing to increased expression of genes associated with proline metabolism. Furthermore, KLU-overexpressing plants showed enhanced drought-stress tolerance and reduced water loss. Collectively, our work illustrates a role for KLU in positively regulating leaf longevity and drought tolerance by synergistically activating cytokinin signaling and promoting proline metabolism. These data promote KLU as a potential ideal genetic target to improve plant fitness.

At-MINI ZINC FINGER2 and Sl-INHIBITOR OF MERISTEM ACTIVITY, a Conserved Missing Link in the Regulation of Floral Meristem Termination in Arabidopsis and Tomato.

In angiosperms, the gynoecium is the last structure to develop within the flower due to the determinate fate of floral meristem (FM) stem cells. The maintenance of stem cell activity before its arrest at the stage called FM termination affects the number of carpels that develop. The necessary inhibition at this stage of WUSCHEL (WUS), which is responsible for stem cell maintenance, involves a two-step mechanism. Direct repression mediated by the MADS domain transcription factor AGAMOUS (AG), followed by indirect repression requiring the C2H2 zinc-finger protein KNUCKLES (KNU), allow for the complete termination of floral stem cell activity. Here, we show that Arabidopsis thaliana MINI ZINC FINGER2 (AtMIF2) and its homolog in tomato (Solanum lycopersicum), INHIBITOR OF MERISTEM ACTIVITY (SlIMA), participate in the FM termination process by functioning as adaptor proteins. AtMIF2 and SlIMA recruit AtKNU and SlKNU, respectively, to form a transcriptional repressor complex together with TOPLESS and HISTONE DEACETYLASE19. AtMIF2 and SlIMA bind to the WUS and SlWUS loci in the respective plants, leading to their repression. These results provide important insights into the molecular mechanisms governing (FM) termination and highlight the essential role of AtMIF2/SlIMA during this developmental step, which determines carpel number and therefore fruit size.

The poly(A) polymerase PAPS1 interacts with the RNA-directed DNA-methylation pathway in sporophyte and pollen development.

RNA-based processes play key roles in the regulation of eukaryotic gene expression. This includes both the processing of pre-mRNAs into mature mRNAs ready for translation and RNA-based silencing processes, such as RNA-directed DNA methylation (RdDM). Polyadenylation of pre-mRNAs is one important step in their processing and is carried out by three functionally specialized canonical nuclear poly(A) polymerases in Arabidopsis thaliana. Null mutations in one of these, termed PAPS1, result in a male gametophytic defect. Using a fluorescence-labelling strategy, we have characterized this defect in more detail using RNA and small-RNA sequencing. In addition to global defects in the expression of pollen-differentiation genes, paps1 null-mutant pollen shows a strong overaccumulation of transposable element (TE) transcripts, yet a depletion of 21- and particularly 24-nucleotide-long short interfering RNAs (siRNAs) and microRNAs (miRNAs) targeting the corresponding TEs. Double-mutant analyses support a specific functional interaction between PAPS1 and components of the RdDM pathway, as evident from strong synergistic phenotypes in mutant combinations involving paps1, but not paps2 paps4, mutations. In particular, the double-mutant of paps1 and rna-dependent rna polymerase 6 (rdr6) shows a synergistic developmental phenotype disrupting the formation of the transmitting tract in the female gynoecium. Thus, our findings in A. thaliana uncover a potentially general link between canonical poly(A) polymerases as components of mRNA processing and RdDM, reflecting an analogous interaction in fission yeast.

Establishment of the Embryonic Shoot Meristem Involves Activation of Two Classes of Genes with Opposing Functions for Meristem Activities.

The shoot meristem, a stem-cell-containing tissue initiated during plant embryogenesis, is responsible for continuous shoot organ production in postembryonic development. Although key regulatory factors including KNOX genes are responsible for stem cell maintenance in the shoot meristem, how the onset of such factors is regulated during embryogenesis is elusive. Here, we present evidence that the two KNOX genes STM and KNAT6 together with the two other regulatory genes BLR and LAS are functionally important downstream genes of CUC1 and CUC2, which are a redundant pair of genes that specify the embryonic shoot organ boundary. Combined expression of STM with any of KNAT6, BLR, and LAS can efficiently rescue the defects of shoot meristem formation and/or separation of cotyledons in cuc1cuc2 double mutants. In addition, CUC1 and CUC2 are also required for the activation of KLU, a cytochrome P450-encoding gene known to restrict organ production, and KLU counteracts STM in the promotion of meristem activity, providing a possible balancing mechanism for shoot meristem maintenance. Together, these results establish the roles for CUC1 and CUC2 in coordinating the activation of two classes of genes with opposite effects on shoot meristem activity.

Fruit shape diversity in the Brassicaceae is generated by varying patterns of anisotropy.

Fruits exhibit a vast array of different 3D shapes, from simple spheres and cylinders to more complex curved forms; however, the mechanism by which growth is oriented and coordinated to generate this diversity of forms is unclear. Here, we compare the growth patterns and orientations for two very different fruit shapes in the Brassicaceae: the heart-shaped Capsella rubella silicle and the near-cylindrical Arabidopsis thaliana silique. We show, through a combination of clonal and morphological analyses, that the different shapes involve different patterns of anisotropic growth during three phases. These experimental data can be accounted for by a tissue-level model in which specified growth rates vary in space and time and are oriented by a proximodistal polarity field. The resulting tissue conflicts lead to deformation of the tissue as it grows. The model allows us to identify tissue-specific and temporally specific activities required to obtain the individual shapes. One such activity may be provided by the valve-identity gene FRUITFULL, which we show through comparative mutant analysis to modulate fruit shape during post-fertilisation growth of both species. Simple modulations of the model presented here can also broadly account for the variety of shapes in other Brassicaceae species, thus providing a simplified framework for fruit development and shape diversity.

Retracing the molecular basis and evolutionary history of the loss of benzaldehyde emission in the genus Capsella.

The transition from pollinator-mediated outbreeding to selfing has occurred many times in angiosperms. This is generally accompanied by a reduction in traits attracting pollinators, including reduced emission of floral scent. In Capsella, emission of benzaldehyde as a main component of floral scent has been lost in selfing C. rubella by mutation of cinnamate-CoA ligase CNL1. However, the biochemical basis and evolutionary history of this loss remain unknown, as does the reason for the absence of benzaldehyde emission in the independently derived selfer Capsella orientalis. We used plant transformation, in vitro enzyme assays, population genetics and quantitative genetics to address these questions. CNL1 has been inactivated twice independently by point mutations in C. rubella, causing a loss of enzymatic activity. Both inactive haplotypes are found within and outside of Greece, the centre of origin of C. rubella, indicating that they arose before its geographical spread. By contrast, the loss of benzaldehyde emission in C. orientalis is not due to an inactivating mutation in CNL1. CNL1 represents a hotspot for mutations that eliminate benzaldehyde emission, potentially reflecting the limited pleiotropy and large effect of its inactivation. Nevertheless, even closely related species have followed different evolutionary routes in reducing floral scent.

Standing genetic variation in a tissue-specific enhancer underlies selfing-syndrome evolution in Capsella.

Mating system shifts recurrently drive specific changes in organ dimensions. The shift in mating system from out-breeding to selfing is one of the most frequent evolutionary transitions in flowering plants and is often associated with an organ-specific reduction in flower size. However, the evolutionary paths along which polygenic traits, such as size, evolve are poorly understood. In particular, it is unclear how natural selection can specifically modulate the size of one organ despite the pleiotropic action of most known growth regulators. Here, we demonstrate that allelic variation in the intron of a general growth regulator contributed to the specific reduction of petal size after the transition to selfing in the genus Capsella Variation within this intron affects an organ-specific enhancer that regulates the level of STERILE APETALA (SAP) protein in the developing petals. The resulting decrease in SAP activity leads to a shortening of the cell proliferation period and reduced number of petal cells. The absence of private polymorphisms at the causal region in the selfing species suggests that the small-petal allele was captured from standing genetic variation in the ancestral out-crossing population. Petal-size variation in the current out-crossing population indicates that several small-effect mutations have contributed to reduce petal-size. These data demonstrate how tissue-specific regulatory elements in pleiotropic genes contribute to organ-specific evolution. In addition, they provide a plausible evolutionary explanation for the rapid evolution of flower size after the out-breeding-to-selfing transition based on additive effects of segregating alleles.

Genome-Wide Analysis of PAPS1-Dependent Polyadenylation Identifies Novel Roles for Functionally Specialized Poly(A) Polymerases in Arabidopsis thaliana.

The poly(A) tail at 3' ends of eukaryotic mRNAs promotes their nuclear export, stability and translational efficiency, and changes in its length can strongly impact gene expression. The Arabidopsis thaliana genome encodes three canonical nuclear poly(A) polymerases, PAPS1, PAPS2 and PAPS4. As shown by their different mutant phenotypes, these three isoforms are functionally specialized, with PAPS1 modifying organ growth and suppressing a constitutive immune response. However, the molecular basis of this specialization is largely unknown. Here, we have estimated poly(A)-tail lengths on a transcriptome-wide scale in wild-type and paps1 mutants. This identified categories of genes as particularly strongly affected in paps1 mutants, including genes encoding ribosomal proteins, cell-division factors and major carbohydrate-metabolic proteins. We experimentally verified two novel functions of PAPS1 in ribosome biogenesis and redox homoeostasis that were predicted based on the analysis of poly(A)-tail length changes in paps1 mutants. When overlaying the PAPS1-dependent effects observed here with coexpression analysis based on independent microarray data, the two clusters of transcripts that are most closely coexpressed with PAPS1 show the strongest change in poly(A)-tail length and transcript abundance in paps1 mutants in our analysis. This suggests that their coexpression reflects at least partly the preferential polyadenylation of these transcripts by PAPS1 versus the other two poly(A)-polymerase isoforms. Thus, transcriptome-wide analysis of poly(A)-tail lengths identifies novel biological functions and likely target transcripts for polyadenylation by PAPS1. Data integration with large-scale co-expression data suggests that changes in the relative activities of the isoforms are used as an endogenous mechanism to co-ordinately modulate plant gene expression.

Antagonistic control of flowering time by functionally specialized poly(A) polymerases in Arabidopsis thaliana.

Polyadenylation is a critical 3'-end processing step during maturation of pre-mRNAs, and the length of the poly(A) tail affects mRNA stability, nuclear export and translation efficiency. The Arabidopsis thaliana genome encodes three canonical nuclear poly(A) polymerase (PAPS) isoforms fulfilling specialized functions, as reflected by their different mutant phenotypes. While PAPS1 affects several processes, such as the immune response, organ growth and male gametophyte development, the roles of PAPS2 and PAPS4 are largely unknown. Here we demonstrate that PAPS2 and PAPS4 promote flowering in a partially redundant manner. The enzymes act antagonistically to PAPS1, which delays the transition to flowering. The opposite flowering-time phenotypes in paps1 and paps2 paps4 mutants are at least partly due to decreased or increased FLC activity, respectively. In contrast to paps2 paps4 mutants, plants with increased PAPS4 activity flower earlier than the wild-type, concomitant with reduced FLC expression. Double mutant analyses suggest that PAPS2 and PAPS4 act independently of the autonomous pathway components FCA, FY and CstF64. The direct polyadenylation targets of the three PAPS isoforms that mediate their effects on flowering time do not include FLC sense mRNA and remain to be identified. Thus, our results uncover a role for canonical PAPS isoforms in flowering-time control, raising the possibility that modulating the balance of the isoform activities could be used to fine tune the transition to flowering.