Nascent Transcripts Research Articles

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics.

High-throughput detection of nascent RNA is critical for studies of transcription and much more challenging than that of mRNA. Recently, several massively parallel nascent RNA sequencing methods were established in eukaryotic cells. Here, we systematically compared 3 classes of methods on the same pure or crude nuclei preparations: GRO-seq for sequence nuclear run-on RNAs, pNET-seq for sequence RNA polymerase II-associated RNAs, and CB RNA-seq for sequence chromatin-bound (CB) RNAs in Arabidopsis (Arabidopsis thaliana). To improve the resolution of CB RNAs, 3'CB RNA-seq was established to sequence the 3' ends of CB RNAs. In addition, we modified pNET-seq to establish the Chromatin Native Elongation Transcript sequencing (ChrNET) method using chromatin as the starting material for RNA immunoprecipitation. Reproducibility, sensitivity and accuracy in detecting nascent transcripts, experimental procedures, and costs were analyzed, which revealed the strengths and weaknesses of each method. We found that pNET and GRO methods best detected active RNA polymerase II. CB RNA-seq is a simple and cost-effective alternative for nascent RNA studies, due to its high correlation with pNET-seq and GRO-seq. Compared with pNET, ChrNET has higher specificity for nascent RNA capture and lower sequencing cost. 3'CB is sensitive to transcription-coupled splicing. Using these methods, we identified 1,404 unknown transcripts, 4,482 unannotated splicing events, and 60 potential recursive splicing events. This comprehensive comparison of different nascent/chromatin RNA sequencing methods highlights the strengths of each method and serves as a guide for researchers aiming to select a method that best meets their study goals.

Binding of the nuclear ribonucleoprotein family member FUS to RNA prevents R-loop RNA:DNA hybrid structures

The protein FUS (FUSed in sarcoma) is a metazoan RNA-binding protein that influences RNA production by all three nuclear polymerases. FUS also binds nascent transcripts, RNA processing factors, RNA polymerases, and transcription machinery. Here, we explored the role of FUS binding interactions for activity during transcription. Invitro run-off transcription assays revealed FUS-enhanced RNA produced by a non-eukaryote polymerase. The activity also reduced the formation of R-loops between RNA products and their DNA template. Analysis by domain mutation and deletion indicated RNA-binding was required for activity. We interpret that FUS binds and sequesters nascent transcripts to prevent R-loops from forming with nearby DNA. DRIP-seq analysis showed that a knockdown of FUS increased R-loop enrichment near expressed genes. Prevention of R-loops by FUS binding to nascent transcripts has the potential to affect transcription by any RNA polymerase, highlighting the broad impact FUS can have on RNA metabolism in cells and disease.

RNA damage: the forgotten target of clinical compounds

Unlike DNA, RNA can be found in every sub-cellular compartment, where it is used to impart the genetic code or perform essential catalytic activities. As a result, damage to RNA is more spatially pervasive than damage to DNA and can have profound effects on gene expression and RNA-dependent activities. The past decade has seen the pathways involved in detecting and responding to damage of specific RNAs defined. These studies largely used high concentrations of tool compounds or deletion of essential factors for the response to RNA damage to study its effects. RNA is damaged by both endogenous and exogenous agents, with the effect of exogenous agents administered as therapeutics the focus of this review. In an effort to formalise studies into clinical RNA damage biology we propose 4 types of RNA damaging drug that we divide into 2 broad classes. Class 1 drugs result from synthesis using non-canonical nucleotides, which are incorporated into RNA in place of the canonical nucleotides. This class is subdivided depending on the outcome of this misincorporation on the nascent transcript. Class 2 drugs result in covalent ligation of moieties that alter RNA structure. This class is subdivided according to the functionality of the covalent ligation—class 2a are monovalent while class 2b are divalent. We discuss the evidence for and mechanisms of RNA damage as well as highlighting the unknown factors that require further investigation to determine the molecular mechanisms of these drugs.

M6A RNA methylation orchestrates transcriptional dormancy during paused pluripotency.

Embryos across metazoan lineages can enter reversible states of developmental pausing, or diapause, in response to adverse environmental conditions. The molecular mechanisms that underlie this remarkable dormant state remain largely unknown. Here we show that N6-methyladenosine (m6A) RNA methylation by Mettl3 is required for developmental pausing in mouse blastocysts and embryonic stem (ES) cells. Mettl3 enforces transcriptional dormancy through two interconnected mechanisms: (1) it promotes global mRNA destabilization and (2) it suppresses global nascent transcription by destabilizing the mRNA of the transcriptional amplifier and oncogene N-Myc, which we identify as a crucial anti-pausing factor. Knockdown of N-Myc rescues pausing in Mettl3-/- ES cells, and forced demethylation and stabilization of Mycn mRNA in paused wild-type ES cells largely recapitulates the transcriptional defects of Mettl3-/- ES cells. These findings uncover Mettl3 as a key orchestrator of the crosstalk between transcriptomic and epitranscriptomic regulation during developmental pausing, with implications for dormancy in adult stem cells and cancer.

Spotter: a single-nucleotide resolution stochastic simulation model of supercoiling-mediated transcription and translation in prokaryotes.

Stochastic simulation models have played an important role in efforts to understand the mechanistic basis of prokaryotic transcription and translation. Despite the fundamental linkage of these processes in bacterial cells, however, most simulation models have been limited to representations of either transcription or translation. In addition, the available simulation models typically either attempt to recapitulate data from single-molecule experiments without considering cellular-scale high-throughput sequencing data or, conversely, seek to reproduce cellular-scale data without paying close attention to many of the mechanistic details. To address these limitations, we here present spotter (Simulation of Prokaryotic Operon Transcription & Translation Elongation Reactions), a flexible, user-friendly simulation model that offers highly-detailed combined representations of prokaryotic transcription, translation, and DNA supercoiling. In incorporating nascent transcript and ribosomal profiling sequencing data, spotter provides a critical bridge between data collected in single-molecule experiments and data collected at the cellular scale. Importantly, in addition to rapidly generating output that can be aggregated for comparison with next-generation sequencing and proteomics data, spotter produces residue-level positional information that can be used to visualize individual simulation trajectories in detail. We anticipate that spotter will be a useful tool in exploring the interplay of processes that are crucially linked in prokaryotes.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

SpinDrop: a droplet microfluidic platform to maximise single-cell sequencing information content

Droplet microfluidic methods have massively increased the throughput of single-cell sequencing campaigns. The benefit of scale-up is, however, accompanied by increased background noise when processing challenging samples and the overall RNA capture efficiency is lower. These drawbacks stem from the lack of strategies to enrich for high-quality material or specific cell types at the moment of cell encapsulation and the absence of implementable multi-step enzymatic processes that increase capture. Here we alleviate both bottlenecks using fluorescence-activated droplet sorting to enrich for droplets that contain single viable cells, intact nuclei, fixed cells or target cell types and use reagent addition to droplets by picoinjection to perform multi-step lysis and reverse transcription. Our methodology increases gene detection rates fivefold, while reducing background noise by up to half. We harness these properties to deliver a high-quality molecular atlas of mouse brain development, despite starting with highly damaged input material, and provide an atlas of nascent RNA transcription during mouse organogenesis. Our method is broadly applicable to other droplet-based workflows to deliver sensitive and accurate single-cell profiling at a reduced cost.

FACT regulates pluripotency through proximal and distal regulation of gene expression in murine embryonic stem cells

BackgroundThe FACT complex is a conserved histone chaperone with critical roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear.ResultsTo explore the mechanism for FACT function, we generated AID-tagged murine embryonic cell lines for FACT subunit SPT16 and paired depletion with nascent transcription and chromatin accessibility analyses. We also analyzed SPT16 occupancy using CUT&RUN and found that SPT16 localizes to both promoter and enhancer elements, with a strong overlap in binding with OCT4, SOX2, and NANOG. Over a timecourse of SPT16 depletion, nucleosomes invade new loci, including promoters, regions bound by SPT16, OCT4, SOX2, and NANOG, and TSS-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (encoding OCT4), Sox2, Nanog, and enhancer RNAs produced from these genes’ associated enhancers are downregulated.ConclusionsWe propose that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency.

SAFB associates with nascent RNAs and can promote gene expression in mouse embryonic stem cells.

Scaffold attachment factor B (SAFB) is a conserved RNA-binding protein that is essential for early mammalian development. However, the functions of SAFB in mouse embryonic stem cells (ESCs) have not been characterized. Using RNA immunoprecipitation followed by RNA-seq (RIP-seq), we examined the RNAs associated with SAFB in wild-type and SAFB/SAFB2 double-knockout ESCs. SAFB predominantly associated with introns of protein-coding genes through purine-rich motifs. The transcript most enriched in SAFB association was the lncRNA Malat1, which also contains a purine-rich region in its 5' end. Knockout of SAFB/SAFB2 led to differential expression of approximately 1000 genes associated with multiple biological processes, including apoptosis, cell division, and cell migration. Knockout of SAFB/SAFB2 also led to splicing changes in a set of genes that were largely distinct from those that exhibited changes in expression level. The spliced and nascent transcripts of many genes whose expression levels were positively regulated by SAFB also associated with high levels of SAFB, implying that SAFB binding promotes their expression. Reintroduction of SAFB into double-knockout cells restored gene expression toward wild-type levels, an effect again observable at the level of spliced and nascent transcripts. Proteomics analysis revealed a significant enrichment of nuclear speckle-associated and RS domain-containing proteins among SAFB interactors. Neither Xist nor Polycomb functions were dramatically altered in SAFB/2 knockout ESCs. Our findings suggest that among other potential functions in ESCs, SAFB promotes the expression of certain genes through its ability to bind nascent RNA.

Domains and residues of the Saccharomyces cerevisiae hnRNP protein Hrp1 important for transcriptional autoregulation and noncoding RNA termination.

Proteins that bind the nascent transcript exiting RNA polymerase II can regulate transcription elongation. The essential Saccharomyces cerevisiae hnRNP protein Hrp1 is one such protein and participates in both cleavage and polyadenylation-coupled and Nrd1-Nab3-Sen1-dependent RNA polymerase II termination. Prior evidence that Hrp1 is a positive RNA polymerase II elongation factor suggests that its release from the elongation complex promotes termination. Here we report the effects of deletions and substitutions in Hrp1 on its autoregulation via an Nrd1-Nab3-Sen1-dependent transcription attenuator in the 5'-UTR of its mRNA and on the function of an Hrp1-dependent Nrd1-Nab3-Sen1 terminator in the SNR82 snoRNA gene. Deletion of either of two central RNA recognition motifs or either of the flanking low-sequence complexity domains is lethal. Smaller, viable deletions in the amino-terminal low-sequence complexity domain cause readthrough of both the HRP1 attenuator and SNR82 terminator. Substitutions that cause readthrough localized mostly to the RNA recognition motifs, although not always to the RNA-binding face. We found that autoregulation of Hrp1 mRNA synthesis is surprisingly robust, overcoming the expected lethal effects of the start codon and frameshift mutations via overexpression of the mRNA up to 40-fold. Our results suggest a model in which binding of attenuator or terminator elements in the nascent transcript by RNA recognition motifs 1 and 2 disrupts interactions between RNA recognition motif 2 and the RNA polymerase II elongation complex, increasing its susceptibility to termination.

A dual-function SNF2 protein drives chromatid resolution and nascent transcripts removal in mitosis.

Mitotic chromatin is largely assumed incompatible with transcription due to changes in the transcription machinery and chromosome architecture. However, the mechanisms of mitotic transcriptional inactivation and their interplay with chromosome assembly remain largely unknown. By monitoring ongoing transcription in Drosophila early embryos, we reveal that eviction of nascent mRNAs from mitotic chromatin occurs after substantial chromosome compaction and is not promoted by condensin I. Instead, we show that the timely removal of transcripts from mitotic chromatin is driven by the SNF2 helicase-like protein Lodestar (Lds), identified here as a modulator of sister chromatid cohesion defects. In addition to the eviction of nascent transcripts, we uncover that Lds cooperates with Topoisomerase 2 to ensure efficient sister chromatid resolution and mitotic fidelity. We conclude that the removal of nascent transcripts upon mitotic entry is not a passive consequence of cell cycle progression and/or chromosome compaction but occurs via dedicated mechanisms with functional parallelisms to sister chromatid resolution.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts.

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

The DUF3715 domain has a conserved role in RNA-directed transposon silencing.

RNA-directed transposon silencing operates in the mammalian soma and germline to safeguard genomic integrity. The piRNA pathway and the HUSH complex identify active transposons through recognition of their nascent transcripts, but mechanistic understanding of how these distinct pathways evolved is lacking. TASOR is an essential component of the HUSH complex. TASOR's DUF3715 domain adopts a pseudo-PARP structure and is required for transposon silencing in a manner independent of complex assembly. TEX15, an essential piRNA pathway factor, also contains the DUF3715 domain. Here, we show that TASOR's and TEX15's DUF3715 domain share extensive structural homology. We found that the DUF3715 domain arose in early eukaryotes and that in vertebrates it is restricted to TEX15, TASOR, and TASORB orthologs. While TASOR-like proteins are found throughout metazoa, TEX15 is vertebrate-specific. The branching of TEX15 and the TASOR-like DUF3715 domain likely occurred in early metazoan evolution. Remarkably, despite this vast evolutionary distance, the DUF3715 domain from divergent TEX15 sequences can functionally substitute the DUF3715 domain of TASOR and mediates transposon silencing. We have thus termed this domain of unknown function as the RNA-directed pseudo-PARP transposon silencing (RDTS) domain. In summary, we show an unexpected functional link between these critical transposon silencing pathways.

RNA polymerase II dynamics shape enhancer-promoter interactions.

How enhancers control target gene expression over long genomic distances remains an important unsolved problem. Here we investigated enhancer-promoter communication by integrating data from nucleosome-resolution genomic contact maps, nascent transcription and perturbations affecting either RNA polymerase II (Pol II) dynamics or the activity of thousands of candidate enhancers. Integration of new Micro-C experiments with published CRISPRi data demonstrated that enhancers spend more time in close proximity to their target promoters in functional enhancer-promoter pairs compared to nonfunctional pairs, which can be attributed in part to factors unrelated to genomic position. Manipulation of the transcription cycle demonstrated a key role for Pol II in enhancer-promoter interactions. Notably, promoter-proximal paused Pol II itself partially stabilized interactions. We propose an updated model in which elements of transcriptional dynamics shape the duration or frequency of interactions to facilitate enhancer-promoter communication.

ARS2 instructs early transcription termination-coupled RNA decay by recruiting ZC3H4 to nascent transcripts

The RNA-binding ARS2 protein is centrally involved in both early RNA polymerase II (RNAPII) transcription termination and transcript decay. Despite its essential nature, the mechanisms by which ARS2 enacts these functions have remained unclear. Here, we show that a conserved basic domain of ARS2 binds a corresponding acidic-rich, short linear motif (SLiM) in the transcription restriction factor ZC3H4. This interaction recruits ZC3H4 to chromatin to elicit RNAPII termination, independent of other early termination pathways defined by the cleavage and polyadenylation (CPA) and Integrator (INT) complexes. We find that ZC3H4, in turn, forms a direct connection to the nuclear exosome targeting (NEXT) complex, hereby facilitating rapid degradation of the nascent RNA. Hence, ARS2 instructs the coupled transcription termination and degradation of the transcript onto which it is bound. This contrasts with ARS2 function at CPA-instructed termination sites where the protein exclusively partakes in RNA suppression via post-transcriptional decay.

Identification of a ΔNp63-Dependent Basal-Like A Subtype-Specific Transcribed Enhancer Program (B-STEP) in Aggressive Pancreatic Ductal Adenocarcinoma.

Subtype-specific enhancer activity analysis via detection of eRNAs on a single-cell level in patient material can be used as a potential tool for treatment stratification.

Abstract A056: AR/FOXA1-mediated active gene repression as a mechanism of prostate carcinogenesis

Abstract The identity of the prostate gland is defined by transcription factors (TF) like the androgen receptor (AR), FOXA1, and NKX3.1 that bind regulatory elements early in development to activate the lineage gene program. For instance, KLK3, is a canonical target gene activated by AR. While the mechanistic interactions between pioneer factors, chromatin remodelers, and epigenetic proteins required for gene activation is well understood, repressive activity of lineage-defining TFs remain largely uncharacterized. Here, we document and mechanistically describe the active, repressive transcriptional functions of AR and FOXA1 in the context of prostate tumorigenesis.We identified a shared repertoire of repressed gene targets of AR and FOXA1 by using transcriptomic profiles (RNA-seq) of PCa cells treated with dihydrotestosterone (DHT) or TF targeting siRNA. Chromatin binding of AR and/or FOXA1 (ChIPseq) at cis-regulatory elements (dCREs) within the gene’s topologically-associated domain was used as primary evidence for direct regulation. Here, the KLK3 and DEPTOR gene loci emerged as candidate-activated or repressed AR/FOXA1 targets in PCa cells. DEPTOR, a tumor suppressor gene, negatively regulates the mTORC1/2 complexes by inhibiting their kinase activity. In contrast to the passive model of repression, where sequestration of co-activators ensues in transcriptional down-regulation of genes, we found ligand-bound AR and FOXA1 bind CREs within DEPTORs’ TAD establishing looping interactions (H3K3me3 HiChIP-seq) with the DEPTOR gene promoter. Functional inactivation of AR or FOXA1 led to a 2-4 fold upregulation of DEPTOR in VCaP, LNCaP, and 22RV1 cells. Instead, DHT-stimulation of AR led to DEPTOR downregulation. To assess the rapidness of this effect, we used a SWI/SNF ATPase PROTAC, AU15330, which triggers instant compaction of CREs and dislodges AR/FOXA1 from the chromatin. Treatment of PCa cells with AU15330 for 4 hours led to an immediate and significant upregulation of DEPTOR, substantiating the direct transcriptional repression by the AR/FOXA1 complex. Consistently, AU-15330 treatment also abolished FOXA1 and AR binding and 3D loop formation into the DEPTOR promoter. Using an orthogonal genetic approach, CRISPR-KRAB-targeted heterochromatinization of the AR/FOXA1 repressive CREs led to increased production of the nascent DEPTOR transcript. Finally, monoclones devoid of the CRE had elevated DEPTOR levels. Altogether, our findings position DEPTOR as an actively repressed target of the AR/FOXA1 transcriptional complex, suggesting the “silencer” duties of driver TF oncogenes to play a critical role in prostate tumorigenesis. Furthermore, the repertoire of activated and repressed gene targets identified in this study serve as exciting loci to study differences in the co-factor machinery and/or other cis or trans-attributes of the encompassing TADs that regulate downstream transcriptional outcomes. We posit the active mode of gene repression plays a driver role in prostate carcinogenesis. Citation Format: Sanjana Eyunni, Abhijit Parolia, Eleanor Young, James George, Mustapha Jaber, Lanbo Xiao, Fengyun Su, Rui Wang, Xuhong Cao, Arul Chinnaiyan. AR/FOXA1-mediated active gene repression as a mechanism of prostate carcinogenesis [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr A056.

Structural basis for control of bacterial RNA polymerase pausing by a riboswitch and its ligand.

Folding of nascent transcripts can be modulated by the RNA polymerase (RNAP) that carries out their transcription, and vice versa. A pause of RNAP during transcription of a preQ1 riboswitch (termed que-PEC) is stabilized by a previously characterized template consensus sequence and the ligand-free conformation of the nascent RNA. Ligand binding to the riboswitch induces RNAP pause release and downstream transcription termination; however, the mechanism by which riboswitch folding modulates pausing is unclear. Here, we report single-particle cryo-electron microscopy reconstructions of que-PEC in ligand-free and ligand-bound states. In the absence of preQ1, the RNA transcript is in an unexpected hyper-translocated state, preventing downstream nucleotide incorporation. Strikingly, on ligand binding, the riboswitch rotates around its helical axis, expanding the surrounding RNAP exit channel and repositioning the transcript for elongation. Our study reveals the tight coupling by which nascent RNA structures and their ligands can functionally regulate the macromolecular transcription machinery.

Spatial single-cell sequencing of meiosis I arrested oocytes indicates acquisition of maternal transcripts from the soma

Maternal RNAs are stored from minutes to decades in oocytes throughout meiosis I arrest in a transcriptionally quiescent state. Recent reports, however, propose a role for nascent transcription in arrested oocytes. Whether arrested oocytes launch nascent transcription in response to environmental or hormonal signals while maintaining the meiosis I arrest remains undetermined. We test this by integrating single-cell RNA sequencing, RNA velocity, and RNA fluorescence in situ hybridization on C.elegans meiosis I arrested oocytes. We identify transcripts that increase as the arrested meiosis I oocyte ages, but rule out extracellular signaling through ERK MAPK and nascent transcription as a mechanism for this increase. We report transcriptacquisition from neighboring somatic cells as a mechanism of transcript increase during meiosis I arrest. These analyses provide a deeper view at single-cell resolution of the RNA landscape of a meiosis I arrested oocyte and as it prepares for oocyte maturation and fertilization.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube.

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.