Search

Male sterility enables hybrid crop breeding to increase yields and has been extensively studied. But thermo-sensitive female sterility, which is an ideal property that may enable full mechanization in hybrid rice breeding, has rarely been investigated due to the absence of such germplasm. Here we identify the spontaneous thermo-sensitive female sterility 1 (tfs1) mutation that confers complete sterility under regular/high temperature and partial fertility under low temperature as a point mutation in ARGONAUTE7 (AGO7). AGO7 associates with miR390 to form an RNA-Induced Silencing Complex (RISC), which triggers the biogenesis of small interfering RNAs (siRNAs) from TRANS-ACTING3 (TAS3) loci by recruiting SUPPRESSOR OF GENE SILENCING (SGS3) and RNA-DEPENDENT RNA POLYMERASE6 (RDR6) to TAS3 transcripts. These siRNAs are known as tasiR-ARFs as they act in trans to repress auxin response factor genes. The mutant TFS1 (mTFS1) protein is compromised in its ability to load the miR390/miR390* duplex and eject miR390* during RISC formation. Furthermore, tasiR-ARF levels are reduced in tfs1 due to the deficiency in RDR6 but not SGS3 recruitment by mTFS1 RISC under regular/high temperature, while low temperature partially restores mTFS1 function in RDR6 recruitment and tasiR-ARF biogenesis. A miR390 mutant also exhibits female sterility, suggesting that female fertility is controlled by the miR390-AGO7 module. Notably, the tfs1 allele introduced into various elite rice cultivars endows thermo-sensitive female sterility. Moreover, field trials confirm the utility of tfs1 as a restorer line in fully mechanized hybrid rice breeding.

Ribosome biogenesis is essential for protein synthesis in gene expression. Yeast eIF5B has been shown biochemically to facilitate 18S ribosomal RNA (rRNA) 3′ end maturation during late-stage 40S ribosomal subunit assembly and gate the transition from translation initiation to elongation. But the genome-wide effects of eIF5B have not been studied at the single-nucleotide resolution in any organism, and 18S rRNA 3′ end maturation is poorly understood in plants. Arabidopsis HOT3/eIF5B1 was found to promote development and heat stress acclimation by translational regulation, but its molecular function remained unknown. Here, we show that HOT3 is a late-stage ribosome biogenesis factor that facilitates 18S rRNA 3′ end processing and is a translation initiation factor that globally impacts the transition from initiation to elongation. By developing and implementing 18S-ENDseq, we revealed previously unknown events in 18S rRNA 3′ end maturation or metabolism. We quantitatively defined processing hotspots and identified adenylation as the prevalent nontemplated RNA addition at the 3′ ends of pre-18S rRNAs. Aberrant 18S rRNA maturation in hot3 further activated RNA interference to generate RDR1- and DCL2/4-dependent risiRNAs mainly from a 3′ portion of 18S rRNA. We further showed that risiRNAs in hot3 were predominantly localized in ribosome-free fractions and were not responsible for the 18S rRNA maturation or translation initiation defects in hot3 . Our study uncovered the molecular function of HOT3/eIF5B1 in 18S rRNA maturation at the late 40S assembly stage and revealed the regulatory crosstalk among ribosome biogenesis, messenger RNA (mRNA) translation initiation, and siRNA biogenesis in plants.

In plants, RNA-directed DNA methylation (RdDM) uses small interfering RNAs (siRNAs) to target transposable elements (TEs) but usually avoids genes. RNA polymerase IV (Pol IV) shapes the landscape of DNA methylation through its pivotal role in siRNA biogenesis. However, how Pol IV is recruited to specific loci, particularly how it avoids genes, is poorly understood. Here, we identified a Pol IV–interacting protein, ZMP (zinc finger, mouse double-minute/switching complex B, Plus-3 protein), which exerts a dual role in regulating siRNA biogenesis and DNA methylation at specific genomic regions. ZMP is required for siRNA biogenesis at some pericentromeric regions and prevents Pol IV from targeting a subset of TEs and genes at euchromatic loci. As a chromatin-associated protein, ZMP prefers regions with depleted histone H3 lysine 4 (H3K4) methylation abutted by regions with H3K4 methylation, probably monitoring changes in local H3K4 methylation status to regulate Pol IV’s chromatin occupancy. Our findings uncover a mechanism governing the specificity of RdDM.

MicroRNAs (miRNAs) play an essential role in plant growth and development, and as such, their biogenesis is fine-tuned via regulation of the core microprocessor components. Here, we report that Arabidopsis AAR2, a homolog of a U5 snRNP assembly factor in yeast and humans, not only acts in splicing but also promotes miRNA biogenesis. AAR2 interacts with the microprocessor component hyponastic leaves 1 (HYL1) in the cytoplasm, nucleus, and dicing bodies. In aar2 mutants, abundance of nonphosphorylated HYL1, the active form of HYL1, and the number of HYL1-labeled dicing bodies are reduced. Primary miRNA (pri-miRNA) accumulation is compromised despite normal promoter activities of MIR genes in aar2 mutants. RNA decay assays show that the aar2-1 mutation leads to faster degradation of pri-miRNAs in a HYL1-dependent manner, which reveals a previously unknown and negative role of HYL1 in miRNA biogenesis. Taken together, our findings reveal a dual role of AAR2 in miRNA biogenesis and pre-messenger RNA splicing.

Meiotic recombination is initiated by the SPORULATION 11 (SPO11)–triggered formation of double-strand breaks (DSBs) that usually occur in open chromatin with active transcriptional features in many eukaryotes. However, gene transcription at DSB sites appears to be detrimental for repair, but the regulatory mechanisms governing transcription at meiotic DSB sites are largely undefined in plants. Here, we demonstrate that the largest DNA polymerase epsilon subunit POL2A interacts with SU(VAR)3 to 9 homologs SUVH2 and SUVH9. N-SIM (structured illumination microscopy) observation shows that the colocalization of SUVH2 with the meiotic DSB marker γ-H2AX is dependent on POL2A. RNA-seq of male meiocytes demonstrates that POL2A and SUVH2 jointly repress the expression of 865 genes, which have several known characteristics associated with meiotic DSB sites. Bisulfite-seq and small RNA-seq of male meiocytes support the idea that the silencing of these genes by POL2A and SUVH2/9 is likely independent of CHH methylation or 24-nt siRNA accumulation. Moreover,pol2a suvh2 suvh9triple mutants have more severe defects in meiotic recombination and fertility compared with eitherpol2aorsuvh2 suvh9. Our results not only identify a epigenetic regulatory mechanism for gene silencing in male meiocytes but also reveal roles for DNA polymerase and SUVH2/9 beyond their classic functions in mitosis.

MicroRNAs (miRNAs) are endogenous 20-to 24-nucleotide (nt) noncoding RNAs that play important regulatory roles in many biological processes in eukaryotes. miRNAs modulate the expression of target genes at the post-transcriptional level by transcript cleavage or translational inhibition. Identification of miRNA target genes have been extensively investigated in Arabidopsis and rice, but an in-depth global analysis of miRNA-mediated target regulation is still lacking in maize. Here, we report a transcriptome-wide identification of miRNA targets by analyzing Parallel Analysis of RNA Ends (PARE) datasets derived from nine different tissues at five developmental stages of the maize (Zea mays L.) B73 cultivar. 246 targets corresponding to 60 miRNAs from 25 families were identified, including transcription factors and other genes. In addition, PARE analysis revealed that miRNAs guide specific target transcript cleavage in a tissue-preferential manner. Interestingly, primary transcripts of MIR159c and MIR169e were found to be cleaved by mature miR159 and miR169, respectively, indicating a positive-feedback regulatory mechanism in miRNA biogenesis. Moreover, several new miRNA-target gene pairs involved in seed germination were identified and experimentally validated. Our PARE analyses generated a wide and detailed miRNA-target interaction atlas, which provides a valuable resource for investigating the roles of miRNAs and their targets in maize.

Scientists have developed many approaches based on PCR or next-generation sequencing to localize and characterize integrated T-DNAs in transgenic plants generated by Agrobacterium tumefaciens-mediated T-DNA transfer. However, none of these methods has the robust ability to handle all transgenic plants with diversified T-DNA patterns. Utilizing the valuable information in the whole-genome sequencing data of transgenic plants, we have developed a comprehensive approach (T-LOC) to localize and characterize T-DNA integration sites (TISs). We evaluated the performance of T-LOC on genome sequencing data from 48 transgenic rice (Oryza sativa) plants that provide real and unbiased resources of T-DNA integration patterns. T-LOC discovered 75 full TISs and reported a diversified pattern of T-DNA integration: the ideal single-copy T-DNA between two borders, multiple-copy of T-DNAs in tandem or inverted repeats, truncated partial T-DNAs with or without the selection hygromycin gene, the inclusion of T-DNA backbone, the integration at the genome repeat region, and the concatenation of multiple ideal or partial T-DNAs. In addition, we reported that DNA fragments from the two A. tumefaciens plasmids can be fused with T-DNA and integrated into the plant genome. Besides, T-LOC characterizes the genomic changes at TISs, including deletion, duplication, accurate repair, and chromosomal rearrangement. Moreover, we validated the robustness of T-LOC using PCR, Sanger sequencing, and Nanopore sequencing. In summary, T-LOC is a robust approach to studying the TISs independent of the integration pattern and can recover all types of TISs in transgenic plants.

Meiosis is a specialized cell division that generates gametes and is essential for sexual reproduction. Studying meiosis in plants, like the model flowering plant Arabidopsis thaliana, contributes to our understanding of the fundamental biology of reproductive biology and has practical implications for improving economically important crop species. In this chapter, we provide a detailed protocol for capillary collection of Arabidopsis male meiocytes followed by total RNA extraction, RNA-Seq, and bioinformatics analysis of small-RNAs (sRNAs) including analysis of sRNA cluster that correlate with genomic features.

The Arabidopsis bHLH010/089/091 (basic helix-loop-helix) genes are functionally redundant and are required for both anther development and normal expression of DYT1-activated anther-related genes. These three genes are conserved in Brassicaceae, suggesting that each of them is under selection pressure; however, little is known about the possible functional differences among these bHLH genes and between the bHLH and DYT1 genes. Here, we compared novel anther transcriptomic data sets from bHLH010/089/091 single and double mutants, with an anther transcriptomic data set from the wild type (WT) and a previously obtained anther transcriptomic data set from the bhlh010 bhlh089 bhlh091 triple mutant. The results revealed molecular phenotypes that support the functional redundancy and divergence of bHLH010, bHLH089, and bHLH091, as well as the functional overlap and difference between them and DYT1. DNA-binding analyses revealed that DYT1 and bHLH089 specifically recognize the TCATGTGC box to activate the expression of target genes, including ATA20, EXL4, and MEE48. In addition, among genes whose expression was affected in the bhlh010 bhlh089 double and bhlh010 bhlh089 bhlh091 triple mutants, genes that are involved in the stress response and cell signaling were enriched, which included 256 genes whose expression was preferentially induced by heat during early flower development. Moreover, the bhlh double mutants exhibited defective pollen development when the plants were grown under elevated temperature, suggesting that bHLH genes are important for anther gene expression under such conditions. These results are consistent with the observation that the heat-induced expression of several genes is less in the bhlh mutants than that in the WT. Therefore, our results provide important insights into the molecular mechanism underlying the activation of direct targets by DYT1-bHLH089 heterodimers and demonstrate the protective roles of bHLH010/089/091 in maintaining fertility upon heat stress.

MicroRNAs (miRNAs) play crucial roles in gene expression regulation through RNA cleavage or translation repression. Although miRNA biogenesis has been extensively studied, new factors involved in miRNA biogenesis are still being found from various genetic screens, suggesting that knowledge of miRNA biogenesis or its regulation is still incomplete. Here, we report the identification of a previously uncharacterized and evolutionarily conserved WD40 domain protein as a player in miRNA biogenesis in Arabidopsis thaliana. A mutation in the REDUCTION IN BLEACHED VEIN AREA (RBV) gene encoding a WD40 domain protein led to the suppression of leaf bleaching caused by an artificial miRNA; the mutation also led to a global reduction in the accumulation of endogenous miRNAs. RBV may act at multiple steps in miRNA biogenesis. RBV promotes the transcription of MIR genes into pri-miRNAs by enhancing the occupancy of RNA polymerase II (Pol II) at MIR gene promoters. The nuclear protein RBV may also participate in pri-miRNA processing, as the mutation of RBV leads to the reduction of dicing body number. Moreover, mutation of RBV leads to a defect of miRNA loading into AGO1.RBV function is required in messenger RNA splicing, as RNA-seq uncovered a global intron retention defect in the mutant. Thus, this previously uncharacterized, evolutionarily conserved, nuclear WD40 domain protein acts in miRNA biogenesis and RNA splicing.

Because of climate change, crops will experience increasing heat stress. However, the ways in which heat stress affects crop growth and yield at the molecular level remain poorly understood. We generated spatiotemporal mRNA and small RNA transcriptome data, spanning seven tissues at three time points, to investigate the effects of heat stress on vegetative and reproductive development in maize (Zea mays). Among the small RNAs significantly induced by heat stress was a plastid-derived 19-nucleotide small RNA, which is possibly the residual footprint of a pentatricopeptide repeat protein. This suggests that heat stress induces the turnover of certain plastid transcripts. Consistently, genes responsible for photosynthesis in chloroplasts were repressed after heat stress. Analysis also revealed that the abundance of 24-nucletide small interfering RNAs from transposable elements was conspicuously reduced by heat stress in tassels and roots; nearby genes showed a similar expression trend. Finally, specific microRNA and passenger microRNA species were identified, which in other plant species have not before been reported as responsive to heat stress. This study generated an atlas of genome-wide transcriptomic responses to heat stress, revealing several key regulators as potential targets for thermotolerance improvement in maize.

As the most common RNA cap in eukaryotes, the 7-methylguanosine (m 7 G) cap impacts nearly all processes that a messenger RNA undergoes, such as splicing, polyadenylation, nuclear export, translation, and degradation. The metabolite and redox agent, nicotinamide adenine diphosphate (NAD + ), can be used as an initiating nucleotide in RNA synthesis to result in NAD + -capped RNAs. Such RNAs have been identified in bacteria, yeast, and human cells, but it is not known whether they exist in plant transcriptomes. The functions of the NAD + cap in RNA metabolism or translation are still poorly understood. Here, through NAD captureSeq, we show that NAD + -capped RNAs are widespread in Arabidopsis thaliana . NAD + -capped RNAs are predominantly messenger RNAs encoded by the nuclear and mitochondrial genomes, but not the chloroplast genome. NAD + -capped transcripts from the nuclear genome appear to be spliced and polyadenylated. Furthermore, although NAD + -capped transcripts constitute a small proportion of the total transcript pool from any gene, they are enriched in the polysomal fraction and associate with translating ribosomes. Our findings implicate the existence of as yet unknown mechanisms whereby the RNA NAD + cap interfaces with RNA metabolic processes as well as translation initiation. More importantly, our findings suggest that cellular metabolic and/or redox states may influence, or be regulated by, mRNA NAD + capping.

PROTEIN PHOSPHATASE4 (PP4) is a highly conserved Ser/Thr protein phosphatase found in yeast, plants, and animals. The composition and functions of PP4 in plants are poorly understood. Here, we uncovered the complexity of PP4 composition and function in Arabidopsis (Arabidopsis thaliana) and identified the composition of one form of PP4 containing the regulatory subunit PP4R3A. We show that PP4R3A, together with one of two redundant catalytic subunit genes, PROTEIN PHOSPHATASE X (PPX)1 and PPX2, promotes the biogenesis of microRNAs (miRNAs). PP4R3A is a chromatin-associated protein that interacts with RNA polymerase II and recruits it to the promoters of miRNA-encoding (MIR) genes to promote their transcription. PP4R3A likely also promotes the cotranscriptional processing of miRNA precursors, because it recruits the microprocessor component HYPONASTIC LEAVES1 to MIR genes and to nuclear dicing bodies. Finally, we show that hundreds of introns exhibit splicing defects in pp4r3a mutants. Together, this study reveals roles for Arabidopsis PP4 in transcription and nuclear RNA metabolism.

Nicotinamide adenine diphosphate (NAD(+)) is a novel messenger RNA 5′ cap in Escherichia coli, yeast, mammals, and Arabidopsis. Transcriptome-wide identification of NAD(+)-capped RNAs (NAD-RNAs) was accomplished through NAD captureSeq, which combines chemoenzymatic RNA enrichment with high-throughput sequencing. NAD-RNAs are enzymatically converted to alkyne-RNAs that are then biotinylated using a copper-catalyzed azide–alkyne cycloaddition (CuAAC) reaction. Originally applied to E. coli RNA, which lacks the m(7)G cap, NAD captureSeq was then applied to eukaryotes without extensive verification of its specificity for NAD-RNAs vs. m(7)G-capped RNAs (m(7)G-RNAs). In addition, the Cu(2+) ion in the CuAAC reaction causes RNA fragmentation, leading to greatly reduced yield and loss of full-length sequence information. We developed an NAD-RNA capture scheme utilizing the copper-free, strain-promoted azide–alkyne cycloaddition reaction (SPAAC). We examined the specificity of CuAAC and SPAAC reactions toward NAD-RNAs and m(7)G-RNAs and found that both prefer the former, but also act on the latter. We demonstrated that SPAAC-NAD sequencing (SPAAC-NAD-seq), when combined with immunodepletion of m(7)G-RNAs, enables NAD-RNA identification with accuracy and sensitivity, leading to the discovery of new NAD-RNA profiles in Arabidopsis. Furthermore, SPAAC-NAD-seq retained full-length sequence information. Therefore, SPAAC-NAD-seq would enable specific and efficient discovery of NAD-RNAs in prokaryotes and, when combined with m(7)G-RNA depletion, in eukaryotes.

Nicotinamide adenine diphosphate (NAD+) is a novel messenger RNA 5' cap in Escherichia coli, yeast, mammals, and Arabidopsis Transcriptome-wide identification of NAD+-capped RNAs (NAD-RNAs) was accomplished through NAD captureSeq, which combines chemoenzymatic RNA enrichment with high-throughput sequencing. NAD-RNAs are enzymatically converted to alkyne-RNAs that are then biotinylated using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) reaction. Originally applied to E. coli RNA, which lacks the m7G cap, NAD captureSeq was then applied to eukaryotes without extensive verification of its specificity for NAD-RNAs vs. m7G-capped RNAs (m7G-RNAs). In addition, the Cu2+ ion in the CuAAC reaction causes RNA fragmentation, leading to greatly reduced yield and loss of full-length sequence information. We developed an NAD-RNA capture scheme utilizing the copper-free, strain-promoted azide-alkyne cycloaddition reaction (SPAAC). We examined the specificity of CuAAC and SPAAC reactions toward NAD-RNAs and m7G-RNAs and found that both prefer the former, but also act on the latter. We demonstrated that SPAAC-NAD sequencing (SPAAC-NAD-seq), when combined with immunodepletion of m7G-RNAs, enables NAD-RNA identification with accuracy and sensitivity, leading to the discovery of new NAD-RNA profiles in Arabidopsis Furthermore, SPAAC-NAD-seq retained full-length sequence information. Therefore, SPAAC-NAD-seq would enable specific and efficient discovery of NAD-RNAs in prokaryotes and, when combined with m7G-RNA depletion, in eukaryotes.

Additional file 1: Figure S1–S17.

Summary Unlike in metazoans, the stepwise biogenesis of microRNAs (miRNAs) occurs within the nucleus in plants. Whether or how the major steps in miRNA biogenesis are coordinated is largely unknown. Here, we show that the plant TREX-2 complex promotes multiple steps in miRNA biogenesis, including transcription, processing and nuclear export. Core subunits of TREX-2, THP1 and SAC3A, interact and co-localize with RNA Polymerase II to promote the transcription of MIR genes in the nucleoplasm. TREX-2 interacts with the microprocessor component SERRATE and promotes the formation of Dicing bodies in the nucleoplasm. THP1 also interacts and co-localizes with the nucleoporin protein NUP1 at the nuclear envelope. NUP1 and THP1 promote the nuclear export of miRNAs and AGO1. These results suggest that TREX-2 coordinates the transcription, processing and export steps in miRNA biogenesis to ensure efficient miRNA production.

The evolution of duplicated genes after polyploidization has been the subject of many evolutionary biology studies. Potato (Solanum tuberosum) and tomato (Solanum lycopersicum) are the first two sequenced genomes of asterids, and share a common polyploidization event. However, the epigenetic role of DNA methylation on the evolution of duplicated genes derived from polyploidization is not fully understood. Here, we explore the role of the DNA methylation in the evolution of duplicated genes in potato and tomato. The overall levels of DNA methylation are different, although patterns of DNA methylation are similar in potato and tomato. Different types of duplicated genes can display different methylation patterns in potato and tomato. In addition, we found that differences in the methylation levels between duplicated genes were associated with gene expression divergence. In particular, for the majority of duplicated gene pairs, one copy is always hyper- or hypo-methylated compared with the other copy across different tomato fruit ripening stages, and these genes are enriched for specific function related to transcription factor (TF) activity. Furthermore, transcription of hundreds of duplicated TFs was shown to be regulated by DNA methylation during fruit ripening stages in tomato, some of which are well-known fruit ripening TFs. Taken together, our results support the notion that DNA methylation may facilitate divergent evolution of duplicated genes and play roles in important biological processes such as tomato fruit ripening.

In Arabidopsis thaliana, the MOS4-ASSOCIATED COMPLEX (MAC) is required for defense and development. The evolutionarily conserved, putative RNA helicase MAC7 is a component of the Arabidopsis MAC, and the human MAC7 homolog, Aquarius, is implicated in pre-mRNA splicing. Here, we show that mac7-1, a partial loss-of-function mutant in MAC7, and two other MAC subunit mutants, mac3a mac3b and prl1 prl2 (pleiotropic regulatory locus), exhibit reduced microRNA (miRNA) levels, indicating that MAC promotes miRNA biogenesis. The mac7-1 mutant shows reduced primary miRNA (pri-miRNA) levels without affecting miRNA gene (MIR) promoter activity or the half-life of pri-miRNA transcripts. As a nuclear protein, MAC7 is not concentrated in dicing bodies, but it affects the localization of HYPONASTIC LEAVES1 (HYL1), a key protein in pri-miRNA processing, to dicing bodies. Immunoprecipitation of HYL1 retrieved 11 known MAC subunits, including MAC7, indicating association between HYL1 and MAC. We propose that MAC7 links MIR transcription to pri-miRNA processing. RNA-seq analysis showed that downregulated genes in MAC subunit mutants are mostly involved in plant defense and stimulus responses, confirming a role of MAC in biotic and abiotic stress responses. We also discovered global intron retention defects in mutants in three subunits of MAC, thus linking MAC function to splicing in Arabidopsis.

The CLAVATA3/ESR-RELATED (CLE) peptide signals are required for cell-cell communication in several plant growth and developmental processes. However, little is known regarding the possible functions of the CLEs in the anther. Here, we show that a T-DNA insertional mutant, and dominant-negative (DN) and overexpression (OX) transgenic plants of the CLE19 gene, exhibited significantly reduced anther size and pollen grain number and abnormal pollen wall formation in Arabidopsis (Arabidopsis thaliana). Interestingly, the DN-CLE19 pollen grains showed a more extensively covered surface, but CLE19-OX pollen exine exhibited clearly missing connections in the network and lacked separation between areas that normally form the lacunae. With a combination of cell biological, genetic, and transcriptomic analyses on cle19, DN-CLE19, and CLE19-OX plants, we demonstrated that CLE19-OX plants produced highly vacuolated and swollen aborted microspores (ams)-like tapetal cells, lacked lipidic tapetosomes and elaioplasts, and had abnormal pollen primexine without obvious accumulation of sporopollenin precursors. Moreover, CLE19 is important for the normal expression of more than 1,000 genes, including the transcription factor gene AMS, 280 AMS-downstream genes, and other genes involved in pollen coat and pollen exine formation, lipid metabolism, pollen germination, and hormone metabolism. In addition, the DN-CLE19(+/+) ams(−/−) plants exhibited the ams anther phenotype and ams(+/−) partially suppressed the DN-CLE19 transgene-induced pollen exine defects. These findings demonstrate that the proper amount of CLE19 signal is essential for the normal expression of AMS and its downstream gene networks in the regulation of anther development and pollen exine formation.

MicroRNAs (miRNAs) are a central player in post-transcriptional regulation of gene expression and are involved in numerous biological processes in eukaryotes. Knowledge of the origins and divergence of miRNAs paves the way for a better understanding of the complexity of the regulatory networks that they participate in. The biogenesis, degradation, and regulatory activities of miRNAs are relatively better understood, but the evolutionary history of miRNAs still needs more exploration. Inverted duplication of target genes, random hairpin sequences and small transposable elements constitute three main models that explain the origination of miRNA genes (MIR). Both inter-species and intra-species divergence of miRNAs exhibits functional adaptation and adaptation to changing environments in evolution. Here we summarize recent progress in studies on the evolution of MIR and related genes.

Background One of the main pathways for cadmium (Cd) transfer from the environment to the human body is through the consumption of leafy vegetables, and Brassica leafy crops tend to be Cd hyper-accumulators. But its response strategies to Cd still lack of systematic study. Results To investigate Brassica response strategies to Cd, we identified two cultivars with different Cd translocation efficiencies and performed mutli-transcriptomic sequencing studies under Cd treatments. Certain transporter families exhibited different temporal expression profiles in the two cultivars and may underlie the different Cd translocation efficiencies. Cd induced a drastic reduction of a 22 nt small RNA, the footprint of a pentatricopeptide repeat (PPR) protein on the chloroplast ndhB transcript and the concomitant down-regulation of the ndhB transcript. A global reduction in the expression of PPR genes was found, revealing previously unknown effects of Cd on organellar gene expression. Analyses of microRNAs (miRNAs) and their target genes by small RNA and degradome sequencing not only revealed Cd-induced changes in miRNAs but also implicated the existence of a regulatory cascade involving bra-miR156, its target gene, and bra-miR397 and bra-miR398 in Cd stress responses. Conclusions The present findings help uncover the impact of Cd stress on the transcriptome of B. parachinensis and provide candidate genes and miRNAs for further investigation.

Plant microRNAs (miRNAs) associate with ARGONAUTE1 (AGO1) to direct post-transcriptional gene silencing and regulate numerous biological processes. Although AGO1 predominantly binds miRNAs in vivo, it also associates with endogenous small interfering RNAs (siRNAs). It is unclear whether the miRNA/siRNA balance affects miRNA activities. Here we report that FIERY1 (FRY1), which is involved in 5′−3′ RNA degradation, regulates miRNA abundance and function by suppressing the biogenesis of ribosomal RNA-derived siRNAs (risiRNAs). In mutants of FRY1 and the nuclear 5′−3′ exonuclease genes XRN2 and XRN3, we find that a large number of 21-nt risiRNAs are generated through an endogenous siRNA biogenesis pathway. The production of risiRNAs correlates with pre-rRNA processing defects in these mutants. We also show that these risiRNAs are loaded into AGO1, causing reduced loading of miRNAs. This study reveals a previously unknown link between rRNA processing and miRNA accumulation., Arabidopsis AGO1 binds both microRNA (miRNA) and endogenous small interfering RNA (siRNA). Here the authors show that 5′−3′ RNA degradation suppresses siRNA biogenesis from ribosomal RNA to prevent overloading of AGO1 with ribosomal siRNA and maintain miRNA-directed post translational gene silencing.

As the most common RNA cap in eukaryotes, the 7-methylguanosine (m(7)G) cap impacts nearly all processes that a messenger RNA undergoes, such as splicing, polyadenylation, nuclear export, translation, and degradation. The metabolite and redox agent, nicotinamide adenine diphosphate (NAD(+)), can be used as an initiating nucleotide in RNA synthesis to result in NAD(+)-capped RNAs. Such RNAs have been identified in bacteria, yeast, and human cells, but it is not known whether they exist in plant transcriptomes. The functions of the NAD(+) cap in RNA metabolism or translation are still poorly understood. Here, through NAD captureSeq, we show that NAD(+)-capped RNAs are widespread in Arabidopsis thaliana. NAD(+)-capped RNAs are predominantly messenger RNAs encoded by the nuclear and mitochondrial genomes, but not the chloroplast genome. NAD(+)-capped transcripts from the nuclear genome appear to be spliced and polyadenylated. Furthermore, although NAD(+)-capped transcripts constitute a small proportion of the total transcript pool from any gene, they are enriched in the polysomal fraction and associate with translating ribosomes. Our findings implicate the existence of as yet unknown mechanisms whereby the RNA NAD(+) cap interfaces with RNA metabolic processes as well as translation initiation. More importantly, our findings suggest that cellular metabolic and/or redox states may influence, or be regulated by, mRNA NAD(+) capping.

One of the main pathways for cadmium (Cd) transfer from the environment to the humans is through the consumption of leafy vegetables, and Brassica leafy crops tend to be Cd hyper-accumulators. To investigate Brassica response strategies to Cd, two cultivars with different Cd tolerance were used to perform transcriptomic studies under Cd treatments. Genes related to vacuolar sequestration (such as ABCCs and YSL family genes) and detoxification (such as GST, SAT and APS) were up-regulated in both cultivars but exhibited different temporal expression profiles in leaves of the two cultivars, which may underlie their different Cd tolerance. Analyses of microRNAs (miRNAs) and their target genes by small RNA and degradome sequencing not only revealed Cd-induced changes in miRNAs but also implicated the existence of a regulatory cascade involving two miRNAs, bra-miR156 and bra-miR397, and their corresponding targets, SPLs and LACs, in Cd stress responses. This cascade down-regulates the expression of several LACs particularly in roots under Cd stress, and it is proposed to be a mechanism to maintain root growth under Cd stress through regulating cell wall lignification. In addition, Cd induced a drastic reduction of a 22-nt small RNA, the footprint of a pentatricopeptide repeat (PPR) protein on the chloroplast ndhB transcript and the concomitant down-regulation of the ndhB transcript. A global reduction in the expression of PPR genes was found, revealing previously unknown effects of Cd on organellar gene expression. The present findings help uncover the impact of Cd stress on the transcriptome of B. parachinensis and reveal the strategies used by B. parachinensis in dealing with Cd stress. It also provides candidate genes and miRNAs for further investigation.

Highlight Identification of over 1600 phosphoproteins and analysis of phosphorylation site motifs provide strong support for extensive protein phosphorylation and the potential regulatory roles for such modifications during anther development., As the male reproductive organ of flowering plants, the stamen consists of the anther and filament. Previous studies on stamen development mainly focused on single gene functions by genetic methods or gene expression changes using comparative transcriptomic approaches, especially in model plants such as Arabidopsis thaliana. However, studies on Arabidopsis anther protein expression and post-translational modifications are still lacking. Here we report proteomic and phosphoproteomic studies on developing Arabidopsis anthers at stages 4–7 and 8–12. We identified 3908 high-confidence phosphorylation sites corresponding to 1637 phosphoproteins. Among the 1637 phosphoproteins, 493 were newly identified, with 952 phosphorylation sites. Phosphopeptide enrichment prior to LC-MS analysis facilitated the identification of low-abundance proteins and regulatory proteins, thereby increasing the coverage of proteomic analysis, and facilitated the analysis of more regulatory proteins. Thirty-nine serine and six threonine phosphorylation motifs were uncovered from the anther phosphoproteome and further analysis supports that phosphorylation of casein kinase II, mitogen-activated protein kinases, and 14-3-3 proteins is a key regulatory mechanism in anther development. Phosphorylated residues were preferentially located in variable protein regions among family members, but they were they were conserved across angiosperms in general. Moreover, phosphorylation might reduce activity of reactive oxygen species scavenging enzymes and hamper brassinosteroid signaling in early anther development. Most of the novel phosphoproteins showed tissue-specific expression in the anther according to previous microarray data. This study provides a community resource with information on the abundance and phosphorylation status of thousands of proteins in developing anthers, contributing to understanding post-translational regulatory mechanisms during anther development.

Plants have evolved a large number of transcription factors (TF), which are enriched among duplicate genes, highlighting their roles in complex regulatory networks. The APETALA2/EREBP-like genes constitute a large plant TF family and participate in development and stress responses. To probe the conservation and divergence of AP2/EREBP genes, we analyzed the duplication patterns of this family in Brassicaceae and identified interacting proteins of representative Arabidopsis AP2/EREBP proteins. We found that many AP2/EREBP duplicates generated early in Brassicaceae history were quickly lost, but many others were retained in all tested Brassicaceae species, suggesting early functional divergence followed by persistent conservation. In addition, the sequences of the AP2 domain and exon numbers were highly conserved in rosids. Furthermore, we used 16 A. thaliana AP2/EREBP proteins as baits in yeast screens and identified 1,970 potential AP2/EREBP-interacting proteins, with a small subset of interactions verified in planta. Many AP2 genes also exhibit reduced expression in an anther-defective mutant, providing a possible link to developmental regulation. The putative AP2-interacting proteins participate in many functions in development and stress responses, including photomorphogenesis, flower development, pathogenesis, drought and cold responses, abscisic acid and auxin signaling. Our results present the AP2/EREBP evolution patterns in Brassicaceae, and support a proposed interaction network of AP2/EREBP proteins and their putative interacting proteins for further study.

A family of DEDDh 3′→5′ exonucleases known as Small RNA Degrading Nucleases (SDNs) initiates the turnover of ARGONAUTE1 (AGO1)-bound microRNAs in Arabidopsis by trimming their 3′ ends. Here, we report the crystal structure of Arabidopsis SDN1 (residues 2-300) in complex with a 9 nucleotide single-stranded RNA substrate, revealing that the DEDDh domain forms rigid interactions with the N-terminal domain and binds 4 nucleotides from the 3′ end of the RNA via its catalytic pocket. Structural and biochemical results suggest that the SDN1 C-terminal domain adopts an RNA Recognition Motif (RRM) fold and is critical for substrate binding and enzymatic processivity of SDN1. In addition, SDN1 interacts with the AGO1 PAZ domain in an RNA-independent manner in vitro, enabling it to act on AGO1-bound microRNAs. These extensive structural and biochemical studies may shed light on a common 3′ end trimming mechanism for 3′→5′ exonucleases in the metabolism of small non-coding RNAs., Small RNA degrading nucleases (SDNs) can degrade short RNAs. Here the authors report the crystal structure of Arabidopsis SDN1 in complex with a single-stranded RNA, and provide new insight into 3′ end trimming mechanism of 3′ to 5′ riboexonucleases in the metabolism of various species of small RNAs.

Anther development, particularly around the time of meiosis, is extremely crucial for plant sexual reproduction. Meanwhile, cell-to-cell communication between somatic (especial tapetum) cells and meiocytes are important for both somatic anther development and meiosis. To investigate possible molecular mechanisms modulating protein activities during anther development, we applied high-resolution mass spectrometry-based proteomic and phosphoproteomic analyses for developing rice (Oryza sativa) anthers around the time of meiosis (RAM). In total, we identified 4984 proteins and 3203 phosphoproteins with 8973 unique phosphorylation sites (p-sites). Among those detected here, 1544 phosphoproteins are currently absent in the Plant Protein Phosphorylation DataBase (P3 DB), substantially enriching plant phosphorylation information. Mapman enrichment analysis showed that 'DNA repair','transcription regulation' and 'signaling' related proteins were overrepresented in the phosphorylated proteins. Ten genetically identified rice meiotic proteins were detected to be phosphorylated at a total of 25 p-sites; moreover more than 400 meiotically expressed proteins were revealed to be phosphorylated and their phosphorylation sites were precisely assigned. 163 putative secretory proteins, possibly functioning in cell-to-cell communication, are also phosphorylated. Furthermore, we showed that DNA synthesis, RNA splicing and RNA-directed DNA methylation pathways are extensively affected by phosphorylation. In addition, our data support 46 kinase-substrate pairs predicted by the rice Kinase-Protein Interaction Map, with SnRK1 substrates highly enriched. Taken together, our data revealed extensive protein phosphorylation during anther development, suggesting an important post-translational modification affecting protein activity.

The anther is the male reproductive organ of flowering plants, and the Arabidopsis bHLH transcription factors encoded by DYSFUNCTIONAL TAPETUM1 (DYT1) and ABORTED MICROSPORE (AMS) are required for control of the complex transcriptional networks regulating anther development. Knowledge of the mechanisms by which the bHLH proteins affect this diverse gene expression is quite limited. We examine here three recently duplicated Arabidopsis bHLH genes, bHLH010, bHLH089 and bHLH091, using evolutionary, genetic, morphological and transcriptomic approaches, and uncover their redundant functions in anther development. These three genes are relatively highly expressed in the tapetum of the Arabidopsis anther; single mutants at each of the bHLH010, bHLH089 and bHLH091 loci are developmentally normal, but the various double and triple combinations progressively exhibit increasingly defective anther phenotypes (abnormal tapetum morphology, delayed callose degeneration, and aborted pollen development), indicating their redundant functions in male fertility. Further transcriptomic and molecular analyses suggest that these three proteins act slightly later than DYT1, and also form protein complexes with DYT1, subsequently affecting the correct expression of many DYT1 target genes in the anther development transcriptional network. This study demonstrated that bHLH010, bHLH089 and bHLH091 together are important for the normal transcriptome of the developing Arabidopsis anther, possibly by forming a feed-forward loop with DYT1.

Flower development is a complex process requiring proper spatiotemporal expression of numerous genes. Accumulating evidence indicates that epigenetic mechanisms, including DNA methylation, play essential roles in modulating gene expression. However, few studies have examined the relationship between DNA methylation and floral gene expression on a genomic scale. Here we present detailed analyses of DNA methylomes at single-base resolution for three Arabidopsis floral periods: meristems, early flowers and late flowers. We detected 1.5 million methylcytosines, and estimated the methylation levels for 24 035 genes. We found that many cytosine sites were methylated de novo from the meristem to the early flower stage, and many sites were demethylated from early to late flowers. A comparison of the transcriptome data of the same three periods revealed that the methylation and demethylation processes were correlated with expression changes of >3000 genes, many of which are important for normal flower development. We also found different methylation patterns for three sequence contexts ((m) CG, (m) CHG and (m) CHH) and in different genic regions, potentially with different roles in gene expression.

Background As key regulators of gene expression in eukaryotes, small RNAs have been characterized in many seed plants, and pathways for their biogenesis, degradation, and action have been defined in model angiosperms. However, both small RNAs themselves and small RNA pathways are not well characterized in other land plants such as lycophytes and ferns, preventing a comprehensive evolutionary perspective on small RNAs in land plants. Results Using 25 representatives from major lineages of lycophytes and ferns, most of which lack sequenced genomes, we characterized small RNAs and small RNA pathways in these plants. We identified homologs of DICER-LIKE (DCL), ARGONAUTE (AGO), and other genes involved in small RNA pathways, predicted over 2600 conserved microRNA (miRNA) candidates, and performed phylogenetic analyses on small RNA pathways as well as miRNAs. Pathways underlying miRNA biogenesis, degradation, and activity were established in the common ancestor of land plants, but the 24-nucleotide siRNA pathway that guides DNA methylation is incomplete in sister species of seed plants, especially lycophytes. We show that the functional diversification of key gene families such as DCL and AGO as observed in angiosperms occurred early in land plants followed by parallel expansion of the AGO family in ferns and angiosperms. We uncovered a conserved AGO subfamily absent in angiosperms. Conclusions Our phylogenetic analyses of miRNAs in bryophytes, lycophytes, ferns, and angiosperms refine the time-of-origin for conserved miRNA families as well as small RNA machinery in land plants. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1291-2) contains supplementary material, which is available to authorized users.

The CLAVATA3/ESR-RELATED (CLE) peptide signals are required for cell-cell communication in several plant growth and developmental processes. However, little is known regarding the possible functions of the CLEs in the anther. Here, we show that a T-DNA insertional mutant, and dominant-negative (DN) and overexpression (OX) transgenic plants of the

Small RNAs are central players in RNA silencing, yet their cytoplasmic compartmentalization and the effects it may have on their activities have not been studied at the genomic scale. Here we report that Arabidopsis microRNAs (miRNAs) and small interfering RNAs (siRNAs) are distinctly partitioned between the endoplasmic reticulum (ER) and cytosol. All miRNAs are associated with membrane-bound polysomes (MBPs) as opposed to polysomes in general. The MBP association is functionally linked to a deeply conserved and tightly regulated activity of miRNAs – production of phased siRNAs (phasiRNAs) from select target RNAs. The phasiRNA precursor RNAs, thought to be noncoding, are on MBPs and are occupied by ribosomes in a manner that supports miRNA-triggered phasiRNA production, suggesting that ribosomes on the rough ER impact siRNA biogenesis. This study reveals global patterns of cytoplasmic partitioning of small RNAs and expands the known functions of ribosomes and ER.

Summary APETALA2 (AP2) is best known for its function in the outer two floral whorls, where it specifies the identities of sepals and petals by restricting the expression of AGAMOUS (AG) to the inner two whorls in Arabidopsis thaliana. Here, we describe a role of AP2 in promoting the maintenance of floral stem cell fate, not by repressing AG transcription, but by antagonizing AG activity in the center of the flower.We performed a genetic screen with ag‐10 plants, which exhibit a weak floral determinacy defect, and isolated a mutant with a strong floral determinacy defect. This mutant was found to harbor another mutation in AG and was named ag‐11. We performed a genetic screen in the ag‐11 background to isolate mutations that suppress the floral determinacy defect. Two suppressor mutants were found to harbor mutations in AP2.While AG is known to shut down the expression of the stem cell maintenance gene WUSCHEL (WUS) to terminate floral stem cell fate, AP2 promotes the expression of WUS. AP2 does not repress the transcription of AG in the inner two whorls, but instead counteracts AG activity.

Drought during rice reproductive development results in yield loss. It is important to understand the functions of drought-responsive genes in reproductive tissues for improving rice yield under water-deficit conditions. We show here that MID1 (MYB Important for Drought Response1), encoding a putative R-R-type MYB-like transcription factor, can improve rice yield under drought. MID1 was primarily expressed in root and leaf vascular tissues, with low level in the tapetum, and was induced by drought and other abiotic stresses. Compared with wild type, MID1-overexpressing plants were more tolerant to drought at both vegetative and reproductive stages and produced more grains under water stress. MID1-overexpressing plants exhibited less severe anther defects such as deformed anther locules, abnormal tapetum, degenerated microspores and expanded middle layer, with improved pollen fertility and higher seed setting rate. MID1 was localized to the nucleus and could activate gene expression in yeast, and its homologs were identified in many other plants with high levels sequence similarity. In addition, candidate MID1-regulated genes were analyzed using RNA-seq and qRT-PCR, including genes crucial for stress responses and anther development, with altered expressions in the florets of MID1-overexpressing plants and RNAi lines. Furthermore, MID1 could bind to the promoters of two drought-related genes (Hsp17.0 and CYP707A5) and one anther developmental gene (KAR) according to ChIP-qPCR data. Our findings suggest that MID1 is a transcriptional regulator that promotes rice male development under drought by modulating the expressions of drought-related and anther developmental genes and provide valuable information for crop improvement.