Haruka Ozaki

CLOCK-Controlled Polyphonic Regulation of Circadian Rhythms through Canonical and Noncanonical E-Boxes

ABSTRACT In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.

CapR: revealing structural specificities of RNA-binding protein target recognition using CLIP-seq data

RNA-binding proteins (RBPs) bind to their target RNA molecules by recognizing specific RNA sequences and structural contexts. The development of CLIP-seq and related protocols has made it possible to exhaustively identify RNA fragments that bind to RBPs. However, no efficient bioinformatics method exists to reveal the structural specificities of RBP–RNA interactions using these data. We present CapR, an efficient algorithm that calculates the probability that each RNA base position is located within each secondary structural context. Using CapR, we demonstrate that several RBPs bind to their target RNA molecules under specific structural contexts. CapR is available at https://sites.google.com/site/fukunagatsu/software/capr.

ADARB1 catalyzes circadian A-to-I editing and regulates RNA rhythm

It has been proposed that the CLOCK–ARNTL (BMAL1) complex drives circadian transcription of thousands of genes, including Per and Cry family genes that encode suppressors of CLOCK–ARNTL-dependent transcription. However, recent studies demonstrated that 70–80% of circadian-oscillating mRNAs have no obvious rhythms in their de novo transcription, indicating the potential importance of post-transcriptional regulation. Our CLOCK-ChIP-seq analysis identified rhythmic expression of adenosine deaminase, RNA-specific, B1 (Adarb1, also known as Adar2), an adenosine-to-inosine (A-to-I) RNA-editing enzyme. RNA-seq showed circadian rhythms of ADARB1-mediated A-to-I editing in a variety of transcripts. In Adarb1-knockout mice, rhythms of large populations of mRNA were attenuated, indicating a profound impact of ADARB1-mediated A-to-I editing on RNA rhythms. Furthermore, Adarb1-knockout mice exhibited short-period rhythms in locomotor activity and gene expression. These phenotypes were associated with abnormal accumulation of CRY2. The present study identifies A-to-I RNA editing as a key mechanism of post-transcriptional regulation in the circadian clockwork.

Parallel selection on gene copy number variations through evolution of three-spined stickleback genomes

BackgroundUnderstanding the genetic basis of adaptive evolution is one of the major goals in evolutionary biology. Recently, it has been revealed that gene copy number variations (GCNVs) constitute significant proportions of genomic diversities within natural populations. However, it has been unclear whether GCNVs are under positive selection and contribute to adaptive evolution. Parallel evolution refers to adaptive evolution of the same trait in related but independent lineages, and three-spined stickleback (Gasterosteus aculeatus) is a well-known model organism. Through identification of genetic variations under parallel selection, i.e., variations shared among related but independent lineages, evidence of positive selection is obtained. In this study, we investigated whole-genome resequencing data from the marine and freshwater groups of three-spined sticklebacks from diverse areas along the Pacific and Atlantic Ocean coastlines, and searched for GCNVs under parallel selection.ResultsWe identified 24 GCNVs that showed significant differences in the numbers of mapped reads between the two groups, and this number was significantly larger than that expected by chance. The derived group, i.e., freshwater group, was typically characterized by larger gene-copy numbers, which implied that gene duplications or multiplications helped with adaptation to the freshwater environment. Some of the identified GCNVs were those of multigenic family genes, which is consistent with the theory that fatal effects due to copy-number changes of multigenic family genes tend to be less than those of single-copy genes.ConclusionThe identification of GCNVs that were likely under parallel selection suggests that contribution of GCNVs should be considered in studies on adaptive evolution.

Flexible selection of diversified Na+/K+-ATPase α-subunit isoforms for osmoregulation in teleosts

Background and methodsMultiple Na+/K+-ATPase (NKA) α-subunit isoforms express differentially in response to salinity transfer in teleosts but we observed that the isoform nomenclature is inconsistent with the phylogenetic relationship of NKA α-genes. We cloned the catalytic NKA α-subunit isoforms in eels and medaka, analyzed the time course of their expressions in osmoregulatory tissues after transfer from freshwater (FW) to seawater (SW), and performed phylogenetic analyses to deduce an evolutionary scenario that illustrates how various duplication events have led to the current genomic arrangement of NKA α-genes in teleosts.Results and discussionFive and six α-subunits were cloned in eels and medaka respectively. In eels, the commonly-reported α1a and α1b isoforms were absent while the α1c isoform was diversified instead (α1c-1, α1c-2, α1c-3, α2, and α3 in eels). Phylogenetic estimation indicated that the α1a and α1b isoforms from salmon, tilapia, and medaka were generated by independent duplication events and thus they are paralogous isoforms. Re-examination of expression changes of known isoforms after salinity challenge revealed that the isoforms selected as predominant SW-types varied among teleost lineages. Diversification of α1 isoforms occurred by various types of gene duplication, or by alternative transcription among tandem genes to form chimeric transcripts, but there is no trend for more α1 copies in euryhaline species. Our data suggest that the isoform switching between FW (α1a predominates) and SW (α1b predominates) that occurs in salmonids is not universal in teleosts. Instead, in eels, α1c-1 was the major α-subunit upregulated gill, intestine, and kidney in SW. Localization of both NKA mRNA and protein showed consistent upregulation in gill and intestine in SW eels, but not in renal distal and collecting tubules, where low transcript expression levels were accompanied by high protein levels, suggesting a tissue-specific translational regulation that determines and fine-tunes the NKA expression. In medaka, α1b was upregulated in SW in anterior intestine while most other α-subunit isoforms were less responsive to salinity changes.ConclusionBy integrating gene expression and phylogenetic results, we propose that the major NKA α-subunits for SW acclimation were not ancestrally selected, but rather were flexibly determined in lineage-specific fashion in teleosts.

MOCCS: Clarifying DNA-binding motif ambiguity using ChIP-Seq data

BACKGROUND As a key mechanism of gene regulation, transcription factors (TFs) bind to DNA by recognizing specific short sequence patterns that are called DNA-binding motifs. A single TF can accept ambiguity within its DNA-binding motifs, which comprise both canonical (typical) and non-canonical motifs. Clarification of such DNA-binding motif ambiguity is crucial for revealing gene regulatory networks and evaluating mutations in cis-regulatory elements. Although chromatin immunoprecipitation sequencing (ChIP-seq) now provides abundant data on the genomic sequences to which a given TF binds, existing motif discovery methods are unable to directly answer whether a given TF can bind to a specific DNA-binding motif. RESULTS Here, we report a method for clarifying the DNA-binding motif ambiguity, MOCCS. Given ChIP-Seq data of any TF, MOCCS comprehensively analyzes and describes every k-mer to which that TF binds. Analysis of simulated datasets revealed that MOCCS is applicable to various ChIP-Seq datasets, requiring only a few minutes per dataset. Application to the ENCODE ChIP-Seq datasets proved that MOCCS directly evaluates whether a given TF binds to each DNA-binding motif, even if known position weight matrix models do not provide sufficient information on DNA-binding motif ambiguity. Furthermore, users are not required to provide numerous parameters or background genomic sequence models that are typically unavailable. MOCCS is implemented in Perl and R and is freely available via https://github.com/yuifu/moccs. CONCLUSIONS By complementing existing motif-discovery software, MOCCS will contribute to the basic understanding of how the genome controls diverse cellular processes via DNA-protein interactions.

A sodium binding system alleviates acute salt stress during seawater acclimation in eels

BackgroundTeleosts transiting from freshwater (FW) to seawater (SW) environments face an immediate osmotic stress from ion influxes and water loss, but some euryhaline species such as eels can maintain a stable plasma osmolality during early SW exposure. The time course changes in the gene expression, protein abundance, and localization of key ion transporters suggested that the reversal of the ion transport systems was gradual, and we investigate how eels utilize a Na-binding strategy to slow down the ion invasion and complement the transporter-mediated osmoregulation.ResultsUsing an electron probe micro-analyzer, we localized bound Na in various eel tissues in response to SW transfer, suggesting that the Na-binding molecules were produced to sequester excess ionic Na+ to negate its osmotic potential, thus preventing acute cellular dehydration. Mucus cells were acutely activated in digestive tract, gill, and skin after SW transfer, producing Na-binding molecule-containing mucus layers that fence off high osmolality of SW. Using gel filtration HPLC, some molecules at 18 kDa were found to bind Na in the luminal secretion of esophagus and intestine, and higher binding was associated with SW transfer. Transcriptome and protein interaction results indicated that downregulation of Notch and β-catenin pathways, and dynamic changes in TGFβ pathways in intestine were involved during early SW transition, supporting the observed histological changes on epithelial desquamation and increased mucus production.ConclusionsThe timing for the activation of the Na-binding mechanism to alleviate the adverse osmotic gradient was temporally complementary to the subsequent remodeling of branchial ionocytes and transporting epithelia of the digestive tract. The strategy to manipulate the osmotic potential of Na+ by specific binding molecules is similar to the osmotically inactive Na described in human skin and muscle. The Na-binding molecules provide a buffer to tolerate the salinity changes, which is advantageous to the estuary and migrating fishes. Our data pave the way to identify this unknown class of molecules and open a new area of vertebrate osmoregulation research.