Cycles of satellite and transposon evolution in Arabidopsis centromeres.

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Standing genetic variation fuels rapid evolution of herbicide resistance in blackgrass.

Repeated herbicide applications in agricultural fields exert strong selection on weeds such as blackgrass (Alopecurus myosuroides), which is a major threat for temperate climate cereal crops. This inadvertent selection pressure provides an opportunity for investigating the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much target-site resistance (TSR) relies on de novo mutations vs. standing variation is important for developing strategies to manage herbicide resistance. We first generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Next, through empirical data in the form of highly accurate long-read amplicons of alleles encoding acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) variants, we showed that most populations with resistance due to TSR mutations-23 out of 27 and six out of nine populations for ACCase and ALS, respectively-contained at least two TSR haplotypes, indicating that soft sweeps are the norm. Finally, through forward-in-time simulations, we inferred that TSR is likely to mainly result from standing genetic variation, with only a minor role for de novo mutations.

Leveraging orthology within maize and Arabidopsis QTL to identify genes affecting natural variation in gravitropism.

Plants typically orient their organs with respect to the Earth's gravity field by a dynamic process called gravitropism. To discover conserved genetic elements affecting seedling root gravitropism, we measured the process in a set of Zea mays (maize) recombinant inbred lines with machine vision and compared the results with those obtained in a similar study of Arabidopsis thaliana. Each of the several quantitative trait loci that we mapped in both species spanned many hundreds of genes, too many to test individually for causality. We reasoned that orthologous genes may be responsible for natural variation in monocot and dicot root gravitropism. If so, pairs of orthologous genes affecting gravitropism may be present within the maize and Arabidopsis QTL intervals. A reciprocal comparison of sequences within the QTL intervals identified seven pairs of such one-to-one orthologs. Analysis of knockout mutants demonstrated a role in gravitropism for four of the seven: CCT2 functions in phosphatidylcholine biosynthesis, ATG5 functions in membrane remodeling during autophagy, UGP2 produces the substrate for cellulose and callose polymer extension, and FAMA is a transcription factor. Automated phenotyping enabled this discovery of four naturally varying components of a conserved process (gravitropism) by making it feasible to conduct the same large-scale experiment in two species.

Plant genetic effects on microbial hubs impact host fitness in repeated field trials.

Although complex interactions between hosts and microbial associates are increasingly well documented, we still know little about how and why hosts shape microbial communities in nature. In addition, host genetic effects on microbial communities vary widely depending on the environment, obscuring conclusions about which microbes are impacted and which plant functions are important. We characterized the leaf microbiota of 200 Arabidopsis thaliana genotypes in eight field experiments and detected consistent host effects on specific, broadly distributed microbial species (operational taxonomic unit [OTUs]). Host genetic effects disproportionately influenced central ecological hubs, with heritability of particular OTUs declining with their distance from the nearest hub within the microbial network. These host effects could reflect either OTUs preferentially associating with specific genotypes or differential microbial success within them. Host genetics associated with microbial hubs explained over 10% of the variation in lifetime seed production among host genotypes across sites and years. We successfully cultured one of these microbial hubs and demonstrated its growth-promoting effects on plants in sterile conditions. Finally, genome-wide association mapping identified many putatively causal genes with small effects on the relative abundance of microbial hubs across sites and years, and these genes were enriched for those involved in the synthesis of specialized metabolites, auxins, and the immune system. Using untargeted metabolomics, we corroborate the consistent association between variation in specialized metabolites and microbial hubs across field sites. Together, our results reveal that host genetic variation impacts the microbial communities in consistent ways across environments and that these effects contribute to fitness variation among host genotypes.

Pushing the limits of HiFi assemblies reveals centromere diversity between two Arabidopsis thaliana genomes.

Although long-read sequencing can often enable chromosome-level reconstruction of genomes, it is still unclear how one can routinely obtain gapless assemblies. In the model plant Arabidopsis thaliana, other than the reference accession Col-0, all other accessions de novo assembled with long-reads until now have used PacBio continuous long reads (CLR). Although these assemblies sometimes achieved chromosome-arm level contigs, they inevitably broke near the centromeres, excluding megabases of DNA from analysis in pan-genome projects. Since PacBio high-fidelity (HiFi) reads circumvent the high error rate of CLR technologies, albeit at the expense of read length, we compared a CLR assembly of accession Eyach15-2 to HiFi assemblies of the same sample. The use of five different assemblers starting from subsampled data allowed us to evaluate the impact of coverage and read length. We found that centromeres and rDNA clusters are responsible for 71% of contig breaks in the CLR scaffolds, while relatively short stretches of GA/TC repeats are at the core of >85% of the unfilled gaps in our best HiFi assemblies. Since the HiFi technology consistently enabled us to reconstruct gapless centromeres and 5S rDNA clusters, we demonstrate the value of the approach by comparing these previously inaccessible regions of the genome between the Eyach15-2 accession and the reference accession Col-0.

Deep haplotype analyses of target-site resistance locus ACCase in blackgrass enabled by pool-based amplicon sequencing.

Rapid adaptation of weeds to herbicide applications in agriculture through resistance development is a widespread phenomenon. In particular, the grass Alopecurus myosuroides is an extremely problematic weed in cereal crops with the potential to manifest resistance in only a few generations. Target-site resistances (TSRs), with their strong phenotypic response, play an important role in this rapid adaptive response. Recently, using PacBio's long-read amplicon sequencing technology in hundreds of individuals, we were able to decipher the genomic context in which TSR mutations occur. However, sequencing individual amplicons are costly and time-consuming, thus impractical to implement for other resistance loci or applications. Alternatively, pool-based approaches overcome these limitations and provide reliable allele frequencies, although at the expense of not preserving haplotype information. In this proof-of-concept study, we sequenced with PacBio High Fidelity (HiFi) reads long-range amplicons (13.2 kb), encompassing the entire ACCase gene in pools of over 100 individuals, and resolved them into haplotypes using the clustering algorithm PacBio amplicon analysis (pbaa), a new application for pools in plants and other organisms. From these amplicon pools, we were able to recover most haplotypes from previously sequenced individuals of the same population. In addition, we analysed new pools from a Germany-wide collection of A. myosuroides populations and found that TSR mutations originating from soft sweeps of independent origin were common. Forward-in-time simulations indicate that TSR haplotypes will persist for decades even at relatively low frequencies and without selection, highlighting the importance of accurate measurement of TSR haplotype prevalence for weed management.

A Truncated Singleton NLR Causes Hybrid Necrosis in Arabidopsis thaliana.

Hybrid necrosis in plants arises from conflict between divergent alleles of immunity genes contributed by different parents, resulting in autoimmunity. We investigate a severe hybrid necrosis case in Arabidopsis thaliana, where the hybrid does not develop past the cotyledon stage and dies 3 weeks after sowing. Massive transcriptional changes take place in the hybrid, including the upregulation of most NLR (nucleotide-binding site leucine-rich repeat) disease-resistance genes. This is due to an incompatible interaction between the singleton TIR-NLR gene DANGEROUS MIX 10 (DM10), which was recently relocated from a larger NLR cluster, and an unlinked locus, DANGEROUS MIX 11 (DM11). There are multiple DM10 allelic variants in the global A. thaliana population, several of which have premature stop codons. One of these, which has a truncated LRR-PL (leucine-rich repeat [LRR]-post-LRR) region, corresponds to the DM10 risk allele. The DM10 locus and the adjacent genomic region in the risk allele carriers are highly differentiated from those in the nonrisk carriers in the global A. thaliana population, suggesting that this allele became geographically widespread only relatively recently. The DM11 risk allele is much rarer and found only in two accessions from southwestern Spain-a region from which the DM10 risk haplotype is absent-indicating that the ranges of DM10 and DM11 risk alleles may be nonoverlapping.

Genetic and epigenetic variation in 5S ribosomal RNA genes reveals genome dynamics in Arabidopsis thaliana.

Organized in tandem repeat arrays in most eukaryotes and transcribed by RNA polymerase III, expression of 5S rRNA genes is under epigenetic control. To unveil mechanisms of transcriptional regulation, we obtained here in depth sequence information on 5S rRNA genes from the Arabidopsis thaliana genome and identified differential enrichment in epigenetic marks between the three 5S rDNA loci situated on chromosomes 3, 4 and 5. We reveal the chromosome 5 locus as the major source of an atypical, long 5S rRNA transcript characteristic of an open chromatin structure. 5S rRNA genes from this locus translocated in the Landsberg erecta ecotype as shown by linkage mapping and chromosome-specific FISH analysis. These variations in 5S rDNA locus organization cause changes in the spatial arrangement of chromosomes in the nucleus. Furthermore, 5S rRNA gene arrangements are highly dynamic with alterations in chromosomal positions through translocations in certain mutants of the RNA-directed DNA methylation pathway and important copy number variations among ecotypes. Finally, variations in 5S rRNA gene sequence, chromatin organization and transcripts indicate differential usage of 5S rDNA loci in distinct ecotypes. We suggest that both the usage of existing and new 5S rDNA loci resulting from translocations may impact neighboring chromatin organization.

Phytohormone sensing in the biotrophic fungus Ustilago maydis - the dual role of the transcription factor Rss1.

The phenolic compound salicylic acid (SA) is a key signalling molecule regulating local and systemic plant defense responses, mainly against biotrophs. Many microbial organisms, including pathogens, share the ability to degrade SA. However, the mechanism by which they perceive SA is unknown. Here we show that Ustilago maydis, the causal agent of corn smut disease, employs a so far uncharacterized SA sensing mechanism. We identified and characterized the novel SA sensing regulator, Rss1, a binuclear zinc cluster protein with dual functions as putative SA receptor and transcriptional activator regulating genes important for SA and tryptophan degradation. Rss1 represents a major component in the identified SA sensing pathway during the fungus' saprophytic stage. However, Rss1 does not have a detectable impact on virulence. The data presented in this work indicate that alternative or redundant sensing cascades exist that regulate the expression of SA-responsive genes in U. maydis during its pathogenic development.

Structural variation and DNA methylation shape the centromere-proximal meiotic crossover landscape in Arabidopsis.

BACKGROUND: Centromeres load kinetochore complexes onto chromosomes, which mediate spindle attachment and allow segregation during cell division. Although centromeres perform a conserved cellular function, their underlying DNA sequences are highly divergent within and between species. Despite variability in DNA sequence, centromeres are also universally suppressed for meiotic crossover recombination, across eukaryotes. However, the genetic and epigenetic factors responsible for suppression of centromeric crossovers remain to be completely defined. RESULTS: To explore the centromere-proximal meiotic recombination landscape, we map 14,397 crossovers against fully assembled Arabidopsis thaliana (A. thaliana) genomes. A. thaliana centromeres comprise megabase satellite repeat arrays that load nucleosomes containing the CENH3 histone variant. Each chromosome contains a structurally polymorphic region of ~3-4 megabases, which lack crossovers and include the satellite arrays. This polymorphic region is flanked by ~1-2 megabase low-recombination zones. These recombination-suppressed regions are enriched for Gypsy/Ty3 retrotransposons, and additionally contain expressed genes with high genetic diversity that initiate meiotic recombination, yet do not crossover. We map crossovers at high-resolution in proximity to CEN3, which resolves punctate centromere-proximal hotspots that overlap gene islands embedded in heterochromatin. Centromeres are densely DNA methylated and the recombination landscape is remodelled in DNA methylation mutants. We observe that the centromeric low-recombining zones decrease and increase crossovers in CG (met1) and non-CG (cmt3) mutants, respectively, whereas the core non-recombining zones remain suppressed. CONCLUSION: Our work relates the genetic and epigenetic organization of A. thaliana centromeres and flanking pericentromeric heterochromatin to the zones of crossover suppression that surround the CENH3-occupied satellite repeat arrays.

RNA-mediated trans-communication can establish paramutation at the b1 locus in maize.

Paramutation is the epigenetic transfer of information between alleles that leads to the heritable change of expression of one allele. Paramutation at the b1 locus in maize requires seven noncoding tandem repeat (b1TR) sequences located approximately 100 kb upstream of the transcription start site of b1, and mutations in several genes required for paramutation implicate an RNA-mediated mechanism. The mediator of paramutation (mop1) gene, which encodes a protein closely related to RNA-dependent RNA polymerases, is absolutely required for paramutation. Herein, we investigate the potential function of mop1 and the siRNAs that are produced from the b1TR sequences. Production of siRNAs from the b1TR sequences depends on a functional mop1 gene, but transcription of the repeats is not dependent on mop1. Further nuclear transcription assays suggest that the b1TR sequences are likely transcribed predominantly by RNA polymerase II. To address whether production of b1TR-siRNAs correlated with paramutation, we examined siRNA production in alleles that cannot undergo paramutation. Alleles that cannot participate in paramutation also produce b1TR-siRNAs, suggesting that b1TR-siRNAs are not sufficient for paramutation in the tissues analyzed. However, when b1TR-siRNAs are produced from a transgene expressing a hairpin RNA, b1 paramutation can be recapitulated. We hypothesize that either the b1TR-siRNAs or the dsRNA template mediates the trans-communication between the alleles that establishes paramutation. In addition, we uncovered a role for mop1 in the biogenesis of a subset of microRNAs (miRNAs) and show that it functions at the level of production of the primary miRNA transcripts.

Extensive sequence duplication in Arabidopsis revealed by pseudo-heterozygosity.

BACKGROUND: It is apparent that genomes harbor much structural variation that is largely undetected for technical reasons. Such variation can cause artifacts when short-read sequencing data are mapped to a reference genome. Spurious SNPs may result from mapping of reads to unrecognized duplicated regions. Calling SNP using the raw reads of the 1001 Arabidopsis Genomes Project we identified 3.3 million (44%) heterozygous SNPs. Given that Arabidopsis thaliana (A. thaliana) is highly selfing, and that extensively heterozygous individuals have been removed, we hypothesize that these SNPs reflected cryptic copy number variation. RESULTS: The heterozygosity we observe consists of particular SNPs being heterozygous across individuals in a manner that strongly suggests it reflects shared segregating duplications rather than random tracts of residual heterozygosity due to occasional outcrossing. Focusing on such pseudo-heterozygosity in annotated genes, we use genome-wide association to map the position of the duplicates. We identify 2500 putatively duplicated genes and validate them using de novo genome assemblies from six lines. Specific examples included an annotated gene and nearby transposon that transpose together. We also demonstrate that cryptic structural variation produces highly inaccurate estimates of DNA methylation polymorphism. CONCLUSIONS: Our study confirms that most heterozygous SNP calls in A. thaliana are artifacts and suggest that great caution is needed when analyzing SNP data from short-read sequencing. The finding that 10% of annotated genes exhibit copy-number variation, and the realization that neither gene- nor transposon-annotation necessarily tells us what is actually mobile in the genome suggests that future analyses based on independently assembled genomes will be very informative.

The integrated LIM-peptidase domain of the CSA1-CHS3/DAR4 paired immune receptor detects changes in DA1 peptidase inhibitors in Arabidopsis.

White blister rust, caused by the oomycete Albugo candida, is a widespread disease of Brassica crops. The Brassica relative Arabidopsis thaliana uses the paired immune receptor complex CSA1-CHS3/DAR4 to resist Albugo infection. The CHS3/DAR4 sensor NLR, which functions together with its partner, the helper NLR CSA1, carries an integrated domain (ID) with homology to DA1 peptidases. Using domain swaps with several DA1 homologs, we show that the LIM-peptidase domain of the family member CHS3/DAR4 functions as an integrated decoy for the family member DAR3, which interacts with and inhibits the peptidase activities of the three closely related peptidases DA1, DAR1, and DAR2. Albugo infection rapidly lowers DAR3 levels and activates DA1 peptidase activity, thereby promoting endoreduplication of host tissues to support pathogen growth. We propose that the paired immune receptor CSA1-CHS3/DAR4 detects the actions of a putative Albugo effector that reduces DAR3 levels, resulting in defense activation.

The Phaseolus vulgaris miR159a precursor encodes a second differentially expressed microRNA.

Plant microRNAs originate from a stem-loop structured single-stranded RNA precursor. Each stem-loop is processed to generate a mature microRNA that is recruited to an ARGONAUTE-containing multi-protein complex to direct silencing of its target mRNA. Here we report that the conserved plant miR159a precursor produces a second 21-nt long RNA with the properties of a microRNA. Its presence in different plant species is supported by its conservation in the stem-loop position and expression as determined by northern blot analysis. We show that successive processing by DCL1 produces this novel microRNA from the same precursor as miR159a. In contrast to the low levels observed in other plant models for the equivalent of miR159.2, in P. vulgaris, the accumulation of miR159.2 is easily detectable and when compared to miR159a, their expression patterns are distinct in different organs and growth conditions. Further evidence of the functionality of miR159.2 comes from its association with silencing complexes as demonstrated by co-immunoprecipitation experiments using an AGO1-specific antibody and processing of an artificial GFP reporter construct containing a complementary target sequence. These results indicate that the second small RNA corresponds to a microRNA, at least partially independent of miR159 activity, and that in plants a miRNA precursor may encode multiple regulatory small RNAs.

Characterization of the Different Levels of Variation in 45S rRNA Genes.

Although commonly regarded as nearly identical in sequence, 45S ribosomal RNA genes represent a massive source of genetic variation at different levels. Unfortunately, due to their repetitive nature and the difficulty to assemble their multiple copies in tandem, these important genomic elements remain largely unexplored in sequencing projects. Here, I describe how to exploit next generation sequencing data to estimate their copy number in an organism and detect true polymorphic sites within and among individuals. Furthermore, for species that carry multiple 45S ribosomal RNA gene clusters, I show how to make use of experimental populations to assign some of these variants to their cluster of origin.

It Is Just a Matter of Time: Balancing Homologous Recombination and Non-homologous End Joining at the rDNA Locus During Meiosis.

Ribosomal RNA genes (rDNAs) are located in large domains of hundreds of rDNA units organized in a head-to-tail manner. The proper and stable inheritance of rDNA clusters is of paramount importance for survival. Yet, these highly repetitive elements pose a potential risk to the genome since they can undergo non-allelic exchanges. Here, we review the current knowledge of the organization of the rDNA clusters in Arabidopsis thaliana and their stability during meiosis. Recent findings suggest that during meiosis, all rDNA loci are embedded within the nucleolus favoring non-homologous end joining (NHEJ) as a repair mechanism, while DNA repair via homologous recombination (HR) appears to be a rare event. We propose a model where (1) frequent meiotic NHEJ events generate abundant single nucleotide polymorphisms and insertions/deletions within the rDNA, resulting in a heterogeneous population of rDNA units and (2) rare HR events dynamically change rDNA unit numbers, only to be observed in large populations over many generations. Based on the latest efforts to delineate the entire rDNA sequence in A. thaliana, we discuss evidence supporting this model. The results compiled so far draw a surprising picture of rDNA sequence heterogeneity between individual units. Furthermore, rDNA cluster sizes have been recognized as relatively stable when observing less than 10 generations, yet emerged as major determinant of genome size variation between different A. thaliana ecotypes. The sequencing efforts also revealed that transcripts from the diverse rDNA units yield heterogenous ribosome populations with potential functional implications. These findings strongly motivate further research to understand the mechanisms that maintain the metastable state of rDNA loci.

Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress.

MicroRNAs (miRNAs) are small RNA molecules recognized as important regulators of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data. We are interested in the identification of miRNAs in Phaseolus vulgaris (common bean) to uncover different plant strategies to cope with adverse conditions and because of its relevance as a crop in developing countries. Here we present the identification of conserved and candidate novel miRNAs in P. vulgaris present in different organs and growth conditions, including drought, abscisic acid treatment, and Rhizobium infection. We also identified cDNA sequences in public databases that represent the corresponding miRNA precursors. In addition, we predicted and validated target mRNAs amongst reported EST and cDNAs for P. vulgaris. We propose that the novel miRNAs present in common bean and other legumes, are involved in regulation of legume-specific processes including adaptation to diverse external cues.

Unstable Inheritance of 45S rRNA Genes in Arabidopsis thaliana.

The considerable genome size variation in Arabidopsis thaliana has been shown largely to be due to copy number variation (CNV) in 45S ribosomal RNA (rRNA) genes. Surprisingly, attempts to map this variation by means of genome-wide association studies (GWAS) failed to identify either of the two likely sources, namely the nucleolus organizer regions (NORs). Instead, GWAS implicated a trans-acting locus, as if rRNA gene CNV was a phenotype rather than a genotype. To explain these results, we investigated the inheritance and stability of rRNA gene copy number using the variety of genetic resources available in A. thaliana - F2 crosses, recombinant inbred lines, the multiparent advanced-generation inter-cross population, and mutation accumulation lines. Our results clearly show that rRNA gene CNV can be mapped to the NORs themselves, with both loci contributing equally to the variation. However, NOR size is unstably inherited, and dramatic copy number changes are visible already within tens of generations, which explains why it is not possible to map the NORs using GWAS. We did not find any evidence of trans-acting loci in crosses, which is also expected since changes due to such loci would take very many generations to manifest themselves. rRNA gene copy number is thus an interesting example of "missing heritability"-a trait that is heritable in pedigrees, but not in the general population.

Verification of Arabidopsis stock collections using SNPmatch, a tool for genotyping high-plexed samples.

Large-scale studies such as the Arabidopsis thaliana '1,001 Genomes' Project require routine genotyping of stocks to avoid sample contamination. To genotype samples efficiently and economically, sequencing must be inexpensive and data processing simple. Here we present SNPmatch, a tool that identifies strains (or inbred lines, or accessions) by matching them to a SNP database. We tested the tool by performing low-coverage resequencing of over 2,000 strains from our lab seed stock collection. SNPmatch correctly genotyped samples from 1-fold coverage sequencing data, and could also identify the parents of F1 or F2 individuals. SNPmatch can be run either on the command line or through AraGeno (https://arageno.gmi.oeaw.ac.at), a web interface that permits sample genotyping from a user-uploaded VCF or BED file.

Epistatic and allelic interactions control expression of ribosomal RNA gene clusters in Arabidopsis thaliana.

BACKGROUND: Ribosomal RNA (rRNA) accounts for the majority of the RNA in eukaryotic cells, and is encoded by hundreds to thousands of nearly identical gene copies, only a subset of which are active at any given time. In Arabidopsis thaliana, 45S rRNA genes are found in two large ribosomal DNA (rDNA) clusters and little is known about the contribution of each to the overall transcription pattern in the species. RESULTS: By taking advantage of genome sequencing data from the 1001 Genomes Consortium, we characterize rRNA gene sequence variation within and among accessions. Notably, variation is not restricted to the pre-rRNA sequences removed during processing, but it is also present within the highly conserved ribosomal subunits. Through linkage mapping we assign these variants to a particular rDNA cluster unambiguously and use them as reporters of rDNA cluster-specific expression. We demonstrate that rDNA cluster-usage varies greatly among accessions and that rDNA cluster-specific expression and silencing is controlled via genetic interactions between entire rDNA cluster haplotypes (alleles). CONCLUSIONS: We show that rRNA gene cluster expression is controlled via complex epistatic and allelic interactions between rDNA haplotypes that apparently regulate the entire rRNA gene cluster. Furthermore, the sequence polymorphism we discovered implies that the pool of rRNA in a cell may be heterogeneous, which could have functional consequences.