Analysis Of Regulatory Sequences Research Articles

Regulatory sequence-based discovery of anti-defense genes in archaeal viruses

In silico identification of viral anti-CRISPR proteins (Acrs) has relied largely on the guilt-by-association method using known Acrs or anti-CRISPR associated proteins (Acas) as the bait. However, the low number and limited spread of the characterized archaeal Acrs and Aca hinders our ability to identify Acrs using guilt-by-association. Here, based on the observation that the few characterized archaeal Acrs and Aca are transcribed immediately post viral infection, we hypothesize that these genes, and many other unidentified anti-defense genes (ADG), are under the control of conserved regulatory sequences including a strong promoter, which can be used to predict anti-defense genes in archaeal viruses. Using this consensus sequence based method, we identify 354 potential ADGs in 57 archaeal viruses and 6 metagenome-assembled genomes. Experimental validation identified a CRISPR subtype I-A inhibitor and the first virally encoded inhibitor of an archaeal toxin-antitoxin based immune system. We also identify regulatory proteins potentially akin to Acas that can facilitate further identification of ADGs combined with the guilt-by-association approach. These results demonstrate the potential of regulatory sequence analysis for extensive identification of ADGs in viruses of archaea and bacteria.

Experimental validation and characterization of putative targets of Escargot and STAT, two master regulators of the intestinal stem cells in Drosophila melanogaster

Many organs contain adult stem cells (ASCs) to replace cells due to damage, disease, or normal tissue turnover. ASCs can divide asymmetrically, giving rise to a new copy of themselves (self-renewal) and a sister that commits to a specific cell type (differentiation). Decades of research have led to the identification of pleiotropic genes whose loss or gain of function affect diverse aspects of normal ASC biology. Genome-wide screens of these so-called genetic “master regulator” (MR) genes, have pointed to hundreds of putative targets that could serve as their downstream effectors. Here, we experimentally validate and characterize the regulation of several putative targets of Escargot (Esg) and the Signal Transducer and Activator of Transcription (Stat92E, a.k.a. STAT), two known MRs in Drosophila intestinal stem cells (ISCs). Our results indicate that regardless of bioinformatic predictions, most experimentally validated targets show a profile of gene expression that is consistent with co-regulation by both Esg and STAT, fitting a rather limited set of co-regulatory modalities. A bioinformatic analysis of proximal regulatory sequences in specific subsets of co-regulated targets identified additional transcription factors that might cooperate with Esg and STAT in modulating their transcription. Lastly, in vivo manipulations of validated targets rarely phenocopied the effects of manipulating Esg and STAT, suggesting the existence of complex genetic interactions among downstream targets of these two MR genes during ISC homeostasis.

Predictive analyses of regulatory sequences with EUGENe

Deep learning has become a popular tool to study cis-regulatory function. Yet efforts to design software for deep-learning analyses in regulatory genomics that are findable, accessible, interoperable and reusable (FAIR) have fallen short of fully meeting these criteria. Here we present elucidating the utility of genomic elements with neural nets (EUGENe), a FAIR toolkit for the analysis of genomic sequences with deep learning. EUGENe consists of a set of modules and subpackages for executing the key functionality of a genomics deep learning workflow: (1) extracting, transforming and loading sequence data from many common file formats; (2) instantiating, initializing and training diverse model architectures; and (3) evaluating and interpreting model behavior. We designed EUGENe as a simple, flexible and extensible interface for streamlining and customizing end-to-end deep-learning sequence analyses, and illustrate these principles through application of the toolkit to three predictive modeling tasks. We hope that EUGENe represents a springboard towards a collaborative ecosystem for deep-learning applications in genomics research.

Exploring the potential role of defensins in differential vector competence of body and head lice for Bartonella quintana

BackgroundThe body and head lice of humans are conspecific, but only the body louse functions as a vector to transmit bacterial pathogens such as Bartonella quintana. Both louse subspecies have only two antimicrobial peptides, defensin 1 and defensin 2. Consequently, any differences in the molecular and functional properties of these two louse subspecies may be responsible for the differential vector competence between them.MethodsTo elucidate the molecular basis of vector competence, we compared differences in the structural properties and transcription factor/microRNA binding sites of the two defensins in body and head lice. Antimicrobial activity spectra were also investigated using recombinant louse defensins expressed via baculovirus.ResultsThe full-length amino acid sequences of defensin 1 were identical in both subspecies, whereas the two amino acid residues in defensin 2 were different between the two subspecies. Recombinant louse defensins showed antimicrobial activities only against the representative Gram-positive Staphylococcus aureus but not against either Gram-negative Escherichia coli or the yeast Candida albicans. However, they did show considerable activity against B. quintana, with body louse defensin 2 being significantly less potent than head louse defensin 2. Regulatory sequence analysis revealed that the gene units of both defensin 1 and defensin 2 in body lice possess decreased numbers of transcription factor-binding sites but increased numbers of microRNA binding sites, suggesting relatively lower transcription activities of body louse defensins.ConclusionsThe significantly lower antibacterial activities of defensin 2 along with the reduced probability of defensin expression in body lice likely contribute to the relaxed immune response to B. quintana proliferation and viability, resulting in higher vector competence of body lice compared to head lice.Graphical

Analysis of regulatory sequences in exosomal DNA of NANOGP8.

Exosomes participate in intercellular communication by transporting functionally active molecules. Such cargo from the original cells comprising proteins, micro-RNA, mRNA, single-stranded (ssDNA) and double-stranded DNA (dsDNA) molecules pleiotropically transforms the target cells. Although cancer cells secrete exosomes carrying a significant level of DNA capable of modulating oncogene expression in a recipient cell, the regulatory mechanism is unknown. We have previously reported that cancer cells produce exosomes containing NANOGP8 DNA. NANOGP8 is an oncogenic paralog of embryonic stem cell transcription factor NANOG and does not express in cells since it is a pseudogene. However, in this study, we evaluated NANOGP8 expression in glioblastoma multiforme (GBM) tissue from a surgically removed brain tumor of a patient. Significantly higher NANOGP8 transcription was observed in GBM cancer stem cells (CSCs) than in GBM cancer cells or neural stem cells (NSCs), despite identical sequences of NANOGP8-upstream genomic region in all the cell lines. This finding suggests that upstream genomic sequences of NANOGP8 may have environment-dependent promoter activity. We also found that the regulatory sequences upstream of exosomal NANOGP8 GBM DNA contain multiple core promoter elements, transcription factor binding sites, and segments of human viruses known for their oncogenic role. The exosomal sequence of NANOGP8-upstream GBM DNA is different from corresponding genomic sequences in CSCs, cancer cells, and NSCs as well as from the sequences reported by NCBI. These sequence dissimilarities suggest that exosomal NANOGP8 GBM DNA may not be a part of the genomic DNA. Exosomes possibly acquire this DNA from other sources where it is synthesized by an unknown mechanism. The significance of exosome-bestowed regulatory elements in the transcription of promoter-less retrogene such as NANOGP8 remains to be determined.

Evolutionary Changes in the Chromatin Landscape Contribute to Reorganization of a Developmental Gene Network During Rapid Life History Evolution in Sea Urchins.

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

RSAT 2022: regulatory sequence analysis tools.

RSAT (Regulatory Sequence Analysis Tools) enables the detection and the analysis of cis-regulatory elements in genomic sequences. This software suite performs (i) de novo motif discovery (including from genome-wide datasets like ChIP-seq/ATAC-seq) (ii) genomic sequences scanning with known motifs, (iii) motif analysis (quality assessment, comparisons and clustering), (iv) analysis of regulatory variations and (v) comparative genomics. RSAT comprises 50 tools. Six public Web servers (including a teaching server) are offered to meet the needs of different biological communities. RSAT philosophy and originality are: (i) a multi-modal access depending on the user needs, through web forms, command-line for local installation and programmatic web services, (ii) a support for virtually any genome (animals, bacteria, plants, totalizing over 10 000 genomes directly accessible). Since the 2018 NAR Web Software Issue, we have developed a large REST API, extended the support for additional genomes and external motif collections, enhanced some tools and Web forms, and developed a novel tool that builds or refine gene regulatory networks using motif scanning (network-interactions). The RSAT website provides extensive documentation, tutorials and published protocols. RSAT code is under open-source license and now hosted in GitHub. RSAT is available at http://www.rsat.eu/.

Construction of Light-Responsive Gene Regulatory Network for Growth, Development and Secondary Metabolite Production in Cordyceps militaris.

Simple SummaryCordyceps militaris is an edible fungus that has been long used in traditional medicine. A large body of research has provided evidence of the medicinal properties of C. militaris extract and its demand has been increasing over the years. This study aims to construct and understand the light-responsive gene regulatory network of the fungi by combining the transcription factor (TF)-target gene interactions with the transcriptomic analysis of C. militaris under a light-programming condition. The study identified several key TFs and their gene targets that regulate growth, development and secondary metabolite production in the fungi under specific light conditions.Cordyceps militaris is an edible fungus that produces many beneficial compounds, including cordycepin and carotenoid. In many fungi, growth, development and secondary metabolite production are controlled by crosstalk between light-signaling pathways and other regulatory cascades. However, little is known about the gene regulation upon light exposure in C. militaris. This study aims to construct a gene regulatory network (GRN) that responds to light in C. militaris. First, a genome-scale GRN was built based on transcription factor (TF)-target gene interactions predicted from the Regulatory Sequence Analysis Tools (RSAT). Then, a light-responsive GRN was extracted by integrating the transcriptomic data onto the genome-scale GRN. The light-responsive network contains 2689 genes and 6837 interactions. From the network, five TFs, Snf21 (CCM_04586), an AT-hook DNA-binding motif TF (CCM_08536), a homeobox TF (CCM_07504), a forkhead box protein L2 (CCM_02646) and a heat shock factor Hsf1 (CCM_05142), were identified as key regulators that co-regulate a large group of growth and developmental genes. The identified regulatory network and expression profiles from our analysis suggested how light may induce the growth and development of C. militaris into a sexual cycle. The light-mediated regulation also couples fungal development with cordycepin and carotenoid production. This study leads to an enhanced understanding of the light-responsive regulation of growth, development and secondary metabolite production in the fungi.

Genetic and molecular characterization of determinant of six-rowed spike of barley carrying vrs1.a4.

Decisive role of reduced vrs1 transcript abundance in six-rowed spike of barley carrying vrs1.a4 was genetically proved and its potential causes were preliminarily analyzed. Six-rowed spike 1 (vrs1) is the major determinant of the six-rowed spike phenotype of barley (Hordeum vulgare L.). Alleles of Vrs1 have been extensively investigated. Allele vrs1.a4 in six-rowed barley is unique in that it has the same coding sequence as Vrs1.b4 in two-rowed barley. The determinant of row-type in vrs1.a4 carriers has not been experimentally identified. Here, we identified Vrs1.b4 in two-rowed accessions and vrs1.a4 in six-rowed accessions from the Qinghai-Tibet Plateau at high frequency. Genetic analyses revealed a single nuclear gene accounting for row-type alteration in these accessions. Physical mapping identified a 0.08-cM (~ 554-kb) target interval on chromosome 2H, wherein Vrs1 was the most likely candidate gene. Further analysis of Vrs1 expression in offspring of the mapping populations or different Vrs1.b4 and vrs1.a4 lines confirmed that downregulated expression of vrs1.a4 causes six-rowed spike. Regulatory sequence analysis found a single 'TA' dinucleotide deletion in vrs1.a4 carriers within a 'TA' tandem-repeat-enriched region ~ 1kb upstream of the coding region. DNA methylation levels did not correspond to the expression difference and therefore did not affect Vrs1 expression. More evidence is needed to verify the causal link between the 'TA' deletion and the downregulated Vrs1 expression and hence the six-rowed spike phenotype.

Electric neurostimulation regulates microglial activation via retinoic acid receptor α signaling

Brain stimulation by electroconvulsive therapy is effective in neuropsychiatric disorders by unknown mechanisms. Microglial toxicity plays key role in neuropsychiatric, neuroinflammatory and degenerative diseases. We examined the mechanism by which electroconvulsive seizures (ECS) regulates microglial phenotype and response to stimuli. Microglial responses were examined by morphological analysis, Iba1 and cytokine expression. ECS did not affect resting microglial phenotype or morphology but regulated their activation by Lipopolysaccharide stimulation. Microglia were isolated after ECS or sham sessions in naïve mice for transcriptome analysis. RNA sequencing identified 141 differentially expressed genes. ECS modulated multiple immune-associated gene families and attenuated neurotoxicity-associated gene expression. Blood brain barrier was examined by injecting Biocytin-TMR tracer. There was no breakdown of the BBB, nor increase in gene-signature of peripheral monocytes, suggesting that ECS effect is mainly on resident microglia. Unbiased analysis of regulatory sequences identified the induction of microglial retinoic acid receptor α (RARα) gene expression and a putative common RARα–binding motif in multiple ECS-upregulated genes. The effects of AM580, a selective RARα agonist on microglial response to LPS was examined in vitro. AM580 prevented LPS-induced cytokine expression and reactive oxygen species production. Chronic murine experimental autoimmune encephalomyelitis (EAE) was utilized to confirm the role RARα signaling as mediator of ECS-induced transcriptional pathway in regulating microglial toxicity. Continuous intracerebroventricular delivery of AM580 attenuated effectively EAE severity. In conclusion, ECS regulates CNS innate immune system responses by activating microglial retinoic acid receptor α pathway, signifying a novel therapeutic approach for chronic neuroinflammatory, neuropsychiatric and neurodegenerative diseases.

Parallel genetics of regulatory sequences using scalable genome editing invivo.

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized invivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds ofmutations in enhancers/promoters and 3' UTRs of 16 genes in C.elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

Massively Parallel Analysis of Regulatory RNA Sequences.

The stability of RNA transcripts is regulated by signals within their sequences, but theidentity ofthose signals still remain elusive in many biological systems. Recently introduced massively parallel tools for the analysis of regulatory RNA sequences provide the ability to detect functional cis-regulatory sequences of post-transcriptional RNA regulation at a much larger scale and resolution than before. Their application formulates theunderlying sequence-based rules and predictsthe impact of genetic variations. Here, we describe the application of UTR-Seq, as a strategy to uncover cis-regulatory signals of RNA stability during early zebrafish embryogenesis. The method combines massively parallel reporter assays (MPRA) with computational regression models. It surveys the effect of tens of thousands of regulatory sequences on RNA stability andanalyzes the results via regression models to identify sequence signals that impact RNA stability and to predict the in vivo effect of sequence variations.

Analysis of Regulatory DNA Sequences by Dual-Luciferase Reporter Assays in Ewing Sarcoma.

Reporter gene assays allow for examining the influence of regulatory DNA sequences on the transcription of target genes. In Ewing sarcoma, the study of these DNA sequences is especially paramount for its main driver mutation is a fusion transcription factor that binds different motifs than its wild-type constituents. Here, we describe the process of analyzing the enhancer activity of regulatory DNA sequences using transfection-based dual-luciferase reporter assays in Ewing sarcoma cell lines. To this end, we provide a protocol for cloning sequences of interest from genomic DNA into a firefly luciferase-containing plasmid, transfecting Ewing sarcoma cells with plasmids and measuring luciferase expression by luminescence. The entire procedure can be completed in 14days.

Novel insights in the production, activity and protective effect of Penicillium expansum antifungal proteins

Antifungal proteins (AFPs) offer a great potential as new biofungicides to control deleterious fungi. The phytopathogenic fungus Penicillium expansum encodes three phylogenetically distinct AFPs, PeAfpA, PeAfpB and PeAfpC. Here, PeAfpA, a potent in vitro self-inhibitory protein, was demonstrated to control the infection caused by P. expansum in Golden apple fruits. We determined the production of the three proteins in different growth media. PeAfpA and PeAfpC were simultaneously produced by P. expansum in three out of the eight media tested as detected by Western blot, whereas PeAfpB was not detected even in those described for class B AFP production. Regardless of the culture medium, the carbon source affected Peafp expression. Notably, the production of PeAfpA was strain-dependent, but analyses of PeafpA regulatory sequences in the three strains studied could not explain differences in protein production. None of the PeAFPs was produced during apple infection, suggesting no relevant role in pathogenesis. PeAfpA together with PeAfpB and also with Penicillium digitatum PdAfpB showed synergistic interaction. The highly active antifungal PeAfpA also showed moderate antibacterial activity. We conclude that there is not a general pattern for Peafp gene expression, protein production or antimicrobial activity and confirm PeAfpA as a promising compound for postharvest conservation.

Just passing through: The auxin gradient of the root meristem.

The root meristem-one of the plant's centers of continuous growth-is a conveyer belt in which cells of different identities are pushed through gradients along the root's longitudinal axis. An auxin gradient has long been implicated in controlling the progression of cell states in the root meristem. Recent work has shown that a PLETHORA (PLT) protein transcription factor gradient, which is under a delayed auxin response, has a dose-dependent effect on the differentiation state of cells. The direct effect of auxin concentration on differential transcriptional outputs remains unclear. Genomic and other analyses of regulatory sequences show that auxin responses are likely controlled by combinatorial inputs from transcription factors outside the core auxin signaling pathway. The passage through the meristem exposes cells to varying positional signals that could help them interpret auxin inputs independent of gradient effects. One open question is whether cells process information from the changes in the gradient over time as they move through the auxin gradient.

GRADED PU.1 LEVELS ACTIVATE GRANULOCYTE VS. MACROPHAGE GENES VIA MULTIPLE (SUPER)ENHANCER ELEMENTS

Precise levels of transcription factor PU.1 are critical for differentiation to distinct blood lineages, and even modest decreases leads to leukemogenesis. However, the mechanism of how target genes sense different PU.1 levels is unknown. To address this, we used myeloid progenitors encoding inducible PU.1 transgene. Intermediate PU.1 activity induced differentiation into granulocyte-like cells, while high PU.1 produced macrophages. Comparative analysis of regulatory sequences (+-50kb) of two genes activated exclusively by intermediate PU.1, Myeloperoxidase (Mpo) and Gelatinase B (Mmp9) identified 14 Mpo and 16 Mmp9 putative enhancers resembling superenhancers. Interestingly, the PU.1 binding motifs within these regions have an intermediate to high affinity (log ods score). To test their functionality we assembled successively these regulatory sequences into luciferase vectors. Specifically, we identified novel enhancer elements in -3.4 and -15 kb of Mpo which responded to intermediate (but not high) PU.1 levels. Interestingly, activity of -3.4 kb enhancer required presence of promoter while the -15kb element required presence of both PP and the -3.4kb element. Similar phenomenon was observed at multiple Mmp9 enhancers. This model is supported by our Chromosome confirmation capture (3C) data. Native activity of tested enhancers was confirmed by expression of bidirectional enhancer RNAs. The mutation of affinity of PU.1 binding sites within enhancers lead to reprogramming of their activity. In conclusion, our data support the model that PU.1 at intermediate concentrations binds to high affinity binding sites in several enhancers of granulocyte genes, causing their successive interaction that leads to transcription activation. Collectively, this suggests that a cooperative assembly of several cell-type specific enhancers is required for optimal Mpo and Mmp9 activation.

A Systematic Bioinformatics Approach to Motif-Based Analysis of Human Locus Control Regions

Locus control regions (LCRs), cis-acting, noncoding regulatory elements with strong transcription-enhancing activity, are conserved in sequence and organization, and exhibit strict gene-specific expression. LCRs have been reported and studied in several mammalian gene systems, signifying that they play an important role in eukaryotic gene expression control. Their highly regulated, stable, and precise levels of expression have made them a strong candidate for use in gene therapy vectors. In this study, we attempted to determine the unique signatures of human LCRs by analyzing a data set of LCR sequences for the presence of motifs through systematic bioinformatics approach. Using web-based regulatory sequence analysis tools (RSAT), motif-based analysis was performed. Detected significant motifs were analyzed further for their identity using Tomtom tool. RSAT analysis revealed that significant motifs are existent within the LCRs. Identity analysis using Tomtom showed that detected significant motifs were comparable with known transcription factor (TF) binding sites and the top scoring motifs belong to zinc finger-containing proteins, an important group of proteins involved in a variety of cellular activities. Correspondence to segment of known motif indicates the biological relevance of the detected motifs. Motif-based analysis is valuable for analyzing the various characteristics of sequences, notably TF binding models in this study. Owning to their unique expression control abilities, LCRs form an important component of integrating vectors, therefore identification of unique signatures present within LCR sequences will be instrumental in the design of new generation of regulatory elements containing LCR sequences.

Bayes Factor-Based Regulatory Gene Network Analysis of Genome-Wide Association Study of Economic Traits in a Purebred Swine Population.

Early stage prediction of economic trait performance is important and directly linked to profitability of farm pig production. Genome-wide association study (GWAS) has been applied to find causative genomic regions of traits. This study established a regulatory gene network using GWAS for critical economic pig characteristics, centered on easily measurable body fat thickness in live animals. We genotyped 2,681 pigs using Illumina Porcine SNP60, followed by GWAS to calculate Bayes factors for 47,697 single nucleotide polymorphisms (SNPs) of seven traits. Using this information, SNPs were annotated with specific genes near genome locations to establish the association weight matrix. The entire network consisted of 226 nodes and 6,921 significant edges. For in silico validation of their interactions, we conducted regulatory sequence analysis of predicted target genes of transcription factors (TFs). Three key regulatory TFs were identified to guarantee maximum coverage: AT-rich interaction domain 3B (ARID3B), glial cell missing homolog 1 (GCM1), and GLI family zinc finger 2 (GLI2). We identified numerous genes targeted by ARID3B, associated with cellular processes. GCM1 and GLI2 were involved in developmental processes, and their shared target genes regulated multicellular organismal process. This system biology-based function analysis might contribute to enhancing understanding of economic pig traits.

Regulation of the cutinases expressed by Aspergillus nidulans and evaluation of their role in cutin degradation.

Four cutinase genes are encoded in the genome of the saprophytic fungus Aspergillus nidulans, but only two of them have proven to codify for active cutinases. However, their overall roles in cutin degradation are unknown, and there is scarce information on the regulatory effectors of their expression. In this work, the expression of the cutinase genes was assayed by multiplex qRT-PCR in cultures grown in media containing both inducer and repressor carbon sources. The genes ancut1 and ancut2 were induced by cutin and its monomers, while ancut3 was constitutively expressed. Besides, cutin induced ancut4 only under oxidative stress conditions. An in silico analysis of the upstream regulatory sequences suggested binding regions for the lipid metabolism transcription factors (TF) FarA for ancut1 and ancut2 while FarB for ancut3. For ancut4, the analysis suggested binding to NapA (the stress response TF). These binding possibilities were experimentally tested by transcriptional analysis using the A. nidulans mutants ANΔfarA, ANΔfarB, and ANΔnapA. Regarding cutin degradation, spectroscopic and chromatographic methods showed similar products from ANCUT1 and ANCUT3. In addition, ANCUT1 produced 9,10-dihydroxy hexadecanoic acid, suggesting an endo-cleavage action of this enzyme. Regarding ANCUT2 and ANCUT4, they produced omega fatty acids. Our results confirmed the cutinolytic activity of the four cutinases, allowed identification of their specific roles in the cutinolytic system and highlighted their differences in the regulatory mechanisms and affinity towards natural substrates. This information is expected to impact the cutinase production processes and broaden their current biotechnological applications.

RSAT variation-tools: An accessible and flexible framework to predict the impact of regulatory variants on transcription factor binding

Gene regulatory regions contain short and degenerated DNA binding sites recognized by transcription factors (TFBS). When TFBS harbor SNPs, the DNA binding site may be affected, thereby altering the transcriptional regulation of the target genes. Such regulatory SNPs have been implicated as causal variants in Genome-Wide Association Study (GWAS) studies. In this study, we describe improved versions of the programs Variation-tools designed to predict regulatory variants, and present four case studies to illustrate their usage and applications. In brief, Variation-tools facilitate i) obtaining variation information, ii) interconversion of variation file formats, iii) retrieval of sequences surrounding variants, and iv) calculating the change on predicted transcription factor affinity scores between alleles, using motif scanning approaches. Notably, the tools support the analysis of haplotypes. The tools are included within the well-maintained suite Regulatory Sequence Analysis Tools (RSAT, http://rsat.eu), and accessible through a web interface that currently enables analysis of five metazoa and ten plant genomes. Variation-tools can also be used in command-line with any locally-installed Ensembl genome. Users can input personal collections of variants and motifs, providing flexibility in the analysis.