Refine
Year of publication
Language
- English (40)
Is part of the Bibliography
- yes (40)
Keywords
- Capsella (4)
- heterostyly (3)
- supergene (3)
- Amplicon sequencing (2)
- Arabidopsis thaliana (2)
- CLE7 (2)
- CYP734A50 (2)
- EMS (2)
- Mixed mating (2)
- MutMap (2)
Institute
Polyadenylation of pre-mRNAs is critical for efficient nuclear export, stability, and translation of the mature mRNAs, and thus for gene expression. The bulk of pre-mRNAs are processed by canonical nuclear poly(A) polymerase (PAPS). Both vertebrate and higher-plant genomes encode more than one isoform of this enzyme, and these are coexpressed in different tissues. However, in neither case is it known whether the isoforms fulfill different functions or polyadenylate distinct subsets of pre-mRNAs. Here we show that the three canonical nuclear PAPS isoforms in Arabidopsis are functionally specialized owing to their evolutionarily divergent C-terminal domains. A strong loss-of-function mutation in PAPS1 causes a male gametophytic defect, whereas a weak allele leads to reduced leaf growth that results in part from a constitutive pathogen response. By contrast, plants lacking both PAPS2 and PAPS4 function are viable with wild-type leaf growth. Polyadenylation of SMALL AUXIN UP RNA (SAUR) mRNAs depends specifically on PAPS1 function. The resulting reduction in SAUR activity in paps1 mutants contributes to their reduced leaf growth, providing a causal link between polyadenylation of specific pre-mRNAs by a particular PAPS isoform and plant growth. This suggests the existence of an additional layer of regulation in plant and possibly vertebrate gene expression, whereby the relative activities of canonical nuclear PAPS isoforms control de novo synthesized poly(A) tail length and hence expression of specific subsets of mRNAs.
Polyadenylation of pre-mRNAs by poly(A) polymerase (PAPS) is a critical process in eukaryotic gene expression. As found in vertebrates, plant genomes encode several isoforms of canonical nuclear PAPS enzymes. In Arabidopsis thaliana these isoforms are functionally specialized, with PAPS1 affecting both organ growth and immune response, at least in part by the preferential polyadenylation of subsets of pre-mRNAs. Here, we demonstrate that the opposite effects of PAPS1 on leaf and flower growth reflect the different identities of these organs, and identify a role for PAPS1 in the elusive connection between organ identity and growth patterns. The overgrowth of paps1 mutant petals is due to increased recruitment of founder cells into early organ primordia, and suggests that PAPS1 activity plays unique roles in influencing organ growth. By contrast, the leaf phenotype of paps1 mutants is dominated by a constitutive immune response that leads to increased resistance to the biotrophic oomycete Hyaloperonospora arabidopsidis and reflects activation of the salicylic acid-independent signalling pathway downstream of ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1)/PHYTOALEXIN DEFICIENT4 (PAD4). These findings provide an insight into the developmental and physiological basis of the functional specialization amongst plant PAPS isoforms.
Induced point mutations are important genetic resources for their ability to create hypo- and hypermorphic alleles that are useful for understanding gene functions and breeding. However, such mutant populations have only been developed for a few temperate maize varieties, mainly B73 and W22, yet no tropical maize inbred lines have been mutagenized and made available to the public to date. We developed a novel Ethyl Methanesulfonate (EMS) induced mutation resource in maize comprising 2050 independent M2 mutant families in the elite tropical maize inbred ML10. By phenotypic screening, we showed that this population is of comparable quality with other mutagenized populations in maize. To illustrate the usefulness of this population for gene discovery, we performed rapid mapping-by-sequencing to clone a fasciated-ear mutant and identify a causal promoter deletion in ZmCLE7 (CLE7). Our mapping procedure does not require crossing to an unrelated parent, thus is suitable for mapping subtle traits and ones affected by heterosis. This first EMS population in tropical maize is expected to be very useful for the maize research community. Also, the EMS mutagenesis and rapid mapping-by-sequencing pipeline described here illustrate the power of performing forward genetics in diverse maize germplasms of choice, which can lead to novel gene discovery due to divergent genetic backgrounds.
Induced point mutations are important genetic resources for their ability to create hypo- and hypermorphic alleles that are useful for understanding gene functions and breeding. However, such mutant populations have only been developed for a few temperate maize varieties, mainly B73 and W22, yet no tropical maize inbred lines have been mutagenized and made available to the public to date. We developed a novel Ethyl Methanesulfonate (EMS) induced mutation resource in maize comprising 2050 independent M2 mutant families in the elite tropical maize inbred ML10. By phenotypic screening, we showed that this population is of comparable quality with other mutagenized populations in maize. To illustrate the usefulness of this population for gene discovery, we performed rapid mapping-by-sequencing to clone a fasciated-ear mutant and identify a causal promoter deletion in ZmCLE7 (CLE7). Our mapping procedure does not require crossing to an unrelated parent, thus is suitable for mapping subtle traits and ones affected by heterosis. This first EMS population in tropical maize is expected to be very useful for the maize research community. Also, the EMS mutagenesis and rapid mapping-by-sequencing pipeline described here illustrate the power of performing forward genetics in diverse maize germplasms of choice, which can lead to novel gene discovery due to divergent genetic backgrounds.
Introduction
The transition from cross-fertilisation (outcrossing) to self-fertilisation (selfing) frequently coincides with changes towards a floral morphology that optimises self-pollination, the selfing syndrome. Population genetic studies have reported the existence of both outcrossing and selfing populations in Arabis alpina (Brassicaceae), which is an emerging model species for studying the molecular basis of perenniality and local adaptation. It is unknown whether its selfing populations have evolved a selfing syndrome.
Methods
Using macro-photography, microscopy and automated cell counting, we compared floral syndromes (size, herkogamy, pollen and ovule numbers) between three outcrossing populations from the Apuan Alps and three selfing populations from the Western and Central Alps (Maritime Alps and Dolomites). In addition, we genotyped the plants for 12 microsatellite loci to confirm previous measures of diversity and inbreeding coefficients based on allozymes, and performed Bayesian clustering.
Results and Discussion
Plants from the three selfing populations had markedly smaller flowers, less herkogamy and lower pollen production than plants from the three outcrossing populations, whereas pistil length and ovule number have remained constant. Compared to allozymes, microsatellite variation was higher, but revealed similar patterns of low diversity and high Fis in selfing populations. Bayesian clustering revealed two clusters. The first cluster contained the three outcrossing populations from the Apuan Alps, the second contained the three selfing populations from the Maritime Alps and Dolomites.
Conclusion
We conclude that in comparison to three outcrossing populations, three populations with high selfing rates are characterised by a flower morphology that is closer to the selfing syndrome. The presence of outcrossing and selfing floral syndromes within a single species will facilitate unravelling the genetic basis of the selfing syndrome, and addressing which selective forces drive its evolution.
Introduction
The transition from cross-fertilisation (outcrossing) to self-fertilisation (selfing) frequently coincides with changes towards a floral morphology that optimises self-pollination, the selfing syndrome. Population genetic studies have reported the existence of both outcrossing and selfing populations in Arabis alpina (Brassicaceae), which is an emerging model species for studying the molecular basis of perenniality and local adaptation. It is unknown whether its selfing populations have evolved a selfing syndrome.
Methods
Using macro-photography, microscopy and automated cell counting, we compared floral syndromes (size, herkogamy, pollen and ovule numbers) between three outcrossing populations from the Apuan Alps and three selfing populations from the Western and Central Alps (Maritime Alps and Dolomites). In addition, we genotyped the plants for 12 microsatellite loci to confirm previous measures of diversity and inbreeding coefficients based on allozymes, and performed Bayesian clustering.
Results and Discussion
Plants from the three selfing populations had markedly smaller flowers, less herkogamy and lower pollen production than plants from the three outcrossing populations, whereas pistil length and ovule number have remained constant. Compared to allozymes, microsatellite variation was higher, but revealed similar patterns of low diversity and high Fis in selfing populations. Bayesian clustering revealed two clusters. The first cluster contained the three outcrossing populations from the Apuan Alps, the second contained the three selfing populations from the Maritime Alps and Dolomites.
Conclusion
We conclude that in comparison to three outcrossing populations, three populations with high selfing rates are characterised by a flower morphology that is closer to the selfing syndrome. The presence of outcrossing and selfing floral syndromes within a single species will facilitate unravelling the genetic basis of the selfing syndrome, and addressing which selective forces drive its evolution.
Gene duplication is a major driver for the increase of biological complexity. The divergence of newly duplicated paralogs may allow novel functions to evolve, while maintaining the ancestral one. Alternatively, partitioning the ancestral function among paralogs may allow parts of that role to follow independent evolutionary trajectories. We studied the REDUCED COMPLEXITY (RCO) locus, which contains three paralogs that have evolved through two independent events of gene duplication, and which underlies repeated events of leaf shape evolution within the Brassicaceae. In particular, we took advantage of the presence of three potentially functional paralogs in Capsella to investigate the extent of functional divergence among them. We demonstrate that the RCO copies control growth in different areas of the leaf. Consequently, the copies that are retained active in the different Brassicaceae lineages contribute to define the leaf dissection pattern. Our results further illustrate how successive gene duplication events and subsequent functional divergence can increase trait evolvability by providing independent evolutionary trajectories to specialized functions that have an additive effect on a given trait.
Elucidating the genetic basis of morphological changes in evolution remains a major challenge in biology [1-3]. Repeated independent trait changes are of particular interest because they can indicate adaptation in different lineages or genetic and developmental constraints on generating morphological variation [4-6]. In animals, changes to "hot spot" genes with minimal pleiotropy and large phenotypic effects underlie many cases of repeated morphological transitions [4-8]. By contrast, only few such genes have been identified from plants [8-11], limiting cross-kingdom comparisons of the principles of morphological evolution. Here, we demonstrate that the REDUCED COMPLEXITY (RCO) locus [12] underlies more than one naturally evolved change in leaf shape in the Brassicaceae. We show that the difference in leaf margin dissection between the sister species Capsella rubella and Capsella grandiflora is caused by cis-regulatory variation in the homeobox gene RCO-A, which alters its activity in the developing lobes of the leaf. Population genetic analyses in the ancestral C. grandiflora indicate that the more-active C. rubella haplotype is derived from a now rare or lost C. grandiflora haplotype via additional mutations. In Arabidopsis thaliana, the deletion of the RCO-A and RCO-B genes has contributed to its evolutionarily derived smooth leaf margin [12], suggesting the RCO locus as a candidate for an evolutionary hot spot. We also find that temperature-responsive expression of RCO-A can explain the phenotypic plasticity of leaf shape to ambient temperature in Capsella, suggesting a molecular basis for the well-known negative correlation between temperature and leaf margin dissection.
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.
In the Bateson–Dobzhansky–Muller model of genetic incompatibilities post-zygotic gene-flow barriers arise by fixation of novel alleles at interacting loci in separated populations. Many such incompatibilities are polymorphic in plants, implying an important role for genetic drift or balancing selection in their origin and evolution. Here we show that NPR1 and RPP5 loci cause a genetic incompatibility between the incipient species Capsella grandiflora and C. rubella, and the more distantly related C. rubella and C. orientalis. The incompatible RPP5 allele results from a mutation in C. rubella, while the incompatible NPR1 allele is frequent in the ancestral C. grandiflora. Compatible and incompatible NPR1 haplotypes are maintained by balancing selection in C. grandiflora, and were divergently sorted into the derived C. rubella and C. orientalis. Thus, by maintaining differentiated alleles at high frequencies, balancing selection on ancestral polymorphisms can facilitate establishing gene-flow barriers between derived populations through lineage sorting of the alternative alleles.
In the Bateson–Dobzhansky–Muller model of genetic incompatibilities post-zygotic gene-flow
barriers arise by fixation of novel alleles at interacting loci in separated populations. Many such incompatibilities are polymorphic in plants, implying an important role for genetic drift or balancing selection in their origin and evolution. Here we show that NPR1 and RPP5 loci cause a genetic incompatibility between the incipient species Capsella grandiflora and C. rubella, and the more distantly related C. rubella and C. orientalis. The incompatible RPP5 allele results from a mutation in C. rubella, while the incompatible NPR1 allele is frequent in the ancestral C. grandiflora. Compatible and incompatible NPR1 haplotypes are maintained by balancing selection in C. grandiflora, and were divergently sorted into the derived C. rubella and
C. orientalis. Thus, by maintaining differentiated alleles at high frequencies, balancing selection
on ancestral polymorphisms can facilitate establishing gene-flow barriers between derived populations through lineage sorting of the alternative alleles.
In the Bateson-Dobzhansky-Muller model of genetic incompatibilities post-zygotic gene-flow barriers arise by fixation of novel alleles at interacting loci in separated populations. Many such incompatibilities are polymorphic in plants, implying an important role for genetic drift or balancing selection in their origin and evolution. Here we show that NPR1 and RPP5 loci cause a genetic incompatibility between the incipient species Capsella grandiflora and C. rubella, and the more distantly related C. rubella and C. orientalis. The incompatible RPP5 allele results from a mutation in C. rubella, while the incompatible NPR1 allele is frequent in the ancestral C. grandiflora. Compatible and incompatible NPR1 haplotypes are maintained by balancing selection in C. grandiflora, and were divergently sorted into the derived C. rubella and C. orientalis. Thus, by maintaining differentiated alleles at high frequencies, balancing selection on ancestral polymorphisms can facilitate establishing gene-flow barriers between derived populations through lineage sorting of the alternative alleles.
The enormous species richness of flowering plants is at least partly due to floral diversification driven by interactions between plants and their animal pollinators [1, 2]. Specific pollinator attraction relies on visual and olfactory floral cues [3-5]; floral scent can not only attract pollinators but also attract or repel herbivorous insects [6-8]. However, despite its central role for plant-animal interactions, the genetic control of floral scent production and its evolutionary modification remain incompletely understood [9-13]. Benzenoids are an important class of floral scent compounds that are generated from phenylalanine via several enzymatic pathways [14-17]. Here we address the genetic basis of the loss of floral scent associated with the transition from outbreeding to selfing in the genus Capsella. While the outbreeding C. grandiflora emits benzaldehyde as a major constituent of its floral scent, this has been lost in the selfing C. rubella. We identify the Capsella CNL1 gene encoding cinnamate: CoA ligase as responsible for this variation. Population genetic analysis indicates that CNL1 has been inactivated twice independently in C. rubella via different novel mutations to its coding sequence. Together with a recent study in Petunia [18], this identifies cinnamate: CoA ligase as an evolutionary hotspot for mutations causing the loss of benzenoid scent compounds in association with a shift in the reproductive strategy of Capsella from pollination by insects to self-fertilization.
The enormous species richness of flowering plants is at least partly due to floral diversification driven by interactions between plants and their animal pollinators [1, 2]. Specific pollinator attraction relies on visual and olfactory floral cues [3-5]; floral scent can not only attract pollinators but also attract or repel herbivorous insects [6-8]. However, despite its central role for plant-animal interactions, the genetic control of floral scent production and its evolutionary modification remain incompletely understood [9-13]. Benzenoids are an important class of floral scent compounds that are generated from phenylalanine via several enzymatic pathways [14-17]. Here we address the genetic basis of the loss of floral scent associated with the transition from outbreeding to selfing in the genus Capsella. While the outbreeding C. grandiflora emits benzaldehyde as a major constituent of its floral scent, this has been lost in the selfing C. rubella. We identify the Capsella CNL1 gene encoding cinnamate: CoA ligase as responsible for this variation. Population genetic analysis indicates that CNL1 has been inactivated twice independently in C. rubella via different novel mutations to its coding sequence. Together with a recent study in Petunia [18], this identifies cinnamate: CoA ligase as an evolutionary hotspot for mutations causing the loss of benzenoid scent compounds in association with a shift in the reproductive strategy of Capsella from pollination by insects to self-fertilization.
In grape (Vitis vinifera), abscisic acid (ABA) accumulates during fruit ripening and is thought to play a pivotal role in this process, but the molecular basis of this control is poorly understood. This work characterizes ABSCISIC ACID RESPONSE ELEMENT-BINDING FACTOR2 (VvABF2), a grape basic leucine zipper transcription factor belonging to a phylogenetic subgroup previously shown to be involved in ABA and abiotic stress signaling in other plant species. VvABF2 transcripts mainly accumulated in the berry, from the onset of ripening to the harvesting stage, and were up-regulated by ABA. Microarray analysis of transgenic grape cells overexpressing VvABF2 showed that this transcription factor up-regulates and/or modifies existing networks related to ABA responses. In addition, grape cells overexpressing VvABF2 exhibited enhanced responses to ABA treatment compared with control cells. Among the VvABF2-mediated responses highlighted in this study, the synthesis of phenolic compounds and cell wall softening were the most strongly affected. VvABF2 overexpression strongly increased the accumulation of stilbenes that play a role in plant defense and human health (resveratrol and piceid). In addition, the firmness of fruits from tomato (Solanum lycopersicum) plants overexpressing VvABF2 was strongly reduced. These data indicate that VvABF2 is an important transcriptional regulator of ABA-dependent grape berry ripening.
Reproductive development of grapevine and berry composition are both strongly influenced by temperature. To date, the molecular mechanisms involved in grapevine berries response to high temperatures are poorly understood. Unlike recent data that addressed the effects on berry development of elevated temperatures applied at the whole plant level, the present work particularly focuses on the fruit responses triggered by direct exposure to heat treatment (HT). In the context of climate change, this work focusing on temperature effect at the microclimate level is of particular interest as it can help to better understand the consequences of leaf removal (a common viticultural practice) on berry development. HT (+8 degrees C) was locally applied to clusters from Cabernet Sauvignon fruiting cuttings at three different developmental stages (middle green, veraison and middle ripening). Samples were collected 1, 7, and 14 days after treatment and used for metabolic and transcriptomic analyses. The results showed dramatic and specific biochemical and transcriptomic changes in heat exposed berries, depending on the developmental stage and the stress duration. When applied at the herbaceous stage, HT delayed the onset of veraison. Heating also strongly altered the berry concentration of amino acids and organic acids (e.g., phenylalanine, raminobutyric acid and malate) and decreased the anthocyanin content at maturity. These physiological alterations could be partly explained by the deep remodeling of transcriptome in heated berries. More than 7000 genes were deregulated in at least one of the nine experimental conditions. The most affected processes belong to the categories "stress responses," protein metabolism" and "secondary metabolism," highlighting the intrinsic capacity of grape berries to perceive HT and to build adaptive responses. Additionally, important changes in processes related to "transport," "hormone" and "cell wall" might contribute to the postponing of veraison. Finally, opposite effects depending on heating duration were observed for genes encoding enzymes of the general phenylpropanoid pathway, suggesting that the HI induced decrease in anthocyanin content may result from a combination of transcript abundance and product degradation.
In addition to their role as a source of reduced carbon, sugars may directly or indirectly control a wide range of activities in plant cells, through transcriptional and post-translational regulation. This control has been studied in detail using Arabidopsis thaliana, where genetic analysis offers many possibilities. Much less is known about perennial woody species. For several years, various aspects of sugar sensing and signalling have been investigated in the grape (Vitis vinifera L.) berry, an organ that accumulates high concentrations of hexoses in the vacuoles of flesh cells. Here we review various aspects of this topic: the molecular basis of sugar transport and its regulation by sugars in grapevine; the functional analysis of several sugar-induced genes; the effects of some biotic and abiotic stresses on the sugar content of the berry; and finally the effects of exogenous sugar supply on the ripening process in field conditions. A picture of complex feedback and multiprocess regulation emerges from these data.
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.
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.
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.