Subgenome-dominant expression and alternative splicing in response to Sclerotinia infection in polyploid Brassica napus and progenitors.

Polyploidy has played an extensive role in the evolution of flowering plants. Allopolyploids, with subgenomes containing duplicated gene pairs called homeologs, can show rapid transcriptome changes including novel alternative splicing (AS) patterns. The extent to which abiotic stress modulates AS of homeologs is a nascent topic in polyploidy research. We subjected both resynthesized and natural lines of polyploid Brassica napus, along with the progenitors Brassica rapa and Brassica oleracea, to infection with the fungal pathogen Sclerotinia sclerotiorum. RNA-sequencing analyses revealed widespread divergence between polyploid subgenomes in both gene expression and AS patterns. Resynthesized B. napus displayed significantly more A and C subgenome biased homeologs under pathogen infection than during uninfected growth. Differential AS (DAS) in response to infection was highest in natural B. napus (12 709 DAS events) and lower in resynthesized B. napus (8863 DAS events). Natural B. napus had more upregulated events and fewer downregulated events. There was a global expression bias towards the B. oleracea-derived (C) subgenome in both resynthesized and natural B. napus, enhanced by widespread non-parental downregulation of the B. rapa-derived (A) homeolog. In the resynthesized B. napus, this resulted in a disproportionate C subgenome contribution to the pathogen defense response, characterized by biases in both transcript expression levels and the proportion of induced genes. Our results elucidate the complex ways in which Sclerotinia infection affects expression and AS of homeologous genes in resynthesized and natural B. napus.

Taxonomically Restricted Genes Are Associated With Responses to Biotic and Abiotic Stresses in Sugarcane (Saccharum spp.).

Orphan genes (OGs) are protein-coding genes that are restricted to particular clades or species and lack homology with genes from other organisms, making their biological functions difficult to predict. OGs can rapidly originate and become functional; consequently, they may support rapid adaptation to environmental changes. Extensive spread of mobile elements and whole-genome duplication occurred in the Saccharum group, which may have contributed to the origin and diversification of OGs in the sugarcane genome. Here, we identified and characterized OGs in sugarcane, examined their expression profiles across tissues and genotypes, and investigated their regulation under varying conditions. We identified 319 OGs in the Saccharum spontaneum genome without detected homology to protein-coding genes in green plants, except those belonging to Saccharinae. Transcriptomic analysis revealed 288 sugarcane OGs with detectable expression levels in at least one tissue or genotype. We observed similar expression patterns of OGs in sugarcane genotypes originating from the closest geographical locations. We also observed tissue-specific expression of some OGs, possibly indicating a complex regulatory process for maintaining diverse functional activity of these genes across sugarcane tissues and genotypes. Sixty-six OGs were differentially expressed under stress conditions, especially cold and osmotic stresses. Gene co-expression network and functional enrichment analyses suggested that sugarcane OGs are involved in several biological mechanisms, including stimulus response and defence mechanisms. These findings provide a valuable genomic resource for sugarcane researchers, especially those interested in selecting stress-responsive genes.

Expression Patterns and Functional Characterization of Arabidopsis Galactose Oxidase-Like Genes Suggest Specialized Roles for Galactose Oxidases in Plants.

Galactose oxidases (GalOxs) are well-known enzymes that have been identified in several fungal species and characterized using structural and enzymatic approaches. However, until very recently, almost no information on their biological functions was available. The Arabidopsis (Arabidopsis thaliana) gene ruby particles in mucilage (RUBY) encodes a putative plant GalOx that is required for pectin cross-linking through modification of galactose (Gal) side chains and promotes cell-cell adhesion between seed coat epidermal cells. RUBY is one member of a family of seven putative GalOxs encoded in the Arabidopsis genome. To examine the function(s) of GalOxs in plants, we studied the remaining six galactose oxidase-like (GOXL) proteins. Like RUBY, four of these proteins (GOXL1, GOXL3, GOXL5 and GOXL6) were found to localize primarily to the apoplast, while GOXL2 and GOXL4 were found primarily in the cytoplasm. Complementation and GalOx assay data suggested that GOXL1, GOXL3 and possibly GOXL6 have similar biochemical activity to RUBY, whereas GOXL5 only weakly complemented and GOXL2 and GOXL4 showed no activity. Members of this protein family separated into four distinct clades prior to the divergence of the angiosperms. There have been recent duplications in Brassicaceae resulting in two closely related pairs of genes that have either retained similarity in expression (GOXL1 and GOXL6) or show expression divergence (GOXL3 and RUBY). Mutant phenotypes were not detected when these genes were disrupted, but their expression patterns suggest that these proteins may function in tissues that require mechanical reinforcements in the absence of lignification.

Global insights into duplicated gene expression and alternative splicing in polyploid Brassica napus under heat, cold, and drought stress.

Polyploidy has been a prevalent process during plant evolution and it has made a major impact on the structure and evolution of plant genomes. Many important crop plants are polyploid. There is considerable interest in expression patterns of duplicated genes in polyploids. Alternative splicing (AS) is a fundamental aspect of gene expression that produces multiple final transcript types from a single type of mRNAs. The effects of abiotic stress conditions on AS in polyploids has received little attention. We conducted a global transcriptome analysis of Brassica napus, an allotetraploid derived from B. rapa (AT ) and B. oleracea (CT ), by RNA-Seq of plants subjected to cold, heat, and drought stress treatments. Analyses of 27,360 pairs of duplicated genes revealed overall AT subgenome biases in gene expression and CT subgenome biases in the extent of alternative splicing under all three stress treatments. More genes increased in expression than decreased in response to the stresses. Negative correlations were found between expression levels and AS frequency for each type of AS. Cold stress produced the greatest changes in gene expression and AS. Cold-induced AS changes were more likely to be shared with those generated by drought than by heat stress. We used homeologs of FLC and CCA1 as case studies to show the dynamics of how duplicates in a polyploid respond to cold stress. Our results suggest that divergence in gene expression and AS patterns between duplicated genes may increase the flexibility of polyploids when responding to abiotic stressors.

Divergence of duplicated genes by repeated partitioning of splice forms and subcellular localization.

Gene duplication is a prominent and recurrent process in plant genomes. Among the possible fates of duplicated genes, subfunctionalization refers to duplicates taking on different parts of the function or expression pattern of the ancestral gene. This partitioning could be accompanied by changes in subcellular localization of the protein products. When alternative splicing of gene products leads to protein products with different subcellular localizations, we propose that after gene duplication there will be partitioning of the alternatively spliced forms such that the products of each duplicate are localized to only one of the original locations, which we refer to as sublocalization. We identified the plastid ascorbate peroxidase (cpAPX) genes across angiosperms and analyzed their duplication history, alternative splicing, and subcellular targeting patterns to identify cases of sublocalization. We found angiosperms typically have one cpAPX gene that generates both thylakoidal APX (tAPX) and stromal APX (sAPX) through alternative splicing. We identified several independent lineage-specific sublocalization cases with specialized paralogues of tAPX and sAPX. We determined that the sublocalization happened through two types of sequence evolution patterns. Our findings suggest that the divergence through sublocalization is key to the retention of paralogous cpAPX genes in angiosperms.

Concerted Divergence after Gene Duplication in Polycomb Repressive Complexes.

Duplicated genes are a major contributor to genome evolution and phenotypic novelty. There are multiple possible evolutionary fates of duplicated genes. Here, we provide an example of concerted divergence of simultaneously duplicated genes whose products function in the same complex. We studied POLYCOMB REPRESSIVE COMPLEX2 (PRC2) in Brassicaceae. The VERNALIZATION (VRN)-PRC2 complex contains VRN2 and SWINGER (SWN), and both genes were duplicated during a whole-genome duplication to generate FERTILIZATION INDEPENDENT SEED2 (FIS2) and MEDEA (MEA), which function in the Brassicaceae-specific FIS-PRC2 complex that regulates seed development. We examined the expression of FIS2, MEA, and their paralogs, compared their cytosine and histone methylation patterns, and analyzed the sequence evolution of the genes. We found that FIS2 and MEA have reproductive-specific expression patterns that are correlated and derived from the broadly expressed VRN2 and SWN in outgroup species. In vegetative tissues of Arabidopsis (Arabidopsis thaliana), repressive methylation marks are enriched in FIS2 and MEA, whereas active marks are associated with their paralogs. We detected comparable accelerated amino acid substitution rates in FIS2 and MEA but not in their paralogs. We also show divergence patterns of the PRC2-associated VERNALIZATION5/VIN3-LIKE2 that are similar to FIS2 and MEA These lines of evidence indicate that FIS2 and MEA have diverged in concert, resulting in functional divergence of the PRC2 complexes in Brassicaceae. This type of concerted divergence is a previously unreported fate of duplicated genes. In addition, the Brassicaceae-specific FIS-PRC2 complex modified the regulatory pathways in female gametophyte and seed development.

Evolution of Cis-Regulatory Elements and Regulatory Networks in Duplicated Genes of Arabidopsis.

Plant genomes contain large numbers of duplicated genes that contribute to the evolution of new functions. Following duplication, genes can exhibit divergence in their coding sequence and their expression patterns. Changes in the cis-regulatory element landscape can result in changes in gene expression patterns. High-throughput methods developed recently can identify potential cis-regulatory elements on a genome-wide scale. Here, we use a recent comprehensive data set of DNase I sequencing-identified cis-regulatory binding sites (footprints) at single-base-pair resolution to compare binding sites and network connectivity in duplicated gene pairs in Arabidopsis (Arabidopsis thaliana). We found that duplicated gene pairs vary greatly in their cis-regulatory element architecture, resulting in changes in regulatory network connectivity. Whole-genome duplicates (WGDs) have approximately twice as many footprints in their promoters left by potential regulatory proteins than do tandem duplicates (TDs). The WGDs have a greater average number of footprint differences between paralogs than TDs. The footprints, in turn, result in more regulatory network connections between WGDs and other genes, forming denser, more complex regulatory networks than shown by TDs. When comparing regulatory connections between duplicates, WGDs had more pairs in which the two genes are either partially or fully diverged in their network connections, but fewer genes with no network connections than the TDs. There is evidence of younger TDs and WGDs having fewer unique connections compared with older duplicates. This study provides insights into cis-regulatory element evolution and network divergence in duplicated genes.

Duplicate gene divergence by changes in microRNA binding sites in Arabidopsis and Brassica.

Gene duplication provides large numbers of new genes that can lead to the evolution of new functions. Duplicated genes can diverge by changes in sequences, expression patterns, and functions. MicroRNAs play an important role in the regulation of gene expression in many eukaryotes. After duplication, two paralogs may diverge in their microRNA binding sites, which might impact their expression and function. Little is known about conservation and divergence of microRNA binding sites in duplicated genes in plants. We analyzed microRNA binding sites in duplicated genes in Arabidopsis thaliana and Brassica rapa. We found that duplicates are more often targeted by microRNAs than singletons. The vast majority of duplicated genes in A. thaliana with microRNA binding sites show divergence in those sites between paralogs. Analysis of microRNA binding sites in genes derived from the ancient whole-genome triplication in B. rapa also revealed extensive divergence. Paralog pairs with divergent microRNA binding sites show more divergence in expression patterns compared with paralog pairs with the same microRNA binding sites in Arabidopsis. Close to half of the cases of binding site divergence are caused by microRNAs that are specific to the Arabidopsis genus, indicating evolutionarily recent gain of binding sites after target gene duplication. We also show rapid evolution of microRNA binding sites in a jacalin gene family. Our analyses reveal a dynamic process of changes in microRNA binding sites after gene duplication in Arabidopsis and highlight the role of microRNA regulation in the divergence and contrasting evolutionary fates of duplicated genes.

Transcriptome analysis indicates considerable divergence in alternative splicing between duplicated genes in Arabidopsis thaliana.

Gene and genome duplication events have created a large number of new genes in plants that can diverge by evolving new expression profiles and functions (neofunctionalization) or dividing extant ones (subfunctionalization). Alternative splicing (AS) generates multiple types of mRNA from a single type of pre-mRNA by differential intron splicing. It can result in new protein isoforms or downregulation of gene expression by transcript decay. Using RNA-seq, we investigated the degree to which alternative splicing patterns are conserved between duplicated genes in Arabidopsis thaliana. Our results revealed that 30% of AS events in α-whole-genome duplicates and 33% of AS events in tandem duplicates are qualitatively conserved within leaf tissue. Loss of ancestral splice forms, as well as asymmetric gain of new splice forms, may account for this divergence. Conserved events had different frequencies, as only 31% of shared AS events in α-whole-genome duplicates and 41% of shared AS events in tandem duplicates had similar frequencies in both paralogs, indicating considerable quantitative divergence. Analysis of published RNA-seq data from nonsense-mediated decay (NMD) mutants indicated that 85% of α-whole-genome duplicates and 89% of tandem duplicates have diverged in their AS-induced NMD. Our results indicate that alternative splicing shows a high degree of divergence between paralogs such that qualitatively conserved alternative splicing events tend to have quantitative divergence. Divergence in AS patterns between duplicates may be a mechanism of regulating expression level divergence.

Protein subcellular relocalization of duplicated genes in Arabidopsis.

Gene duplications during eukaroytic evolution, by successive rounds of polyploidy and by smaller scale duplications, have provided an enormous reservoir of new genes for the evolution of new functions. Preservation of many duplicated genes can be ascribed to changes in sequences, expression patterns, and functions. Protein subcellular relocalization (protein targeting to a new location within the cell) is another way that duplicated genes can diverge. We studied subcellular relocalization of gene pairs duplicated during the evolution of the Brassicaceae including gene pairs from the alpha whole genome duplication that occurred at the base of the family. We analyzed experimental localization data from green fluorescent protein experiments for 128 duplicate pairs in Arabidopsis thaliana, revealing 19 pairs with subcellular relocalization. Many more of the duplicate pairs with relocalization than with the same localization showed an accelerated rate of amino acid sequence evolution in one duplicate, and one gene showed evidence for positive selection. We studied six duplicate gene pairs in more detail. We used gene family analysis with several pairs to infer which gene shows relocalization. We identified potential sequence mutations through comparative analysis that likely result in relocalization of two duplicated gene products. We show that four cases of relocalization have new expression patterns, compared with orthologs in outgroup species, including two with novel expression in pollen. This study provides insights into subcellular relocalization of evolutionarily recent gene duplicates and features of genes whose products have been relocalized.

Frequent changes in expression profile and accelerated sequence evolution of duplicated imprinted genes in arabidopsis.

Eukaryotic genomes have large numbers of duplicated genes that can evolve new functions or expression patterns by changes in coding and regulatory sequences, referred to as neofunctionalization. In flowering plants, some duplicated genes are imprinted in the endosperm, where only one allele is expressed depending on its parental origin. We found that 125 imprinted genes in Arabidopsis arose from gene duplication events during the evolution of the Brassicales. Analyses of 46 gene pairs duplicated by an ancient whole-genome duplication (alpha WGD) indicated that many imprinted genes show an accelerated rate of amino acid changes compared with their paralogs. Analyses of microarray expression data from 63 organ types and developmental stages indicated that many imprinted genes have expression patterns restricted to flowers and/or seeds in contrast to their broadly expressed paralogs. Assays of expression in orthologs from outgroup species revealed that some imprinted genes have acquired an organ-specific expression pattern restricted to flowers and/or seeds. The changes in expression pattern and the accelerated sequence evolution in the imprinted genes suggest that some of them may have undergone neofunctionalization. The imprinted genes MPC, HOMEODOMAIN GLABROUS6 (HDG6), and HDG3 are particularly interesting cases that have different functions from their paralogs. This study indicates that a large number of imprinted genes in Arabidopsis are evolutionarily recent duplicates and that many of them show changes in expression profiles and accelerated sequence evolution. Acquisition of imprinting is a mode of duplicate gene divergence in plants that is more common than previously thought.

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids compared with their parents from natural populations.

Hybridization is a prominent process among natural plant populations that can result in phenotypic novelty, heterosis, and changes in gene expression. The effects of intraspecific hybridization on F1 hybrid gene expression were investigated using parents from divergent, natural populations of Cirsium arvense, an invasive Compositae weed. Using an RNA-seq approach, the expression of 68,746 unigenes was quantified in parents and hybrids. The expression levels of 51% of transcripts differed between parents, a majority of which had less than 1.25× fold-changes. More unigenes had higher expression in the invasive parent (P1) than the noninvasive parent (P2). Of those that were divergently expressed between parents, 10% showed additive and 81% showed nonadditive (transgressive or dominant) modes of gene action in the hybrids. A majority of the dominant cases had P2-like expression patterns in the hybrids. Comparisons of allele-specific expression also enabled a survey of cis- and trans-regulatory effects. Cis- and trans-regulatory divergence was found at 70% and 68% of 62,281 informative single-nucleotide polymorphism sites, respectively. Of the 17% of sites exhibiting both cis- and trans-effects, a majority (70%) had antagonistic regulatory interactions (cis x trans); trans-divergence tended to drive higher expression of the P1 allele, whereas cis-divergence tended to increase P2 transcript abundance. Trans-effects correlated more highly than cis with parental expression divergence and accounted for a greater proportion of the regulatory divergence at sites with additive compared with nonadditive inheritance patterns. This study explores the nature of, and types of mechanisms underlying, expression changes that occur in upon intraspecific hybridization in natural populations.

Features of evolutionarily conserved alternative splicing events between Brassica and Arabidopsis.

Alternative splicing (AS) generates multiple types of mRNA from a single type of pre-mRNA by differential intron splicing. It can result in new protein isoforms or down-regulation of gene expression by transcript decay. The evolutionary conservation of AS events in plants is largely unexplored and only a small number of AS events have been identified as conserved between divergent species. We performed a large-scale analysis of cDNA data from Brassica and Arabidopsis to identify and further characterize conserved AS events. We identified 537 conserved AS events in 485 genes. Alternative donor and acceptor events are significantly overrepresented among conserved events, whereas intron retention and exon skipping events are underrepresented. Conserved AS events are significantly shorter, less likely to be in the 3'UTR, and they are enriched for genes whose products function in the chloroplast. AS modified a functional domain for about half of the genes with conserved events. We further characterized three genes with conserved AS events. This study identifies many AS events that are conserved between Brassica and Arabidopsis, revealing features of conserved AS events. Many of the conserved AS events may have important, but uncharacterized, functions.

Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants.

The importance of gene gain through duplication has long been appreciated. In contrast, the importance of gene loss has only recently attracted attention. Indeed, studies in organisms ranging from plants to worms and humans suggest that duplication of some genes might be better tolerated than that of others. Here we have undertaken a large-scale study to investigate the existence of duplication-resistant genes in the sequenced genomes of 20 flowering plants. We demonstrate that there is a large set of genes that is convergently restored to single-copy status following multiple genome-wide and smaller scale duplication events. We rule out the possibility that such a pattern could be explained by random gene loss only and therefore propose that there is selection pressure to preserve such genes as singletons. This is further substantiated by the observation that angiosperm single-copy genes do not comprise a random fraction of the genome, but instead are often involved in essential housekeeping functions that are highly conserved across all eukaryotes. Furthermore, single-copy genes are generally expressed more highly and in more tissues than non-single-copy genes, and they exhibit higher sequence conservation. Finally, we propose different hypotheses to explain their resistance against duplication.

Genomics of Compositae weeds: EST libraries, microarrays, and evidence of introgression.

PREMISE OF STUDY: Weeds cause considerable environmental and economic damage. However, genomic characterization of weeds has lagged behind that of model plants and crop species. Here we describe the development of genomic tools and resources for 11 weeds from the Compositae family that will serve as a basis for subsequent population and comparative genomic analyses. Because hybridization has been suggested as a stimulus for the evolution of invasiveness, we also analyze these genomic data for evidence of hybridization. METHODS: We generated 22 expressed sequence tag (EST) libraries for the 11 targeted weeds using Sanger, 454, and Illumina sequencing, compared the coverage and quality of sequence assemblies, and developed NimbleGen microarrays for expression analyses in five taxa. When possible, we also compared the distributions of Ks values between orthologs of congeneric taxa to detect and quantify hybridization and introgression. RESULTS: Gene discovery was enhanced by sequencing from multiple tissues, normalization of cDNA libraries, and especially greater sequencing depth. However, assemblies from short sequence reads sometimes failed to resolve close paralogs. Substantial introgression was detected in Centaurea and Helianthus, but not in Ambrosia and Lactuca. CONCLUSIONS: Transcriptome sequencing using next-generation platforms has greatly reduced the cost of genomic studies of nonmodel organisms, and the ESTs and microarrays reported here will accelerate evolutionary and molecular investigations of Compositae weeds. Our study also shows how ortholog comparisons can be used to approximately estimate the genome-wide extent of introgression and to identify genes that have been exchanged between hybridizing taxa.

Organ and cell type-specific complementary expression patterns and regulatory neofunctionalization between duplicated genes in Arabidopsis thaliana.

Duplicated genes can contribute to the evolution of new functions and they are common in eukaryotic genomes. After duplication, genes can show divergence in their sequence and/or expression patterns. Qualitative complementary expression, or reciprocal expression, is when only one copy is expressed in some organ or tissue types and only the other copy is expressed in others, indicative of regulatory subfunctionalization or neofunctionalization. From analyses of two microarray data sets with 83 different organ types, developmental stages, and cell types in Arabidopsis thaliana, we determined that 30% of whole-genome duplicate pairs and 38% of tandem duplicate pairs show reciprocal expression patterns. We reconstructed the ancestral state of expression patterns to infer that considerably more cases of reciprocal expression resulted from gain of a new expression pattern (regulatory neofunctionalization) than from partitioning of ancestral expression patterns (regulatory subfunctionalization). Pollen was an especially common organ type for expression gain, resulting in contrasting expression of some duplicates in pollen. Many of the gene pairs with reciprocal expression showed asymmetric sequence rate evolution, consistent with neofunctionalization, and the more rapidly evolving copy often showed a more restricted expression pattern. A gene with reciprocal expression in pollen, involved in brassinosteroid signal transduction, has evolved more rapidly than its paralog, and it shows evidence for a new function in pollen. This study indicates the evolutionary importance of reciprocal expression patterns between gene duplicates, showing that they are common, often associated with regulatory neofunctionalization, and may be a factor allowing for retention and divergence of duplicated genes.

Extensive changes to alternative splicing patterns following allopolyploidy in natural and resynthesized polyploids.

Polyploidy has been a common process during the evolution of eukaryotes, especially plants, leading to speciation and the evolution of new gene functions. Gene expression levels and patterns can change, and gene silencing can occur in allopolyploids--phenomena sometimes referred to as "transcriptome shock." Alternative splicing (AS) creates multiple mature mRNAs from a single type of precursor mRNA. Here we examined the evolution of AS patterns after polyploidy, with natural and two resynthesized allotetraploid Brassica napus lines, using RT-PCR and sequencing assays of 82 AS events in duplicated gene pairs (homeologs). Comparing the AS patterns between the two homeologs in natural B. napus revealed that many of the gene pairs show different AS patterns, with a few showing variation that was organ specific or induced by abiotic stress treatments. In the resynthesized allotetraploids, 26-30% of the duplicated genes showed changes in AS compared with the parents, including many cases of AS event loss after polyploidy. Parallel losses of many AS events after allopolyploidy were detected in the two independently resynthesized lines. More changes occurred in parallel between the two lines than changes specific to each line. The PASTICCINO gene showed partitioning of two AS events between the two homeologs in the resynthesized allopolyploids. AS changes after allopolyploidy were much more common than homeolog silencing. Our findings indicate that AS patterns can change rapidly after polyploidy, that many genes are affected, and that AS changes are an important component of the transcriptome shock experienced by new allopolyploids.

Differential contributions to the transcriptome of duplicated genes in response to abiotic stresses in natural and synthetic polyploids.

Polyploidy has occurred throughout plant evolution and can result in considerable changes to gene expression when it takes place and over evolutionary time. Little is known about the effects of abiotic stress conditions on duplicate gene expression patterns in polyploid plants. We examined the expression patterns of 60 duplicated genes in leaves, roots and cotyledons of allotetraploid Gossypium hirsutum in response to five abiotic stress treatments (heat, cold, drought, high salt and water submersion) using single-strand conformation polymorphism assays, and 20 genes in a synthetic allotetraploid. Over 70% of the genes showed stress-induced changes in the relative expression levels of the duplicates under one or more stress treatments with frequent variability among treatments. Twelve pairs showed opposite changes in expression levels in response to different abiotic stress treatments. Stress-induced expression changes occurred in the synthetic allopolyploid, but there was little correspondence in patterns between the natural and synthetic polyploids. Our results indicate that abiotic stress conditions can have considerable effects on duplicate gene expression in a polyploid, with the effects varying by gene, stress and organ type. Differential expression in response to environmental stresses may be a factor in the preservation of some duplicated genes in polyploids.

Dramatic change in function and expression pattern of a gene duplicated by polyploidy created a paternal effect gene in the Brassicaceae.

New gene formation by polyploidy has been an ongoing process during the evolution of various eukaryotes that has contributed greatly to the large number of genes in their genomes. After duplication, some genes that are retained can acquire new functions or expression patterns, or subdivide their functions or expression patterns between duplicates. Here, we show that SHORT SUSPENSOR (SSP) and Brassinosteroid Kinase 1 (BSK1) are paralogs duplicated by a polyploidy event that occurred in the Brassicaceae family about 23 Ma. SSP is involved in paternal control of zygote elongation in Arabidopsis thaliana by transcription in the sperm cells of pollen and then translation in the zygote, whereas BSK1 is involved in brassinosteroid signal transduction. Comparative analysis of expression in 63 different organs and developmental stages revealed that BSK1 and SSP have opposite expression patterns in pollen compared with all other parts of the plant. We determined that BSK1 retains the ancestral expression pattern and function. Thus, SSP has diverged in function after duplication from a component of the brassinosteroid signaling pathway to a paternal regulator of the timing of zygote elongation. The ancestral function of SSP was lost by deletions in the kinase domain. Our sequence rate analysis revealed that SSP but not BSK1 has experienced a greatly accelerated rate of amino acid sequence changes and relaxation of purifying selection. In addition, SSP has been duplicated to create a new gene (SSP-like1) with a completely different expression pattern, a shorter coding sequence that has lost a critical functional domain, and a greatly accelerated rate of amino acid sequence evolution along with evidence for positive selection, together indicative of neofunctionalization. This study illustrates two dramatic examples of neofunctionalization following gene duplication by complete changes in expression pattern and function. In addition, our findings indicate that paternal control of zygote elongation by SSP is an evolutionarily recent innovation in the Brassicaceae family.