Mating system is associated with seed phenotypes upon loss of RNA-directed DNA methylation in Brassicaceae.

In plants, de novo DNA methylation is guided by 24-nt short interfering (si)RNAs in a process called RNA-directed DNA methylation (RdDM). Primarily targeted at transposons, RdDM causes transcriptional silencing and can indirectly influence expression of neighboring genes. During reproduction, a small number of siRNA loci are dramatically upregulated in the maternally derived seed coat, suggesting that RdDM might have a special function during reproduction. However, the developmental consequence of RdDM has been difficult to dissect because disruption of RdDM does not result in overt phenotypes in Arabidopsis (Arabidopsis thaliana), where the pathway has been most thoroughly studied. In contrast, Brassica rapa mutants lacking RdDM have a severe seed production defect, which is determined by the maternal sporophytic genotype. To explore the factors that underlie the different phenotypes of these species, we produced RdDM mutations in 3 additional members of the Brassicaceae family: Camelina sativa, Capsella rubella, and Capsella grandiflora. Among these 3 species, only mutations in the obligate outcrosser, C. grandiflora, displayed a seed production defect similar to Brassica rapa mutants, suggesting that mating system is a key determinant for reproductive phenotypes in RdDM mutants.

Evolutionary analysis of the LORELEI gene family in plants reveals regulatory subfunctionalization.

A signaling complex comprising members of the LORELEI (LRE)-LIKE GPI-anchored protein (LLG) and Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) families perceive RAPID ALKALINIZATION FACTOR (RALF) peptides and regulate growth, reproduction, immunity, and stress responses in Arabidopsis (Arabidopsis thaliana). Genes encoding these proteins are members of multigene families in most angiosperms and could generate thousands of signaling complex variants. However, the links between expansion of these gene families and the functional diversification of this critical signaling complex as well as the evolutionary factors underlying the maintenance of gene duplicates remain unknown. Here, we investigated LLG gene family evolution by sampling land plant genomes and explored the function and expression of angiosperm LLGs. We found that LLG diversity within major land plant lineages is primarily due to lineage-specific duplication events, and that these duplications occurred both early in the history of these lineages and more recently. Our complementation and expression analyses showed that expression divergence (i.e. regulatory subfunctionalization), rather than functional divergence, explains the retention of LLG paralogs. Interestingly, all but one monocot and all eudicot species examined had an LLG copy with preferential expression in male reproductive tissues, while the other duplicate copies showed highest levels of expression in female or vegetative tissues. The single LLG copy in Amborella trichopoda is expressed vastly higher in male compared to in female reproductive or vegetative tissues. We propose that expression divergence plays an important role in retention of LLG duplicates in angiosperms.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

A Conserved Long Intergenic Non-coding RNA Containing snoRNA Sequences, lncCOBRA1, Affects Arabidopsis Germination and Development.

Long non-coding RNAs (lncRNAs) are an increasingly studied group of non-protein coding transcripts with a wide variety of molecular functions gaining attention for their roles in numerous biological processes. Nearly 6,000 lncRNAs have been identified in Arabidopsis thaliana but many have yet to be studied. Here, we examine a class of previously uncharacterized lncRNAs termed CONSERVED IN BRASSICA RAPA (lncCOBRA) transcripts that were previously identified for their high level of sequence conservation in the related crop species Brassica rapa, their nuclear-localization and protein-bound nature. In particular, we focus on lncCOBRA1 and demonstrate that its abundance is highly tissue and developmental specific, with particularly high levels early in germination. lncCOBRA1 contains two snoRNAs domains within it, making it the first sno-lincRNA example in a non-mammalian system. However, we find that it is processed differently than its mammalian counterparts. We further show that plants lacking lncCOBRA1 display patterns of delayed germination and are overall smaller than wild-type plants. Lastly, we identify the proteins that interact with lncCOBRA1 and propose a novel mechanism of lincRNA action in which it may act as a scaffold with the RACK1A protein to regulate germination and development, possibly through a role in ribosome biogenesis.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships.

Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.

N6-methyladenosine and RNA secondary structure affect transcript stability and protein abundance during systemic salt stress in Arabidopsis.

After transcription, a messenger RNA (mRNA) is further post-transcriptionally regulated by several features including RNA secondary structure and covalent RNA modifications (specifically N6-methyladenosine, m6A). Both RNA secondary structure and m6A have been demonstrated to regulate mRNA stability and translation and have been independently linked to plant responses to soil salinity levels. However, the effect of m6A on regulating RNA secondary structure and the combinatorial interplay between these two RNA features during salt stress response has yet to be studied. Here, we globally identify RNA-protein interactions and RNA secondary structure during systemic salt stress. This analysis reveals that RNA secondary structure changes significantly during salt stress, and that it is independent of global changes in RNA-protein interactions. Conversely, we find that m6A is anti-correlated with RNA secondary structure in a condition-dependent manner, with salt-specific m6A correlated with a decrease in mRNA secondary structure during salt stress. Taken together, we suggest that salt-specific m6A deposition and the associated loss of RNA secondary structure results in increases in mRNA stability for transcripts encoding abiotic stress response proteins and ultimately increases in protein levels from these stabilized transcripts. In total, our comprehensive analyses reveal important post-transcriptional regulatory mechanisms involved in plant long-term salt stress response and adaptation.

Evolutionary and biochemical analyses reveal conservation of the Brassicaceae telomerase ribonucleoprotein complex.

The telomerase ribonucleoprotein complex (RNP) is essential for genome stability and performs this role through the addition of repetitive DNA to the ends of chromosomes. The telomerase enzyme is composed of a reverse transcriptase (TERT), which utilizes a template domain in an RNA subunit (TER) to reiteratively add telomeric DNA at the ends of chromosomes. Multiple TERs have been identified in the model plant Arabidopsis thaliana. Here we combine a phylogenetic and biochemical approach to understand how the telomerase RNP has evolved in Brassicaceae, the family that includes A. thaliana. Because of the complex phylogenetic pattern of template domain loss and alteration at the previously characterized A. thaliana TER loci, TER1 and TER2, across the plant family Brassicaceae, we bred double mutants from plants with a template deletion at AtTER1 and T-DNA insertion at AtTER2. These double mutants exhibited no telomere length deficiency, a definitive indication that neither of these loci encode a functional telomerase RNA. Moreover, we determined that the telomerase components TERT, Dyskerin, and the KU heterodimer are under strong purifying selection, consistent with the idea that the TER with which they interact is also conserved. To test this hypothesis further, we analyzed the substrate specificity of telomerase from species across Brassicaceae and determined that telomerase from close relatives bind and extend substrates in a similar manner, supporting the idea that TERs in different species are highly similar to one another and are likely encoded from an orthologous locus. Lastly, TERT proteins from across Brassicaceae were able to complement loss of function tert mutants in vivo, indicating TERTs from other species have the ability to recognize the native TER of A. thaliana. Finally, we immunoprecipitated the telomerase complex and identified associated RNAs via RNA-seq. Using our evolutionary data we constrained our analyses to conserved RNAs within Brassicaceae that contained a template domain. These analyses revealed a highly expressed locus whose disruption by a T-DNA resulted in a telomeric phenotype similar to the loss of other telomerase core proteins, indicating that the RNA has an important function in telomere maintenance.

Divergence-Based Introgression Polarization.

Introgressive hybridization results in the transfer of genetic material between species, often with fitness implications for the recipient species. The development of statistical methods for detecting the signatures of historical introgression in whole-genome data has been a major area of focus. Although existing techniques are able to identify the taxa that exchanged genes during introgression using a four-taxon system, most methods do not explicitly distinguish which taxon served as donor and which as recipient during introgression (i.e., polarization of introgression directionality). Existing methods that do polarize introgression are often only able to do so when there is a fifth taxon available and that taxon is sister to one of the taxa involved in introgression. Here, we present divergence-based introgression polarization (DIP), a method for polarizing introgression using patterns of sequence divergence across whole genomes, which operates in a four-taxon context. Thus, DIP can be applied to infer the directionality of introgression when additional taxa are not available. We use simulations to show that DIP can polarize introgression and identify potential sources of bias in the assignment of directionality, and we apply DIP to a well-described hominin introgression event.

Duplication and functional divergence of a calcium sensor in the Brassicaceae.

The presence of varied numbers of CALCINEURIN B-LIKE10 (CBL10) calcium sensor genes in species across the Brassicaceae and the demonstrated role of CBL10 in salt tolerance in Arabidopsis thaliana and Eutrema salsugineum provided a unique opportunity to determine if CBL10 function is modified in different species and linked to salt tolerance. Salinity effects on species growth and cross-species complementation were used to determine the extent of conservation and divergence of CBL10 function in four species representing major lineages within the core Brassicaceae (A. thaliana, E. salsugineum, Schrenkiella parvula, and Sisymbrium irio) as well as the first diverging lineage (Aethionema arabicum). Evolutionary and functional analyses indicate that CBL10 duplicated within expanded lineage II of the Brassicaceae and that, while portions of CBL10 function are conserved across the family, there are species-specific variations in CBL10 function. Paralogous CBL10 genes within a species diverged in expression and function probably contributing to the maintenance of the duplicated gene pairs. Orthologous CBL10 genes diverged in function in a species-specific manner, suggesting that functions arose post-speciation. Multiple CBL10 genes and their functional divergence may have expanded calcium-mediated signaling responses and contributed to the ability of certain members of the Brassicaceae to maintain growth in salt-affected soils.

N6-Methyladenosine Inhibits Local Ribonucleolytic Cleavage to Stabilize mRNAs in Arabidopsis.

N6-methyladenosine (m6A) is a dynamic, reversible, covalently modified ribonucleotide that occurs predominantly toward 3' ends of eukaryotic mRNAs and is essential for their proper function and regulation. In Arabidopsis thaliana, many RNAs contain at least one m6A site, yet the transcriptome-wide function of m6A remains mostly unknown. Here, we show that many m6A-modified mRNAs in Arabidopsis have reduced abundance in the absence of this mark. The decrease in abundance is due to transcript destabilization caused by cleavage occurring 4 or 5 nt directly upstream of unmodified m6A sites. Importantly, we also find that, upon agriculturally relevant salt treatment, m6A is dynamically deposited on and stabilizes transcripts encoding proteins required for salt and osmotic stress response. Overall, our findings reveal that m6A generally acts as a stabilizing mark through inhibition of site-specific cleavage in plant transcriptomes, and this mechanism is required for proper regulation of the salt-stress-responsive transcriptome.

Evidence for a Unique DNA-Dependent RNA Polymerase in Cereal Crops.

Gene duplication is an important driver for the evolution of new genes and protein functions. Duplication of DNA-dependent RNA polymerase (Pol) II subunits within plants led to the emergence of RNA Pol IV and V complexes, each of which possess unique functions necessary for RNA-directed DNA Methylation. Comprehensive identification of Pol V subunit orthologs across the monocot radiation revealed a duplication of the largest two subunits within the grasses (Poaceae), including critical cereal crops. These paralogous Pol subunits display sequence conservation within catalytic domains, but their carboxy terminal domains differ in length and character of the Ago-binding platform, suggesting unique functional interactions. Phylogenetic analysis of the catalytic region indicates positive selection on one paralog following duplication, consistent with retention via neofunctionalization. Positive selection on residue pairs that are predicted to interact between subunits suggests that paralogous subunits have evolved specific assembly partners. Additional Pol subunits as well as Pol-interacting proteins also possess grass-specific paralogs, supporting the hypothesis that a novel Pol complex with distinct function has evolved in the grass family, Poaceae.

Phylogenetics of Camelina Crantz. (Brassicaceae) and insights on the origin of gold-of-pleasure (Camelina sativa).

Camelina sativa (false flax or gold-of-pleasure) is an Old World oilseed crop that fell out of use in the mid 20th Century but has recently gained renewed interest as a biofuel source. The crop is hexaploid, and its relationship to its diploid and polyploid congeners has remained unresolved. Using 54 accessions representing five species sampled across Camelina's center of diversity in Turkey and the Caucasus, we performed phylogenetic and genetic diversity analyses using RADseq genotyping and ITS sequencing. Flow cytometry was performed to assess relationships between genome size and phylogenetic groupings. Accessions fell into distinct, highly-supported clades that accord with named species, indicating that morphological characters can reliably distinguish members of the genus. A phylogenetically distinct lineage from Turkey may represent a currently unrecognized diploid species. In most analyses, C. sativa accessions nest within those of C. microcarpa, suggesting that the crop is descended from this wild hexaploid species. This inference is further supported by their similar genome size, and by lower genetic diversity in C. sativa, which is consistent with a domestication bottleneck. These analyses provide the first definitive phylogeny of C. sativa and its wild relatives, and they point to C. microcarpa as the crop's wild ancestor.

Evolinc: A Tool for the Identification and Evolutionary Comparison of Long Intergenic Non-coding RNAs.

Long intergenic non-coding RNAs (lincRNAs) are an abundant and functionally diverse class of eukaryotic transcripts. Reported lincRNA repertoires in mammals vary, but are commonly in the thousands to tens of thousands of transcripts, covering ~90% of the genome. In addition to elucidating function, there is particular interest in understanding the origin and evolution of lincRNAs. Aside from mammals, lincRNA populations have been sparsely sampled, precluding evolutionary analyses focused on their emergence and persistence. Here we present Evolinc, a two-module pipeline designed to facilitate lincRNA discovery and characterize aspects of lincRNA evolution. The first module (Evolinc-I) is a lincRNA identification workflow that also facilitates downstream differential expression analysis and genome browser visualization of identified lincRNAs. The second module (Evolinc-II) is a genomic and transcriptomic comparative analysis workflow that determines the phylogenetic depth to which a lincRNA locus is conserved within a user-defined group of related species. Here we validate lincRNA catalogs generated with Evolinc-I against previously annotated Arabidopsis and human lincRNA data. Evolinc-I recapitulated earlier findings and uncovered an additional 70 Arabidopsis and 43 human lincRNAs. We demonstrate the usefulness of Evolinc-II by examining the evolutionary histories of a public dataset of 5,361 Arabidopsis lincRNAs. We used Evolinc-II to winnow this dataset to 40 lincRNAs conserved across species in Brassicaceae. Finally, we show how Evolinc-II can be used to recover the evolutionary history of a known lincRNA, the human telomerase RNA (TERC). These latter analyses revealed unexpected duplication events as well as the loss and subsequent acquisition of a novel TERC locus in the lineage leading to mice and rats. The Evolinc pipeline is currently integrated in CyVerse's Discovery Environment and is free for use by researchers.

Evolutionary interplay between sister cytochrome P450 genes shapes plasticity in plant metabolism.

Expansion of the cytochrome P450 gene family is often proposed to have a critical role in the evolution of metabolic complexity, in particular in microorganisms, insects and plants. However, the molecular mechanisms underlying the evolution of this complexity are poorly understood. Here we describe the evolutionary history of a plant P450 retrogene, which emerged and underwent fixation in the common ancestor of Brassicales, before undergoing tandem duplication in the ancestor of Brassicaceae. Duplication leads first to gain of dual functions in one of the copies. Both sister genes are retained through subsequent speciation but eventually return to a single copy in two of three diverging lineages. In the lineage in which both copies are maintained, the ancestral functions are split between paralogs and a novel function arises in the copy under relaxed selection. Our work illustrates how retrotransposition and gene duplication can favour the emergence of novel metabolic functions.

Evolution of Arabidopsis protection of telomeres 1 alters nucleic acid recognition and telomerase regulation.

Protection of telomeres (POT1) binds chromosome ends, recognizing single-strand telomeric DNA via two oligonucleotide/oligosaccharide binding folds (OB-folds). The Arabidopsis thaliana POT1a and POT1b paralogs are atypical: they do not exhibit telomeric DNA binding, and they have opposing roles in regulating telomerase activity. AtPOT1a stimulates repeat addition processivity of the canonical telomerase enzyme, while AtPOT1b interacts with a regulatory lncRNA that represses telomerase activity. Here, we show that OB1 of POT1a, but not POT1b, has an intrinsic affinity for telomeric DNA. DNA binding was dependent upon a highly conserved Phe residue (F65) that in human POT1 directly contacts telomeric DNA. F65A mutation of POT1aOB1 abolished DNA binding and diminished telomerase repeat addition processivity. Conversely, AtPOT1b and other POT1b homologs from Brassicaceae and its sister family, Cleomaceae, naturally bear a non-aromatic amino acid at this position. By swapping Val (V63) with Phe, AtPOT1bOB1 gained the capacity to bind telomeric DNA and to stimulate telomerase repeat addition processivity. We conclude that, in the context of DNA binding, variation at a single amino acid position promotes divergence of the AtPOT1b paralog from the ancestral POT1 protein.

A Genomic Analysis of Factors Driving lincRNA Diversification: Lessons from Plants.

Transcriptomic analyses from across eukaryotes indicate that most of the genome is transcribed at some point in the developmental trajectory of an organism. One class of these transcripts is termed long intergenic noncoding RNAs (lincRNAs). Recently, attention has focused on understanding the evolutionary dynamics of lincRNAs, particularly their conservation within genomes. Here, we take a comparative genomic and phylogenetic approach to uncover factors influencing lincRNA emergence and persistence in the plant family Brassicaceae, to which Arabidopsis thaliana belongs. We searched 10 genomes across the family for evidence of > 5000 lincRNA loci from A. thaliana From loci conserved in the genomes of multiple species, we built alignments and inferred phylogeny. We then used gene tree/species tree reconciliation to examine the duplication history and timing of emergence of these loci. Emergence of lincRNA loci appears to be linked to local duplication events, but, surprisingly, not whole genome duplication events (WGD), or transposable elements. Interestingly, WGD events are associated with the loss of loci for species having undergone relatively recent polyploidy. Lastly, we identify 1180 loci of the 6480 previously annotated A. thaliana lincRNAs (18%) with elevated levels of conservation. These conserved lincRNAs show higher expression, and are enriched for stress-responsiveness and cis-regulatory motifs known as conserved noncoding sequences (CNSs). These data highlight potential functional pathways and suggest that CNSs may regulate neighboring genes at both the genomic and transcriptomic level. In sum, we provide insight into processes that may influence lincRNA diversification by providing an evolutionary context for previously annotated lincRNAs.

The Argonaute-binding platform of NRPE1 evolves through modulation of intrinsically disordered repeats.

Argonaute (Ago) proteins are important effectors in RNA silencing pathways, but they must interact with other machinery to trigger silencing. Ago hooks have emerged as a conserved motif responsible for interaction with Ago proteins, but little is known about the sequence surrounding Ago hooks that must restrict or enable interaction with specific Argonautes. Here we investigated the evolutionary dynamics of an Ago-binding platform in NRPE1, the largest subunit of RNA polymerase V. We compared NRPE1 sequences from > 50 species, including dense sampling of two plant lineages. This study demonstrates that the Ago-binding platform of NRPE1 retains Ago hooks, intrinsic disorder, and repetitive character while being highly labile at the sequence level. We reveal that loss of sequence conservation is the result of relaxed selection and frequent expansions and contractions of tandem repeat arrays. These factors allow a complete restructuring of the Ago-binding platform over 50-60 million yr. This evolutionary pattern is also detected in a second Ago-binding platform, suggesting it is a general mechanism. The presence of labile repeat arrays in all analyzed NRPE1 Ago-binding platforms indicates that selection maintains repetitive character, potentially to retain the ability to rapidly restructure the Ago-binding platform.

Ancient Origin and Recent Innovations of RNA Polymerase IV and V.

Small RNA-mediated chromatin modification is a conserved feature of eukaryotes. In flowering plants, the short interfering (si)RNAs that direct transcriptional silencing are abundant and subfunctionalization has led to specialized machinery responsible for synthesis and action of these small RNAs. In particular, plants possess polymerase (Pol) IV and Pol V, multi-subunit homologs of the canonical DNA-dependent RNA Pol II, as well as specialized members of the RNA-dependent RNA Polymerase (RDR), Dicer-like (DCL), and Argonaute (AGO) families. Together these enzymes are required for production and activity of Pol IV-dependent (p4-)siRNAs, which trigger RNA-directed DNA methylation (RdDM) at homologous sequences. p4-siRNAs accumulate highly in developing endosperm, a specialized tissue found only in flowering plants, and are rare in nonflowering plants, suggesting that the evolution of flowers might coincide with the emergence of specialized RdDM machinery. Through comprehensive identification of RdDM genes from species representing the breadth of the land plant phylogeny, we describe the ancient origin of Pol IV and Pol V, suggesting that a nearly complete and functional RdDM pathway could have existed in the earliest land plants. We also uncover innovations in these enzymes that are coincident with the emergence of seed plants and flowering plants, and recent duplications that might indicate additional subfunctionalization. Phylogenetic analysis reveals rapid evolution of Pol IV and Pol V subunits relative to their Pol II counterparts and suggests that duplicates were retained and subfunctionalized through Escape from Adaptive Conflict. Evolution within the carboxy-terminal domain of the Pol V largest subunit is particularly striking, where illegitimate recombination facilitated extreme sequence divergence.

Evolution of the Telomere-Associated Protein POT1a in Arabidopsis thaliana Is Characterized by Positive Selection to Reinforce Protein-Protein Interaction.

Gene duplication is a major driving force in genome evolution. Here, we explore the nature and origin of the POT1 gene duplication in Arabidopsis thaliana. Protection of Telomeres (POT1) is a conserved multifunctional protein that modulates telomerase activity and its engagement with telomeres. Arabidopsis thaliana encodes two divergent POT1 paralogs termed AtPOT1a and AtPOT1b. AtPOT1a positively regulates telomerase activity, whereas AtPOT1b is proposed to negatively regulate telomerase and promote chromosome end protection. Phylogenetic analysis uncovered two independent POT1 duplication events in the plant kingdom, including one at the base of Brassicaceae. Tests for positive selection implemented in PAML revealed that the Brassicaceae POT1a lineage experienced positive selection postduplication and identified three amino acid residues with signatures of positive selection. A sensitive and quantitative genetic complementation assay was developed to assess POT1a function in A. thaliana. The assay showed that AtPOT1a is functionally distinct from single-copy POT1 genes in other plants. Moreover, for two of the sites with a strong signature of positive selection, substitutions that swap the amino acids in AtPOT1a for residues found in AtPOT1b dramatically compromised AtPOT1a function in vivo. In vitro-binding studies demonstrated that all three sites under positive selection specifically enhance the AtPOT1a interaction with CTC1, a core component of the highly conserved CST (CTC1/STN1/TEN1) telomere protein complex. Our results reveal a molecular mechanism for the role of these positively selected sites in AtPOT1a. The data also provide an important empirical example to refine theories of duplicate gene retention, as the outcome of positive selection here appears to be reinforcement of an ancestral function, rather than neofunctionalization. We propose that this outcome may not be unusual when the duplicated protein is a component of a multisubunit complex whose function is in part specified by other members.

Global analysis of the RNA-protein interaction and RNA secondary structure landscapes of the Arabidopsis nucleus.

Posttranscriptional regulation in eukaryotes requires cis- and trans-acting features and factors including RNA secondary structure and RNA-binding proteins (RBPs). However, a comprehensive view of the structural and RBP interaction landscape of nuclear RNAs has yet to be compiled for any organism. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approaches to globally profile these features in Arabidopsis seedling nuclei in vivo. We reveal anticorrelated patterns of secondary structure and RBP binding throughout nuclear mRNAs that demarcate sites of alternative splicing and polyadenylation. We also uncover a collection of protein-bound sequence motifs, and identify their structural contexts, co-occurrences in transcripts encoding functionally related proteins, and interactions with putative RBPs. Finally, using these motifs, we find that the chloroplast RBP CP29A also interacts with nuclear mRNAs. In total, we provide a simultaneous view of the RNA secondary structure and RBP interaction landscapes in a eukaryotic nucleus.