Siwaret Arikit

Spatiotemporally dynamic, cell-type–dependent premeiotic and meiotic phasiRNAs in maize anthers

Significance By RNA profiling of 10 stages of maize anthers plus mature pollen, we found two distinct classes of phased small-interfering RNAs (phasiRNAs): 21-nt premeiotic phasiRNAs, after germinal and somatic cell specification, and 24-nt meiotic phasiRNAs coordinately accumulated during meiosis and persist into pollen. Sequencing of RNA from five male-sterile, anther developmental mutants—ocl4, mac1, ms23, msca1, and ameiotic1—demonstrated the involvement of specific somatic layers. Premeiotic phasiRNAs require a functional epidermis, whereas meiotic phasiRNAs require a normal tapetum. Mammalian germ cells express “prepachytene” or “pachytene” PIWI-interacting RNAs (piRNAs). Whereas differences in biogenesis indicate independent origins, grass phasiRNAs and mammalian piRNAs share developmental timing, a lack of obvious targets, and an impact on male fertility, suggesting a possible evolutionary convergence. Maize anthers, the male reproductive floral organs, express two classes of phased small-interfering RNAs (phasiRNAs). PhasiRNA precursors are transcribed by RNA polymerase II and map to low-copy, intergenic regions similar to PIWI-interacting RNAs (piRNAs) in mammalian testis. From 10 sequential cohorts of staged maize anthers plus mature pollen we find that 21-nt phased siRNAs from 463 loci appear abruptly after germinal and initial somatic cell fate specification and then diminish, whereas 24-nt phasiRNAs from 176 loci coordinately accumulate during meiosis and persist as anther somatic cells mature and haploid gametophytes differentiate into pollen. Male-sterile ocl4 anthers defective in epidermal signaling lack 21-nt phasiRNAs. Male-sterile mutants with subepidermal defects—mac1 (excess meiocytes), ms23 (defective pretapetal cells), and msca1 (no normal soma or meiocytes)—lack 24-nt phasiRNAs. ameiotic1 mutants (normal soma, no meiosis) accumulate both 21-nt and 24-nt phasiRNAs, ruling out meiotic cells as a source or regulator of phasiRNA biogenesis. By in situ hybridization, miR2118 triggers of 21-nt phasiRNA biogenesis localize to epidermis; however, 21-PHAS precursors and 21-nt phasiRNAs are abundant subepidermally. The miR2275 trigger, 24-PHAS precursors, and 24-nt phasiRNAs all accumulate preferentially in tapetum and meiocytes. Therefore, each phasiRNA type exhibits independent spatiotemporal regulation with 21-nt premeiotic phasiRNAs dependent on epidermal and 24-nt meiotic phasiRNAs dependent on tapetal cell differentiation. Maize phasiRNAs and mammalian piRNAs illustrate putative convergent evolution of small RNAs in male reproduction.

An Atlas of Soybean Small RNAs Identifies Phased siRNAs from Hundreds of Coding Genes

An extensive analysis of small RNAs in soybean identified many miRNAs and phased, secondary siRNA (phasiRNA) loci; some of these miRNAs were the triggers of the phasiRNA loci. Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

Plant MicroRNAs Display Differential 3′ Truncation and Tailing Modifications That Are ARGONAUTE1 Dependent and Conserved Across Species

Analysis of hen1 methyltransferase mutants of Arabidopsis, rice, and maize shows widespread 3′ truncation and tailing of microRNAs in patterns that are different among microRNA families but highly conserved across plant species, suggesting an important endogenous function in wild-type plants. Plant small RNAs are 3′ methylated by the methyltransferase HUA1 ENHANCER1 (HEN1). In plant hen1 mutants, 3′ modifications of small RNAs, including oligo-uridylation (tailing), are associated with accelerated degradation of microRNAs (miRNAs). By sequencing small RNAs of the wild type and hen1 mutants from Arabidopsis thaliana, rice (Oryza sativa), and maize (Zea mays), we found 3′ truncation prior to tailing is widespread in these mutants. Moreover, the patterns of miRNA truncation and tailing differ substantially among miRNA families but are conserved across species. The same patterns are also observable in wild-type libraries from a broad range of species, only at lower abundances. ARGONAUTE (AGO1), even with defective slicer activity, can bind these truncated and tailed variants of miRNAs. An ago1 mutation in hen1 suppressed such 3′ modifications, indicating that they occur while miRNAs are in association with AGO1, either during or after RNA-induced silencing complex assembly. Our results showed AGO1-bound miRNAs are actively 3′ truncated and tailed, possibly reflecting the activity of cofactors acting in conserved patterns in miRNA degradation.

Biogenesis and function of rice small RNAs from non-coding RNA precursors.

Non-coding RNAs, especially small RNAs, play important roles in many biological processes. Several small RNA types, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are well-described in rice (Oryza sativa), although much remains to be learned about their function. Many small RNAs along with their targets have been characterized with deep sequencing technologies. Some special classes of these small RNAs have been found to be unique to rice or within the larger group of grasses. The functional and biological roles of numerous plants small RNAs have been described in detail, including functions as varied as the regulation of tissue development, phase transition, or abiotic and biotic stress resistance. Mutant analysis has proven useful in the genetic identification of components involved in small RNA biogenesis and also in identification of regulatory functions of small RNAs. Although many small RNAs have been identified by deep sequencing in rice, their precise regulatory functions and cell-type specificity are in many cases still unknown.

Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce

Lettuce (Lactuca sativa) is a major crop and a member of the large, highly successful Compositae family of flowering plants. Here we present a reference assembly for the species and family. This was generated using whole-genome shotgun Illumina reads plus in vitro proximity ligation data to create large superscaffolds; it was validated genetically and superscaffolds were oriented in genetic bins ordered along nine chromosomal pseudomolecules. We identify several genomic features that may have contributed to the success of the family, including genes encoding Cycloidea-like transcription factors, kinases, enzymes involved in rubber biosynthesis and disease resistance proteins that are expanded in the genome. We characterize 21 novel microRNAs, one of which may trigger phasiRNAs from numerous kinase transcripts. We provide evidence for a whole-genome triplication event specific but basal to the Compositae. We detect 26% of the genome in triplicated regions containing 30% of all genes that are enriched for regulatory sequences and depleted for genes involved in defence.

Extensive Families of miRNAs and PHAS Loci in Norway Spruce Demonstrate the Origins of Complex phasiRNA Networks in Seed Plants

In eudicot plants, the miR482/miR2118 superfamily regulates and instigates the production of phased secondary small interfering RNAs (siRNAs) from NB-LRR (nucleotide binding leucine-rich repeat) genes that encode disease resistance proteins. In grasses, this miRNA family triggers siRNA production specifically in reproductive tissues from long noncoding RNAs. To understand this functional divergence, we examined the small RNA population in the ancient gymnosperm Norway spruce (Picea abies). As many as 41 miRNA families in spruce were found to trigger phasiRNA (phased, secondary siRNAs) production from diverse PHAS loci, with a remarkable 19 miRNA families capable of targeting over 750 NB-LRR genes to generate phasiRNAs. miR482/miR2118, encoded in spruce by at least 24 precursor loci, targets not only NB-LRR genes to trigger phasiRNA production (as in eudicots) but also noncoding PHAS loci, generating phasiRNAs preferentially in male or female cones, reminiscent of its role in the grasses. These data suggest a dual function of miR482/miR2118 present in gymnosperms that was selectively yet divergently retained in flowering plants. A few MIR482/MIR2118 precursors possess an extremely long stem-loop structure, one arm of which shows significant sequence similarity to spruce NB-LRR genes, suggestive of an evolutionary origin from NB-LRR genes through gene duplication. We also characterized an expanded miR390-TAS3 (TRANS-ACTING SIRNA GENE 3)-ARF (AUXIN RESPONSIVE FACTOR) pathway, comprising 18 TAS3 genes of diverse features. Finally, we annotated spruce miRNAs and their targets. Taken together, these data expand our understanding of phasiRNA network in plants and the evolution of plant miRNAs, particularly miR482/miR2118 and its functional diversification.

Composition and Expression of Conserved MicroRNA Genes in Diploid Cotton (Gossypium) Species

MicroRNAs are ubiquitous in plant genomes but vary greatly in their abundance within and conservation among plant lineages. To gain insight into the evolutionary birth/death dynamics of microRNA families, we sequenced small RNA and 5′-end PARE libraries generated from two closely related species of Gossypium. Here, we demonstrate that 33 microRNA families, with similar copy numbers and average evolutionary rates, are conserved in the two congeneric cottons. Analysis of the presence/absence of these microRNA families in other land plants sheds light on their depth of phylogenetic origin and lineage-specific loss/gain. Conserved microRNA families in Gossypium exhibit a striking interspecific asymmetry in expression, potentially connected to relative proximity to neighboring transposable elements. A complex correlated expression pattern of microRNA target genes with their controlling microRNAs indicates that possible functional divergence of conserved microRNA families can also exist even within a single plant genus.

Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j‐3p in soybean nodulation

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation.

A PCR-based marker for a locus conferring the aroma in Myanmar rice (Oryza sativa L.)

Aromatic rice is an important commodity for international trade, which has encouraged the interest of rice breeders to identify the genetic control of rice aroma. The recessive Os2AP gene, which is located on chromosome 8, has been reported to be associated with rice aroma. The 8-bp deletion in exon 7 is an aromatic allele that is present in most aromatic accessions, including the most popular aromatic rice varieties, Jasmine and Basmati. However, other mutations associated with aroma have been detected, but the other mutations are less frequent. In this study, we report an aromatic allele, a 3-bp insertion in exon 13 of Os2AP, as a major allele found in aromatic rice varieties from Myanmar. The insertion is in frame and causes an additional tyrosine (Y) in the amino acid sequence. However, the mutation does not affect the expression of the Os2AP gene. A functional marker for detecting this allele was developed and tested in an aroma-segregating F2 population. The aroma phenotypes and genotypes showed perfect co-segregation of this population. The marker was also used for screening a collection of aromatic rice varieties collected from different geographical sites of Myanmar. Twice as many aromatic Myanmar rice varieties containing the 3-bp insertion allele were found as the varieties containing the 8-bp deletion allele, which suggested that the 3-bp insertion allele originated in regions of Myanmar.

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.Several models have been proposed to explain the emergence of sex chromosomes. Here, through comparative genomics and mutant analysis, Harkess et al. show that linked but separate genes on the Y chromosome are responsible for sex determination in Asparagus, supporting a two-gene model for sex chromosome evolution.