Cultivated strawberry emerged from the hybridization of two wild octoploid species, both descendants from the merger of four diploid progenitor species into a single nucleus more than 1 million years ago. Here we report a near-complete chromosome-scale assembly for cultivated octoploid strawberry (Fragaria × ananassa) and uncovered the origin and evolutionary processes that shaped this complex allopolyploid. We identified the extant relatives of each diploid progenitor species and provide support for the North American origin of octoploid strawberry. We examined the dynamics among the four subgenomes in octoploid strawberry and uncovered the presence of a single dominant subgenome with significantly greater gene content, gene expression abundance, and biased exchanges between homoeologous chromosomes, as compared with the other subgenomes. Pathway analysis showed that certain metabolomic and disease-resistance traits are largely controlled by the dominant subgenome. These findings and the reference genome should serve as a powerful platform for future evolutionary studies and enable molecular breeding in strawberry.

Highlights d m 6 A sites are highly maintained at different Arabidopsis developmental stages d Many m 6 A sites stabilize Arabidopsis adult leaf transcripts d m 6 A stabilizes mRNAs through inhibition of local ribonucleolytic cleavage d Dynamic, stress-specific m 6 A sites stabilize stress response protein transcripts

SUMMARY 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.

Telomerase replenishes telomere tracts by reiteratively copying its RNA template, TER. Unlike other model organisms, Arabidopsis thaliana harbors two divergent TER genes. However, only TER1 is required for telomere maintenance. Here we examine the function of TER2. We show that TER2 is spliced and its 39 end is truncated in vivo to generate a third TER isoform, TER2 S . TERT preferentially associates with TER2 > TER1 > TER2 S . Moreover, TER2 and TER2 S assemble with Ku and POT1b (protection of telomeres), forming RNP (ribonucleoprotein) complexes distinct from TER1 RNP. Plants null for TER2 display increased telomerase enzyme activity, while TER2 overexpression inhibits telomere synthesis from TER1 and leads to telomere shortening. These findings argue that TER2 negatively regulates telomerase by sequestering TERT in a nonproductive RNP complex. Introduction of DNA double-strand breaks by zeocin leads to an immediate and specific spike in TER2 and a concomitant decrease in telomerase enzyme activity. This response is not triggered by replication stress or telomere dysfunction and is abrogated in ter2 mutants. We conclude that Arabidopsis telomerase is modulated by TER2, a novel DNA damageinduced noncoding RNA that works in concert with the canonical TER to promote genome integrity.

Dyskerin binds the H/ACA box of human telomerase RNA and is a core telomerase subunit required for RNP biogenesis and enzyme function in vivo. Missense mutations in dyskerin result in dyskeratosis congenita, a complex syndrome characterized by bone marrow failure, telomerase enzyme deficiency, and progressive telomere shortening. Here we demonstrate that dyskerin also contributes to telomere maintenance in Arabidopsis thaliana. We report that both AtNAP57, the Arabidopsis dyskerin homolog, and AtTERT, the telomerase catalytic subunit, accumulate in the plant nucleolus, and AtNAP57 associates with active telomerase RNP particles in an RNA-dependent manner. Furthermore, AtNAP57 interacts in vitro with AtPOT1a, a novel component of Arabidopsis telomerase. Although a null mutation in AtNAP57 is lethal, AtNAP57, like AtTERT, is not haploinsufficient for telomere maintenance in Arabidopsis. However, introduction of an AtNAP57 allele containing a T66A mutation decreased telomerase activity in vitro, disrupted telomere length regulation on individual chromosome ends in vivo, and established a new, shorter telomere length set point. These results imply that T66A NAP57 behaves as a dominant-negative inhibitor of telomerase. We conclude that dyskerin is a conserved component of the telomerase RNP complex in higher eukaryotes that is required for maximal enzyme activity in vivo.

Conversion of a double-strand break into a telomere is a dangerous, potentially lethal event. However, little is known about the mechanism and control of de novo telomere formation (DNTF). DNTF can be instigated by the insertion of a telomere repeat array (TRA) into the host genome, which seeds the formation of a new telomere, resulting in chromosome truncation. Such events are rare and concentrated at chromosome ends. Here, we introduce tetraploid Arabidopsis thaliana as a robust genetic model for DNTF. Transformation of a 2.6-kb TRA into tetraploid plants resulted in a DNTF efficiency of 56%, fivefold higher than in diploid plants and 50-fold higher than in human cells. DNTF events were recovered across the entire genome, indicating that genetic redundancy facilitates recovery of DNTF events. Although TRAs as short as 100 bp seeded new telomeres, these tracts were unstable unless they were extended above a 1-kb size threshold. Unexpectedly, DNTF efficiency increased in plants lacking telomerase, and DNTF rates were lower in plants null for Ku70 or Lig4, components of the nonhomologous end-joining repair pathway. We conclude that multiple competing pathways modulate DNTF, and that tetraploid Arabidopsis will be a powerful model for elucidating the molecular details of these processes.

After transcription, a messenger RNA (mRNA) is further post‐transcriptionally regulated by several features including RNA secondary structure and covalent RNA modifications (specifically N 6 ‐methyladenosine, m 6 A). Both RNA secondary structure and m 6 A 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 m 6 A 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 m 6 A is anti‐correlated with RNA secondary structure in a condition‐dependent manner, with salt‐specific m 6 A correlated with a decrease in mRNA secondary structure during salt stress. Taken together, we suggest that salt‐specific m 6 A 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.

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.