24-nt siRNAs Research Articles

Dicer-2 mutations in Aedes aegypti cells lead to a diminished antiviral function against Rift Valley fever virus and Bunyamwera virus infection.

Mosquitoes are known to transmit different arthropod-borne viruses belonging to various virus families. The exogenous small interfering RNA pathway plays an important role in the mosquito defence against such virus infections, with Dicer-2 (Dcr2) as one of the key proteins that initiates the cleavage of viral dsRNAs into 21 nt long virus-derived small interfering RNAs. Previous data identified the importance of various motifs in Dcr2 for its small interfering RNA (siRNA)-mediated antiviral activity. However, all these data focus on positive-strand RNA viruses, although negative-strand RNA viruses, like Bunyaviricetes, include several important mosquito-borne viruses. Here, we aim to investigate the importance of different domains of Dcr2 for antiviral activity against viruses of the Bunyaviricetes. For this, we used the Aedes aegypti-derived Dcr2 knock-out cell line Aag2-AF319 to study the importance of the helicase, RNase III and PIWI-Argonaute-Zwille domains of Dcr2 on the antiviral activity of two viruses belonging to different families of the Bunyaviricetes: the Rift Valley fever virus (RVFV) vaccine strain MP12 (Phenuiviridae, Phlebovirus) and the Bunyamwera orthobunyavirus (BUNV; Peribunyaviridae, Orthobunyavirus). All three domains were determined to be critical for the antiviral activity against both RVFV and BUNV. Interestingly, one specific mutation in the helicase domain (KN) did not result in a loss of antiviral activity for RVFV, but for BUNV, despite losing the ability to produce 21 nt siRNAs.

Integrated multi-omics reveals contrasting epigenetic patterns in leaf and root morphogenesis in Aegilops speltoides

Aegilops speltoides, the closest ancestor of the wheat B subgenome, has been well studied genomically. However, the epigenetic landscape of Ae. speltoides and the effects of epigenetics on its growth and development remain poorly understood. Here, we present a comprehensive multi-omics atlas of leaves and roots in Ae. speltoides, encompassing transcriptome, DNA methylation, histone modifications, and small RNA profiling. Divergent DNA methylation levels were detected between leaves and roots, and were associated with differences in accumulated 24-nt siRNAs. DNA methylation changes in promoters and gene bodies showed strong connections with altered expression between leaves and roots. Transcriptional regulatory networks (TRN) reconstructed between leaves and roots were driven by tissue-specific TF families. DNA methylation and histone modification act together as switches that shape root and leaf morphogenesis by modulating the binding of tissue-specific TFs to their target genes. The TRNs in leaves and roots reshaped during wheat polyploidization were associated with alterations in epigenetic modifications. Collectively, these results not only shed light on the critical contribution of epigenetic regulation in the morphogenesis of leaves and roots in Ae. speltoides but also provide new insights for future investigations into the complex interplay of genetic and epigenetic factors in the developmental biology of common wheat.

A novel RAV transcription factor from pear interacts with viral RNA-silencing suppressors to inhibit viral infection.

In plants, RNA silencing constitutes a strong defense against viral infection, which viruses counteract with RNA-silencing suppressors (RSSs). Understanding the interactions between viral RSSs and host factors is crucial for elucidating the molecular arms race between viruses and host plants. We report that the helicase motif (Hel) of the replicase encoded by apple stem grooving virus (ASGV)-the main virus affecting pear trees in China-is an RSS that can inhibit both local and systemic RNA silencing, possibly by binding double-stranded (ds) siRNA. The transcription factor related to ABSCISIC ACID INSENSITIVE3/VIVIPAROUS1 from pear (PbRAV1) enters the cytoplasm and binds Hel through its C terminus, thereby attenuating its RSS activity by reducing its binding affinity to 21- and 24-nt ds siRNA, and suppressing ASGV infection. PbRAV1 can also target p24, an RSS encoded by grapevine leafroll-associated virus 2 (GLRaV-2), with similar negative effects on p24's suppressive function and inhibition of GLRaV-2 infection. Moreover, like the positive role of the PbRAV1 homolog from grapevine (VvRAV1) in p24's previously reported RSS activity, ASGV Hel can also hijack VvRAV1 and employ the protein to sequester 21-nt ds siRNA, thereby enhancing its own RSS activity and promoting ASGV infection. Furthermore, PbRAV1 neither interacts with CP, an RSS encoded by grapevine inner necrosis virus, nor has any obvious effect on CP's RSS activity. Our results identify an RSS encoded by ASGV and demonstrate that PbRAV1, representing a novel type of RAV transcription factor, plays a defensive role against viral infection by targeting viral RSSs.

Endogenous viral elements are targeted by RNA silencing pathways in banana.

Endogenous banana streak virus (eBSV) integrants derived from three distinct species, present in Musa balbisiana (B) but not Musa acuminata (A) banana genomes are able to reconstitute functional episomal viruses causing banana streak disease in interspecific triploid AAB banana hybrids but not in the diploid (BB) parent line, which harbours identical eBSV loci. Here, we investigated the regulation of these eBSV. In-depth characterization of siRNAs, transcripts and methylation derived from eBSV using Illumina and bisulfite sequencing were carried out on eBSV-free Musa acuminata AAA plants and BB or AAB banana plants with eBSV. eBSV loci produce low-abundance transcripts covering most of the viral sequence and generate predominantly 24-nt siRNAs. siRNA accumulation is restricted to duplicated and inverted viral sequences present in eBSV. Both siRNA-accumulating and nonaccumulating sequences of eBSV in BB plants are heavily methylated in all three CG, CHG and CHH contexts. Our data suggest that eBSVs are controlled at the epigenetic level in BB diploids. This regulation not only prevents their awakening and systemic infection of the plant but is also probably involved in the inherent resistance of the BB plants to mealybug-transmitted viral infection. These findings are thus of relevance to other plant resources hosting integrated viruses.

Asymmetric bulges within hairpin RNA transgenes influence small RNA size, secondary siRNA production and viral defence.

Small RNAs (sRNAs) are essential for normal plant development and range in size classes of 21-24 nucleotides. The 22nt small interfering RNAs (siRNAs) and miRNAs are processed by Dicer-like 2 (DCL2) and DCL1 respectively and can initiate secondary siRNA production from the target transcript. 22nt siRNAs are under-represented due to competition between DCL2 and DCL4, while only a small number of 22nt miRNAs exist. Here we produce abundant 22nt siRNAs and other siRNA size classes using long hairpin RNA (hpRNA) transgenes. By introducing asymmetric bulges into the antisense strand of hpRNA, we shifted the dominant siRNA size class from 21nt of the traditional hpRNA to 22, 23 and 24nt of the asymmetric hpRNAs. The asymmetric hpRNAs effectively silenced a β-glucuronidase (GUS) reporter transgene and the endogenous ethylene insensitive-2 (EIN2) and chalcone synthase (CHS) genes. Furthermore, plants containing the asymmetric hpRNA transgenes showed increased amounts of 21nt siRNAs downstream of the hpRNA target site compared to plants with the traditional hpRNA transgenes. This indicates that these asymmetric hpRNAs are more effective at inducing secondary siRNA production to amplify silencing signals. The 22nt asymmetric hpRNA constructs enhanced virus resistance in plants compared to the traditional hpRNA constructs.

ARGONAUTE1-binding Tudor domain proteins function in small interfering RNA production for RNA-directed DNA methylation.

Transposable elements (TEs) contribute to plant evolution, development, and adaptation to environmental changes, but the regulatory mechanisms are largely unknown. RNA-directed DNA methylation (RdDM) is 1 TE regulatory mechanism in plants. Here, we identified that novel ARGONAUTE 1 (AGO1)-binding Tudor domain proteins Precocious dissociation of sisters C/E (PDS5C/E) are involved in 24-nt siRNA production to establish RdDM on TEs in Arabidopsis thaliana. PDS5 family proteins are subunits of the eukaryote-conserved cohesin complex. However, the double mutant lacking angiosperm-specific subfamily PDS5C and PDS5E (pds5c/e) exhibited different developmental phenotypes and transcriptome compared with those of the double mutant lacking eukaryote-conserved subfamily PDS5A and PDS5B (pds5a/b), suggesting that the angiosperm-specific PDS5C/E subfamily has a unique function in angiosperm plants. Proteome and imaging analyses revealed that PDS5C/E interact with AGO1. The pds5c/e double mutant had defects in 24-nt siRNA accumulation and CHH DNA methylation on TEs. In addition, some lncRNAs that accumulated in the pds5c/e mutant were targeted by AGO1-loading 21-nt miRNAs and 21-nt siRNAs. These results indicate that PDS5C/E and AGO1 participate in 24-nt siRNA production for RdDM in the cytoplasm. These findings indicate that angiosperm plants evolved a new regulator, the PDS5C/E subfamily, to control the increase in TEs during angiosperm evolution.

Transcriptome and small RNAome profiling uncovers how a recombinant begomovirus evades RDRγ-mediated silencing of viral genes and outcompetes its parental virus in mixed infection.

Tomato yellow leaf curl virus (TYLCV, genus Begomovirus, family Geminiviridae) causes severe disease of cultivated tomatoes. Geminiviruses replicate circular single-stranded genomic DNA via rolling-circle and recombination-dependent mechanisms, frequently generating recombinants in mixed infections. Circular double-stranded intermediates of replication also serve as templates for Pol II bidirectional transcription. IS76, a recombinant derivative of TYLCV with a short sequence in the bidirectional promoter/origin-of-replication region acquired from a related begomovirus, outcompetes TYLCV in mixed infection and breaks disease resistance in tomato Ty-1 cultivars. Ty-1 encodes a γ-clade RNA-dependent RNA polymerase (RDRγ) implicated in Dicer-like (DCL)-mediated biogenesis of small interfering (si)RNAs directing gene silencing. Here, we profiled transcriptome and small RNAome of Ty-1 resistant and control susceptible plants infected with TYLCV, IS76 or their combination at early and late infection stages. We found that RDRγ boosts production rates of 21, 22 and 24 nt siRNAs from entire genomes of both viruses and modulates DCL activities in favour of 22 and 24 nt siRNAs. Compared to parental TYLCV, IS76 undergoes faster transition to the infection stage favouring rightward transcription of silencing suppressor and coat protein genes, thereby evading RDRγ activity and facilitating its DNA accumulation in both single and mixed infections. In coinfected Ty-1 plants, IS76 efficiently competes for host replication and transcription machineries, thereby impairing TYLCV replication and transcription and forcing its elimination associated with further increased siRNA production. RDRγ is constitutively overexpressed in Ty-1 plants, which correlates with begomovirus resistance, while siRNA-generating DCLs (DCL2b/d, DCL3, DCL4) and genes implicated in siRNA amplification (α-clade RDR1) and function (Argonaute2) are upregulated to similar levels in TYLCV- and IS76-infected susceptible plants. Collectively, IS76 recombination facilitates replication and promotes expression of silencing suppressor and coat proteins, which allows the recombinant virus to evade the negative impact of RDRγ-boosted production of viral siRNAs directing transcriptional and posttranscriptional silencing.

Pigmentation of soybean seed coats via a mutation that abolishes production of multiple-phased siRNAs of chalcone synthase genes.

Lack of pigmentation in seed coats of soybean is caused by natural RNA silencing of chalcone synthase (CHS) genes. This phenomenon is an evolutionary consequence of structural changes in DNA that resulted in the production of double-stranded RNAs (dsRNAs) that trigger RNA degradation. Here we determined that a mutant with pigmented seed coats derived from a cultivar that lacked the pigmentation had a deletion between DNA regions ICHS1 and a cytochrome P450 gene; the deletion included GmIRCHS, a candidate gene that triggers CHS RNA silencing via production of CHS dsRNAs. We also characterized CHS short interfering RNAs (siRNAs) produced in the wild-type seed coats that had CHS RNA silencing. Phased 21-nt CHS siRNAs were detected in all 21 phases and were widely distributed in exon 2 of CHS7, which indicates commonality in the pattern of RNA degradation in natural CHS RNA silencing between distantly related species. These results with the similarities in the rearrangements found in spontaneous mutants suggest that the structural organization that generates dsRNAs that trigger phased siRNA production is vulnerable to further structural changes, which eventually abolish the induction of RNA silencing.

Identification and Expressional Analysis of siRNAs Responsive to Fusarium graminearum Infection in Wheat

The outbreak of Fusarium head blight (FHB) poses a serious threat to wheat production as it leads to both significant yield losses and accumulation of several mycotoxins including deoxynivalenol (DON) in the grains, which are harmful to human and livestock. To date, hundreds of FHB-resistance-related quantitative trait loci (QTLs) have been reported, but only a few of them have been cloned and used for breeding. Small interfering RNAs (siRNA) have been reported in plants to mediate host defense against pathogens, but they have rarely been reported in wheat-FHB interaction. In order to identify the key siRNAs that can potentially be used in the improvement of resistance to FHB, siRNAs from the spikes of an FHB-resistant variety Sumai 3 and an FHB-susceptible variety of Chinese Spring (CS) were sequenced after F. graminearum infection and mock inoculation, respectively. The expression patterns of the siRNAs of interest were analyzed. A total of 4019 siRNAs of high-confidence were identified, with 131 being CS-specific, 309 Sumai 3-specific and 3071 being common in both varieties. More than 87% of these siRNAs were 24 nt in length. An overall down-regulation trend was found for siRNAs in the spikes of both varieties after being infected with F. graminearum. The expression patterns for Triticum aestivum Dicer-like 3 (TaDCL3) that synthesizes 24 nt siRNAs were validated by qRT-PCR, which were positively correlated with those of the siRNAs. A total of 85% of the differentially expressed genes putatively targeted by the siRNAs were significantly up-regulated after infection, showing a negative correlation with the overall down-regulated expression of siRNAs. Interestingly, the majority of the up-regulated genes are annotated as disease resistance. These results suggested that the inhibition of siRNA by F. graminearum up-regulated the disease resistance genes, which were putatively suppressed by siRNAs through RNA-directed DNA methylation (RdDM). Consequently, the resistant capability to F. graminearum infection was enhanced. This study provides novel clues for investigating the function of siRNA in wheat-F. graminearum interaction.

A transposable element-derived siRNAs involve DNA hypermethylation at the promoter of OsGSTZ4 for cadmium tolerance in rice

Environmental contaminants such as cadmium (Cd) pose high risks to crop production and human health. The genetic basis for regulation of Cd stress-responsive genes for plant adaptation to adverse environments remains poorly understood. In this study, we characterized a rice Zeta family glutathione-S-transferase (OsGSTZ4) gene for Cd detoxification. Heterologous expression of OsGSTZ4 in a yeast (Saccharomyces cerevisiae) conferred cellular Cd tolerance. Transgenic rice overexpressing OsGSTZ4 improved plant growth, attenuated Cd-induced toxicity, and accumulated more Cd in roots. OsGSTZ4 transcription was rapidly induced 3 h after Cd exposure and then declined to the basal level. This was followed by (days after Cd treatment) by CHH hypermethylation (by 41.2 %) at a MITE (Miniature Inverted-repeat Transposable Element) transposable element (TE) inserted in the 5′-untranscribed region (UTR) (-1,722 ∼ -1,392 bp) of OsGSTZ4. Meanwhile, three 24-nt siRNAs derived from the TE (-1,722 ∼ -1,471 bp) were detected and was also rapidly enriched under Cd stress. To validate the possibility that Cd-induced change in OsGSTZ4 expression correlates with the siRNAs-involved CHH methylation through an RdDM (RNA-directed DNA methylation) pathway, genetic analyses were performed. We found that the CHH methylation at the promoter and transcript level of OsGSTZ4 were compromised in the osdrm2 (loss of function for CHH methylation) and osrdr2i (defective in RNA-dependent RNA polymerase 2) but did not change in other types of methyltransferases such as osmet1, ossdg714 or osros1. Promoter deletion analyses confirmed that the siRNA target sequences were essential for the proper expression of OsGSTZ4. Our studies reveal an unusual feedback mechanism by which the Cd-induced rapid OsGSTZ4 expression for Cd tolerance would interplay with the late CHH hypermethylation to silence the TE through the 24-nt siRNAs- and Osdrm2-mediated RdDM pathway, and help understand the diversity of gene regulation via an epigenetic mechanism for rice adaptation to the environmental stress.

Trans-grafting plum pox virus resistance from transgenic plum rootstocks to apricot scions.

Trans-grafting could be a strategy to transfer virus resistance from a transgenic rootstock to a wild type scion. However contradictory results have been obtained in herbaceous and woody plants. This work was intended to determine if the resistance to sharka could be transferred from transgenic plum rootstocks to wild-type apricot scions grafted onto them. To this end, we conducted grafting experiments of wild- type apricots onto plum plants transformed with a construction codifying a hairpin RNA designed to silence the PPV virus and studied if the resistance was transmitted from the rootstock to the scion. Our data support that the RNA-silencing-based PPV resistance can be transmitted from PPV-resistant plum rootstocks to non-transgenic apricot scions and that its efficiency is augmented after successive growth cycles. PPV resistance conferred by the rootstocks was robust, already occurring within the same growing cycle and maintained in successive evaluation cycles. The RNA silencing mechanism reduces the relative accumulation of the virus progressively eliminating the virus from the wild type scions grafted on the transgenic resistant PPV plants. There was a preferential accumulation of the 24nt siRNAs in the scions grafted onto resistant rootstocks that was not found in the scions grafted on the susceptible rootstock. This matched with a significantly lower relative accumulation of hpRNA in the resistant rootstocks compared with the susceptible or the tolerant ones. Using transgenic rootstocks should mitigate public concerns about transgenes dispersion and eating transgenic food and allow conferring virus resistance to recalcitrant to transformation cultivars or species.

RETINOBLASTOMA-RELATED interactions with key factors of the RNA-directed DNA methylation (RdDM) pathway and its influence on root development.

Our study presents evidence for a novel mechanism for RBR function in transcriptional gene silencing by interacting with key players of the RdDM pathway in Arabidopsis and several plant clades. Transposable elements and other repetitive elements are silenced by the RNA-directed DNA methylation pathway (RdDM). In RdDM, POLIV-derived transcripts are converted into double-stranded RNA (dsRNA) by the activity of RDR2 and subsequently processed into 24 nucleotide short interfering RNAs (24-nt siRNAs) by DCL3. 24-nt siRNAs serve as guides to direct AGO4-siRNA complexes to chromatin-bound POLV-derived transcripts generated from the template/target DNA. The interaction between POLV, AGO4, DMS3, DRD1, RDM1 and DRM2 promotes DRM2-mediated de novo DNA methylation. The Arabidopsis Retinoblastoma protein homolog (RBR) is a master regulator of the cell cycle, stem cell maintenance, and development. We in silico predicted and explored experimentally the protein-protein interactions (PPIs) between RBR and members of the RdDM pathway. We found that the largest subunits of POLIV and POLV (NRPD1 and NRPE1), the shared second largest subunit of POLIV and POLV (NRPD/E2), RDR1, RDR2, DCL3, DRM2, and SUVR2 contain canonical and non-canonical RBR binding motifs and several of them are conserved since algae and bryophytes. We validated experimentally PPIs between Arabidopsis RBR and several of the RdDM pathway proteins. Moreover, seedlings from loss-of-function mutants in RdDM and RBR show similar phenotypes in the root apical meristem. We show that RdDM and SUVR2 targets are up-regulated in the 35S:AmiGO-RBR background.

Dicer-like 3a mediates intergenerational resistance against root-knot nematodes in rice via hormone responses.

In a continuously changing and challenging environment, passing down the memory of encountered stress factors to offspring could provide an evolutionary advantage. In this study, we demonstrate the existence of 'intergenerational acquired resistance' in the progeny of rice (Oryza sativa) plants attacked by the belowground parasitic nematode Meloidogyne graminicola. Transcriptome analyses revealed that genes involved in defense pathways are generally downregulated in progeny of nematode-infected plants under uninfected conditions but show a stronger induction upon nematode infection. This phenomenon was termed "spring loading" and depends on initial downregulation by the 24nt siRNA biogenesis gene dicer-like 3a (dcl3a) involved in the RNA-directed DNA methylation pathway. Knock-down of dcl3a led to increased nematode susceptibility and abolished intergenerational acquired resistance, as well as jasmonic acid/ethylene spring loading in the offspring of infected plants. The importance of ethylene signaling in intergenerational resistance was confirmed by experiments on a knock-down line of ethylene insensitive 2 (ein2b), which lacks intergenerational acquired resistance. Taken together, these data indicate a role for DCL3a in regulating plant defense pathways during both within-generation and intergenerational resistance against nematodes in rice.

Antiviral immune response reveals host-specific virus infections in natural ant populations.

Hosts can carry many viruses in their bodies, but not all of them cause disease. We studied ants as a social host to determine both their overall viral repertoire and the subset of actively infecting viruses across natural populations of three subfamilies: the Argentine ant (Linepithema humile, Dolichoderinae), the invasive garden ant (Lasius neglectus, Formicinae) and the red ant (Myrmica rubra, Myrmicinae). We used a dual sequencing strategy to reconstruct complete virus genomes by RNA-seq and to simultaneously determine the small interfering RNAs (siRNAs) by small RNA sequencing (sRNA-seq), which constitute the host antiviral RNAi immune response. This approach led to the discovery of 41 novel viruses in ants and revealed a host ant-specific RNAi response (21 vs. 22 nt siRNAs) in the different ant species. The efficiency of the RNAi response (sRNA/RNA read count ratio) depended on the virus and the respective ant species, but not its population. Overall, we found the highest virus abundance and diversity per population in Li. humile, followed by La. neglectus and M. rubra. Argentine ants also shared a high proportion of viruses between populations, whilst overlap was nearly absent in M. rubra. Only one of the 59 viruses was found to infect two of the ant species as hosts, revealing high host-specificity in active infections. In contrast, six viruses actively infected one ant species, but were found as contaminants only in the others. Disentangling spillover of disease-causing infection from non-infecting contamination across species is providing relevant information for disease ecology and ecosystem management.

Enzymatic reactions of AGO4 in RNA-directed DNA methylation: siRNA duplex loading, passenger strand elimination, target RNA slicing, and sliced target retention.

RNA-directed DNA methylation in plants is guided by 24-nt siRNAs generated in parallel with 23-nt RNAs of unknown function. We show that 23-nt RNAs function as passenger strands during 24-nt siRNA incorporation into AGO4. The 23-nt RNAs are then sliced into 11- and 12-nt fragments, with 12-nt fragments remaining associated with AGO4. Slicing recapitulated with recombinant AGO4 and synthetic RNAs reveals that siRNAs of 21-24 nt, with any 5'-terminal nucleotide, can guide slicing, with sliced RNAs then retained by AGO4. In vivo, RdDM target locus RNAs that copurify with AGO4 also display a sequence signature of slicing. Comparing plants expressing slicing-competent versus slicing-defective AGO4 shows that slicing elevates cytosine methylation levels at virtually all RdDM loci. We propose that siRNA passenger strand elimination and AGO4 tethering to sliced target RNAs are distinct modes by which AGO4 slicing enhances RNA-directed DNA methylation.

Semaphorin 3A-Neuropilin-1 Signaling Modulates MMP13 Expression in Human Osteoarthritic Chondrocytes.

Osteoarthritis (OA) is a complex disorder of diarthrodial joints caused by multiple risk factors and is characterized by articular cartilage destruction as well as changes in other articular tissues. Semaphorin 3A (Sema3A), known to be a chemo-repellent for sensory nerve fibers, has recently been implicated in cartilage OA pathophysiology. We demonstrated that the expression of SEMA3A and its receptor neuropilin-1 (NRP1) are synchronously upregulated in chondrocytes isolated from knee cartilage of OA patients compared to non-OA control chondrocytes. In addition, we observed that during in vitro passaging of OA chondrocytes, the Nrp-1 level increases, whereas the Sema3A level decreases. In this study, we aimed to uncover how Sema3A-Nrp-1 signaling affects metabolism and viability of OA chondrocytes via siRNA-mediated inhibition of Nrp-1 expression. We observed a decreased proliferation rate and an increase in adhesion and senescence after Nrp-1 silencing. Moreover, MMP13 gene expression was reduced by approximately 75% in NRP1 knockdown OA chondrocytes, whereas MMP13 expression was induced by Sema3A treatment in control (nt siRNA) OA chondrocytes, accompanied by an impaired AKT phosphorylation. These findings suggest a potential catabolic function of Sema3A signaling in OA chondrocytes by inducing MMP13 expression and by compromising pro-survival AKT activation. We propose that targeting the Sema3A-Nrp-1 signaling axis might be an opportunity to interfere with OA pathogenesis and progression.

The dicing activity of DCL3 and DCL4 is negatively affected by flavonoids.

The dicing activities of DCL3 and DCL4 are inhibited by accumulated metabolites in soybean leaves. Epicatechin and 7,4'-dihydroxyflavone inhibited Arabidopsis DCL3 and DCL4 in vitro. Flavonoids are major secondary metabolites in plants, and soybean (Glycine max L.) is a representative plant that accumulates flavonoids, including isoflavonoids, to high levels. Naturally-occurring RNA interference (RNAi) against the chalcone synthase (CHS) gene represses flavonoid (anthocyanin) biosynthesis in an organ-specific manner, resulting in a colorless (yellow) seed coat in many soybean cultivars. To better understand seed coat-specific naturally-occurring RNAi in soybean, we characterized soybean Dicer-like (DCL) 3 and 4, which play critical roles in RNAi. Using a previously established dicing assay, two dicing activities producing 24- and 21-nt siRNAs, corresponding to DCL3 and DCL4, respectively, were detected in soybean. Dicing activity was detected in colorless seed coats where RNAi against CHS genes was found, but no dicing activity was detected in leaves where CHS expression was prevalent. Biochemical analysis revealed that soybean leaves contained two types of inhibitors effective for Arabidopsis Dicers (AtDCL3 and AtDCL4), one of which was a heat-labile high molecular weight compound of 50 to 100 kD while another was a low molecular weight substance. We found that some flavonoids, such as epicatechin and 7,4'-dihydroxyflavone, inhibited both AtDCL3 and AtDCL4, but AtDCL4 was more sensitive to these flavonoids than AtDCL3. These results suggest that flavonoids inhibit the dicing activity of DCL4 and thereby attenuate RNAi in soybean leaves.

DNA polymerase epsilon interacts with SUVH2/9 to repress the expression of genes associated with meiotic DSB hotspot in Arabidopsis

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 suvh9 triple mutants have more severe defects in meiotic recombination and fertility compared with either pol2a or suvh2 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.

Cheminformatics Modeling of Gene Silencing for Both Natural and Chemically Modified siRNAs.

In designing effective siRNAs for a specific mRNA target, it is critically important to have predictive models for the potency of siRNAs. None of the published methods characterized the chemical structures of individual nucleotides constituting a siRNA molecule; therefore, they cannot predict the potency of gene silencing by chemically modified siRNAs (cm-siRNA). We propose a new approach that can predict the potency of gene silencing by cm-siRNAs, which characterizes each nucleotide (NT) using 12 BCUT cheminformatics descriptors describing its charge distribution, hydrophobic and polar properties. Thus, a 21-NT siRNA molecule is described by 252 descriptors resulting from concatenating all the BCUT values of its composing nucleotides. Partial Least Square is employed to develop statistical models. The Huesken data (2431 natural siRNA molecules) were used to perform model building and evaluation for natural siRNAs. Our results were comparable with or superior to those from Huesken’s algorithm. The Bramsen dataset (48 cm-siRNAs) was used to build and test the models for cm-siRNAs. The predictive r2 of the resulting models reached 0.65 (or Pearson r values of 0.82). Thus, this new method can be used to successfully model gene silencing potency by both natural and chemically modified siRNA molecules.

Low cell number ChIP-seq reveals chromatin state-based regulation of gene transcription in the rice male meiocytes.

Low cell number ChIP-seq reveals chromatin state-based regulation of gene transcription in the rice male meiocytes.