Stem-loop Region Research Articles

A Novel SARS-CoV-2 Viral Sequence Bioinformatic Pipeline Has Found Genetic Evidence That the Viral 3' Untranslated Region (UTR) Is Evolving and Generating Increased Viral Diversity.

An unprecedented amount of SARS-CoV-2 sequencing has been performed, however, novel bioinformatic tools to cope with and process these large datasets is needed. Here, we have devised a bioinformatic pipeline that inputs SARS-CoV-2 genome sequencing in FASTA/FASTQ format and outputs a single Variant Calling Format file that can be processed to obtain variant annotations and perform downstream population genetic testing. As proof of concept, we have analyzed over 229,000 SARS-CoV-2 viral sequences up until November 30, 2020. We have identified over 39,000 variants worldwide with increased polymorphisms, spanning the ORF3a gene as well as the 3′ untranslated (UTR) regions, specifically in the conserved stem loop region of SARS-CoV-2 which is accumulating greater observed viral diversity relative to chance variation. Our analysis pipeline has also discovered the existence of SARS-CoV-2 hypermutation with low frequency (less than in 2% of genomes) likely arising through host immune responses and not due to sequencing errors. Among annotated non-sense variants with a population frequency over 1%, recurrent inactivation of the ORF8 gene was found. This was found to be present in the newly identified B.1.1.7 SARS-CoV-2 lineage that originated in the United Kingdom. Almost all VOC-containing genomes possess one stop codon in ORF8 gene (Q27∗), however, 13% of these genomes also contains another stop codon (K68∗), suggesting that ORF8 loss does not interfere with SARS-CoV-2 spread and may play a role in its increased virulence. We have developed this computational pipeline to assist researchers in the rapid analysis and characterization of SARS-CoV-2 variation.

Crystal structure of human RPP20-RPP25 proteins in complex with the P3 domain of lncRNA RMRP

Human RNase MRP ribonucleoprotein complex is an essential endoribonuclease involved in the processing of ribosomal RNAs, mitochondrial RNAs and certain messenger RNAs. Its RNA subunit RMRP catalyzes the cleavage of substrate RNAs, and the protein components of RNase MRP are required for activity. RMRP mutations are associated with several types of inherited developmental disorders, but the pathogenic mechanism is largely unknown. Recent structural studies shed lights on the catalytic mechanism of yeast RNase MRP and the closely related RNase P; however, the structural and catalytic mechanism of RMRP in human RNase MRP complex remains unclear. Here we report the crystal structure of the P3 domain of RMRP in complex with the RPP20 and RPP25 proteins of human RNase MRP, which shows that the P3 RNA binds to a conserved positively-charged surface of the RPP20-RPP25 heterodimer through its distal stem and internal loop regions. The disease-related mutations of RMRPP3 are mostly located at the protein-RNA interface and are likely to weaken the binding of P3 to RPP20-RPP25. Moreover, the structure reveals a homodimeric organization of the entire RPP20-RPP25-RMRPP3 complex, which might mediate the dimerization of human RNase MRP complex in cells. These findings provide structural clues to the assembly and pathogenesis of human RNase MRP complex and also reveal a tetrameric feature of RPP20-RPP25 evolutionarily conserved with that of the archaeal Alba proteins.

Post-transcriptional regulation of redox homeostasis by the small RNA SHOxi in haloarchaea

ABSTRACT While haloarchaea are highly resistant to oxidative stress, a comprehensive understanding of the processes regulating this remarkable response is lacking. Oxidative stress-responsive small non-coding RNAs (sRNAs) have been reported in the model archaeon, Haloferax volcanii, but targets and mechanisms have not been elucidated. Using a combination of high throughput and reverse molecular genetic approaches, we elucidated the functional role of the most up-regulated intergenic sRNA during oxidative stress in H. volcanii, named S mall RNA in H aloferax Oxi dative Stress (SHOxi). SHOxi was predicted to form a stable secondary structure with a conserved stem-loop region as the potential binding site for trans-targets. NAD-dependent malic enzyme mRNA, identified as a putative target of SHOxi, interacted directly with a putative ‘seed’ region within the predicted stem loop of SHOxi. Malic enzyme catalyzes the oxidative decarboxylation of malate into pyruvate using NAD+ as a cofactor. The destabilization of malic enzyme mRNA, and the decrease in the NAD+/NADH ratio, resulting from the direct RNA–RNA interaction between SHOxi and its trans-target was essential for the survival of H. volcanii to oxidative stress. These findings indicate that SHOxi likely regulates redox homoeostasis during oxidative stress by the post-transcriptional destabilization of malic enzyme mRNA. SHOxi-mediated regulation provides evidence that the fine-tuning of metabolic cofactors could be a core strategy to mitigate damage from oxidative stress and confer resistance. This study is the first to establish the regulatory effects of sRNAs on mRNAs during the oxidative stress response in Archaea.

Nucleocapsid protein preferentially binds the stem-loop of duplex/quadruplex hybrid that unfolds the quadruplex structure

NCp7 protein binds the duplex/quadruplex hybrid structure, which decreases the thermal stability of DNA and unfolds the G-quadruplex structure. Interestingly, the duplex in the stem-loop region is the more favorable binding site of NCp7. The NCp7 binding twists the top G-tetrad, weakens hydrogen bonding and causes K+ ejection, hence disrupting the G4 structure.

Characterization of UVA-Induced Alterations to Transfer RNA Sequences.

Ultraviolet radiation (UVR) adversely affects the integrity of DNA, RNA, and their nucleoside modifications. By employing liquid chromatography–tandem mass spectrometry (LC–MS/MS)-based RNA modification mapping approaches, we identified the transfer RNA (tRNA) regions most vulnerable to photooxidation. Photooxidative damage to the anticodon and variable loop regions was consistently observed in both modified and unmodified sequences of tRNA upon UVA (λ 370 nm) exposure. The extent of oxidative damage measured in terms of oxidized guanosine, however, was higher in unmodified RNA compared to its modified version, suggesting an auxiliary role for nucleoside modifications. The type of oxidation product formed in the anticodon stem–loop region varied with the modification type, status, and whether the tRNA was inside or outside the cell during exposure. Oligonucleotide-based characterization of tRNA following UVA exposure also revealed the presence of novel photoproducts and stable intermediates not observed by nucleoside analysis alone. This approach provides sequence-specific information revealing potential hotspots for UVA-induced damage in tRNAs.

Visualizing a protonated RNA state that modulates microRNA-21 maturation.

MicroRNAs are evolutionarily conserved small, noncoding RNAs that regulate diverse biological processes. Due to their essential regulatory roles, microRNA biogenesis is tightly regulated, where protein factors are often found to interact with specific primary and precursor microRNAs for regulation. Here, using NMR relaxation dispersion spectroscopy and mutagenesis, we reveal that the precursor of oncogenic microRNA-21 exists as a pH-dependent ensemble that spontaneously reshuffles the secondary structure of the entire apical stem-loop region, including the Dicer cleavage site. We show that the alternative excited conformation transiently sequesters the bulged adenine into a noncanonical protonated A+-G mismatch, conferring a substantial enhancement in Dicer processing over its ground conformational state. These results indicate that microRNA maturation efficiency may be encoded in the intrinsic dynamic ensemble of primary and precursor microRNAs, providing a potential means of regulating microRNA biogenesis in response to environmental and cellular stimuli.

Real-time audio and visual display of the Coronavirus genome

BackgroundThis paper describes a web based tool that uses a combination of sonification and an animated display to inquire into the SARS-CoV-2 genome. The audio data is generated in real time from a variety of RNA motifs that are known to be important in the functioning of RNA. Additionally, metadata relating to RNA translation and transcription has been used to shape the auditory and visual displays. Together these tools provide a unique approach to further understand the metabolism of the viral RNA genome. This audio provides a further means to represent the function of the RNA in addition to traditional written and visual approaches.ResultsSonification of the SARS-CoV-2 genomic RNA sequence results in a complex auditory stream composed of up to 12 individual audio tracks. Each auditory motive is derived from the actual RNA sequence or from metadata. This approach has been used to represent transcription or translation of the viral RNA genome. The display highlights the real-time interaction of functional RNA elements. The sonification of codons derived from all three reading frames of the viral RNA sequence in combination with sonified metadata provide the framework for this display. Functional RNA motifs such as transcription regulatory sequences and stem loop regions have also been sonified. Using the tool, audio can be generated in real-time from either genomic or sub-genomic representations of the RNA. Given the large size of the viral genome, a collection of interactive buttons has been provided to navigate to regions of interest, such as cleavage regions in the polyprotein, untranslated regions or each gene. These tools are available through an internet browser and the user can interact with the data display in real time.ConclusionThe auditory display in combination with real-time animation of the process of translation and transcription provide a unique insight into the large body of evidence describing the metabolism of the RNA genome. Furthermore, the tool has been used as an algorithmic based audio generator. These audio tracks can be listened to by the general community without reference to the visual display to encourage further inquiry into the science.

MAC5, an RNA-binding protein, protects pri-miRNAs from SERRATE-dependent exoribonuclease activities

MAC5 is a component of the conserved MOS4-associated complex. It plays critical roles in development and immunity. Here we report that MAC5 is required for microRNA (miRNA) biogenesis. MAC5 interacts with Serrate (SE), which is a core component of the microprocessor that processes primary miRNA transcripts (pri-miRNAs) into miRNAs and binds the stem-loop region of pri-miRNAs. MAC5 is essential for both the efficient processing and the stability of pri-miRNAs. Interestingly, the reduction of pri-miRNA levels in mac5 is partially caused by XRN2/XRN3, the nuclear-localized 5'-to-3' exoribonucleases, and depends on SE. These results reveal that MAC5 plays a dual role in promoting pri-miRNA processing and stability through its interaction with SE and/or pri-miRNAs. This study also uncovers that pri-miRNAs need to be protected from nuclear RNA decay machinery, which is connected to the microprocessor.

MRNA adenosine methylase (MTA) deposits m6A on pri-miRNAs to modulate miRNA biogenesis in Arabidopsis thaliana

In Arabidopsis thaliana, the METTL3 homolog, mRNA adenosine methylase (MTA) introduces N6-methyladenosine (m6A) into various coding and noncoding RNAs of the plant transcriptome. Here, we show that an MTA-deficient mutant (mta) has decreased levels of microRNAs (miRNAs) but accumulates primary miRNA transcripts (pri-miRNAs). Moreover, pri-miRNAs are methylated by MTA, and RNA structure probing analysis reveals a decrease in secondary structure within stem-loop regions of these transcripts in mta mutant plants. We demonstrate interaction between MTA and both RNA Polymerase II and TOUGH (TGH), a plant protein needed for early steps of miRNA biogenesis. Both MTA and TGH are necessary for efficient colocalization of the Microprocessor components Dicer-like 1 (DCL1) and Hyponastic Leaves 1 (HYL1) with RNA Polymerase II. We propose that secondary structure of miRNA precursors induced by their MTA-dependent m6A methylation status, together with direct interactions between MTA and TGH, influence the recruitment of Microprocessor to plant pri-miRNAs. Therefore, the lack of MTA in mta mutant plants disturbs pri-miRNA processing and leads to the decrease in miRNA accumulation. Furthermore, our findings reveal that reduced miR393b levels likely contributes to the impaired auxin response phenotypes of mta mutant plants.

TRBP-Dicer interaction may enhance HIV-1 TAR RNA translation via TAR RNA processing, repressing host-cell apoptosis.

ABSTRACTThe transactivating response (TAR) RNA-binding protein (TRBP) has been identified as a double-stranded RNA (dsRNA)-binding protein, which associates with a stem-loop region known as the TAR element in human immunodeficiency virus-1 (HIV-1). However, TRBP is also known to be an enhancer of RNA silencing, interacting with Dicer, an enzyme that belongs to the RNase III family. Dicer cleaves long dsRNA into small dsRNA fragments called small interfering RNA or microRNA (miRNA) to mediate RNA silencing. During HIV-1 infection, TAR RNA-mediated translation is suppressed by the secondary structure of 5′UTR TAR RNA. However, TRBP binding to TAR RNA relieves its inhibitory action of translation and Dicer processes HIV-1 TAR RNA to generate TAR miRNA. However, whether the interaction between TRBP and Dicer is necessary for TAR RNA translation or TAR miRNA processing remains unclear. In this study, we constructed TRBP mutants that were unable to interact with Dicer by introducing mutations into amino acid residues necessary for the interaction. Furthermore, we established cell lines expressing such TRBP mutants. Then, we revealed that the TRBP–Dicer interaction is essential for both the TAR-containing RNA translation and the TAR miRNA processing in HIV-1.

Ligand screening to pre-miRNA 149 G-quadruplex investigated by molecular dynamics

Using a molecular dynamics approach, the study of the interaction between six different known ligands and a predicted pre-miRNA 149 RNA G-quadruplex (rG4) structure is reported. The stabilization of rG4 structures formed within the pre-miRNA stem-loop regions using small ligands is an attractive anticancer strategy. Particularly, miRNA-149 is upregulated in a variety of cancers such as prostate cancer and is therefore a potential target for drug development. The results show that ligands C8 and PhenDC3 interact with the rG4 structure via stacking interactions with the end G-quartets. Ligands [16]phenN2, [32]phen2N4 and pyridostatin on the other hand bind the loops/groove interface of the rG4 being H-bonding and electrostatic interactions the driving force of the interaction. The C8 precursor, C8-NH2, emphasizes the structural nuances of the rG4 short loops as the lack of a large terminal aromatic moiety produced a mixed stacking-groove binding mode. Overall, this study may help the design of specific ligands for pre-miRNA rG4 towards anticancer therapeutics development.Communicated by Ramaswamy H. Sarma

Beacons Contribute Valuable Empirical Information to Theoretical 3-D Aptamer-Peptide Binding.

DNA aptamers were developed against five different peptides from the known binding regions of anti-Cytomegalovirus and anti-Herpes Simplex Virus-2 antibodies and the aptamers were ranked by relative affinity based on an ELISA-like (ELASA) microplate assay. The secondary structures of the top five highest affinity aptamers were studied for stem-loop commonalities and the most probable peptide binding sites. Two of these stem-loop structures were converted into beacons by addition of TYE 665 dye on the 5' end and Iowa Black quencher on the 3' end. When competed against increasing concentrations of each of the five peptides, only three of the possible ten interactions demonstrated "lights on" fluorescence beacon responses. When modeled by generation of PDB files, after passage through PATCHDOCK and YASARA, two of the aptamer beacon-peptide interactions showed no theoretical evidence of separating the G-C stem-loop region, despite clear empirical evidence of separation of the fluorophore and quencher beyond the Förster distance leading to abundant fluorescence. And in the second beacon's case, YASARA modeling suggested that the beacon was always open despite clear empirical evidence that it was not (no fluorescence response) and only opened in the presence of one of the five peptides. These results are interpreted as a demonstration that 3-dimensional docking software such as PATCHDOCK and YASARA, which are based on rigid receptor-ligand shape complementarity may not reflect the "induced-fit" interactions between aptamers and their cognate targets. Therefore, for the most complete and accurate picture of aptamer-peptide binding, several theoretical and empirical (e.g., beacon fluorescence) analysis methods may be needed.

Phylogenetic patterns and molecular evolution among ‘True citrus fruit trees’ group (Rutaceae family and Aurantioideae subfamily)

Cytoplasmic chloroplast DNAs was explored to establish phylogenetic relationships among ‘True citrus’ group and elucidate its molecular evolution. Conserved group I intron (PQRS) sequence motifs were detected in Citrus genus and showed differences of nucleotide variation between elements and some species. The P8 stem-loop region was the only variable region in trnL intron of Citrus species. The deletion of 6 bp has been detected between sites (40–45) in the partial alignment of P8 stem-loop resulted in a change in conformation of the secondary structures in the ancestral species and other derived species and this could be explained by the evolution and speciation of this genus. The consistency index (CI=1) and the retention index (RI=1) values suggest that the changes are mostly apomorphic despite the absence of homoplasy (HI=0). The phylogenetic tree confirms and shows a branching of C. limon, C. sinensis, C. paradisi, C. tangelo and C. aurantium in a monophyletic clade. This result so far to C. reticulata ruled out the maternal contribution of mandarin gene pool in the elaboration of these secondary Citrus species. The Maximum Likelihood tree indicated that the genus Citrus is quite distant from the related genera Poncirus, Microcitrus, Eremocitrus and Fortunella except only Fortunella margarita species. Here new insight is given about mechanism involving in the evolution of the trnL intron in ‘True Citrus’. From the study it appears that an explicit rejection of the null hypothesis and a scenario of a demography history and recent expansion of ‘True Citrus’ genotypes seem to be retained that may be an important factor influencing the power of genetic diversity in ‘True Citrus’.

The complete mitochondrial genome of Budorcas taxicolor tibetana (Artiodactyla: Bovidae) and comparison with other Caprinae species: Insight into the phylogeny of the genus Budorcas

Takin (Budorcas taxicolor) is one of the most endangered species. However, the taxonomy of takin is still in dispute. Here, the complete mitochondrial genome of Budorcas taxicolor tibetana was reported for the first time, which featured a typical circular molecule of 16,665 bp in length, including 13 PCGs, 22 tRNAs, two rRNAs, and one control region. A + T content was higher than G + C content. All of the genes were encoded on the heavy strand, except for eight tRNAs and ND6 gene. The OL region was 49 bp in length and highly conserved in the synthesis and stem-loop regions, and all of the substitutions and indels were found only in the loop structure. Three types of tandem repeat units, six pairs of hairpin loop structure (TACAT, ATGTA) and six CSBs were predicted in the control region. Our results clearly revealed the systematic status of Budorcas species, and the phylogenetic analyses indicated that Budorcas was closer to Capra and Pseudois, rather than to Ovis, which should be merged into the subfamily of Caprinae.

Discovery of a Branched Peptide That Recognizes the Rev Response Element (RRE) RNA and Blocks HIV-1 Replication.

We synthesized and screened a unique 46 656-member library composed of unnatural amino acids that revealed several hits against RRE IIB RNA. Among the hit peptides identified, peptide 4A5 was found to be selective against competitor RNAs and inhibited HIV-1 Rev-RRE RNA interaction in cell culture in a p24 ELISA assay. Biophysical characterization in a ribonuclease protection assay suggested that 4A5 bound to the stem-loop region in RRE IIB while SHAPE MaP probing with 234 nt RRE RNA indicated additional interaction with secondary Rev binding sites. Taken together, our investigation suggests that HIV replication is inhibited by 4A5 blocking binding of Rev and subsequent multimerization.

Complex Analyses of Short Inverted Repeats in All Sequenced Chloroplast DNAs.

Chloroplasts are key organelles in the management of oxygen in algae and plants and are therefore crucial for all living beings that consume oxygen. Chloroplasts typically contain a circular DNA molecule with nucleus-independent replication and heredity. Using “palindrome analyser” we performed complete analyses of short inverted repeats (S-IRs) in all chloroplast DNAs (cpDNAs) available from the NCBI genome database. Our results provide basic parameters of cpDNAs including comparative information on localization, frequency, and differences in S-IR presence. In a total of 2,565 cpDNA sequences available, the average frequency of S-IRs in cpDNA genomes is 45 S-IRs/per kbp, significantly higher than that found in mitochondrial DNA sequences. The frequency of S-IRs in cpDNAs generally decreased with S-IR length, but not for S-IRs 15, 22, 24, or 27 bp long, which are significantly more abundant than S-IRs with other lengths. These results point to the importance of specific S-IRs in cpDNA genomes. Moreover, comparison by Levenshtein distance of S-IR similarities showed that a limited number of S-IR sequences are shared in the majority of cpDNAs. S-IRs are not located randomly in cpDNAs, but are length-dependently enriched in specific locations, including the repeat region, stem, introns, and tRNA regions. The highest enrichment was found for 12 bp and longer S-IRs in the stem-loop region followed by 12 bp and longer S-IRs located before the repeat region. On the other hand, S-IRs are relatively rare in rRNA sequences and around introns. These data show nonrandom and conserved arrangements of S-IRs in chloroplast genomes.

The use of the hypervariable P8 region of trnL(UAA) intron for identification of orchid species: Evidence from restriction site polymorphism analysis.

The P8 stem-loop region of the trnL intron, which is known to be hypervariable in size with multiple repeat motifs and created difficulties in alignment, is always excluded in phylogenetic as well as barcode analyses. This region was investigated for species discrimination in 98 taxa of orchids belonging to the tribe Vandeae using in silico mapping of restriction site polymorphism. The length of the P8 regions varied from 200 nucleotides in Aerides rosea to 669 nucleotides in Dendrophylax sallei. Forty two taxa had unique lengths, while as many as eight shared a common length of 521 nucleotides. Of the 35 restriction endonucleases producing digestions in the P8 regions, three, viz., AgsI, ApoI and TspDTI turned out to have recognition sites across all the 98 taxa being studied. When their restriction data were combined, 92 taxa could be discriminated leaving three taxon pairs. However, Acampe papillosa and Aeranthes arachnites despite having similar restriction sites differed in their P8 lengths. This is the first report on thorough investigation of the P8 region of trnL intron for search of species specific restriction sites and hence its use as a potential plant DNA barcode.

Fluorometric determination of lead(II) and mercury(II) based on their interaction with a complex formed between graphene oxide and a DNAzyme.

The authors have designed a DNAzyme where graphene oxide (GO) interacts with the ssDNA stem loop region. The DNAzyme strand and substrate strand are hybridized and bind to the surface of GO which act as a signal reporter, while GO act as a strong quencher. The presence of Pb(II) ion disturbs the GO-DNAzyme complex and causes internal cleavage of the DNAzyme complex. On addition of Thioflavin T (ThT) as a quadruplex inducer, fluorescence intensity (best measured at excitation/emission peaks of 425/490nm) is strongly enhanced. Subsequent addition of Hg(II) to ThT/G-quadruplex complex decreases fluorescence because the G-quadruplex is unwinding to form a T-Hg(II)-T dsDNA system. Therefore, the change in fluorescence intensity of ThT is directly correlated to the concentration of Pb(II) and Hg(II). As a result, the assay is highly selective and sensitive. The limits of detection are 96 pM for Pb(II) and 356 pM for Hg(II). Moreover, the method was applied to the detection of the two ions in spiked real samples and gave satisfactory results. Graphical abstract A label free sensitive and selective "on-off" fluorescent assay for detection of Pb(II) and Hg(II) based on graphene oxide -DNAzyme complex with fluorogenic dye thioflavin T. The limits of detection are 96 pM (Pb2+) and 356 pM (Hg2+).

Regulation of a muralytic enzyme-encoding gene by two non-coding RNAs

ABSTRACTNon-coding regulatory RNAs fine-tune gene expression post-transcriptionally. In the streptomycetes, rpfA – encoding a muralytic enzyme required for establishing and exiting dormancy – is flanked by non-coding regulatory RNA elements both upstream (riboswitch) and downstream [antisense small RNA (sRNA)]. In Streptomyces coelicolor, the upstream riboswitch decreases rpfA transcript abundance in response to the second messenger cyclic di-AMP, itself involved in cell wall metabolism and dormancy. There is, however, no obvious expression platform associated with this riboswitch and consequently, its mechanism of action is entirely unknown. Using in vitro transcription assays, we discovered that the rpfA riboswitch promoted premature transcription termination in response to cyclic di-AMP. Through an extensive mutational analysis, we determined that attenuation required ligand binding and involved an unusual extended stem-loop region unique to a subset of rpfA riboswitches in the actinobacteria. At the other end of the rpfA gene, an antisense sRNA, termed Scr3097, is expressed opposite the predicted rpfA terminator. Using northern blotting, we found that Scr3097 accumulation mirrored that of the rpfA mRNA. In liquid culture, we detected Scr3097 exclusively in exponential-phase cells, and in plate-grown culture, we observed the sRNA primarily in differentiating cultures. Using mutational analyses, we found that the sRNA increased rpfA mRNA abundance in cells. Taken together, our work revealed multiple regulatory RNAs controlling rpfA expression in the streptomycetes.

The RNA binding protein EWS is broadly involved in the regulation of pri-miRNA processing in mammalian cells.

The Ewing Sarcoma protein (EWS) is a multifaceted RNA binding protein (RBP) with established roles in transcription, pre-mRNA processing and DNA damage response. By generating high quality EWS–RNA interactome, we uncovered its specific and prevalent interaction with a large subset of primary microRNAs (pri-miRNAs) in mammalian cells. Knockdown of EWS reduced, whereas overexpression enhanced, the expression of its target miRNAs. Biochemical analysis revealed that multiple elements in target pri-miRNAs, including the sequences flanking the stem–loop region, contributed to high affinity EWS binding and sequence swap experiments between target and non-target demonstrated that the flanking sequences provided the specificity for enhanced pri-miRNA processing by the Microprocessor Drosha/DGCR8. Interestingly, while repressing Drosha expression, as reported earlier, we found that EWS was able to enhance the recruitment of Drosha to chromatin. Together, these findings suggest that EWS may positively and negatively regulate miRNA biogenesis via distinct mechanisms, thus providing a new foundation to understand the function of EWS in development and disease.