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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Distinct roles for SOS1 in the convergent evolution of salt tolerance in Eutrema salsugineum and Schrenkiella parvula.

Eutrema salsugineum and Schrenkiella parvula are salt-tolerant relatives of the salt-sensitive species Arabidopsis thaliana. An important component of salt tolerance is the regulation of Na(+) ion homeostasis, which occurs in part through proteins encoded by the Cation/Proton Antiporter-1 (CPA1) gene family. We used a combination of evolutionary and functional analyses to examine the role of CPA1 genes in the salt tolerance of E. salsugineum and Sc. parvula, and found evidence that changes in CPA1-mediated Na(+) extrusion may contribute to the salt tolerance of both species. Specifically, we found that a member of the CPA1 family, the Na(+)/H(+) antiporter gene Salt Overly Sensitive 1 (SOS1), evolved under positive selection in E. salsugineum. In the absence of activation by the SOS2 kinase/SOS3 calcium-binding protein complex, SOS1 from E. salsugineum (EsSOS1) confers greater salt tolerance than SOS1 from Sc. parvula (SpSOS1) and Ar. thaliana (AtSOS1) when expressed in a salt-sensitive strain of Saccharomyces cerevisiae. A single amino acid change in the putative autoinhibitory domain is required but not sufficient for the enhanced salt tolerance conferred by EsSOS1. When activated by SOS2 and SOS3, both EsSOS1 and SpSOS1 confer greater salt tolerance than AtSOS1. Enhanced SOS1-mediated Na(+) extrusion therefore appears to contribute to the salt tolerance of both E. salsugineum and Sc. parvula, although through apparently different mechanisms.

Extending the model of Arabidopsis telomere length and composition across Brassicaceae.

Telomeres are repetitive TG-rich DNA elements essential for maintaining the stability of genomes and replicative capacity of cells in almost all eukaryotes. Most of what is known about telomeres in plants comes from the angiosperm Arabidopsis thaliana, which has become an important comparative model for telomere biology. Arabidopsis tolerates numerous insults to its genome, many of which are catastrophic or lethal in other eukaryotic systems such as yeast and vertebrates. Despite the importance of Arabidopsis in establishing a model for the structure and regulation of plant telomeres, only a handful of studies have used this information to assay components of telomeres from across land plants, or even among the closest relatives of Arabidopsis in the plant family Brassicaceae. Here, we determined how well Arabidopsis represents Brassicaceae by comparing multiple aspects of telomere biology in species that represent major clades in the family tree. Specifically, we determined the telomeric repeat sequence, measured bulk telomere length, and analyzed variation in telomere length on syntenic chromosome arms. In addition, we used a phylogenetic approach to infer the evolutionary history of putative telomere-binding proteins, CTC1, STN1, TEN1 (CST), telomere repeat-binding factor like (TRFL), and single Myb histone (SMH). Our analyses revealed conservation of the telomeric DNA repeat sequence, but considerable variation in telomere length among the sampled species, even in comparisons of syntenic chromosome arms. We also found that the single-stranded and double-stranded telomeric DNA-binding complexes CST and TRFL, respectively, differ in their pattern of gene duplication and loss. The TRFL and SMH gene families have undergone numerous duplication events, and these duplicate copies are often retained in the genome. In contrast, CST components occur as single-copy genes in all sampled genomes, even in species that experienced recent whole genome duplication events. Taken together, our results place the Arabidopsis model in the context of other species in Brassicaceae, making the family the best characterized plant group in regard to telomere architecture.

Diversification times among Brassica (Brassicaceae) crops suggest hybrid formation after 20 million years of divergence.

PREMISE OF THE STUDY: Cruciferous vegetables, many of which are in the genus Brassica (Brassicaceae), are prized for their nutritive value and have been cultivated for thousands of years. There are numerous wild northwestern Mediterranean species in the tribe Brassiceae, and it is therefore assumed this center of diversity is also the region of origin. Within the tribe, the Nigra and Oleracea clades contain the three diploid Brassica crops, B. oleracea, B. rapa, and B. nigra. These three species hybridized in the past to form the tetraploid crop species B. juncea, B. carinata, and B. napus. Collectively, these crop Brassicas have been thought to be closely related because they can still hybridize. METHODS: Using a combination of molecular phylogenetics, diversification analysis, and historical biogeography, we evaluated the relationships and origins of four nested clades: the tribe Brassiceae, the Nigra-Oleracea clade, the core Oleracea (includes B. oleracea + B. rapa and their respective wild relatives), and Brassica oleracea and relatives. KEY RESULTS: We found evidence that the tribe originated around the intersection forming between the Arabian Peninsula and Saharan Africa approximately 24 million years ago (Mya). Our data also suggest that the maternal genomes of the three diploid crop Brassicas are not closely related and that the Nigra-Oleracea clade diverged 20 Mya. Finally, our analyses indicate that the core Oleracea lineage giving rise to B. oleracea + B. rapa originated ≈3 Mya in the northeastern Mediterranean, from where ancestors of B. oleracea spread through Europe and B. rapa to Asia. CONCLUSIONS: These results challenge previous hypotheses about the biogeographic origins of the tribe Brassiceae and the crop Brassica species and appear to be correlated with major geological and climatic events in the Mediterranean basin.

Dated molecular phylogenies indicate a Miocene origin for Arabidopsis thaliana.

Dated molecular phylogenies are the basis for understanding species diversity and for linking changes in rates of diversification with historical events such as restructuring in developmental pathways, genome doubling, or dispersal onto a new continent. Valid fossil calibration points are essential to the accurate estimation of divergence dates, but for many groups of flowering plants fossil evidence is unavailable or limited. Arabidopsis thaliana, the primary genetic model in plant biology and the first plant to have its entire genome sequenced, belongs to one such group, the plant family Brassicaceae. Thus, the timing of A. thaliana evolution and the history of its genome have been controversial. We bring previously overlooked fossil evidence to bear on these questions and find the split between A. thaliana and Arabidopsis lyrata occurred about 13 Mya, and that the split between Arabidopsis and the Brassica complex (broccoli, cabbage, canola) occurred about 43 Mya. These estimates, which are two- to threefold older than previous estimates, indicate that gene, genomic, and developmental evolution occurred much more slowly than previously hypothesized and that Arabidopsis evolved during a period of warming rather than of cooling. We detected a 2- to 10-fold shift in species diversification rates on the branch uniting Brassicaceae with its sister families. The timing of this shift suggests a possible impact of the Cretaceous-Paleogene mass extinction on their radiation and that Brassicales codiversified with pierid butterflies that specialize on mustard-oil-producing plants.

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.

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.

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.

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.