Natural variation in biogenesis efficiency of individual Arabidopsis thaliana microRNAs

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

After silencing suppression: miRNA targets strike back

Within the continuous tug-of-war between plants and microbes, RNA silencing stands out as a key battleground. Pathogens, in their quest to colonize host plants, have evolved a diverse arsenal of silencing suppressors as a common strategy to undermine the host's RNA silencing-based defenses. When RNA silencing malfunctions in the host, genes that are usually targeted and silenced by microRNAs (miRNAs) become active and can contribute to the reprogramming of host cells, providing an additional defense mechanism. A growing body of evidence suggests that miRNAs may act as intracellular sensors to enable a rapid response to pathogen threats. Herein we review how plant miRNA targets play a crucial role in immune responses against different pathogens.

Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na⁺ -K⁺ /H⁺ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

The PtoTCP20-miR396d-PtoGRF15 module regulates secondary vascular development in Populus

Secondary vascular development is a key biological characteristic of woody plants and the basis of wood formation. Our understanding of gene expression regulation and dynamic changes in microRNAs (miRNAs) during secondary vascular development is still limited. Here we present an integrated analysis of the miRNA and mRNA transcriptome of six phase-specific tissues-the shoot apex, procambium, primary vascular tissue, cambium, secondary phloem, and secondary xylem-in Populus tomentosa. Several novel regulatory modules, including the PtoTCP20-miR396d-PtoGRF15 module, were identified during secondary vascular development in Populus. A series of biochemical and molecular experiments confirmed that PtoTCP20 activated transcription of the miR396d precursor gene and that miR396d targeted PtoGRF15 to downregulate its expression. Plants overexpressing miR396d (35S:miR396d) showed enhanced secondary growth and increased xylem production. Conversely, during the transition from primary to secondary vascular development, plants with downregulated PtoTCP20expression (PtoTCP20-SRDX), downregulated miR396 expression (35S:STTM396), and PtoGRF15 overexpression (35S:PtoGRF15) showed delayed secondary growth. Novel regulatory modules were identified by integrated analysis of the miRNA and mRNA transcriptome, and the regulatory role of the PtoTCP20-miR396d-PtoGRF15 signaling cascade in secondary vascular development was validated in Populus, providing information to support improvements in forest cultivation and wood properties.

Genetic variants in microRNA biogenesis genes as novel indicators for secondary growth in Populus

MicroRNAs (miRNAs) function as key regulators of complex traits, but how genetic alterations in miRNA biogenesis genes (miRBGs) affect quantitative variation has not been elucidated. We conducted transcript analyses and association genetics to investigate how miRBGs, miRNA genes (MIRNAs) and their respective targets contribute to secondary growth in a natural population of 435 Populus tomentosa individuals. This analysis identified 29 843 common single-nucleotide polymorphisms (SNPs; frequency > 0.10) within 682 genes (80 miRBGs, 152 MIRNAs, and 457 miRNA targets). Single-SNP association analysis found SNPs in 234 candidate genes exhibited significant additive/dominant effects on phenotypes. Among these, specific candidates that associated with the same traits produced 791 miRBG-MIRNA-target combinations, suggesting possible genetic miRBG-MIRNA and MIRNA-target interactions, providing an important clue for the regulatory mechanisms of miRBGs. Multi-SNP association found 4672 epistatic pairs involving 578 genes that showed significant associations with traits and identified 106 miRBG-MIRNA-target combinations. Two multi-hierarchical networks were constructed based on correlations of miRBG-miRNA and miRNA-target expression to further probe the mechanisms of trait diversity underlying changes in miRBGs. Our study opens avenues for the investigation of miRNA function in perennial plants and underscored miRBGs as potentially modulating quantitative variation in traits.

Extracting genotype information of Arabidopsis thaliana recombinant inbred lines from transcript profiles established with high-density oligonucleotide arrays

Polymorphic probes identified via a sequence-based approach are suitable to infer the genotypes of recombinant inbred lines from hybridisation intensities of GeneChip ^® transcript profiling experiments. The sequences of the probes of the ATH1 GeneChip^® exactly match transcript sequences of the Arabidopsis thaliana reference genome Col-0, whereas nucleotide differences and/or insertions/deletions may be observed for transcripts of other A. thaliana accessions. Individual probes of the GeneChip^® that show sequence polymorphisms between different A. thaliana accessions may serve as single-feature polymorphism (SFP) markers, provided that the sequence changes cause differences in hybridisation intensity for the accessions of interest. A sequence-based approach identified features on the high-density oligonucleotide array that showed sequence polymorphisms between A. thaliana accessions Col-0 and C24. Hybridisation intensities of polymorphic probes were extracted from genome-wide transcript profiles of Col-0/C24 and C24/Col-0 recombinant inbred lines and assessed after standardisation via sliding window analyses to identify SFP markers. The genotypes of the recombinant inbred lines were determined with the SFP markers and the resulting data were integrated with information, which had been established previously with single nucleotide polymorphism and insertion/deletion markers, to enrich the linkage map of the Col-0/C24 and C24/Col-0 recombinant inbred populations. Congruence between the molecular marker map and the sequence maps of the A. thaliana Col-0 chromosomes proved the reliability of the genotype information which was deduced from the transcript profiles of the Col-0/C24 and C24/Col-0 recombinant inbred lines.

microRNA/microRNA* complementarity is important for the regulation pattern of NFYA5 by miR169 under dehydration shock in Arabidopsis

MicroRNAs regulate gene expression at the mRNA and translational levels. Although our previous research showed that expression of miR169 and one of its targets, NFYA5, is down- and up-regulated by drought stress, respectively, the current study shows that expression of both miR169 and NFYA5 are induced by dehydration shock. Unlike overexpression of MIR169a/b, overexpression of MIR169i/l did not decrease NFYA5 transcripts but increased NFYA5 protein levels. NFYA5 protein abundance also increased in mir169a knock-out mutants and 35S::MIR169l mir169a double mutants. When bulge #11 and bulge #16 in the miR169/miR169* duplex were mutated, both NFYA5 transcripts and the corresponding protein levels were lower in 35S::MIR169l*-mut transgenic plants than in wild-type (WT) plants, and the 35S::MIR169l*-mut transgenic plants were as sensitive to drought stress as the 35S::MIR169a plants. The mRNA and protein levels of NFYA5 did not differ substantially between WT and 35S::MIR169a*-mut transgenic plants when the two bulges were introduced in the miR169a/miR169a* duplex. Both bulge #11 and bulge #16 in the miR169/miR169* duplex were essential for different regulation patterns of NFYA5 by miR169a and miR169l. These results increase the understanding of regulatory specialization in one MIR family, and also increase our understanding of the importance of microRNA/microRNA* secondary structure.

Adaptive evolution and functional innovation of Populus-specific recently evolved microRNAs

Lineage-specific microRNAs (miRNAs) undergo rapid turnover during evolution; however, their origin and functional importance have remained controversial. Here, we examine the origin, evolution, and potential roles in local adaptation of Populus-specific miRNAs, which originated after the recent salicoid-specific, whole-genome duplication. RNA sequencing was used to generate extensive, comparable miRNA and gene expression data for six tissues. A natural population of Populus trichocarpa and closely related species were used to study the divergence rates, evolution, and adaptive variation of miRNAs. MiRNAs that originated in 5' untranslated regions had higher expression levels and their expression showed high correlation with their host genes. Compared with conserved miRNAs, a significantly higher proportion of Populus-specific miRNAs appear to target genes that were duplicated in salicoids. Examination of single nucleotide polymorphisms in Populus-specific miRNA precursors showed high amounts of population differentiation. We also characterized the newly emerged MIR6445 family, which could trigger the production of phased small interfering RNAs from NAC mRNAs, which encode a transcription factor with primary roles in a variety of plant developmental processes. Together, these observations provide evolutionary insights into the birth and potential roles of Populus-specific miRNAs in genome maintenance, local adaptation, and functional innovation.

Genetic variations and miRNA-target interactions contribute to natural phenotypic variations in Populus

Variation in regulatory factors, including microRNAs (miRNAs), contributes to variation in quantitative and complex traits. However, in plants, variants in miRNAs and their target genes that contribute to natural phenotypic variation, and the underlying regulatory networks, remain poorly characterized. We investigated the associations and interactions of single-nucleotide polymorphisms (SNPs) in miRNAs and their target genes with phenotypes in 435 individuals from a natural population of Populus. We used RNA-seq to identify 217 miRNAs differentially expressed in a tension wood system, and identified 1196 candidate target genes; degradome sequencing confirmed 60 of the target sites. In addition, 72 miRNA-target pairs showed significant co-expression. Gene ontology (GO) term analysis showed that most of the genes in the co-regulated pairs participate in biological regulation. Genome resequencing found 5383 common SNPs (frequency ≥ 0.05) in 139 miRNAs and 31 037 SNPs in 819 target genes. Single-SNP association analyses identified 232 significant associations between wood traits (P ≤ 0.05) and SNPs in 102 miRNAs and 1387 associations with 478 target genes. Among these, 102 miRNA-target pairs associated with the same traits. Multi-SNP associations found 102 epistatic pairs associated with traits. Furthermore, a reconstructed regulatory network contained 12 significantly co-expressed pairs, including eight miRNAs and nine targets associated with traits. Lastly, both expression and genetic association showed that miR156i, miR156j, miR396a and miR6445b were involved in the formation of tension wood. This study shows that variants in miRNAs and target genes contribute to natural phenotypic variation and annotated roles and interactions of miRNAs and their target genes by genetic association analysis.

Evolutionarily Dynamic, but Robust, Targeting of Resistance Genes by the miR482/2118 Gene Family in the Solanaceae

Plants are exposed to pathogens around the clock. A common resistance response in plants upon pathogen detection is localized cell death. Given the irreversible nature of this response, multiple layers of negative regulation are present to prevent the untimely or misexpression of resistance genes. One layer of negative regulation is provided by a recently discovered microRNA (miRNA) gene family, miR482/2118. This family targets the transcripts of resistance genes in plants. We investigated the evolutionary history and specificity of this miRNA gene family within the Solanaceae. This plant family includes many important crop species, providing a set of well-defined resistance gene repertoires. Across 14 species from the Solanaceae, we identified eight distinct miR482/2118 gene family members. Our studies show conservation of miRNA type and number in the group of wild tomatoes and, to a lesser extent, throughout the Solanaceae. The eight orthologous miRNA gene clusters evolved under different evolutionary constraints, allowing for individual subfunctionalization of the miRNAs. Despite differences in the predicted targeting behavior of each miRNA, the miRNA-R-gene network is robust due to its high degree of interconnectivity and redundant targeting. Our data suggest that the miR482/2118 gene family acts as an evolutionary buffer for R-gene sequence diversity.

Heterochrony underpins natural variation in Cardamine hirsuta leaf form

A key problem in biology is whether the same processes underlie morphological variation between and within species. Here, by using plant leaves as an example, we show that the causes of diversity at these two evolutionary scales can be divergent. Some species like the model plant Arabidopsis thaliana have simple leaves, whereas others like the A. thaliana relative Cardamine hirsuta bear complex leaves comprising leaflets. Previous work has shown that these interspecific differences result mostly from variation in local tissue growth and patterning. Now, by cloning and characterizing a quantitative trait locus (QTL) for C. hirsuta leaf shape, we find that a different process, age-dependent progression of leaf form, underlies variation in this trait within species. This QTL effect is caused by cis-regulatory variation in the floral repressor ChFLC, such that genotypes with low-expressing ChFLC alleles show both early flowering and accelerated age-dependent changes in leaf form, including faster leaflet production. We provide evidence that this mechanism coordinates leaf development with reproductive timing and may help to optimize resource allocation to the next generation.

Natural Variation Identifies ICARUS1, a Universal Gene Required for Cell Proliferation and Growth at High Temperatures in Arabidopsis thaliana

Plants are highly sensitive to environmental changes and even small variations in ambient temperature have severe consequences on their growth and development. Temperature affects multiple aspects of plant development, but the processes and mechanisms underlying thermo-sensitive growth responses are mostly unknown. Here we exploit natural variation in Arabidopsis thaliana to identify and characterize novel components and processes mediating thermo-sensitive growth responses in plants. Phenotypic screening of wild accessions identified several strains displaying pleiotropic growth defects, at cellular and organism levels, specifically at high ambient temperatures. Positional cloning and characterization of the underlying gene revealed that ICARUS1 (ICA1), which encodes a protein of the tRNA^His guanylyl transferase (Thg1) superfamily, is required for plant growth at high temperatures. Transcriptome and gene marker analyses together with DNA content measurements show that ICA1 loss-of-function results in down regulation of cell cycle associated genes at high temperatures, which is linked with a block in G2/M transition and endoreduplication. In addition, plants with mutations in ICA1 show enhanced sensitivity to DNA damage. Characterization of additional strains that carry lesions in ICA1, but display normal growth, shows that alternative splicing is likely to alleviate the deleterious effects of some natural mutations. Furthermore, analyses of worldwide and regional collections of natural accessions indicate that ICA1 loss-of-function has arisen several times independently, and that these occur at high frequency in some local populations. Overall our results suggest that ICA1-mediated-modulation of fundamental processes such as tRNA^His maturation, modify plant growth responses to temperature changes in a quantitative and reversible manner, in natural populations.

The increase in average temperatures across the globe has been predicted to have negative impacts on agricultural productivity. Therefore, there is a need to understand the molecular mechanisms that underlie plant growth responses to varying temperature regimes. At present, very little is known about the genes and pathways that modulate thermo-sensory growth responses in plants. In this article, the authors exploit natural variation in the commonly occurring weed thale cress (Arabidopsis thaliana) and identify a gene referred to as ICARUS1 to be required for plant growth at higher ambient temperatures. Plants carrying lesions in this gene stop growing at high temperatures and revert to growth when temperatures reduce. Using a combination of computational, molecular and cell biological approaches, the authors demonstrate that allelic variation at ICARUS1, which encodes an enzyme required for the fundamental biochemical process of tRNA^His maturation, underlies variation in thermo-sensory growth responses of A. thaliana. Furthermore, the authors discover that the deleterious impact of a natural mutation in ICARUS1 is suppressed through alternative splicing, thus suggesting the potential for alternative splicing to buffer the impacts of some natural mutations. These results support that modulation of fundamental processes, in addition to transcriptional regulation, mediate thermo-sensory growth responses in plants.

Natural variation in epigenetic pathways affects the specification of female gamete precursors in Arabidopsis

In angiosperms, the transition to the female gametophytic phase relies on the specification of premeiotic gamete precursors from sporophytic cells in the ovule. In Arabidopsis thaliana, a single diploid cell is specified as the premeiotic female gamete precursor. Here, we show that ecotypes of Arabidopsis exhibit differences in megasporogenesis leading to phenotypes reminiscent of defects in dominant mutations that epigenetically affect the specification of female gamete precursors. Intraspecific hybridization and polyploidy exacerbate these defects, which segregate quantitatively in F2 populations derived from ecotypic hybrids, suggesting that multiple loci control cell specification at the onset of female meiosis. This variation in cell differentiation is influenced by the activity of ARGONAUTE9 (AGO9) and RNA-DEPENDENT RNA POLYMERASE6 (RDR6), two genes involved in epigenetic silencing that control the specification of female gamete precursors. The pattern of transcriptional regulation and localization of AGO9 varies among ecotypes, and abnormal gamete precursors in ovules defective for RDR6 share identity with ectopic gamete precursors found in selected ecotypes. Our results indicate that differences in the epigenetic control of cell specification lead to natural phenotypic variation during megasporogenesis. We propose that this mechanism could be implicated in the emergence and evolution of the reproductive alternatives that prevail in flowering plants.

Post-transcriptional and post-translational regulations of drought and heat response in plants: a spider's web of mechanisms

Drought and heat tolerance are complex quantitative traits. Moreover, the adaptive significance of some stress-related traits is more related to plant survival than to agronomic performance. A web of regulatory mechanisms fine-tunes the expression of stress-related traits and integrates both environmental and developmental signals. Both post-transcriptional and post-translational modifications contribute substantially to this network with a pivotal regulatory function of the transcriptional changes related to cellular and plant stress response. Alternative splicing and RNA-mediated silencing control the amount of specific transcripts, while ubiquitin and SUMO modify activity, sub-cellular localization and half-life of proteins. Interactions across these modification mechanisms ensure temporally and spatially appropriate patterns of downstream-gene expression. For key molecular components of these regulatory mechanisms, natural genetic diversity exists among genotypes with different behavior in terms of stress tolerance, with effects upon the expression of adaptive morphological and/or physiological target traits.

Rapid divergence and high diversity of miRNAs and miRNA targets in the Camelineae

MicroRNAs (miRNAs) are short RNAs involved in gene regulation through translational inhibition and transcript cleavage. After processing from imperfect fold-back structures, miRNAs are incorporated into RNA-induced silencing complexes (RISCs) before targeting transcripts with varying degrees of complementarity. Some miRNAs are evolutionarily deep-rooted, and sequence complementarity with their targets is maintained through purifying selection. Both Arabidopsis and Capsella belong to the tribe Camelineae in the Brassicaceae, with Capsella rubella serving as an outgroup to the genus Arabidopsis. The genome sequence of C. rubella has recently been released, which allows characterization of its miRNA complement in comparison with Arabidopsis thaliana and Arabidopsis lyrata. Through next-generation sequencing, we identify high-confidence miRNA candidates specific to the C. rubella lineage. Only a few lineage-specific miRNAs have been studied for evolutionary constraints, and there have been no systematic studies of miRNA target diversity within or divergence between closely related plant species. Therefore we contrast sequence variation in miRNAs and their targets within A. thaliana, and between A. thaliana, A. lyrata and C. rubella. We document a surprising amount of small-scale variation in miRNA-target pairs, where many miRNAs are predicted to have species-specific targets in addition to ones that are shared between species. Our results emphasize that the transitive nature of many miRNA-target pairs can be observed even on a relatively short evolutionary time-scale, with non-random occurrences of differences in miRNAs and their complements in the miRNA precursors, the miRNA* sequences.

Natural loss-of-function mutation of EDR1 conferring resistance to tomato powdery mildew in Arabidopsis thaliana accession C24

To screen for potentially novel types of resistance to tomato powdery mildew Oidium neolycopersici, a disease assay was performed on 123 Arabidopsis thaliana accessions. Forty accessions were fully resistant, and one, C24, was analysed in detail. By quantitative trait locus (QTL) analysis of an F2 population derived from C24 × Sha (susceptible accession), two QTLs associated with resistance were identified in C24. Fine mapping of QTL-1 on chromosome 1 delimited the region to an interval of 58 kb encompassing 15 candidate genes. One of these was Enhanced Disease Resistance 1 (EDR1). Evaluation of the previously obtained edr1 mutant of Arabidopsis accession Col-0, which was identified because of its resistance to powdery mildew Golovinomyces cichoracearum, showed that it also displayed resistance to O. neolycopersici. Sequencing of EDR1 in our C24 germplasm (referred to as C24-W) revealed two missing nucleotides in the second exon of EDR1 resulting in a premature stop codon. Remarkably, C24 obtained from other laboratories does not contain the EDR1 mutation. To verify the identity of C24-W, a DNA region containing a single nucleotide polymorphism (SNP) unique to C24 was sequenced showing that C24-W contains the C24-specific nucleotide. C24-W showed enhanced resistance to O. neolycopersici compared with C24 not containing the edr1 mutation. Furthermore, C24-W displayed a dwarf phenotype, which was not associated with the mutation in EDR1 and was not caused by the differential accumulation of pathogenesis-related genes. In conclusion, we identified a natural edr1 mutant in the background of C24.

The miR9863 family regulates distinct Mla alleles in barley to attenuate NLR receptor-triggered disease resistance and cell-death signaling

Barley (Hordeum vulgare L.) Mla alleles encode coiled-coil (CC), nucleotide binding, leucine-rich repeat (NB-LRR) receptors that trigger isolate-specific immune responses against the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). How Mla or NB-LRR genes in grass species are regulated at post-transcriptional level is not clear. The microRNA family, miR9863, comprises four members that differentially regulate distinct Mla alleles in barley. We show that miR9863 members guide the cleavage of Mla1 transcripts in barley, and block or reduce the accumulation of MLA1 protein in the heterologous Nicotiana benthamiana expression system. Regulation specificity is determined by variation in a unique single-nucleotide-polymorphism (SNP) in mature miR9863 family members and two SNPs in the Mla miR9863-binding site that separates these alleles into three groups. Further, we demonstrate that 22-nt miR9863s trigger the biogenesis of 21-nt phased siRNAs (phasiRNAs) and together these sRNAs form a feed-forward regulation network for repressing the expression of group I Mla alleles. Overexpression of miR9863 members specifically attenuates MLA1, but not MLA10-triggered disease resistance and cell-death signaling. We propose a key role of the miR9863 family in dampening immune response signaling triggered by a group of MLA immune receptors in barley.

Plants rely on cell-surface and intracellular immune receptors to sense pathogen invasion and to mediate defense responses. However, uncontrolled activation of immune responses is harmful to plant growth and development. Small RNAs have recently been shown to fine-tune the expression of intracellular immune receptors and contribute to the regulation of defense signaling in dicot plants, while similar processes have not been well documented in monocot grain crops, such as barley and wheat. Here, we show that, in barley, some members of the miR9863 family target a subset of Mla alleles that confer race-specific disease resistance to the powdery mildew fungus. These miRNAs act on Mla transcripts by cleavage and translational repression. Production of a type of trans-acting small RNAs, designated as phasiRNAs, enhances the effects of miRNA regulation on Mla targets. We propose that Mla-mediated immune signaling is fine-tuned by the miRNAs at later stage of MLA activation to avoid overloading of immune responses in barley cells.

The small RNA diversity from Medicago truncatula roots under biotic interactions evidences the environmental plasticity of the miRNAome

BACKGROUND

Legume roots show a remarkable plasticity to adapt their architecture to biotic and abiotic constraints, including symbiotic interactions. However, global analysis of miRNA regulation in roots is limited, and a global view of the evolution of miRNA-mediated diversification in different ecotypes is lacking.

RESULTS

In the model legume Medicago truncatula, we analyze the small RNA transcriptome of roots submitted to symbiotic and pathogenic interactions. Genome mapping and a computational pipeline identify 416 miRNA candidates, including known and novel variants of 78 miRNA families present in miRBase. Stringent criteria of pre-miRNA prediction yield 52 new mtr-miRNAs, including 27 miRtrons. Analyzing miRNA precursor polymorphisms in 26 M. truncatula ecotypes identifies higher sequence polymorphism in conserved rather than Medicago-specific miRNA precursors. An average of 19 targets, mainly involved in environmental responses and signalling, is predicted per novel miRNA. We identify miRNAs responsive to bacterial and fungal pathogens or symbionts as well as their related Nod and Myc-LCO symbiotic signals. Network analyses reveal modules of new and conserved co-expressed miRNAs that regulate distinct sets of targets, highlighting potential miRNA-regulated biological pathways relevant to pathogenic and symbiotic interactions.

CONCLUSIONS

We identify 52 novel genuine miRNAs and large plasticity of the root miRNAome in response to the environment, and also in response to purified Myc/Nod signaling molecules. The new miRNAs identified and their sequence variation across M. truncatula ecotypes may be crucial to understand the adaptation of root growth to the soil environment, notably in the agriculturally important legume crops.

Developmental patterning by gradients of mobile small RNAs

Development of multicellular organisms depends on intercellular communication via mobile signals that provide positional information to coordinate cell fate decisions. In addition to peptide ligands, transcription factors, and hormones, plants use small RNAs as positional instructive signals. The unique patterning properties of small RNA gradients resulting from regulated mobility suggest conceptual similarities to the function of animal morphogens, and provide robustness and precision to the formation of cell fate boundaries. While common principles may underlie the formation, stability, and interpretation of both plant small RNA and animal morphogen gradients, the unique nature of small RNAs with respect to their biogenesis and target regulation imply important differences as well. In this review, we discuss the patterning properties of mobile small RNAs and highlight recent studies that have advanced our understanding of how small RNAs move, and how the graded accumulation that underlies their patterning activity could be created, maintained, and interpreted.

Genetic architecture of naturally occurring quantitative traits in plants: an updated synthesis

Deciphering the genetic and molecular bases of quantitative variation is a long-standing challenge in plant biology because it is essential for understanding evolution and for accelerating plant breeding. Recent multi-trait analyses at different phenotypic levels are uncovering the pleiotropy and the genetic regulation underlying high-level complex traits. Thus, the number of known causal loci, genes and nucleotide polymorphisms is expanding. Current plant causal catalogs contain ∼400 genes and natural polymorphisms revealing several dysfunctional allelic series that involve multiple mutations. In addition, repeated evolution of quantitative traits mediated by large effect alleles is found across plant phylogeny. Finally, systematic analyses of genetic and environmental interactions are beginning to elucidate the molecular mechanisms of relevant interactions.

Determinants beyond both complementarity and cleavage govern microR159 efficacy in Arabidopsis

Plant microRNAs (miRNAs) are critical regulators of gene expression, however little attention has been given to the principles governing miRNA silencing efficacy. Here, we utilize the highly conserved Arabidopsis miR159-MYB33/MYB65 regulatory module to explore these principles. Firstly, we show that perfect central complementarity is not required for strong silencing. Artificial miR159 variants with two cleavage site mismatches can potently silence MYB33/MYB65, fully complementing a loss-of-function mir159 mutant. Moreover, these miR159 variants can cleave MYB33/MYB65 mRNA, however cleavage appears attenuated, as the ratio of cleavage products to full length transcripts decreases with increasing central mismatches. Nevertheless, high levels of un-cleaved MYB33/MYB65 transcripts are strongly silenced by a non-cleavage mechanism. Contrary to MIR159a variants that strongly silenced endogenous MYB33/MYB65, artificial MYB33 variants with central mismatches to miR159 are not efficiently silenced. We demonstrate that differences in the miRNA:target mRNA stoichiometry underlie this paradox. Increasing miR159 abundance in the MYB33 variants results in a strong silencing outcome, whereas increasing MYB33 transcript levels in the MIR159a variants results in a poor silencing outcome. Finally, we identify highly conserved nucleotides that flank the miR159 binding site in MYB33, and demonstrate that they are critical for efficient silencing, as mutation of these flanking nucleotides attenuates silencing at a level similar to that of central mismatches. This implies that the context in which the miRNA binding site resides is a key determinant in controlling the degree of silencing and that a miRNA "target site" encompasses sequences that extend beyond the miRNA binding site. In conclusion, our findings dismiss the notion that miRNA:target complementarity, underpinned by central matches, is the sole dictator of the silencing outcome.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

BACKGROUND

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

RESULTS

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

CONCLUSIONS

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

Can we always sweep the details of RNA-processing under the carpet?

RNA molecules follow a succession of enzyme-mediated processing steps from transcription to maturation. The participating enzymes, for example the spliceosome for mRNAs and Drosha and Dicer for microRNAs, are also produced in the cell and their copy-numbers fluctuate over time. Enzyme copy-number changes affect the processing rate of the substrate molecules; high enzyme numbers increase the processing rate, while low enzyme numbers decrease it. We study different RNA-processing cascades where enzyme copy-numbers are either fixed or fluctuate. We find that for the fixed enzyme copy-numbers, the substrates at steady-state are Poisson-distributed, and the whole RNA cascade dynamics can be understood as a single birth-death process of the mature RNA product. In this case, solely fluctuations in the timing of RNA processing lead to variation in the number of RNA molecules. However, we show analytically and numerically that when enzyme copy-numbers fluctuate, the strength of RNA fluctuations increases linearly with the RNA transcription rate. This linear effect becomes stronger as the speed of enzyme dynamics decreases relative to the speed of RNA dynamics. Interestingly, we find that under certain conditions, the RNA cascade can reduce the strength of fluctuations in the expression level of the mature RNA product. Finally, by investigating the effects of processing polymorphisms, we show that it is possible for the effects of transcriptional polymorphisms to be enhanced, reduced or even reversed. Our results provide a framework to understand the dynamics of RNA processing.

Morphogenesis of simple leaves: regulation of leaf size and shape

Plants produce new organs throughout their life span. Leaves first initiate as rod-like structures protruding from the shoot apical meristem, while they need to pass through different developmental stages to become the flat organ specialized in photosynthesis. Leaf morphogenesis is an active process regulated by many genes and pathways that can generate organs with a wide variety of sizes and shapes. Important differences in leaf architecture can be seen among different species, but also in single individuals. A key aspect of leaf morphogenesis is the precise control of cell proliferation. Modification or manipulation of this process may lead to leaves with different sizes and shapes, and changes in the organ margins and curvature. Many genes required for leaf development have been identified in Arabidopsis thaliana, and the mechanisms underlying leaf morphogenesis are starting to be unraveled at the molecular level.

RNA sequencing reveals small RNAs differentially expressed between incipient Japanese threespine sticklebacks

Background

Non-coding small RNAs, ranging from 20 to 30 nucleotides in length, mediate the regulation of gene expression and play important roles in many biological processes. One class of small RNAs, microRNAs (miRNAs), are highly conserved across taxa and mediate the regulation of the chromatin state and the post-transcriptional regulation of messenger RNA (mRNA). Another class of small RNAs is the Piwi-interacting RNAs, which play important roles in the silencing of transposons and other functional genes. Although the biological functions of the different small RNAs have been elucidated in several laboratory animals, little is known regarding naturally occurring variation in small RNA transcriptomes among closely related species.

Results

We employed next-generation sequencing technology to compare the expression profiles of brain small RNAs between sympatric species of the Japanese threespine stickleback (Gasterosteus aculeatus). We identified several small RNAs that were differentially expressed between sympatric Pacific Ocean and Japan Sea sticklebacks. Potential targets of several small RNAs were identified as repetitive sequences. Female-biased miRNA expression from the old X chromosome was also observed, and it was attributed to the degeneration of the Y chromosome.

Conclusions

Our results suggest that expression patterns of small RNA can differ between incipient species and may be a potential mechanism underlying differential mRNA expression and transposon activity.

Classification and comparison of small RNAs from plants

Regulatory small RNAs, which range in size from 20 to 24 nucleotides, are ubiquitous components of endogenous plant transcriptomes, as well as common responses to exogenous viral infections and introduced double-stranded RNA (dsRNA). Endogenous small RNAs derive from the processing of helical RNA precursors and can be categorized into several groups based on differences in biogenesis and function. A major distinction can be observed between small RNAs derived from single-stranded precursors with a hairpin structure [referred to here as hairpin RNAs (hpRNAs)] and those derived from dsRNA precursors [small interfering RNAs (siRNAs)]. hpRNAs in plants can be divided into two secondary groups: microRNAs and those that are not microRNAs. The currently known siRNAs fall mostly into one of three secondary groups: heterochromatic siRNAs, secondary siRNAs, and natural antisense transcript siRNAs. Tertiary subdivisions can be identified within many of the secondary classifications as well. Comparisons between the different classes of plant small RNAs help to illuminate key goals for future research.

STA1, an Arabidopsis pre-mRNA processing factor 6 homolog, is a new player involved in miRNA biogenesis

MicroRNAs (miRNAs) are small regulatory RNAs that have important regulatory roles in numerous developmental and metabolic processes in most eukaryotes. In Arabidopsis, DICER-LIKE1 (DCL1), HYPONASTIC LEAVES 1, SERRATE, HUA ENHANCER1 and HASTY are involved in processing of primary miRNAs (pri-miRNAs) to yield precursor miRNAs (pre-miRNAs) and eventually miRNAs. In addition to these components, mRNA cap-binding proteins, CBP80/ABA HYPERSENSITIVE1 and CBP20, also participate in miRNA biogenesis. Here, we show that STABILIZED1 (STA1), an Arabidopsis pre-mRNA processing factor 6 homolog, is also involved in the biogenesis of miRNAs. Similar to other miRNA biogenesis-defective mutants, sta1-1 accumulated significantly lower levels of mature miRNAs and concurrently higher levels of pri-miRNAs than wild type. The dramatic reductions of mature miRNAs were associated with the accumulation of their target gene transcripts and developmental defects. Furthermore, sta1-1 impaired splicing of intron containing pri-miRNAs and decreased transcript levels of DCL1. These results suggest that STA1 is involved in miRNA biogenesis directly by functioning in pri-miRNA splicing and indirectly by modulating the DCL1 transcript level.

Regulation of apomixis: learning from sexual experience

Apomixis is a natural form of asexual reproduction through seeds that leads to viable offspring genetically identical to the mother plant. New evidence from sexual model species indicates that the regulation of female gametogenesis and seed formation is also directed by epigenetic mechanisms that are crucial to control events that distinguish sexuality from apomixis, with important implications for our understanding of the evolutionary forces that shape structural variation and diversity in plant reproduction.

Transgenerational defense induction and epigenetic inheritance in plants

Rapidly accumulating evidence shows that herbivore and pathogen attack of plants can generate particular defense phenotypes across generations. What was once thought to be an oddity of plant defense induction now appears to be a taxonomically widespread phenomenon with strong potential to impact the ecology and evolution of species interactions. DNA methylation, histone modifications, and small RNAs each contribute to transgenerational defense initiation; examples in several species demonstrate that this induction can last for multiple generations. Priming of the offspring generation for more rapid induction following subsequent attack has also been reported. The extent to which transgenerational induction is predictable, detectable in nature, and subject to manipulation will determine the ability of researchers to decipher its role in plant-herbivore and plant-pathogen interactions.