Variable Motifs Research Articles

Epigenetic dissection of human blood group genes reveals regulatory elements and detailed characteristics of KEL and four other loci.

Blood typing is essential for safe transfusions and is performed serologically or genetically. Genotyping predominantly focuses on coding regions, but non-coding variants may affect gene regulation, as demonstrated in the ABO, FY and XG systems. To uncover regulatory loci, we expanded a recently developed bioinformatics pipeline for discovery of non-coding variants by including additional epigenetic datasets. Multiple datasets including ChIP-seq with erythroid transcription factors (TFs), histone modifications (H3K27ac, H3K4me1), and chromatin accessibility (ATAC-seq) were analyzed. Candidate regulatory regions were investigated for activity (luciferase assays) and TF binding (electrophoretic mobility shift assay, EMSA, and mass spectrometry, MS). In total, 814 potential regulatory sites in 47 blood-group-related genes were identified where one or more erythroid TFs bound. Enhancer candidates in CR1, EMP3, ABCB6, and ABCC4 indicated by ATAC-seq, histone markers, and co-occupancy of 4 TFs (GATA1/KLF1/RUNX1/NFE2) were investigated but only CR1 and ABCC4 showed increased transcription. Co-occupancy of GATA1 and KLF1 was observed in the KEL promoter, previously reported to contain GATA1 and Sp1 sites. TF binding energy scores decreased when three naturally occurring variants were introduced into GATA1 and KLF1 motifs. Two of three GATA1 sites and the KLF1 site were confirmed functionally. EMSA and MS demonstrated increased GATA1 and KLF1 binding to the wild-type compared to variant motifs. This combined bioinformatics and experimental approach revealed multiple candidate regulatory regions and predicted TF co-occupancy sites. The KEL promoter was characterized in detail, indicating that two adjacent GATA1 and KLF1 motifs are most crucial for transcription.

Abstract B057: Oncogenic Foxl2 directs enhancer reprogramming and changes in 3D genome structure in ovarian granulosa cell tumors

Abstract Background: Adult-type granulosa cell tumors (aGCTs) are a rare monogenic disease. Nearly all aGCTs carry a missense point mutation in the Forkhead domain containing FOXL2 (Foxl2 p.C134W) transcription factor, but the oncogenic mechanism of this mutation is unknown. Other Forkhead family transcription factors have well-described “pioneer” activity, binding to compacted, nucleosome-bound DNA and increasing accessibility for other regulatory proteins. Objectives: To investigate if oncogenic Foxl2-C134W has pioneering activity and if mutant Foxl2 operates its oncogenic activity through enhancer reprogramming. Methods: CRISPR/Cas9 editing was used to generate isogenic granulosa cell tumor cells lacking either the FOXL2 wild-type allele (single knock-out; SKO) or both the mutant and wild-type FOXL2 alleles (double knock-out; DKO). ATAC-Seq and endogenous Foxl2 ChIP-Seq were performed on these isogenic lines to determine the differential chromatin accessibility at Foxl2-bound regulatory regions across genotypes. Promoter capture HiC was performed using a bespoke genome-wide hybrid capture probe set, CHiCAGO was used to identify significant HiC interactions, and an activity-by-contact model was constructed to identify enhancer elements bound by endogenous Foxl2-C134W and determine their target promoter contact frequency (promoter capture HiC) and enhancer activity (geometric mean of H3K27ac and ATACseq signal. Results: Endogenous ChIP-seq of Foxl2-C134W from the SKO cell line identified 1147 high-confidence peaks. De novo motif discovery identified the canonical Foxl2 binding motif (TGTTTACATT, P = 1.4 × 10−7) in 44.9% of peaks, as well as a novel variant motif (TGTTTTGTCT, P = 6.7 × 10−15) in 68.5% of peaks. Median chromatin accessibility at Foxl2-C134W peak regions, as measured by ATAC-seq, was significantly decreased in the DKO cells compared to either SKO or parental cell lines (P < 2 × 10-16 for both comparisons). Decreased chromatin accessibility at Foxl2-C134W ChIP-seq peaks in the DKO cells was driven by peaks containing the novel variant Foxl2 binding motif. Promoter capture HiC through activity-by-contact model identified promoter-enhancer interactions lost in the DKO cell line involving several genes, including ADAMTS5, FOXF2, and PRICKLE1 (adjusted P < 0.05). Conclusion: Foxl2-C134W exhibits “pioneering” activity, increasing chromatin accessibility at key gene regulatory elements with associated re-programming of promoter-enhancer interactions. The oncogenic mechanism of Foxl2-C134W in granulosa cell tumors may involve changes in DNA binding specificity, re-directing this pioneering function to sites containing a novel variant Foxl2 binding motif. Citation Format: Veena K. Vuttaradhi, Eleonora Khlebus, Thomas Welte, Barrett Lawson, Robert Tyler Hillman. Oncogenic Foxl2 directs enhancer reprogramming and changes in 3D genome structure in ovarian granulosa cell tumors [abstract]. In: Proceedings of the AACR Special Conference on Ovarian Cancer; 2023 Oct 5-7; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_2):Abstract nr B057.

Single cell multiomic profiling of high-grade serous carcinoma.

e17535 Background: High-grade serous ovarian carcinoma (HGSOC) is the most frequent and deadliest type of ovarian cancer (OC). With the recent advances in transcriptomic and epigenomic profiling of cancer at single cell resolution, it has become clear that intra-tumor heterogeneity driven by genetic and epigenetic factors may be the basis of drug resistance. Current treatments are not homogeneously effective against all the cancer cell subpopulations, thus enabling resistance. With single cell sequencing technology, characterization of cell signatures and regulatory mechanisms may translate to prognostic markers and novel drug targets. Methods: The transcriptomic and epigenomic landscape of HGSOC was mapped for five tumor samples from debulking surgeries and one tissue patient-derived xenograft (PDX) using the 10X genomics sequencing platform. This strategy permits the combination of single cell ATAC-seq and single cell RNA-seq within the same nucleus (Multiome). Isolation of nuclei, library preparation and sequencing were performed following 10X genomics protocols. Sample quality and library quality were assessed by nuclear morphology and fragment size analysis, respectively. Sequenced reads were mapped to the human genome and downstream analysis including copy number variant prediction, pathway enrichment, motif enrichment and gene regulatory network analysis were performed. The sampled patients continue to be followed for clinical outcomes, such as response to therapy. Clinical data including patient age, grade, and stage of cancer, debulking status, CA-125 levels, and neoadjuvant status, were used in conjunction with the genomics data to characterize patient specific molecular regulatory signatures. Results: Cells from the six samples (total N= 26,421) were projected via UMAP. Immune and stromal cells were shown to cluster by cell-type while cancer cells clustered by patient. Cancer cells were identified as cycling cells, ciliated cells, cells with enrichment of the JAK-STAT signaling pathway, and cells exhibiting cancer stemness signatures. Transcription factor motifs and binding-site enrichment in open chromatin regions reveal transcription regulation-based subpopulations. Using the combined gene expression and open chromatin information from these cells, there is potential to uncover the genetic regulatory network(s) that drives treatment resistance. Conclusions: With single cell technology, specific clusters of cancer and tumor microenvironment cells can be classified. Beyond characterizing patient specific signatures, multiome enables the discovery of genomic, transcriptomic and epigenomic signatures that provide insight in tumor progression and treatment resistance.

Abstract 1462: Oncogenic Foxl2 is a chromatin-remodeling pioneer transcription factor in adult-type ovarian granulosa cell tumors

Abstract Background: A missense mutation in the Forkhead domain-containing FOXL2 (Foxl2 p.C134W) transcription factor is found in nearly all adult-type granulosa cell tumors, but the oncogenic mechanism of this mutation is not known. Other Forkhead family transcription factors have well-described “pioneer” activity, binding to compacted, nucleosome-bound DNA and increasing accessibility for other regulatory proteins. Objectives: To develop novel cell culture model systems and determine whether oncogenic Foxl2-C134W has pioneering activity. Methods: CRISPR/Cas9 editing was used to generate isogenic granulosa cell tumor cells lacking either the FOXL2 wild-type allele (single knock-out; SKO) or both the mutant and wild-type FOXL2 alleles (double knock-out; DKO). ATAC-Seq and endogenous Foxl2 ChIP-Seq were performed on these isogenic lines to determine the differential chromatin accessibility at Foxl2-bound regulatory regions across genotypes. ENCODE pipelines and data standards were used for analysis, and the irreproducible discovery rate was used to identify high-reliability ATAC-seq and ChIP-seq peaks. De novo motifs were identified with STREME. Promoter capture HiC was performed using a bespoke genome-wide hybrid capture probe set, and CHiCAGO was used to identify significant HiC interactions. Results: Endogenous ChIP-seq of Foxl2-C134W from the SKO cell line identified 1147 high-confidence peaks. De novo motif discovery identified the canonical Foxl2 binding motif (TGTTTACATT, P = 1.4 x 10-7) in 44.9% of peaks, as well as a novel variant motif (TGTTTTGTCT, P = 6.7 x 10-15) in 68.5% of peaks. Median chromatin accessibility at Foxl2-C134W peak regions, as measured by ATAC-seq, was significantly decreased in the DKO cells compared to either SKO or parental cell lines (P < 2 x 10-16 for both comparisons). No difference was observed between SKO and parental cell lines (P = 0.19). Decreased chromatin accessibility at Foxl2-C134W ChIP-seq peaks in the DKO cells was driven by peaks containing the novel variant Foxl2 binding motif. Promoter capture HiC identified promoter-enhancer interactions lost in the DKO cell line involving several genes, including HUNK, NOVA1, andSFRP2 (adjusted P < 0.05). Conclusion: Foxl2-C134W exhibits “pioneering” activity, increasing chromatin accessibility at key gene regulatory elements with associated re-programming of promoter-enhancer interactions. The oncogenic mechanism of Foxl2-C134W in granulosa cell tumors may involve changes in DNA binding specificity, re-directing this pioneering function to sites containing a novel variant Foxl2 binding motif. Citation Format: Veena Vuttaradhi, Eleonora Khlebus, Thomas Welte, Namrata Khurana, Barrett Lawson, Robert Tyler Hillman. Oncogenic Foxl2 is a chromatin-remodeling pioneer transcription factor in adult-type ovarian granulosa cell tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1462.

A variant of guanidine-IV riboswitches exhibits evidence of a distinct ligand specificity

ABSTRACT Riboswitches are regulatory RNAs that specifically bind a small molecule or ion. Like metabolite-binding proteins, riboswitches can evolve new ligand specificities, and some examples of this phenomenon have been validated. As part of work based on comparative genomics to discover novel riboswitches, we encountered a candidate riboswitch with striking similarities to the recently identified guanidine-IV riboswitch. This candidate riboswitch, the Gd4v motif, is predicted in four distinct bacterial phyla, thus almost as widespread as the guanidine-IV riboswitch. Bioinformatic and experimental analysis suggest that the Gd4v motif is a riboswitch that binds a ligand other than guanidine. It is found associated with gene classes that differ from genes regulated by confirmed guanidine riboswitches. In inline-probing assays, we showed that free guanidine binds only weakly to one of the tested sequences of the variant. Further tested compounds did not show binding, attenuation of transcription termination, or activation of a genetic reporter construct. We characterized an N-acetyltransferase frequently associated with the Gd4v motif and compared its substrate preference to an N-acetyltransferase that occurs under control of guanidine-IV riboswitches. The substrates of this Gd4v-motif-associated enzyme did not show activity for Gd4v RNA binding or transcription termination. Hence, the ligand of the candidate riboswitch motif remains unidentified. The variant RNA motif is predominantly found in gut metagenome sequences, hinting at a ligand that is highly relevant in this environment. This finding is a first step to determining the identity of this unknown ligand, and understanding how guanidine-IV-riboswitch-like structures can evolve to bind different ligands.

Abstract 2356: The epigenomic characterization of clear cell sarcoma

Abstract Clear cell sarcoma (CCS) is an aggressive malignancy that affects adolescents and young adults. CCS is a translocation-driven sarcoma subtype meaning a genetic fusion event defines and drives the malignancy. The chromosomal translocation occurring between EWSR1 on chromosome 22 and ATF1 on chromosome 12 results in the EWSR1-ATF1 fusion oncogene (hereafter abbreviated EA1) which drives more than 90% of CCS cases. The EA1 oncoprotein is a transcription factor that reprograms gene expression to drive CCS. However, the target genes of EA1 and how EA1 affects target gene expression levels are largely unknown. Elucidation of these effectors is critical in order to understand how the EA1 oncoprotein drives the epigenetic reprogramming in CCS. We utilized a robust mouse model that recapitulates CCS in which mice express TATCre-inducible EA1. Using CCS tumors that form in these mice, ChIP-sequencing identified where EA1 is binding along the genome. Analysis in combination with RNA-sequencing data revealed target genes that are bound by EA1 and either upregulated or downregulated compared to control tissue. Target gene expression levels were corroborated in a human context by RNA-sequencing of human tumors and a CCS human cell line. RNAi knockdown of EA1 in the human cell line confirmed functionality of the fusion protein at these target genes. We show that EA1 can bind at promoter regions to directly regulate transcription of target genes. These binding sites include the canonical ATF binding motif. The EA1 target genes overlap with known wildtype ATF1 target genes. Many upregulated genes are involved in cell cycle and proliferation as expected. Alternatively, EA1 can bind intergenic regions including H3K27ac-defined super enhancers to distally regulate target genes. These binding sites contain unique DNA motifs that are different from the canonical ATF binding motif by 1 base pair. One highly recurrent motif identified matches the motif of the AP1 complex, a putative co-factor that may be functioning with EA1. These so-called variant motifs appear to define enhancer regions where EA1 can bind along with AP1. In summary, EA1 binds the expected ATF motif at promoter regions to directly regulate target gene expression. Secondly, EA1 binds intergenic enhancer regions defined by variant motifs including the AP1 motif. Target genes can either be upregulated or downregulated by the fusion, suggesting that there are secondary epigenetic mechanisms. One such model involves differential co-factor composition, such as AP1 co-binding, that dictates where EA1 binds or subsequent transcriptional effects. This research sets up the foundation for how EA1 functions as a transcription factor including the distribution of EA1 across the genome, the motifs that EA1 binds, and the key genes that are regulated by EA1. These epigenetic findings are crucial to understanding how a powerful fusion oncoprotein drives this aggressive, under-studied malignancy. Citation Format: Benjamin B. Ozenberger, Li Li, Kevin B. Jones. The epigenomic characterization of clear cell sarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2356.

A Horse Is a Horse

A Horse Is a Horse

The Signaling Pathway That cGAMP Riboswitches Found: Analysis and Application of Riboswitches to Study cGAMP Signaling in Geobacter sulfurreducens.

The Hypr cGAMP signaling pathway was discovered via the function of the riboswitch. In this study, we show the development of a method for affinity capture followed by sequencing to identify non-coding RNA regions that bind nucleotide signals such as cGAMP. The RNAseq of affinity-captured cGAMP riboswitches from the Geobacter sulfurreducens transcriptome highlights general challenges that remain for this technique. Furthermore, by applying riboswitch reporters in vivo, we identify new growth conditions and transposon mutations that affect cGAMP levels in G. sulfurreducens. This work reveals an extensive regulatory network and supports a second functional cGAMP synthase gene in G. sulfurreducens. The activity of the second synthase was validated using riboswitch-based fluorescent biosensors, and is the first known example of an active enzyme with a variant GGDDF motif.

GGGCTA repeats can fold into hairpins poorly unfolded by replication protein A: a possible origin of the length-dependent instability of GGGCTA variant repeats in human telomeres.

Human telomeres are composed of GGGTTA repeats and interspersed with variant repeats. The GGGCTA variant motif was identified in the proximal regions of human telomeres about 10 years ago and was shown to display a length-dependent instability. In parallel, a structural study showed that four GGGCTA repeats folded into a non-canonical G-quadruplex (G4) comprising a Watson–Crick GCGC tetrad. It was proposed that this non-canonical G4 might be an additional obstacle for telomere replication. In the present study, we demonstrate that longer GGGCTA arrays fold into G4 and into hairpins. We also demonstrate that replication protein A (RPA) efficiently binds to GGGCTA repeats structured into G4 but poorly binds to GGGCTA repeats structured into hairpins. Our results (along with results obtained with a more stable variant motif) suggest that GGGCTA hairpins are at the origin of GGGCTA length-dependent instability. They also suggest, as working hypothesis, that failure of efficient binding of RPA to GGGCTA structured into hairpins might be involved in the mechanism of GGGCTA array instability. On the basis of our present and past studies about telomeric G4 and their interaction with RPA, we propose an original point of view about telomeric G4 and the evolution of telomeric motifs.

Fra-1 regulates its target genes via binding to remote enhancers without exerting major control on chromatin architecture in triple negative breast cancers.

The ubiquitous family of dimeric transcription factors AP-1 is made up of Fos and Jun family proteins. It has long been thought to operate principally at gene promoters and how it controls transcription is still ill-understood. The Fos family protein Fra-1 is overexpressed in triple negative breast cancers (TNBCs) where it contributes to tumor aggressiveness. To address its transcriptional actions in TNBCs, we combined transcriptomics, ChIP-seqs, machine learning and NG Capture-C. Additionally, we studied its Fos family kin Fra-2 also expressed in TNBCs, albeit much less. Consistently with their pleiotropic effects, Fra-1 and Fra-2 up- and downregulate individually, together or redundantly many genes associated with a wide range of biological processes. Target gene regulation is principally due to binding of Fra-1 and Fra-2 at regulatory elements located distantly from cognate promoters where Fra-1 modulates the recruitment of the transcriptional co-regulator p300/CBP and where differences in AP-1 variant motif recognition can underlie preferential Fra-1- or Fra-2 bindings. Our work also shows no major role for Fra-1 in chromatin architecture control at target gene loci, but suggests collaboration between Fra-1-bound and -unbound enhancers within chromatin hubs sometimes including promoters for other Fra-1-regulated genes. Our work impacts our view of AP-1.

Transcriptome-wide identification of WRKY family genes and their expression profiling toward salicylic acid in Camellia japonica

ABSTRACT The ornamental plant Camellia japonica is widely distributed worldwide and is susceptible to various environmental stresses. The WRKY transcription factor (TF) is an important node of plant tolerance. However, WRKY TFs from C. japonica have not been reported yet. In this study, 48 CjWRKYs, namely, CjWRKY1 to CjWRKY48, were identified. Protein structure analysis revealed that CjWRKY proteins contain a highly conserved motif (WRKYGQK) and two variant motifs (WRKYGKK and WRKYGRK). Phylogenetic analysis indicated that the 48 CjWRKYs can be divided into three groups, which are further classified into six subgroups, namely, I-C, II-a, II-b, II-c, II-e, and III, which contain 10, 6, 8, 13, 7, and 4 members, respectively. The expression patterns of 15 CjWRKYs under salicylic acid (SA) treatment were investigated by real-time quantitative PCR (qRT-PCR). Results showed that the 15 CjWRKYs could be induced by SA treatment. This study is the first to screen CjWRKYs and identify the expression profile of CjWRKYs under SA treatment and provides a theoretical basis for analyzing the function of CjWRKY genes to SA stress tolerance in C. japonica.

Crystal Structure of a Variant PAM2 Motif of LARP4B Bound to the MLLE Domain of PABPC1.

Eukaryotic cells determine the protein output of their genetic program by regulating mRNA transcription, localization, translation and turnover rates. This regulation is accomplished by an ensemble of RNA-binding proteins (RBPs) that bind to any given mRNA, thus forming mRNPs. Poly(A) binding proteins (PABPs) are prominent members of virtually all mRNPs that possess poly(A) tails. They serve as multifunctional scaffolds, allowing the recruitment of diverse factors containing a poly(A)-interacting motif (PAM) into mRNPs. We present the crystal structure of the variant PAM motif (termed PAM2w) in the N-terminal part of the positive translation factor LARP4B, which binds to the MLLE domain of the poly(A) binding protein C1 cytoplasmic 1 (PABPC1). The structural analysis, along with mutational studies in vitro and in vivo, uncovered a new mode of interaction between PAM2 motifs and MLLE domains.

Dynamics of Influenza A (H5N1) virus protein sequence diversity.

BackgroundInfluenza A (H5N1) virus is a global concern with potential as a pandemic threat. High sequence variability of influenza A viruses is a major challenge for effective vaccine design. A continuing goal towards this is a greater understanding of influenza A (H5N1) proteome sequence diversity in the context of the immune system (antigenic diversity), the dynamics of mutation, and effective strategies to overcome the diversity for vaccine design.MethodsHerein, we report a comprehensive study of the dynamics of H5N1 mutations by analysis of the aligned overlapping nonamer positions (1–9, 2–10, etc.) of more than 13,000 protein sequences of avian and human influenza A (H5N1) viruses, reported over at least 50 years. Entropy calculations were performed on 9,408 overlapping nonamer position of the proteome to study the diversity in the context of immune system. The nonamers represent the predominant length of the binding cores for peptides recognized by the cellular immune system. To further dissect the sequence diversity, each overlapping nonamer position was quantitatively analyzed for four patterns of sequence diversity motifs: index, major, minor and unique.ResultsAlmost all of the aligned overlapping nonamer positions of each viral proteome exhibited variants (major, minor, and unique) to the predominant index sequence. Each variant motif displayed a characteristic pattern of incidence change in relation to increased total variants. The major variant exhibited a restrictive pyramidal incidence pattern, with peak incidence at 50% total variants. Post this peak incidence, the minor variants became the predominant motif for majority of the positions. Unique variants, each sequence observed only once, were present at nearly all of the nonamer positions. The diversity motifs (index and variants) demonstrated complex inter-relationships, with motif switching being a common phenomenon. Additionally, 25 highly conserved sequences were identified to be shared across viruses of both hosts, with half conserved to several other influenza A subtypes.DiscussionThe presence of distinct sequences (nonatypes) at nearly all nonamer positions represents a large repertoire of reported viral variants in the proteome, which influence the variability dynamics of the viral population. This work elucidated and provided important insights on the components that make up the viral diversity, delineating inherent patterns in the organization of sequence changes that function in the viral fitness-selection. Additionally, it provides a catalogue of all the mutational changes involved in the dynamics of H5N1 viral diversity for both avian and human host populations. This work provides data relevant for the design of prophylactics and therapeutics that overcome the diversity of the virus, and can aid in the surveillance of existing and future strains of influenza viruses.

Molecular characterization of insulin-like peptides in the brown planthopper, Nilaparvata lugens (Hemiptera: Delphacidae).

Insulin-like peptides (ILPs) including insulin, insulin-like growth factor (IGF) and relaxin are evolutionarily conserved hormones in metazoans, and they are involved in diverse physiological processes. The migratory brown planthopper (BPH), Nilaparvata lugens, encodes four ILP genes (Nlilp1, Nlilp2, Nlilp3 and Nlilp4) but their physiological roles are largely unknown. Sequence analysis showed that NlILP1 contained a relaxin-specific G protein-coupled receptor-binding motif and a variant motif of cysteine residues, and NlILP2 and NlILP4 resembled vertebrate IGFs. RNA interference (RNAi)-mediated gene silencing showed that depletion of each of Nlilp1, 2 and 3 significantly delayed the developmental duration of nymphs, and this effect could be exacerbated by double or triple gene depletion. Depletion of Nlilp1, Nlilp2 or Nlilp3 induces the accumulation of glucose, trehalose and glycogen, which is contradictory to depletion of the insulin receptor (NlInR1) in the BPH. Depletion of Nlilp1 significantly enhanced starvation resistance in both females and males although its extent was smaller than NlInR1 depletion. A parental RNAi assay showed that depletion of each of Nlilp1-4 dramatically impaired female fecundity. These findings indicate that NlILP1-4 have redundant and distinct roles in physiological processes in the BPH, thereby enhancing our understanding of the contribution of each NlILP to the ecological success of this species in natural habitats.

The KLDpT activation loop motif is critical for MARK kinase activity.

MAP/microtubule-affinity regulating kinases (MARK1-4) are members of the AMPK family of Ser/Thr-specific kinases, which phosphorylate substrates at consensus LXRXXSXXXL motifs. Within microtubule-associated proteins, MARKs also mediate phosphorylation of variant KXGS or ζXKXGSXXNΨ motifs, interfering with the ability of tau and MAP2/4 to bind to microtubules. Here we show that, although MARKs and the closely related salt-inducible kinases (SIKs) phosphorylate substrates with consensus AMPK motifs comparably, MARKs are more potent in recognizing variant ζXKXGSXXNΨ motifs on cellular tau. In studies to identify regions of MARKs that confer catalytic activity towards variant sites, we found that the C-terminal kinase associated-1 (KA1) domain in MARK1-3 mediates binding to microtubule-associated proteins CLASP1/2; but this interaction is dispensable for ζXKXGSXXNΨ phosphorylation. Mutational analysis of MARK2 revealed that the N-terminal kinase domain of MARK2 is sufficient for phosphorylation of both consensus and variant ζXKXGSXXNΨ sites. Within this domain, the KLDpT activation loop motif promotes MARK2 activity both intracellularly and in vitro, but has no effect on SIK2 activity. As KLDpT is conserved in all vertebrates MARKs, we conclude that this sequence is crucial for MARK-dependent regulation of cellular polarity.

LARP4A recognizes polyA RNA via a novel binding mechanism mediated by disordered regions and involving the PAM2w motif, revealing interplay between PABP, LARP4A and mRNA.

LARP4A belongs to the ancient RNA-binding protein superfamily of La-related proteins (LARPs). In humans, it acts mainly by stabilizing mRNAs, enhancing translation and controlling polyA lengths of heterologous mRNAs. These activities are known to implicate its association with mRNA, protein partners and translating ribosomes, albeit molecular details are missing. Here, we characterize the direct interaction between LARP4A, oligoA RNA and the MLLE domain of the PolyA-binding protein (PABP). Our study shows that LARP4A–oligoA association entails novel RNA recognition features involving the N-terminal region of the protein that exists in a semi-disordered state and lacks any recognizable RNA-binding motif. Against expectations, we show that the La module, the conserved RNA-binding unit across LARPs, is not the principal determinant for oligoA interaction, only contributing to binding to a limited degree. Furthermore, the variant PABP-interacting motif 2 (PAM2w) featured in the N-terminal region of LARP4A was found to be important for both RNA and PABP recognition, revealing a new role for this protein–protein binding motif. Our analysis demonstrates the mutual exclusive nature of the PAM2w-mediated interactions, thereby unveiling a tantalizing interplay between LARP4A, polyA and PABP.

Genome wide analysis of W-box element in Arabidopsis thaliana reveals TGAC motif with genes down regulated by heat and salinity

To design, synthetic promoters leading to stress-specific induction of a transgene, the study of cis-regulatory elements is of great importance. Cis-regulatory elements play a major role in regulating the gene expression spatially and temporally at the transcriptional level. The present work focuses on one of the important cis-regulatory element, W-box having TGAC as a core motif which serves as a binding site for the members of the WRKY transcription factor family. In the present study, we have analyzed the occurrence frequency of TGAC core motifs for varying spacer lengths (ranging from 0 to 30 base pairs) across the Arabidopsis thaliana genome in order to determine the biological and functional significance of these conserved sequences. Further, the available microarray data was used to determine the role of TGAC motif in abiotic stresses namely salinity, osmolarity and heat. It was observed that TGAC motifs with spacer sequences like TGACCCATTTTGAC and TGACCCATGAATTTTGAC had a significant deviation in frequency and were thought to be favored for transcriptional bindings. The microarray data analysis revealed the involvement of TGAC motif mainly with genes down-regulated under abiotic stress conditions. These results were further confirmed by the transient expression studies with promoter-reporter cassettes carrying TGAC and TGAC-ACGT variant motifs with spacer lengths of 5 and 10.

Structural basis for brassinosteroid response by BIL1/BZR1

BRZ-INSENSITIVE-LONG HYPOCOTYL 1 (BIL1)/BRASSINAZOLE-RESISTANT 1 (BZR1) is a master transcription factor of brassinosteroid (BR) signalling. The varieties of nucleobase recognition of the NN-BRRE-core motif (NNCGTG), one of variant G-box motifs, distinguish BIL1/BZR1 from basic helix-loop-helix transcription factors, underlying the specific regulation of BR-responsive genes. Here, we show the non-canonical bHLH dimer formation of BIL1/BZR1 to optimize the interaction network with DNA and the orientation of a key residue for NN-BRRE-core motif recognition.

Identification of the WRKY gene family and functional analysis of two genes in Caragana intermedia

BackgroundWRKY transcription factors, one of the largest families of transcriptional regulators in plants, play important roles in plant development and various stress responses. The WRKYs of Caragana intermedia are still not well characterized, although many WRKYs have been identified in various plant species.ResultsWe identified 53 CiWRKY genes from C. intermedia transcriptome data, 28 of which exhibited complete open reading frames (ORFs). These CiWRKYs were divided into three groups via phylogenetic analysis according to their WRKY domains and zinc finger motifs. Conserved domain analysis showed that the CiWRKY proteins contain a highly conserved WRKYGQK motif and two variant motifs (WRKYGKK and WKKYEEK). The subcellular localization of CiWRKY26 and CiWRKY28–1 indicated that these two proteins localized exclusively to nuclei, supporting their role as transcription factors. The expression patterns of the 28 CiWRKYs with complete ORFs were examined through quantitative real-time PCR (qRT-PCR) in various tissues and under different abiotic stresses (drought, cold, salt, high-pH and abscisic acid (ABA)). The results showed that each CiWRKY responded to at least one stress treatment. Furthermore, overexpression of CiWRKY75–1 and CiWRKY40–4 in Arabidopsis thaliana suppressed the drought stress tolerance of the plants and delayed leaf senescence, respectively.ConclusionsFifty-three CiWRKY genes from the C. intermedia transcriptome were identified and divided into three groups via phylogenetic analysis. The expression patterns of the 28 CiWRKYs under different abiotic stresses suggested that each CiWRKY responded to at least one stress treatment. Overexpression of CiWRKY75–1 and CiWRKY40–4 suppressed the drought stress tolerance of Arabidopsis and delayed leaf senescence, respectively. These results provide a basis for the molecular mechanism through which CiWRKYs mediate stress tolerance.

THiCweed: fast, sensitive detection of sequence features by clustering big datasets.

We present THiCweed, a new approach to analyzing transcription factor binding data from high-throughput chromatin immunoprecipitation-sequencing (ChIP-seq) experiments. THiCweed clusters bound regions based on sequence similarity using a divisive hierarchical clustering approach based on sequence similarity within sliding windows, while exploring both strands. ThiCweed is specially geared toward data containing mixtures of motifs, which present a challenge to traditional motif-finders. Our implementation is significantly faster than standard motif-finding programs, able to process 30 000 peaks in 1–2 h, on a single CPU core of a desktop computer. On synthetic data containing mixtures of motifs it is as accurate or more accurate than all other tested programs. THiCweed performs best with large ‘window’ sizes (≥50 bp), much longer than typical binding sites (7–15 bp). On real data it successfully recovers literature motifs, but also uncovers complex sequence characteristics in flanking DNA, variant motifs and secondary motifs even when they occur in <5% of the input, all of which appear biologically relevant. We also find recurring sequence patterns across diverse ChIP-seq datasets, possibly related to chromatin architecture and looping. THiCweed thus goes beyond traditional motif finding to give new insights into genomic transcription factor-binding complexity.