Identification of MIR390a precursor processing-defective mutants in Arabidopsis by direct genome sequencing

Principles of miRNA/miRNA* function in plant MIRNA processing

MicroRNAs (miRNAs) are essential regulators of gene expression, defined by their unique biogenesis, which requires the precise excision of the small RNA from an imperfect fold-back precursor. Unlike their animal counterparts, plant miRNA precursors exhibit variations in sizes and shapes. Plant MIRNAs can undergo processing in a base-to-loop or loop-to-base direction, with DICER-LIKE1 (DCL1) releasing the miRNA after two cuts (two-step MIRNAs) or more (sequential MIRNAs). In this study, we demonstrate the critical role of the miRNA/miRNA* duplex region in the processing of miRNA precursors. We observed that endogenous MIRNAs frequently experience suboptimal processing in vivo due to mismatches in the miRNA/miRNA* duplex, a key region that fine-tunes miRNA levels. Enhancing the interaction energy of the miRNA/miRNA* duplex in two-step MIRNAs results in a substantial increase in miRNA levels. Conversely, sequential MIRNAs display distinct and specific requirements for the miRNA/miRNA* duplexes along their foldback structure. Our work establishes a connection between the miRNA/miRNA* structure and precursor processing mechanisms. Furthermore, we reveal a link between the biological function of miRNAs and the processing mechanism of their precursors with the evolution of plant miRNA/miRNA* duplex structures.

Evidence for an RNAi-independent role of Arabidopsis DICER-LIKE2 in growth inhibition and basal antiviral resistance

Flowering plant genomes encode four or five DICER-LIKE (DCL) enzymes that produce small interfering RNAs (siRNAs) and microRNAs, which function in RNA interference (RNAi). Different RNAi pathways in plants effect transposon silencing, antiviral defense, and endogenous gene regulation. DCL2 acts genetically redundantly with DCL4 to confer basal antiviral defense. However, DCL2 may also counteract DCL4 since knockout of DCL4 causes growth defects that are suppressed by DCL2 inactivation. Current models maintain that RNAi via DCL2-dependent siRNAs is the biochemical basis of both effects. Here, we report that DCL2-mediated antiviral resistance and growth defects cannot be explained by the silencing effects of DCL2-dependent siRNAs. Both functions are defective in genetic backgrounds that maintain high levels of DCL2-dependent siRNAs, either with specific point mutations in DCL2 or with reduced DCL2 dosage because of heterozygosity for dcl2 knockout alleles. Intriguingly, all DCL2 functions require its catalytic activity, and the penetrance of DCL2-dependent growth phenotypes in dcl4 mutants correlates with DCL2 protein levels but not with levels of major DCL2-dependent siRNAs. We discuss this requirement and correlation with catalytic activity but not with resulting siRNAs, in light of other findings that reveal a DCL2 function in innate immunity activation triggered by cytoplasmic double-stranded RNA.

Transgene-free, virus-based gene silencing in plants by artificial microRNAs derived from minimal precursors

Artificial microRNAs (amiRNAs) are highly specific, 21-nucleotide (nt) small RNAs designed to silence target transcripts. In plants, their application as biotechnological tools for functional genomics or crop improvement is limited by the need of transgenically expressing long primary miRNA (pri-miRNA) precursors to produce the amiRNAs in vivo. Here, we analyzed the minimal structural and sequence requirements for producing effective amiRNAs from the widely used, 521-nt long AtMIR390a pri-miRNA from Arabidopsis thaliana. We functionally screened in Nicotiana benthamiana a large collection of constructs transiently expressing amiRNAs against endogenous genes and from artificially shortened MIR390-based precursors and concluded that highly effective and accurately processed amiRNAs can be produced from a chimeric precursor of only 89 nt. This minimal precursor was further validated in A. thaliana transgenic plants expressing amiRNAs against endogenous genes. Remarkably, minimal but not full-length precursors produce authentic amiRNAs and induce widespread gene silencing in N. benthamiana when expressed from an RNA virus, which can be applied into leaves by spraying infectious crude extracts. Our results reveal that the length of amiRNA precursors can be shortened without affecting silencing efficacy, and that viral vectors including minimal amiRNA precursors can be applied in a transgene-free manner to induce whole-plant gene silencing.

Argonaute7 (AGO7) optimizes arbuscular mycorrhizal fungal associations and enhances competitive growth in Nicotiana attenuata

Plants interact with arbuscular mycorrhizal fungi (AMF) and in doing so, change transcript levels of many miRNAs and their targets. However, the identity of an Argonaute (AGO) that modulates this interaction remains unknown, including in Nicotiana attenuata. We examined how the silencing of NaAGO1/2/4/7/and 10 by RNAi influenced plant-competitive ability under low-P conditions when they interact with AMF. Furthermore, the roles of seven miRNAs, predicted to regulate signaling and phosphate homeostasis, were evaluated by transient overexpression. Only NaAGO7 silencing by RNAi (irAGO7) significantly reduced the competitive ability under P-limited conditions, without changes in leaf or root development, or juvenile-to-adult phase transitions. In plants growing competitively in the glasshouse, irAGO7 roots were over-colonized with AMF, but they accumulated significantly less phosphate and the expression of their AMF-specific transporters was deregulated. Furthermore, the AMF-induced miRNA levels were inversely regulated with the abundance of their target transcripts. miRNA overexpression consistently decreased plant fitness, with four of seven-tested miRNAs reducing mycorrhization rates, and two increasing mycorrhization rates. Overexpression of Na-miR473 and Na-miRNA-PN59 downregulated targets in GA, ethylene, and fatty acid metabolism pathways. We infer that AGO7 optimizes competitive ability and colonization by regulating miRNA levels and signaling pathways during a plant's interaction with AMF.

AGO1 and HSP90 buffer different genetic variants in Arabidopsis thaliana

Argonaute 1 (AGO1), the principal protein component of microRNA-mediated regulation, plays a key role in plant growth and development. AGO1 physically interacts with the chaperone HSP90, which buffers cryptic genetic variation in plants and animals. We sought to determine whether genetic perturbation of AGO1 in Arabidopsis thaliana would also reveal cryptic genetic variation, and if so, whether AGO1-dependent loci overlap with those dependent on HSP90. To address these questions, we introgressed a hypomorphic mutant allele of AGO1 into a set of mapping lines derived from the commonly used Arabidopsis strains Col-0 and Ler. Although we identified several cases in which AGO1 buffered genetic variation, none of the AGO1-dependent loci overlapped with those buffered by HSP90 for the traits assayed. We focused on 1 buffered locus where AGO1 perturbation uncoupled the traits days to flowering and rosette leaf number, which are otherwise closely correlated. Using a bulk segregant approach, we identified a nonfunctional Ler hua2 mutant allele as the causal AGO1-buffered polymorphism. Introduction of a nonfunctional hua2 allele into a Col-0 ago1 mutant background recapitulated the Ler-dependent ago1 phenotype, implying that coupling of these traits involves different molecular players in these closely related strains. Taken together, our findings demonstrate that even though AGO1 and HSP90 buffer genetic variation in the same traits, these robustness regulators interact epistatically with different genetic loci, suggesting that higher-order epistasis is uncommon. Plain Language Summary Argonaute 1 (AGO1), a key player in plant development, interacts with the chaperone HSP90, which buffers environmental and genetic variation. We found that AGO1 buffers environmental and genetic variation in the same traits; however, AGO1-dependent and HSP90-dependent loci do not overlap. Detailed analysis of a buffered locus found that a nonfunctional HUA2 allele decouples days to flowering and rosette leaf number in an AGO1-dependent manner, suggesting that the AGO1-dependent buffering acts at the network level.

A conserved SNP variation in the pre-miR396c flanking region in Oryza sativa indica landraces correlates with mature miRNA abundance

Plant precursor miRNAs (pre-miRNA) have conserved evolutionary footprints that correlate with mode of miRNA biogenesis. In plants, base to loop and loop to base modes of biogenesis have been reported. Conserved structural element(s) in pre-miRNA play a major role in turn over and abundance of mature miRNA. Pre-miR396c sequences and secondary structural characteristics across Oryza species are presented. Based on secondary structure, twelve Oryza pre-miR396c sequences are divided into three groups, with the precursor from halophytic Oryza coarctata forming a distinct group. The miRNA-miRNA* duplex region is completely conserved across eleven Oryza species as are other structural elements in the pre-miRNA, suggestive of an evolutionarily conserved base-to-loop mode of miRNA biogenesis. SNPs within O. coarctata mature miR396c sequence and miRNA* region have the potential to alter target specificity and association with the RNA-induced silencing complex. A conserved SNP variation, rs10234287911 (G/A), identified in O. sativa pre-miR396c sequences alters base pairing above the miRNA-miRNA* duplex. The more stable structure conferred by the 'A_10234287911' allele may promote better processing vis-à-vis the structure conferred by 'G_10234287911' allele. We also examine pri- and pre-miR396c expression in cultivated rice under heat and salinity and their correlation with miR396c expression.

MIR390 Is Involved in Regulating Anthracnose Resistance in Apple

As an important cash crop in China, apple has a good flavor and is rich in nutrients. Fungal attacks have become a major obstacle in apple cultivation. Colletotrichum gloeosporioides is one of the most devastating fungal pathogens in apple. Thus, discovering resistance genes in response to C. gloeosporioides may aid in designing safer control strategies and facilitate the development of apple resistance breeding. A previous study reported that 'Hanfu' autotetraploid apple displayed higher C. gloeosporioides resistance than 'Hanfu' apple, and the expression level of mdm-MIR390b was significantly upregulated in autotetraploid plants compared to that in 'Hanfu' plants, as demonstrated by digital gene expression (DGE) analysis. It is still unclear, however, whether mdm-MIR390b regulates apple anthracnose resistance. Apple MIR390b was transformed into apple 'GL-3' plants to identify the functions of mdm-MIR390b in anthracnose resistance. C. gloeosporioides treatment analysis indicated that the overexpression of mdm-MIR390b reduced fungal damage to apple leaves and fruit. Physiology analysis showed that mdm-MIR390b increased C. gloeosporioides resistance by improving superoxide dismutase (SOD) and peroxidase (POD) activity to alleviate the damage caused by O₂^- and H₂O₂. Our results demonstrate that mdm-MIR390b can improve apple plants' anthracnose resistance.

Biogenesis, conservation, and function of miRNA in liverworts

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules, 18-24 nucleotides long, that control multiple gene regulatory pathways via post-transcriptional gene silencing in eukaryotes. To develop a comprehensive picture of the evolutionary history of miRNA biogenesis and action in land plants, studies on bryophyte representatives are needed. Here, we review current understanding of liverwort MIR gene structure, miRNA biogenesis, and function, focusing on the simple thalloid Pellia endiviifolia and the complex thalloid Marchantia polymorpha. We review what is known about conserved and non-conserved miRNAs, their targets, and the functional implications of miRNA action in M. polymorpha and P. endiviifolia. We note that most M. polymorpha miRNAs are encoded within protein-coding genes and provide data for 23 MIR gene structures recognized as independent transcriptional units. We identify M. polymorpha genes involved in miRNA biogenesis that are homologous to those identified in higher plants, including those encoding core microprocessor components and other auxiliary and regulatory proteins that influence the stability, folding, and processing of pri-miRNAs. We analyzed miRNA biogenesis proteins and found similar domain architecture in most cases. Our data support the hypothesis that almost all miRNA biogenesis factors in higher plants are also present in liverworts, suggesting that they emerged early during land plant evolution.

CRISPR/Cas9-mediated fine-tuning of miRNA expression in tetraploid potato

MicroRNAs (miRNAs) are small noncoding RNAs, which modulate the abundance and spatiotemporal accumulation of target mRNAs at transcriptional and post-transcriptional levels and through that play important roles in several biological processes in plants. Here we show that in polyploid species, CRISPR/Cas9 system can be used for fine-tuning of miRNA expression, which can have broader range of applications compared to knock-out mutants. We established the complete pipeline for CRISPR-Cas9-mediated modulation of miRNA expression in potato. It consists of (1) design and assembly of dual sgRNA CRISPR/Cas9 constructs, (2) transient transfection of protoplasts following fast and efficient screening by high resolution melting analysis to select functional sgRNAs, and (3) stable transformation of potato explants with functional sgRNAs and selection of regenerated transgenic lines with desired mutations and desired miRNA abundance based on sequencing and RT-qPCR. We show that miRNA-editing using dual sgRNA approach results in different types of mutations among transgenic lines but also in different alleles of the same plant, which are target site-dependent. The most frequent were short deletions, but we also detected 1-nt insertions (T or G), deletions between two sgRNAs and larger deletions. miRNA abundance correlates with the frequency and type of introduced mutations, as more extensive mutations in more alleles result in lower miRNA abundance. Interestingly, some mutated loci can generate alternative miRNAs, now novel targets were however predicted for those. In all transgenic lines with Cas9 expression, we detected mutations, suggesting high efficiency of Cas9-editing. We confirmed the miRNA-editing efficiency of our optimised approach in two different potato genotypes and three different loci.

MicroRNA candidate miRcand137 in apple is induced by Botryosphaeria dothidea for impairing host defense

MicroRNA (miRNA)-mediated gene silencing is a master gene regulatory pathway in plant-pathogen interactions. The differential accumulation of miRNAs among plant varieties alters the expression of target genes, affecting plant defense responses and causing resistance differences among varieties. Botryosphaeria dothidea is an important phytopathogenic fungus of apple (Malus domestica). Malus hupehensis (Pamp.) Rehder, a wild apple species, is highly resistant, whereas the apple cultivar "Fuji" is highly susceptible. Here, we identified a 22-nt miRNA candidate named miRcand137 that compromises host resistance to B. dothidea infection and whose processing was affected by precursor sequence variation between M. hupehensis and "Fuji." miRcand137 guides the direct cleavage of and produced target-derived secondary siRNA against Ethylene response factor 14 (ERF14), a transcriptional activator of pathogenesis-related homologs that confers disease resistance to apple. We showed that miRcand137 acts as an inhibitor of apple immunity by compromising ERF14-mediated anti-fungal defense and revealed a negative association between miRcand137 expression and B. dothidea sensitivity in both resistant and susceptible apples. Furthermore, MIRCAND137 was transcriptionally activated by the invading fungi but not by the fungal elicitor, implying B. dothidea induced host miRcand137 as an infection strategy. We propose that the inefficient miRcand137 processing in M. hupehensis decreased pathogen-initiated miRcand137 accumulation, leading to higher resistance against B. dothidea.

Phased small RNA-mediated systemic signaling in plants

Systemic acquired resistance (SAR) involves the generation of systemically transported signal that arms distal plant parts against secondary infections. We show that two phased 21-nucleotide (nt) trans-acting small interfering RNA3a RNAs (tasi-RNA) derived from TAS3a and synthesized within 3 hours of pathogen infection are the early mobile signal in SAR. TAS3a undergoes alternate polyadenylation, resulting in the generation of 555- and 367-nt transcripts. The 555-nt transcripts likely serves as the sole precursor for tasi-RNAs D7 and D8, which cleave Auxin response factors (ARF) 2, 3, and 4 to induce SAR. Conversely, increased expression of ARF3 represses SAR. Knockout mutations in TAS3a or RNA silencing components required for tasi-RNA biogenesis compromise SAR without altering levels of known SAR-inducing chemicals. Both tasi-ARFs and the 367-nt transcripts are mobile and transported via plasmodesmata. Together, we show that tasi-ARFs are the early mobile signal in SAR.

Nested miRNA Secondary Structure Is a Unique Determinant of miR159 Efficacy in Arabidopsis

MicroRNAs (miRNAs) are 20- to 24-nucleotide small RNAs, and whenever a pri-miRNA precursor includes another miRNA precursor, and both of these precursors may generate independent non overlapping mature miRNAs, we called them nested miRNAs. However, the functional and regulatory roles of nested miRNA structures in plants are still unknown. In this study, the Arabidopsis nested miR159a structure, which consists of two nested miRNAs, miR159a.1, and miR159a.2, was used as a model to determine miRNA-mediated gene silencing in plants. Complementation analysis of nested miR159a structures revealed that the miR159a structure can differentially complement the mir159ab phenotype, and a duplex nested structure in the tail end region of the pre-miR159a fold back may have a possible dominant function, indicating the importance of the flanking sequence of the stem in the cleavage of the mature miRNA. Furthermore, continuously higher expression of the miR159a.2 duplex in the severe leaf curl phenotype indicates that miR159a.2 is functional in Arabidopsis and suggests that in plants, a miRNA precursor may encode multiple regulatory small RNAs. Taken together, our study demonstrates that the nested miR159a structure regulated by duplex mutations of miR159a has a unique pattern and provides novel insight into silencing efficacy of Arabidopsis miR159a.

A combinatorial strategy to identify various types of QTLs for quantitative traits using extreme phenotype individuals in an F2 population

Theoretical and applied studies demonstrate the difficulty of detecting extremely over-dominant and small-effect genes for quantitative traits via bulked segregant analysis (BSA) in an F₂ population. To address this issue, we proposed an integrated strategy for mapping various types of quantitative trait loci (QTLs) for quantitative traits via a combination of BSA and whole-genome sequencing. In this strategy, the numbers of read counts of marker alleles in two extreme pools were used to predict the numbers of read counts of marker genotypes. These observed and predicted numbers were used to construct a new statistic, G_w, for detecting quantitative trait genes (QTGs), and the method was named dQTG-seq1. This method was significantly better than existing BSA methods. If the goal was to identify extremely over-dominant and small-effect genes, another reserved DNA/RNA sample from each extreme phenotype F₂ plant was sequenced, and the observed numbers of marker alleles and genotypes were used to calculate G_w to detect QTGs; this method was named dQTG-seq2. In simulated and real rice dataset analyses, dQTG-seq2 could identify many more extremely over-dominant and small-effect genes than BSA and QTL mapping methods. dQTG-seq2 may be extended to other heterogeneous mapping populations. The significance threshold of G_w in this study was determined by permutation experiments. In addition, a handbook for the R software dQTG.seq, which is available at https://cran.r-project.org/web/packages/dQTG.seq/index.html, has been provided in the supplemental materials for the users' convenience. This study provides a new strategy for identifying all types of QTLs for quantitative traits in an F₂ population.

miR160 Interacts in vivo With Pinus pinaster AUXIN RESPONSE FACTOR 18 Target Site and Negatively Regulates Its Expression During Conifer Somatic Embryo Development

MicroRNAs (miRNAs) are key regulators of several plant developmental processes including embryogenesis. Most miRNA families are conserved across major groups of plant species, but their regulatory roles have been studied mainly in model species like Arabidopsis and other angiosperms. In gymnosperms, miRNA-dependent regulation has been less studied since functional approaches in these species are often difficult to establish. Given the fundamental roles of auxin signaling in somatic embryogenesis (SE) induction and embryo development, we investigated a previously predicted interaction between miR160 and a putative target encoding AUXIN RESPONSE FACTOR 18 in Pinus pinaster (PpARF18) embryonic tissues. Phylogenetic analysis of AUXIN RESPONSE FACTOR 18 (ARF18) from Pinus pinaster and Picea abies, used here as a model system of conifer embryogenesis, showed their close relatedness to AUXIN RESPONSE FACTOR (ARF) genes known to be targeted by miR160 in other species, including Arabidopsis ARF10 and ARF16. By using a luciferase (LUC) reporter system for miRNA activity in Arabidopsis protoplasts, we have confirmed that P. pinaster miR160 (ppi-miR160) interacts in vivo with PpARF18 target site. When the primary miR160 from P. pinaster was overexpressed in protoplasts under non-limiting levels of ARGONAUTE1, a significant increase of miR160 target cleavage activity was observed. In contrast, co-expression of the primary miRNA and the target mimic MIM160 led to a decrease of miR160 activity. Our results further support that this interaction is functional during consecutive stages of SE in the conifer model P. abies. Expression analyses conducted in five stages of development, from proembryogenic masses (PEMs) to the mature embryo, show that conifer ARF18 is negatively regulated by miR160 toward the fully developed mature embryo when miR160 reached its highest expression level. This study reports the first in vivo validation of a predicted target site of a conifer miRNA supporting the conservation of miR160 interaction with ARF targets in gymnosperms. The approach used here should be useful for future characterization of miRNA functions in conifer embryogenesis.

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA-target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA-target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA-target interactions.

Systematic Characterization of MicroRNA Processing Modes in Plants With Parallel Amplification of RNA Ends

In plants, the RNase III-type enzyme Dicer-like 1 (DCL1) processes most microRNAs (miRNAs) from their primary transcripts called pri-miRNAs. Four distinct processing modes (i.e., short base to loop, sequential base to loop, short loop to base, and sequential loop to base) have been characterized in Arabidopsis, mainly by the Specific Parallel Amplification of RNA Ends (SPARE) approach. However, SPARE is a targeted cloning method which requires optimization of cloning efficiency and specificity for each target. PARE (Parallel Amplification of RNA Ends) is an untargeted method per se and is widely used to identify miRNA mediated target slicing events. A major concern with PARE in characterizing miRNA processing modes is the potential contamination of mature miRNAs. Here, we provide a method to estimate miRNA contamination levels and showed that most publicly available PARE libraries have negligible miRNA contamination. Both the numbers and processing modes detected by PARE were similar to those identified by SPARE in Arabidopsis. PARE also determined the processing modes of 36 Arabidopsis miRNAs that were unexplored by SPARE, suggesting that it can complement the SPARE approach. Using publicly available PARE datasets, we identified the processing modes of 36, 91, 90, and 54 miRNAs in maize, rice, soybean, and tomato, respectively, and demonstrated that the processing mode was conserved overall within each miRNA family. Through its power of tracking miRNA processing remnants, PARE also facilitated miRNA characterization and annotation.

E183K Mutation in Chalcone Synthase C2 Causes Protein Aggregation and Maize Colorless

Flavonoids give plants their rich colors and play roles in a number of physiological processes. In this study, we identified a novel colorless maize mutant showing reduced pigmentation throughout the whole life cycle by EMS mutagenesis. E183K mutation in maize chalcone synthase C2 (ZmC2) was mapped using MutMap strategy as the causal for colorless, which was further validated by transformation in Arabidopsis. We evaluated transcriptomic and metabolic changes in maize first sheaths caused by the mutation. The downstream biosynthesis was blocked while very few genes changed their expression pattern. ZmC2-E183 site is highly conserved in chalcone synthase among Plantae kingdom and within species' different varieties. Through prokaryotic expression, transient expression in maize leaf protoplasts and stable expression in Arabidopsis, we observed that E183K and other mutations on E183 could cause almost complete protein aggregation of chalcone synthase. Our findings will benefit the characterization of flavonoid biosynthesis and contribute to the body of knowledge on protein aggregation in plants.

Anno genominis XX: 20 years of Arabidopsis genomics

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

Recent advances in the regulation of plant miRNA biogenesis

MicroRNAs (miRNAs) are essential non-coding riboregulators of gene expression in plants and animals. In plants, miRNAs guide their effector protein named ARGONAUTE (AGO) to find target RNAs for gene silencing through target RNA cleavage or translational inhibition. miRNAs are derived from primary miRNA transcripts (pri-miRNAs), most of which are transcribed by the DNA-dependent RNA polymerase II. In plants, an RNase III enzyme DICER-LIKE1-containing complex processes pri-miRNAs in the nucleus into miRNAs. To ensure proper function of miRNAs, plants use multiple mechanisms to control miRNA accumulation. On one hand, pri-miRNA levels are controlled through transcription and stability. On the other hand, the activities of the DCL1 complex are regulated by many protein factors at transcriptional, post-transcriptional and post-translational levels. Notably, recent studies reveal that pri-miRNA structure/sequence features and modifications also play important roles in miRNA biogenesis. In this review, we summarize recent progresses on the mechanisms regulating miRNA biogenesis.

Effective Mapping by Sequencing to Isolate Causal Mutations in the Tomato Genome

Forward genetic analysis remains as one of the most powerful tools for assessing gene functions, although the identification of the causal mutation responsible for a given phenotype has been a tedious and time-consuming task until recently. Advances in deep sequencing technologies have provided new approaches for the exploitation of natural and artificially induced genetic diversity, thus accelerating the discovery of novel allelic variants. In this chapter, a mapping-by-sequencing forward genetics approach is described to identify causal mutations in tomato (Solanum lycopersicum L.), a major crop species that is also a model species for plant biology and breeding.

Identification of key sequence features required for microRNA biogenesis in plants

MicroRNAs (miRNAs) are endogenous small RNAs of ∼21 nt that regulate multiple biological pathways in multicellular organisms. They derive from longer transcripts that harbor an imperfect stem-loop structure. In plants, the ribonuclease type III DICER-LIKE1 assisted by accessory proteins cleaves the precursor to release the mature miRNA. Numerous studies highlight the role of the precursor secondary structure during plant miRNA biogenesis; however, little is known about the relevance of the precursor sequence. Here, we analyzed the sequence composition of plant miRNA primary transcripts and found specifically located sequence biases. We show that changes in the identity of specific nucleotides can increase or abolish miRNA biogenesis. Most conspicuously, our analysis revealed that the identity of the nucleotides at unpaired positions of the precursor plays a crucial role during miRNA biogenesis in Arabidopsis.

The secondary structure of miRNA precursor sequences is known to affect processing by DICER-like proteins. Here Rojas et al. show that additional sequence features also play a regulatory role in plants with nucleotide identity at unpaired positions substantially impacting processing efficiency.

Next Generation Sequencing Based Forward Genetic Approaches for Identification and Mapping of Causal Mutations in Crop Plants: A Comprehensive Review

The recent advancements in forward genetics have expanded the applications of mutation techniques in advanced genetics and genomics, ahead of direct use in breeding programs. The advent of next-generation sequencing (NGS) has enabled easy identification and mapping of causal mutations within a short period and at relatively low cost. Identifying the genetic mutations and genes that underlie phenotypic changes is essential for understanding a wide variety of biological functions. To accelerate the mutation mapping for crop improvement, several high-throughput and novel NGS based forward genetic approaches have been developed and applied in various crops. These techniques are highly efficient in crop plants, as it is relatively easy to grow and screen thousands of individuals. These approaches have improved the resolution in quantitative trait loci (QTL) position/point mutations and assisted in determining the functional causative variations in genes. To be successful in the interpretation of NGS data, bioinformatics computational methods are critical elements in delivering accurate assembly, alignment, and variant detection. Numerous bioinformatics tools/pipelines have been developed for such analysis. This article intends to review the recent advances in NGS based forward genetic approaches to identify and map the causal mutations in the crop genomes. The article also highlights the available bioinformatics tools/pipelines for reducing the complexity of NGS data and delivering the concluding outcomes.

PhasiRNAs in Plants: Their Biogenesis, Genic Sources, and Roles in Stress Responses, Development, and Reproduction

Phased secondary small interfering RNAs (phasiRNAs) constitute a major category of small RNAs in plants, but most of their functions are still poorly defined. Some phasiRNAs, known as trans-acting siRNAs, are known to target complementary mRNAs for degradation and to function in development. However, the targets or biological roles of other phasiRNAs remain speculative. New insights into phasiRNA biogenesis, their conservation, and their variation across the flowering plants continue to emerge due to the increased availability of plant genomic sequences, deeper and more sophisticated sequencing approaches, and improvements in computational biology and biochemical/molecular/genetic analyses. In this review, we survey recent progress in phasiRNA biology, with a particular focus on two classes associated with male reproduction: 21-nucleotide (accumulate early in anther ontogeny) and 24-nucloetide (produced in somatic cells during meiosis) phasiRNAs. We describe phasiRNA biogenesis, function, and evolution and define the unanswered questions that represent topics for future research.

Regulation of Female Germline Specification via Small RNA Mobility in Arabidopsis

In the ovules of most sexually reproducing plants, one hypodermal cell differentiates into a megaspore mother cell (MMC), which gives rise to the female germline. Trans-acting small interfering RNAs known as tasiR-ARFs have been suggested to act non-cell-autonomously to prevent the formation of multiple MMCs by repressing AUXIN RESPONSE FACTOR3 (ARF3) expression in Arabidopsis (Arabidopsis thaliana), but the underlying mechanisms are unknown. Here, we examined tasiR-ARF-related intercellular regulatory mechanisms. Expression analysis revealed that components of the tasiR-ARF biogenesis pathway are restricted to distinct ovule cell types, thus limiting tasiR-ARF production to the nucellar epidermis. We also provide data suggesting tasiR-ARF movement along the mediolateral axis into the hypodermal cells and basipetally into the chalaza. Furthermore, we used cell type-specific promoters to express ARF3m, which is resistant to tasiR-ARF regulation, in different ovule cell layers. ARF3m expression in hypodermal cells surrounding the MMC, but not in epidermal cells, led to a multiple-MMC phenotype, suggesting that tasiR-ARFs repress ARF3 in these hypodermal cells to suppress ectopic MMC fate. RNA sequencing analyses in plants with hypodermally expressed ARF3m showed that ARF3 potentially regulates MMC specification through phytohormone pathways. Our findings uncover intricate spatial restriction of tasiR-ARF biogenesis, which together with tasiR-ARF mobility enables cell-cell communication in MMC differentiation.

A conserved sequence signature is essential for robust plant miRNA biogenesis

Micro (mi)RNAs are 20–22nt long non-coding RNA molecules involved in post-transcriptional silencing of targets having high base-pair complementarity. Plant miRNAs are processed from long Pol II-transcripts with specific stem-loop structures by Dicer-like (DCL) 1 protein. Although there were reports indicating how a specific region is selected for miRNA biogenesis, molecular details were unclear. Here, we show that the presence of specific GC-rich sequence signature within miRNA/miRNA* region is required for the precise miRNA biogenesis. The involvement of GC-rich signatures in precise processing and abundance of miRNAs was confirmed through detailed molecular and functional analysis. Consistent with the presence of the miRNA-specific GC signature, target RNAs of miRNAs also possess conserved complementary sequence signatures in their miRNA binding motifs. The selection of these GC signatures was dependent on an RNA binding protein partner of DCL1 named HYL1. Finally, we demonstrate a direct application of this discovery for enhancing the abundance and efficiency of artificial miRNAs that are popular in plant functional genomic studies.

Role of phasiRNAs from two distinct phasing frames of GhMYB2 loci in cis- gene regulation in the cotton genome

BACKGROUND

Phased small interfering RNA (phasiRNA) is primarily derived from the 22-nt miRNA targeting loci. GhMYB2, a gene with potential roles in cotton fiber cell fate determination, is a target gene of miR828 and miR858 in the generation of phasiRNAs.

RESULTS

In the presented work, through the evaluation of phasing scores and phasiRNA distribution pattern, we found that phasiRNAs from GhMYB2 were derived from the 3' cleavage fragments of 22-nt miR828 and 21-nt miR858 respectively. These two miRNA targeting sites initiated two phasing frames on transcripts of one locus. By means of RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE), we further demonstrated that phasiRNAs derived from the two phasing frames played a role in cis-regulation of GhMYB2. The phasiRNAs derived from GhMYB2 were expressed in the somatic tissues, especially in anther and hypocotyl. We further employed our previous small RNA sequencing data as well as the degradome data of cotton fiber bearing ovules, anthers, hypocotyls and embryogenic calli tissues published in public databases, to validate the expression, phasing pattern and functions of phasiRNAs.

CONCLUSIONS

The presenting research provide insights of the molecular mechanism of phasiRNAs in regulation of GhMYB2 loci.

Next-generation forward genetic screens: using simulated data to improve the design of mapping-by-sequencing experiments in Arabidopsis

Forward genetic screens have successfully identified many genes and continue to be powerful tools for dissecting biological processes in Arabidopsis and other model species. Next-generation sequencing technologies have revolutionized the time-consuming process of identifying the mutations that cause a phenotype of interest. However, due to the cost of such mapping-by-sequencing experiments, special attention should be paid to experimental design and technical decisions so that the read data allows to map the desired mutation. Here, we simulated different mapping-by-sequencing scenarios. We first evaluated which short-read technology was best suited for analyzing gene-rich genomic regions in Arabidopsis and determined the minimum sequencing depth required to confidently call single nucleotide variants. We also designed ways to discriminate mutagenesis-induced mutations from background Single Nucleotide Polymorphisms in mutants isolated in Arabidopsis non-reference lines. In addition, we simulated bulked segregant mapping populations for identifying point mutations and monitored how the size of the mapping population and the sequencing depth affect mapping precision. Finally, we provide the computational basis of a protocol that we already used to map T-DNA insertions with paired-end Illumina-like reads, using very low sequencing depths and pooling several mutants together; this approach can also be used with single-end reads as well as to map any other insertional mutagen. All these simulations proved useful for designing experiments that allowed us to map several mutations in Arabidopsis.

Conserved and novel roles of miR164-CUC2 regulatory module in specifying leaf and floral organ morphology in strawberry

MicroRNAs (miRNAs) are a kind of short noncoding RNA (20-24 nt), playing versatile roles in plant growth and development. Strawberry generates leaves and flowers with unique features. However, few miRNAs have been functionally characterised in strawberry, especially for their developmental regulation. Here, we identified one ethyl methanesulfonate (EMS) mutant, deeply serrated (des), in the woodland strawberry Fragaria vesca that has wrinkled leaves with deeper serrations, serrated petals and deformed carpels. The causative mutation occurs in the 19^th nucleotide of the FvemiR164a mature sequence. Overexpressing FveMIR164A rescued the phenotypes of des/fvemir164a except the petal serrations. Furthermore, we identified two allelic mutants of FveCUC2a, one target of FvemiR164a, which developed leaves with smooth margins and fused leaflets. Phenotypes of the double mutant fvemir164a fvecuc2a indicated that the two genes act linearly in leaf and carpel development, but synergistically in the development of other floral organs and inflorescence architecture. This work demonstrates the conserved and novel roles of the miR164-CUC2 module in leaf and flower development in different plant species, and reveals that the 19^th nucleotide of FvemiR164a is important for its processing.

Efficiency and precision of microRNA biogenesis modes in plants

Many evolutionarily conserved microRNAs (miRNAs) in plants regulate transcription factors with key functions in development. Hence, mutations in the core components of the miRNA biogenesis machinery cause strong growth defects. An essential aspect of miRNA biogenesis is the precise excision of the small RNA from its precursor. In plants, miRNA precursors are largely variable in size and shape and can be processed by different modes. Here, we optimized an approach to detect processing intermediates during miRNA biogenesis. We characterized a miRNA whose processing is triggered by a terminal branched loop. Plant miRNA processing can be initiated by internal bubbles, small terminal loops or branched loops followed by dsRNA segments of 15–17 bp. Interestingly, precision and efficiency vary with the processing modes. Despite the various potential structural determinants present in a single a miRNA precursor, DCL1 is mostly guided by a predominant structural region in each precursor in wild-type plants. However, our studies in fiery1, hyl1 and se mutants revealed the existence of cleavage signatures consistent with the recognition of alternative processing determinants. The results provide a general view of the mechanisms underlying the specificity of miRNA biogenesis in plants.

Directions for research and training in plant omics: Big Questions and Big Data

A key remit of the NSF-funded "Arabidopsis Research and Training for the 21st Century" (ART-21) Research Coordination Network has been to convene a series of workshops with community members to explore issues concerning research and training in plant biology, including the role that research using Arabidopsis thaliana can play in addressing those issues. A first workshop focused on training needs for bioinformatic and computational approaches in plant biology was held in 2016, and recommendations from that workshop have been published (Friesner et al., Plant Physiology, 175, 2017, 1499). In this white paper, we provide a summary of the discussions and insights arising from the second ART-21 workshop. The second workshop focused on experimental aspects of omics data acquisition and analysis and involved a broad spectrum of participants from academics and industry, ranging from graduate students through post-doctorates, early career and established investigators. Our hope is that this article will inspire beginning and established scientists, corporations, and funding agencies to pursue directions in research and training identified by this workshop, capitalizing on the reference species Arabidopsis thaliana and other valuable plant systems.

Identification and profiling of narrow-leafed lupin (Lupinus angustifolius) microRNAs during seed development

BACKGROUND

Whilst information regarding small RNAs within agricultural crops is increasing, the miRNA composition of the nutritionally valuable pulse narrow-leafed lupin (Lupinus angustifolius) remains unknown.

RESULTS

By conducting a genome- and transcriptome-wide survey we identified 7 Dicer-like and 16 Argonaute narrow-leafed lupin genes, which were highly homologous to their legume counterparts. We identified 43 conserved miRNAs belonging to 16 families, and 13 novel narrow-leafed lupin-specific miRNAs using high-throughput sequencing of small RNAs from foliar and root and five seed development stages. We observed up-regulation of members of the miRNA families miR167, miR399, miR156, miR319 and miR164 in narrow-leafed lupin seeds, and confirmed expression of miR156, miR166, miR164, miR1507 and miR396 using quantitative RT-PCR during five narrow-leafed lupin seed development stages. We identified potential targets for the conserved and novel miRNAs and were able to validate targets of miR399 and miR159 using 5' RLM-RACE. The conserved miRNAs are predicted to predominately target transcription factors and 93% of the conserved miRNAs originate from intergenic regions. In contrast, only 43% of the novel miRNAs originate from intergenic regions and their predicted targets were more functionally diverse.

CONCLUSION

This study provides important insights into the miRNA gene regulatory networks during narrow-leafed lupin seed development.

The role of a class III gibberellin 2-oxidase in tomato internode elongation

A network of environmental inputs and internal signaling controls plant growth, development and organ elongation. In particular, the growth-promoting hormone gibberellin (GA) has been shown to play a significant role in organ elongation. The use of tomato as a model organism to study elongation presents an opportunity to study the genetic control of internode-specific elongation in a eudicot species with a sympodial growth habit and substantial internodes that can and do respond to external stimuli. To investigate internode elongation, a mutant with an elongated hypocotyl and internodes but wild-type petioles was identified through a forward genetic screen. In addition to stem-specific elongation, this mutant, named tomato internode elongated -1 (tie-1) is more sensitive to the GA biosynthetic inhibitor paclobutrazol and has altered levels of intermediate and bioactive GAs compared with wild-type plants. The mutation responsible for the internode elongation phenotype was mapped to GA2oxidase 7, a class III GA 2-oxidase in the GA biosynthetic pathway, through a bulked segregant analysis and bioinformatic pipeline, and confirmed by transgenic complementation. Furthermore, bacterially expressed recombinant TIE protein was shown to have bona fide GA 2-oxidase activity. These results define a critical role for this gene in internode elongation and are significant because they further the understanding of the role of GA biosynthetic genes in organ-specific elongation.

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.

Consensus Coexpression Network Analysis Identifies Key Regulators of Flower and Fruit Development in Wild Strawberry

The diploid strawberry, Fragaria vesca, is a developing model system for the economically important Rosaceae family. Strawberry fleshy fruit develops from the floral receptacle and its ripening is nonclimacteric. The external seed configuration of strawberry fruit facilitates the study of seed-to-fruit cross tissue communication, particularly phytohormone biosynthesis and transport. To investigate strawberry fruit development, we previously generated spatial and temporal transcriptome data profiling F. vesca flower and fruit development pre- and postfertilization. In this study, we combined 46 of our existing RNA-seq libraries to generate coexpression networks using the Weighted Gene Co-Expression Network Analysis package in R. We then applied a post-hoc consensus clustering approach and used bootstrapping to demonstrate consensus clustering's ability to produce robust and reproducible clusters. Further, we experimentally tested hypotheses based on the networks, including increased iron transport from the receptacle to the seed postfertilization and characterized a F. vesca floral mutant and its candidate gene. To increase their utility, the networks are presented in a web interface (www.fv.rosaceaefruits.org) for easy exploration and identification of coexpressed genes. Together, the work reported here illustrates ways to generate robust networks optimized for the mining of large transcriptome data sets, thereby providing a useful resource for hypothesis generation and experimental design in strawberry and related Rosaceae fruit crops.

Characterization of DCL4 missense alleles provides insights into its ability to process distinct classes of dsRNA substrates

In the model plant Arabidopsis thaliana, four Dicer-like proteins (DCL1-4) mediate the production of various classes of small RNAs (sRNAs). Among these four proteins, DCL4 is by far the most versatile RNaseIII-like enzyme, and previously identified dcl4 missense alleles were shown to uncouple the production of the various classes of DCL4-dependent sRNAs. Yet little is known about the molecular mechanism behind this uncoupling. Here, by studying the subcellular localization, interactome and binding to the sRNA precursors of three distinct dcl4 missense alleles, we simultaneously highlight the absolute requirement of a specific residue in the helicase domain for the efficient production of all DCL4-dependent sRNAs, and identify, within the PAZ domain, an important determinant of DCL4 versatility that is mandatory for the efficient processing of intramolecular fold-back double-stranded RNA (dsRNA) precursors, but that is dispensable for the production of small interfering RNAs (siRNAs) from RDR-dependent dsRNA susbtrates. This study not only provides insights into the DCL4 mode of action, but also delineates interesting tools to further study the complexity of RNA silencing pathways in plants, and possibly other organisms.

Mapping causal mutations by exome sequencing in a wheat TILLING population: a tall mutant case study

Forward genetic screens of induced mutant plant populations are powerful tools to identify genes underlying phenotypes of interest. Using traditional techniques, mapping causative mutations from forward screens is a lengthy, multi-step process, requiring the identification of a broad genetic region followed by candidate gene sequencing to characterize the causal variant. Mapping by whole genome sequencing accelerates the identification of causal mutations by simultaneously defining a mapping region and providing information on the induced genetic variants. In wheat, although the availability of a high-quality draft genome assembly facilitates mapping and mutation calling, whole genome resequencing remains prohibitively expensive due to its large genome. In the current study, we used exome sequencing as a complexity reduction strategy to detect mutations associated with a target phenotype. In a segregating wheat EMS population, we identified a clear peak region on chromosome arm 4BS associated with increased plant height. Although none of the significant SNPs seemed causative for the mutant phenotype, they were sufficient to identify a linked ~ 1.9 Mb deletion encompassing nine genes. These genes included Rht-B1, which is known to have a strong effect on plant height and is a strong candidate for the observed phenotype. We performed simulation experiments to determine the impacts of sequencing depth and bulk size and discuss the importance of considering each factor when designing mapping-by-sequencing experiments in wheat. This approach can accelerate the identification of candidate causal point mutations or linked deletions underlying important phenotypes.

Mapping causal mutations by exome sequencing in a wheat TILLING population: a tall mutant case study

Forward genetic screens of induced mutant plant populations are powerful tools to identify genes underlying phenotypes of interest. Using traditional techniques, mapping causative mutations from forward screens is a lengthy, multi-step process, requiring the identification of a broad genetic region followed by candidate gene sequencing to characterize the causal variant. Mapping by whole genome sequencing accelerates the identification of causal mutations by simultaneously defining a mapping region and providing information on the induced genetic variants. In wheat, although the availability of a high-quality draft genome assembly facilitates mapping and mutation calling, whole genome resequencing remains prohibitively expensive due to its large genome. In the current study, we used exome sequencing as a complexity reduction strategy to detect mutations associated with a target phenotype. In a segregating wheat EMS population, we identified a clear peak region on chromosome arm 4BS associated with increased plant height. Although none of the significant SNPs seemed causative for the mutant phenotype, they were sufficient to identify a linked ~ 1.9 Mb deletion encompassing nine genes. These genes included Rht-B1, which is known to have a strong effect on plant height and is a strong candidate for the observed phenotype. We performed simulation experiments to determine the impacts of sequencing depth and bulk size and discuss the importance of considering each factor when designing mapping-by-sequencing experiments in wheat. This approach can accelerate the identification of candidate causal point mutations or linked deletions underlying important phenotypes.

Rice stripe virus NS3 protein regulates primary miRNA processing through association with the miRNA biogenesis factor OsDRB1 and facilitates virus infection in rice

MicroRNAs (miRNAs) are small regulatory RNAs processed from primary miRNA transcripts, and plant miRNAs play important roles in plant growth, development, and response to infection by microbes. Microbial infections broadly alter miRNA biogenesis, but the underlying mechanisms remain poorly understood. In this study, we report that the Rice stripe virus (RSV)-encoded nonstructural protein 3 (NS3) interacts with OsDRB1, an indispensable component of the rice (Oryza sativa) miRNA-processing complex. Moreover, the NS3-OsDRB1 interaction occurs at the sites required for OsDRB1 self-interaction, which is essential for miRNA biogenesis. Further analysis revealed that NS3 acts as a scaffold between OsDRB1 and pri-miRNAs to regulate their association and aids in vivo processing of pri-miRNAs. Genetic evidence in Arabidopsis showed that NS3 can partially substitute for the function of double-stranded RNA binding domain (dsRBD) of AtDRB1/AtHYL1 during miRNA biogenesis. As a result, NS3 induces the accumulation of several miRNAs, most of which target pivotal genes associated with development or pathogen resistance. In contrast, a mutant version of NS3 (mNS3), which still associated with OsDRB1 but has defects in pri-miRNA binding, reduces accumulation of these miRNAs. Transgenic rice lines expressing NS3 exhibited significantly higher susceptibility to RSV infection compared with non-transgenic wild-type plants, whereas the transgenic lines expressing mNS3 showed a less-sensitive response. Our findings revealed a previously unknown mechanism in which a viral protein hijacks OsDRB1, a key component of the processing complex, for miRNA biogenesis and enhances viral infection and pathogenesis in rice.

MicroRNAs (miRNAs) regulate gene expression at the transcriptional or post-transcriptional level and have emerged as key players in regulating plant growth, development and response to biotic and abiotic stresses. Accumulating evidences suggest that miRNAs are pivotal modulators of host–virus interactions, but how virus regulates miRNA accumulation remains poorly understood. Here, we report that NS3 protein encoded by Rice stripe virus (RSV) regulates the processing of several primary miRNA transcripts (pri-miRNAs) by acting as an intermediary to modulate the association of pri-miRNAs and OsDRB1, a key factor of the pri-miRNA processing complex. NS3 increases recruitment of pri-miRNA to the processing complex by its association with OsDRB1 at the sites required for OsDRB1 dimer formation and induces several miRNAs accumulations as well as target genes repression, promoting the sensitivity of rice to RSV infection. Together these findings reveal a novel mechanism by which RSV regulates pri-miRNA processing, leading to enhanced viral infection.

Sequence-Based Mapping and Genome Editing Reveal Mutations in Stickleback Hps5 Cause Oculocutaneous Albinism and the casper Phenotype

Here, we present and characterize the spontaneous X-linked recessive mutation casper, which causes oculocutaneous albinism in threespine sticklebacks (Gasterosteus aculeatus). In humans, Hermansky-Pudlak syndrome results in pigmentation defects due to disrupted formation of the melanin-containing lysosomal-related organelle (LRO), the melanosome. casper mutants display not only reduced pigmentation of melanosomes in melanophores, but also reductions in the iridescent silver color from iridophores, while the yellow pigmentation from xanthophores appears unaffected. We mapped casper using high-throughput sequencing of genomic DNA from bulked casper mutants to a region of the stickleback X chromosome (chromosome 19) near the stickleback ortholog of Hermansky-Pudlak syndrome 5 (Hps5). casper mutants have an insertion of a single nucleotide in the sixth exon of Hps5, predicted to generate an early frameshift. Genome editing using CRISPR/Cas9 induced lesions in Hps5 and phenocopied the casper mutation. Injecting single or paired Hps5 guide RNAs revealed higher incidences of genomic deletions from paired guide RNAs compared to single gRNAs. Stickleback Hps5 provides a genetic system where a hemizygous locus in XY males and a diploid locus in XX females can be used to generate an easily scored visible phenotype, facilitating quantitative studies of different genome editing approaches. Lastly, we show the ability to better visualize patterns of fluorescent transgenic reporters in Hps5 mutant fish. Thus, Hps5 mutations present an opportunity to study pigmented LROs in the emerging stickleback model system, as well as a tool to aid in assaying genome editing and visualizing enhancer activity in transgenic fish.

Evolutionary Footprints Reveal Insights into Plant MicroRNA Biogenesis

MicroRNAs (miRNAs) are endogenous small RNAs that recognize target sequences by base complementarity and play a role in the regulation of target gene expression. They are processed from longer precursor molecules that harbor a fold-back structure. Plant miRNA precursors are quite variable in size and shape, and are recognized by the processing machinery in different ways. However, ancient miRNAs and their binding sites in target genes are conserved during evolution. Here, we designed a strategy to systematically analyze MIRNAs from different species generating a graphical representation of the conservation of the primary sequence and secondary structure. We found that plant MIRNAs have evolutionary footprints that go beyond the small RNA sequence itself, yet their location along the precursor depends on the specific MIRNA We show that these conserved regions correspond to structural determinants recognized during the biogenesis of plant miRNAs. Furthermore, we found that the members of the miR166 family have unusual conservation patterns and demonstrated that the recognition of these precursors in vivo differs from other known miRNAs. Our results describe a link between the evolutionary conservation of plant MIRNAs and the mechanisms underlying the biogenesis of these small RNAs and show that the MIRNA pattern of conservation can be used to infer the mode of miRNA biogenesis.

Tetrahymena as a Unicellular Model Eukaryote: Genetic and Genomic Tools

Tetrahymena thermophila is a ciliate model organism whose study has led to important discoveries and insights into both conserved and divergent biological processes. In this review, we describe the tools for the use of Tetrahymena as a model eukaryote, including an overview of its life cycle, orientation to its evolutionary roots, and methodological approaches to forward and reverse genetics. Recent genomic tools have expanded Tetrahymena's utility as a genetic model system. With the unique advantages that Tetrahymena provide, we argue that it will continue to be a model organism of choice.

Bioinformatics Study of Structural Patterns in Plant MicroRNA Precursors

According to the RNA world theory, RNAs which stored genetic information and catalyzed chemical reactions had their contribution in the formation of current living organisms. In recent years, researchers studied this molecule diversity, i.a. focusing on small non-coding regulatory RNAs. Among them, of particular interest is evolutionarily ancient, 19–24 nt molecule of microRNA (miRNA). It has been already recognized as a regulator of gene expression in eukaryotes. In plants, miRNA plays a key role in the response to stress conditions and it participates in the process of growth and development. MicroRNAs originate from primary transcripts (pri-miRNA) encoded in the nuclear genome. They are processed from single-stranded stem-loop RNA precursors containing hairpin structures. While the mechanism of mature miRNA production in animals is better understood, its biogenesis in plants remains less clear. Herein, we present the results of bioinformatics analysis aimed at discovering how plant microRNAs are recognized within their precursors (pre-miRNAs). The study has been focused on sequential and structural motif identification in the neighbourhood of microRNA.

Base-pair opening dynamics of the microRNA precursor pri-miR156a affect temperature-responsive flowering in Arabidopsis

Internal and environmental cues, including ambient temperature changes, regulate the timing of flowering in plants. Arabidopsis miR156 represses flowering and plays an important role in the regulation of temperature-responsive flowering. However, the molecular basis of miR156 processing at lower temperatures remains largely unknown. Here, we performed nuclear magnetic resonance studies to investigate the base-pair opening dynamics of model RNAs at 16 °C and investigated the in vivo effects of the mutant RNAs on temperature-responsive flowering. The A9C and A10CG mutations in the B5 bulge of the lower stem of pri-miR156a stabilized the C15∙G98 and U16∙A97 base-pairs at the cleavage site of pri-miR156a at 16 °C. Consistent with this, production of mature miR156 was severely affected in plants overexpressing the A9C and A10CG constructs and these plants exhibited almost no delay in flowering at 16 °C. The A10G and A9AC mutations did not strongly affect C15∙G98 and U16∙A97 base-pairs at 16 °C, and plants overexpressing A10G and A9AC mutants of miR156 produced more mature miR156 than plants overexpressing the A9C and A10CG mutants and showed a strong delay in flowering at 16 °C. Interestingly, the A9AC mutation had distinct effects on the opening dynamics of the C15∙G98 and U16∙A97 base-pairs between 16 °C and 23 °C, and plants expressing the A9AC mutant miR156 showed only a moderate delay in flowering at 16 °C. Based on these results, we propose that fine-tuning of the base-pair stability at the cleavage site is essential for efficient processing of pri-miR156a at a low temperature and for reduced flowering sensitivity to ambient temperature changes.

Origin and evolution of MIR1444 genes in Salicaceae

miR1444s are functionally significant miRNAs targeting polyphenol oxidase (PPO) genes for cleavage. MIR1444 genes were reported only in Populus trichocarpa. Through the computational analysis of 215 RNA-seq data, four whole genome sequences of Salicaceae species and deep sequencing of six P. trichocarpa small RNA libraries, we investigated the origin and evolution history of MIR1444s. A total of 23 MIR1444s were identified. Populus and Idesia species contain two MIR1444 genes, while Salix includes only one. Populus and Idesia MIR1444b genes and Salix MIR1444s were phylogenetically separated from Populus and Idesia MIR1444a genes. Ptr-miR1444a and ptr-miR1444b showed sequence divergence. Compared with ptr-miR1444b, ptr-miR1444a started 2 nt upstream of precursor, resulting in differential regulation of PPO targets. Sequence alignments showed that MIR1444 genes exhibited extensive similarity to their PPO targets, the characteristics of MIRs originated from targets through an inverted gene duplication event. Genome sequence comparison showed that MIR1444 genes in Populus and Idesia were expanded through the Salicoid genome duplication event. A copy of MIR1444 gene was lost in Salix through DNA segment deletion during chromosome rearrangements. The results provide significant information for the origin of plant miRNAs and the mechanism of Salicaceae gene evolution and divergence.

Sensitive whole mount in situ localization of small RNAs in plants

Small RNAs, such as microRNAs (miRNAs), regulate gene expression and play important roles in many plant processes. Although our knowledge of their biogenesis and mode of action has significantly progressed, we still have comparatively little information about their biological functions. In particular, knowledge about their spatio-temporal expression patterns rely on either indirect detection by use of reporter constructs or labor-intensive direct detection by in situ hybridization on sectioned material. None of the current approaches allows a systematic investigation of small RNA expression patterns. Here, we present a sensitive method for in situ detection of miRNAs and siRNAs in intact plant tissues that utilizes both double-labeled probes and a specific cross-linker. We determined the expression patterns of several small RNAs in diverse plant tissues.

Proximal disruption of base pairing of the second stem in the upper stem of pri-miR156a caused ambient temperature-sensitive flowering in Arabidopsis

MicroRNAs are generated from primary transcripts (pri-miRNAs) that form hairpin structures. Plant miRNAs play an important role in regulating flowering; however, little is known about the role of their structures in ambient temperature-responsive flowering. We recently showed that disruption of base pairing in the second stem (S2) in the upper stem of pri-miR156a caused hypersensitive flowering in response to ambient temperature changes. To further substantiate our findings on the role of S2 of pri-miR156a, we analyzed the effects of serial disruption (from the proximal or distal sides) of base-pairing in S2 of pri-miR156a on temperature-dependent flowering. We found that flowering time was gradually delayed with increasing size of the proximal disruption of S2 at 16°C. Particularly, disrupting base pairing of 5 nucleotides from the proximal side caused flowering to be hypersensitive to ambient temperature changes, which is similar to the phenotype of plants overexpressing pri-miR156a with a disruption of S2 (156-DBP-S2). However, disrupting base pairing from the distal side did not cause late flowering at 16°C and thus did not cause temperature-sensitive flowering. These results supported our notion that the second stem (S2) in the upper stem of pri-miR156a plays a role in the regulation of ambient temperature-responsive flowering.

Base-pair opening dynamics of primary miR156a using NMR elucidates structural determinants important for its processing level and leaf number phenotype in Arabidopsis

MicroRNAs originate from primary transcripts containing hairpin structures. The levels of mature miR156 influence the leaf number prior to flowering in the life cycle of plants. To understand the molecular mechanism of biogenesis of primary miR156a (pri-miR156a) to mature miR156, a base-pair opening dynamics study was performed using model RNAs mimicking the cleavage site of wild type and B5 bulge-stabilizing mutant pri-miR156a constructs. We also determined the mature miR156 levels and measured leaf numbers at flowering of plants overexpressing the wild type and mutant constructs. Our results suggest that the stabilities and/or opening dynamics of the C15·G98 and U16·A97 base-pairs at the cleavage site are essential for formation of the active conformation and for efficient processing of pri-miR156a, and that mutations of the B5 bulge can modulate mature miR156 levels as well as miR156-driven leaf number phenotypes via changes in the base-pair stability of the cleavage site.

Structural determinants of miR156a precursor processing in temperature-responsive flowering in Arabidopsis

MicroRNAs originate from primary transcripts (pri-miRNAs) containing hairpin structures. Plant pri-miRNAs have highly variable structures and little is known about the information encoded in their secondary structures. Arabidopsis miR156 is an ambient temperature-responsive miRNA and plays an important role in regulating flowering time. To identify the structural determinants for miR156 processing, we analyzed the effects of mutations introduced in the upper stem of pri-miR156a on its temperature-dependent processing and flowering time. The levels of pri-miR156a and mature miR156 were opposite at different temperatures. Mutations in the upper stem, especially the region closer to the miR156a/miR156a* duplex, reduced miR156 processing at 23 °C and 16 °C and caused a less severe phenotype compared with the un-mutated construct. Mutation in the second stem near the first cleavage site of pri-miR156a affected miR156 processing at 23 °C, but not at 16 °C. This was also seen in pri-miR172a, another ambient temperature-responsive miRNA. Replacement of the upper stem of pri-miR156a with that of pri-miR172a severely affected miR156 processing and flowering time. These results suggested that the upper stem of pri-miR156a is important for miR156 processing at different temperatures. In particular, the second stem adjacent to the first cleavage site plays a role in the regulation of ambient temperature-responsive flowering.