Reproductive phasiRNA loci and DICER-LIKE5, but not microRNA loci, diversified in monocotyledonous plants.

In monocots other than maize (Zea mays) and rice (Oryza sativa), the repertoire and diversity of microRNAs (miRNAs) and the populations of phased, secondary, small interfering RNAs (phasiRNAs) are poorly characterized. To remedy this, we sequenced small RNAs (sRNA) from vegetative and dissected inflorescence tissue in 28 phylogenetically diverse monocots and from several early-diverging angiosperm lineages, as well as publicly available data from 10 additional monocot species. We annotated miRNAs, small interfering RNAs (siRNAs) and phasiRNAs across the monocot phylogeny, identifying miRNAs apparently lost or gained in the grasses relative to other monocot families, as well as a number of transfer RNA fragments misannotated as miRNAs. Using our miRNA database cleaned of these misannotations, we identified conservation at the 8th, 9th, 19th, and 3'-end positions that we hypothesize are signatures of selection for processing, targeting, or Argonaute sorting. We show that 21-nucleotide (nt) reproductive phasiRNAs are far more numerous in grass genomes than other monocots. Based on sequenced monocot genomes and transcriptomes, DICER-LIKE5, important to 24-nt phasiRNA biogenesis, likely originated via gene duplication before the diversification of the grasses. This curated database of phylogenetically diverse monocot miRNAs, siRNAs, and phasiRNAs represents a large collection of data that should facilitate continued exploration of sRNA diversification in flowering plants.

Genome-Wide Reinforcement of DNA Methylation Occurs during Somatic Embryogenesis in Soybean.

Somatic embryogenesis is an important tissue culture technique that sometimes leads to phenotypic variation via genetic and/or epigenetic changes. To understand the genomic and epigenomic impacts of somatic embryogenesis, we characterized soybean (Glycine max) epigenomes sampled from embryos at 10 different stages ranging from 6 weeks to 13 years of continuous culture. We identified genome-wide increases in DNA methylation from cultured samples, especially at CHH sites. The hypermethylation almost exclusively occurred in regions previously possessing non-CG methylation and was accompanied by increases in the expression of genes encoding the RNA-directed DNA methylation (RdDM) machinery. The epigenomic changes were similar between somatic and zygotic embryogenesis. Following the initial global wave of hypermethylation, rare decay events of maintenance methylation were observed, and the extent of the decay increased with time in culture. These losses in DNA methylation were accompanied by downregulation of genes encoding the RdDM machinery and transcriptome reprogramming reminiscent of transcriptomes during late-stage seed development. These results reveal a process for reinforcing already silenced regions to maintain genome integrity during somatic embryogenesis over the short term, which eventually decays at certain loci over longer time scales.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

PMS1T, producing phased small-interfering RNAs, regulates photoperiod-sensitive male sterility in rice.

Phased small-interfering RNAs (phasiRNAs) are a special class of small RNAs, which are generated in 21- or 24-nt intervals from transcripts of precursor RNAs. Although phasiRNAs have been found in a range of organisms, their biological functions in plants have yet to be uncovered. Here we show that phasiRNAs generated by the photopheriod-sensetive genic male sterility 1 (Pms1) locus were associated with photoperiod-sensitive male sterility (PSMS) in rice, a germplasm that started the two-line hybrid rice breeding. The Pms1 locus encodes a long-noncoding RNA PMS1T that was preferentially expressed in young panicles. PMS1T was targeted by miR2118 to produce 21-nt phasiRNAs that preferentially accumulated in the PSMS line under long-day conditions. A single nucleotide polymorphism in PMS1T nearby the miR2118 recognition site was critical for fertility change, likely leading to differential accumulation of the phasiRNAs. This result suggested possible roles of phasiRNAs in reproductive development of rice, demonstrating the potential importance of this RNA class as regulators in biological processes.

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses.

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

Physiological stressors and invasive plant infections alter the small RNA transcriptome of the rice blast fungus, Magnaporthe oryzae.

BACKGROUND: The rice blast fungus, Magnaporthe oryzae is a destructive pathogen of rice and other related crops, causing significant yield losses worldwide. Endogenous small RNAs (sRNAs), including small interfering RNAs (siRNAs) and microRNAs (miRNAs) are critical components of gene regulation in many eukaryotic organisms. Recently several new species of sRNAs have been identified in fungi. This fact along with the availability of genome sequence makes M. oryzae a compelling target for sRNA profiling. We have examined sRNA species and their biosynthetic genes in M. oryzae, and the degree to which these elements regulate fungal stress responses. To this end, we have characterized sRNAs under different physiological stress conditions, which had not yet been examined in this fungus. RESULTS: The resulting libraries are composed of more than 37 million total genome matched reads mapping to intergenic regions, coding sequences, retrotransposons, inverted, tandem, and other repeated regions of the genome with more than half of the small RNAs arising from intergenic regions. The 24 nucleotide (nt) size class of sRNAs was predominant. A comparison to transcriptional data of M. oryzae undergoing the same physiological stresses indicates that sRNAs play a role in transcriptional regulation for a small subset of genes. Support for this idea comes from generation and characterization of mutants putatively involved in sRNAs biogenesis; our results indicate that the deletion of Dicer-like genes and an RNA-Dependent RNA Polymerase gene increases the transcriptional regulation of this subset of genes, including one involved in virulence. CONCLUSIONS: Various physiological stressors and in planta conditions alter the small RNA profile of the rice blast fungus. Characterization of sRNA biosynthetic mutants helps to clarify the role of sRNAs in transcriptional control.

Transcriptome profiling of the rice blast fungus during invasive plant infection and in vitro stresses.

BACKGROUND: Rice blast is the most threatening disease to cultivated rice. Magnaporthe oryzae, its causal agent, is likely to encounter environmental challenges during invasive growth in its host plants that require shifts in gene expression to establish a compatible interaction. Here, we tested the hypothesis that gene expression patterns during in planta invasive growth are similar to in vitro stress conditions, such as nutrient limitation, temperature up shift and oxidative stress, and determined which condition most closely mimicked that of in planta invasive growth. Gene expression data were collected from these in vitro experiments and compared to fungal gene expression during the invasive growth phase at 72 hours post-inoculation in compatible interactions on two grass hosts, rice and barley. RESULTS: We identified 4,973 genes that were differentially expressed in at least one of the in planta and in vitro stress conditions when compared to fungal mycelia grown in complete medium, which was used as reference. From those genes, 1,909 showed similar expression patterns between at least one of the in vitro stresses and rice and/or barley. Hierarchical clustering of these 1,909 genes showed three major clusters in which in planta conditions closely grouped with the nutrient starvation conditions. Out of these 1,909 genes, 55 genes and 129 genes were induced and repressed in all treatments, respectively. Functional categorization of the 55 induced genes revealed that most were either related to carbon metabolism, membrane proteins, or were involved in oxidoreduction reactions. The 129 repressed genes showed putative roles in vesicle trafficking, signal transduction, nitrogen metabolism, or molecular transport. CONCLUSIONS: These findings suggest that M. oryzae is likely primarily coping with nutrient-limited environments at the invasive growth stage 72 hours post-inoculation, and not with oxidative or temperature stresses.

Construction and characterization of a knock-down RNA interference line of OsNRPD1 in rice (Oryza sativa ssp japonica cv Nipponbare).

In plants, RNA-directed DNA methylation (RdDM) is a silencing mechanism relying on the production of 24-nt small interfering RNAs (siRNAs) by RNA POLYMERASE IV (Pol IV) to trigger methylation and inactivation of transposable elements (TEs). We present the construction and characterization of osnrpd1, a knock-down RNA interference line of OsNRPD1 gene that encodes the largest subunit of Pol IV in rice (Oryza sativa ssp japonica cv Nipponbare). We show that osnrpd1 displays a lower accumulation of OsNRPD1 transcripts, associated with an overall reduction of 24-nt siRNAs and DNA methylation level in all three contexts, CG, CHG and CHH. We uncovered new insertions of known active TEs, the LTR retrotransposons Tos17 and Lullaby and the long interspersed nuclear element-type retrotransposon Karma. However, we did not observe any clear developmental phenotype, contrary to what was expected for a mutant severely affected in RdDM. In addition, despite the presence of many putatively functional TEs in the rice genome, we found no evidence of in planta global reactivation of transposition. This knock-down of OsNRPD1 likely led to a weakly affected line, with no effect on development and a limited effect on transposition. We discuss the possibility that a knock-out mutation of OsNRPD1 would cause sterility in rice. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Characterization of Plant Small RNAs by Next Generation Sequencing.

Plant small RNAs are ∼20 to 24 nucleotide noncoding RNAs that typically have repressive regulatory roles in gene expression, functioning at the transcriptional or post-transcriptional level. This influence on regulation of developmental and physiological processes has direct effects on phenotype. High-throughput sequencing technologies have enabled the sequencing of millions of small RNAs. Along with decreased sequencing costs, recent improvements in small RNA library construction have facilitated the ability to use minimal amounts of input RNA for analysis. This unit describes steps to isolate total RNA from limited amounts of plant tissue to construct small RNA libraries and perform small RNA data processing. © 2017 by John Wiley & Sons, Inc.

The asparagus genome sheds light on the origin and evolution of a young Y chromosome.

Sex chromosomes evolved from autosomes many times across the eukaryote phylogeny. Several models have been proposed to explain this transition, some involving male and female sterility mutations linked in a region of suppressed recombination between X and Y (or Z/W, U/V) chromosomes. Comparative and experimental analysis of a reference genome assembly for a double haploid YY male garden asparagus (Asparagus officinalis L.) individual implicates separate but linked genes as responsible for sex determination. Dioecy has evolved recently within Asparagus and sex chromosomes are cytogenetically identical with the Y, harboring a megabase segment that is missing from the X. We show that deletion of this entire region results in a male-to-female conversion, whereas loss of a single suppressor of female development drives male-to-hermaphrodite conversion. A single copy anther-specific gene with a male sterile Arabidopsis knockout phenotype is also in the Y-specific region, supporting a two-gene model for sex chromosome evolution.

MicroRNA Maturation and MicroRNA Target Gene Expression Regulation Are Severely Disrupted in Soybean dicer-like1 Double Mutants.

Small nonprotein-coding microRNAs (miRNAs) are present in most eukaryotes and are central effectors of RNA silencing-mediated mechanisms for gene expression regulation. In plants, DICER-LIKE1 (DCL1) is the founding member of a highly conserved family of RNase III-like endonucleases that function as core machinery proteins to process hairpin-like precursor transcripts into mature miRNAs, small regulatory RNAs, 21-22 nucleotides in length. Zinc finger nucleases (ZFNs) were used to generate single and double-mutants of putative soybean DCL1 homologs, DCL1a and DCL1b, to confirm their functional role(s) in the soybean miRNA pathway. Neither DCL1 single mutant, dcl1a or dcl1b plants, exhibited a pronounced morphological or molecular phenotype. However, the dcl1a/dcl1b double mutant expressed a strong morphological phenotype, characterized by reduced seed size and aborted seedling development, in addition to defective miRNA precursor transcript processing efficiency and deregulated miRNA target gene expression. Together, these findings indicate that the two soybean DCL1 paralogs, DCL1a and DCL1b, largely play functionally redundant roles in the miRNA pathway and are essential for normal plant development.

Transcriptomics of the rice blast fungus Magnaporthe oryzae in response to the bacterial antagonist Lysobacter enzymogenes reveals candidate fungal defense response genes.

Plants and animals have evolved a first line of defense response to pathogens called innate or basal immunity. While basal defenses in these organisms are well studied, there is almost a complete lack of understanding of such systems in fungal species, and more specifically, how they are able to detect and mount a defense response upon pathogen attack. Hence, the goal of the present study was to understand how fungi respond to biotic stress by assessing the transcriptional profile of the rice blast pathogen, Magnaporthe oryzae, when challenged with the bacterial antagonist Lysobacter enzymogenes. Based on microscopic observations of interactions between M. oryzae and wild-type L. enzymogenes strain C3, we selected early and intermediate stages represented by time-points of 3 and 9 hours post-inoculation, respectively, to evaluate the fungal transcriptome using RNA-seq. For comparative purposes, we also challenged the fungus with L. enzymogenes mutant strain DCA, previously demonstrated to be devoid of antifungal activity. A comparison of transcriptional data from fungal interactions with the wild-type bacterial strain C3 and the mutant strain DCA revealed 463 fungal genes that were down-regulated during attack by C3; of these genes, 100 were also found to be up-regulated during the interaction with DCA. Functional categorization of genes in this suite included those with roles in carbohydrate metabolism, cellular transport and stress response. One gene in this suite belongs to the CFEM-domain class of fungal proteins. Another CFEM class protein called PTH11 has been previously characterized, and we found that a deletion in this gene caused advanced lesion development by C3 compared to its growth on the wild-type fungus. We discuss the characterization of this suite of 100 genes with respect to their role in the fungal defense response.

Getting the most out of your fungal microarray data: two cost- and time-effective methods.

Advances in genome sequencing technologies have facilitated production of a wealth of fungal data; within the last 5 years, experimental costs and labor have diminished, shifting the production bottleneck from genomic data generation to data analysis. Genome sequences and microarrays now exist for many fungi, and transcriptional profiling has been shown to be an efficient way to examine how the entire genome changes in response to many different environments or treatments. Multiple platforms, programs, and protocols exist for analyzing such data, making this task daunting for the bench-based scientist. Furthermore, many existing programs are expensive and require license renewals on a yearly basis for each user in the laboratory. Costs may be prohibitively high for bench-based scientists in academia. Our combined experiences with this kind of analysis have favored two programs, depending upon whether the scientist is working with single- or dual-channel hybridization data. Our protocols are aimed toward helping the bench-based PI get the most possible information from their data, without the need for expensive software or an experienced bioinformaticist.