Search

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter., Histone variant H2A.Z has been suggested to contribute to the regulation of promoter accessibility. Here, the authors present high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for human Pol II promoters and provide evidence that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

Ageing is a degenerative process that leads to tissue dysfunction and death. A proposed cause of ageing is the accumulation of epigenetic noise that disrupts gene expression patterns, leading to decreases in tissue function and regenerative capacity1–3. Changes to DNA methylation patterns over time form the basis of ageing clocks4, but whether older individuals retain the information needed to restore these patterns—and, if so, whether this could improve tissue function—is not known. Over time, the central nervous system (CNS) loses function and regenerative capacity5–7. Using the eye as a model CNS tissue, here we show that ectopic expression of Oct4 (also known as Pou5f1), Sox2 and Klf4 genes (OSK) in mouse retinal ganglion cells restores youthful DNA methylation patterns and transcriptomes, promotes axon regeneration after injury, and reverses vision loss in a mouse model of glaucoma and in aged mice. The beneficial effects of OSK-induced reprogramming in axon regeneration and vision require the DNA demethylases TET1 and TET2. These data indicate that mammalian tissues retain a record of youthful epigenetic information—encoded in part by DNA methylation—that can be accessed to improve tissue function and promote regeneration in vivo. Expression of three Yamanaka transcription factors in mouse retinal ganglion cells restores youthful DNA methylation patterns, promotes axon regeneration after injury, and reverses vision loss in a mouse model of glaucoma and in aged mice, suggesting that mammalian tissues retain a record of youthful epigenetic information that can be accessed to improve tissue function.

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2′-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

Epigenetic “clocks” based on DNA methylation (DNAme) have emerged as the most robust and widely employed aging biomarkers, but conventional methods for applying them are expensive and laborious. Here, we developTagmentation-based Indexing forMethylationSequencing (TIME-Seq), a highly multiplexed and scalable method for low-cost epigenetic clocks. Using TIME-Seq, we applied multi-tissue and tissue-specific epigenetic clocks to over 1,600 mouse DNA samples. We also discovered a novel approach for age prediction from shallow sequencing (e.g., 10,000 reads) by adaptingscAgefor bulk measurements. In benchmarking experiments, TIME-Seq performed favorably against prevailing methods and could quantify the effects of interventions thought to accelerate, slow, and reverse aging in mice. Finally, we built and validated a highly accurate human blood clock from 1,056 demographically representative individuals. Our methods increase the scalability and reduce the cost of epigenetic age predictions by more than 100-fold, enabling accurate aging biomarkers to be applied in more large-scale animal and human studies.

Elucidating the transcriptional regulatory networks that underlie growth and development requires robust ways to define the complete set of transcription factor (TF) binding sites. Although TF-binding sites are known to be generally located within accessible chromatin regions (ACRs), pinpointing these DNA regulatory elements globally remains challenging. Current approaches primarily identify binding sites for a single TF (e.g. ChIP-seq), or globally detect ACRs but lack the resolution to consistently define TF-binding sites (e.g. DNAse-seq, ATAC-seq). To address this challenge, we developed MNase-defined cistrome-Occupancy Analysis (MOA-seq), a high-resolution (< 30 bp), high-throughput, and genome-wide strategy to globally identify putative TF-binding sites within ACRs. We used MOA-seq on developing maize ears as a proof of concept, able to define a cistrome of 145,000 MOA footprints (MFs). While a substantial majority (76%) of the known ATAC-seq ACRs intersected with the MFs, only a minority of MFs overlapped with the ATAC peaks, indicating that the majority of MFs were novel and not detected by ATAC-seq. MFs were associated with promoters and significantly enriched for TF-binding and long-range chromatin interaction sites, including for the well-characterized FASCIATED EAR4, KNOTTED1, and TEOSINTE BRANCHED1. Importantly, the MOA-seq strategy improved the spatial resolution of TF-binding prediction and allowed us to identify 215 motif families collectively distributed over more than 100,000 non-overlapping, putatively-occupied binding sites across the genome. Our study presents a simple, efficient, and high-resolution approach to identify putative TF footprints and binding motifs genome-wide, to ultimately define a native cistrome atlas., Author summary Understanding gene regulation remains a central goal of modern biology. Delineating the full set of regulatory DNA elements that orchestrate this regulation requires information at two scales; the broad landscape of accessible chromatin, and the site-specific binding of transcription factors (TFs) at discrete cis-regulatory DNA elements. Here we describe a single assay that uses micrococcal nuclease (MNase) as a structural probe to simultaneously reveal regions of accessible chromatin in addition to high-resolution footprints with signatures of TF-occupied cis-elements. We have used maize developing ear tissue as proof of concept, showing the method detects known TF-binding sites. This genome-wide assay not only defines chromatin landscapes, but crucially enables global discovery and mapping of sequence motifs underlying small footprints of ~30 bp to produce an atlas of candidate TF occupancy.

All living things experience entropy, manifested as a loss of inherited genetic and epigenetic information over time. In yeast, epigenetic changes result in a loss of cell identity and sterility, both hallmarks of yeast aging. In mammals, epigenetic information is also lost over time, but what causes it to be lost and whether it is a cause or a consequence of aging is not known. Using a transgenic mouse system called "ICE" (for Inducible Changes to the Epigenome), we show that the process of repairing non-mutagenic DNA breaks accelerates age-related physiological, cognitive, and molecular changes, including the erosion of the epigenetic landscape, a loss of cellular identity, cellular senescence and advancement of the epigenetic clock. Epigenetic reprogramming through ectopic expression of Oct4, Sox2 and Klf4 (OSK) restores patterns of youthful gene expression. These data support a model in which a loss of epigenetic information is a cause of aging in mammals.

In this issue of EMBO Molecular Medicine, Stebbing et al (2020b) validate an artificial intelligence‐assisted prediction that a drug used to treat rheumatoid arthritis could be a potent weapon against COVID‐19. Using liver organoids infected with SARS‐CoV‐2, they confirm dual antiviral and anti‐inflammatory activities and show that its administration in four COVID‐19 patients is correlated with disease improvement, paving the way for more rigorous placebo‐controlled trials., As the COVID‐19 pandemic is demonstrating, adversity fosters innovation. Here, the marriage of machine learning and rapid clinical trials provides hope for progress.

Background The functional genome of agronomically important plant species remains largely unexplored, yet presents a virtually untapped resource for targeted crop improvement. Functional elements of regulatory DNA revealed through profiles of chromatin accessibility can be harnessed for fine-tuning gene expression to optimal phenotypes in specific environments. Result Here, we investigate the non-coding regulatory space in the maize (Zea mays) genome during early reproductive development of pollen- and grain-bearing inflorescences. Using an assay for differential sensitivity of chromatin to micrococcal nuclease (MNase) digestion, we profile accessible chromatin and nucleosome occupancy in these largely undifferentiated tissues and classify at least 1.6% of the genome as accessible, with the majority of MNase hypersensitive sites marking proximal promoters, but also 3′ ends of maize genes. This approach maps regulatory elements to footprint-level resolution. Integration of complementary transcriptome profiles and transcription factor occupancy data are used to annotate regulatory factors, such as combinatorial transcription factor binding motifs and long non-coding RNAs, that potentially contribute to organogenesis, including tissue-specific regulation between male and female inflorescence structures. Finally, genome-wide association studies for inflorescence architecture traits based solely on functional regions delineated by MNase hypersensitivity reveals new SNP-trait associations in known regulators of inflorescence development as well as new candidates. Conclusions These analyses provide a comprehensive look into the cis-regulatory landscape during inflorescence differentiation in a major cereal crop, which ultimately shapes architecture and influences yield potential.

Whole-genome sequencing and annotation have allowed planetary protection engineers to assess the functional capabilities of microorganisms isolated from spacecraft hardware and associated surfaces. Here, we report draft genomes of six strains isolated from the InSight mission, determined using Oxford Nanopore- and Illumina-based sequencing.

Summary Zika virus (ZIKV) and dengue virus (DENV) are two closely related flaviviruses that lead to different clinical outcomes. The mechanism for the distinct pathogenesis of ZIKV and DENV is poorly understood. Here, we investigate ZIKV and DENV infection of macrophages using a human pluripotent stem cell (hPSC)-derived macrophage model and discover key virus-specific responses. ZIKV and DENV productively infect hPSC-derived macrophages. DENV, but not ZIKV, infection of macrophages strongly activates macrophage migration inhibitory factor (MIF) secretion and decreases macrophage migration. Neutralization of MIF leads to improved migratory ability of DENV-infected macrophages. In contrast, ZIKV-infected macrophages exhibit prolonged migration and express low levels of pro-inflammatory cytokines and chemokines. Mechanistically, ZIKV disrupts the nuclear factor κB (NF-κB)-MIF positive feedback loop by inhibiting the NF-κB signaling pathway. Our results demonstrate the utility of hPSC-derived macrophages in infectious disease modeling and suggest that the distinct impact of ZIKV and DENV on macrophage immune response may underlie different pathogenesis of Zika and dengue diseases., Highlights • An hPSC-derived tissue-resident macrophage model for ZIKV and DENV infection • ZIKV-, but not DENV-, infected macrophages maintain migratory capacity • ZIKV, but not DENV, inhibits pro-inflammatory cytokines and chemokines expression • ZIKV disrupts NF-κB-MIF positive feedback loop by inhibiting NF-κB pathway, In this article, Tang and colleagues demonstrate the utility of hPSC-derived tissue-resident macrophages in infectious disease modeling and show differential responses of macrophages to ZIKV and DENV infection. ZIKV-, but not DENV-, infected macrophages exhibit prolonged migration and express low levels of pro-inflammatory cytokines and chemokines by disrupting the NF-κB-MIF positive feedback loop via inhibition of the NF-κB signaling pathway.

This protocol is an extension to: Nat. Protoc. 6, 870-895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.

The prairie vole (Microtus ochrogaster) is an important model for the study of social monogamy and dual parental care of offspring. Characterization of specific host species-microbe strain interactions is critical for understanding the effects of the microbiota on mood and behavior. The five metagenome-assembled genome sequences reported here represent an important step in defining the prairie vole microbiome.

SUMMARYThere are numerous hallmarks of aging in mammals, but no unifying cause has been identified. In budding yeast, aging is associated with a loss of epigenetic information that occurs in response to genome instability, particularly DNA double-strand breaks (DSBs). Mammals also undergo predictable epigenetic changes with age, including alterations to DNA methylation patterns that serve as epigenetic “age” clocks, but what drives these changes is not known. Using a transgenic mouse system called “ICE” (for induciblechanges to theepigenome), we show that a tissue’s response to non-mutagenic DSBs reorganizes the epigenome and accelerates physiological, cognitive, and molecular changes normally seen in older mice, including advancement of the epigenetic clock. These findings implicate DSB-induced epigenetic drift as a conserved cause of aging from yeast to mammals.One Sentence SummaryDNA breaks induce epigenomic changes that accelerate the aging clock in mammals

Ageing is a degenerative process leading to tissue dysfunction and death. A proposed cause of ageing is the accumulation of epigenetic noise, which disrupts youthful gene expression patterns that are required for cells to function optimally and recover from damage1–3. Changes to DNA methylation patterns over time form the basis of an ‘ageing clock’4, 5, but whether old individuals retain information to reset the clock and, if so, whether this would improve tissue function is not known. Of all the tissues in the body, the central nervous system (CNS) is one of the first to lose regenerative capacity6, 7. Using the eye as a model tissue, we show that expression of Oct4, Sox2, and Klf4 genes (OSK) in mice resets youthful gene expression patterns and the DNA methylation age of retinal ganglion cells, promotes axon regeneration after optic nerve crush injury, and restores vision in a mouse model of glaucoma and in normal old mice. This process, which we call recovery of information via epigenetic reprogramming or REVIVER, requires the DNA demethylases Tet1 and Tet2, indicating that DNA methylation patterns don’t just indicate age, they participate in ageing. Thus, old tissues retain a faithful record of youthful epigenetic information that can be accessed for functional age reversal.

Summary The development of dengue antivirals and vaccine has been hampered by the incomplete understanding of molecular mechanisms of dengue virus (DENV) infection and pathology, partly due to the limited suitable cell culture or animal models that can capture the comprehensive cellular changes induced by DENV. In this study, we differentiated human pluripotent stem cells (hPSCs) into hepatocytes, one of the target cells of DENV, to investigate various aspects of DENV-hepatocyte interaction. hPSC-derived hepatocyte-like cells (HLCs) supported persistent and productive DENV infection. The activation of interferon pathways by DENV protected bystander cells from infection and protected the infected cells from massive apoptosis. Furthermore, DENV infection activated the NF-κB pathway, which led to production of proinflammatory cytokines and downregulated many liver-specific genes such as albumin and coagulation factor V. Our study demonstrates the utility of hPSC-derived hepatocytes as an in vitro model for DENV infection and reveals important aspects of DENV-host interactions., Graphical Abstract, Highlights • Hepatocyte-like cells (HLCs) derived from hPSCs support productive DENV infection • DENV infection of HLCs leads to cellular responses that are relevant to pathology • HLCs represent a new model system to study dengue virus infection of liver, Infection/disease models modeling the comprehensive cellular changes induced by DENV infection are limited. Here, Tang and colleagues use pluripotent stem cells to produce a renewable cell model that captures in vitro phenotypes relevant for dengue pathogenesis. Combining targeted differentiation and viral infection represents a new model to study the interface of host genetics and dengue infection and pathogenesis.

SUMMARYAll living things experience entropy, manifested as a loss of inherited genetic and epigenetic information over time. As budding yeast cells age, epigenetic changes result in a loss of cell identity and sterility, both hallmarks of yeast aging. In mammals, epigenetic information is also lost over time, but what causes it to be lost and whether it is a cause or a consequence of aging is not known. Here we show that the transient induction of genomic instability, in the form of a low number of non-mutagenic DNA breaks, accelerates many of the chromatin and tissue changes seen during aging, including the erosion of the epigenetic landscape, a loss of cellular identity, advancement of the DNA methylation clock and cellular senescence. These data support a model in which a loss of epigenetic information is a cause of aging in mammals.One Sentence SummaryThe act of repairing DNA breaks induces chromatin reorganization and a loss of cell identity that may contribute to mammalian aging

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program. RT is regulated during development in discrete chromosomal units, coordinated with transcriptional activity and 3D genome organization. Here, we derived distinct cell types from F1 hybrid musculus × castaneus mouse crosses and exploited the high single-nucleotide polymorphism (SNP) density to characterize allelic differences in RT (Repli-seq), genome organization (Hi-C and promoter-capture Hi-C), gene expression (total nuclear RNA-seq), and chromatin accessibility (ATAC-seq). We also present HARP, a new computational tool for sorting SNPs in phased genomes to efficiently measure allele-specific genome-wide data. Analysis of six different hybrid mESC clones with different genomes (C57BL/6, 129/sv, and CAST/Ei), parental configurations, and gender revealed significant RT asynchrony between alleles across ∼12% of the autosomal genome linked to subspecies genomes but not to parental origin, growth conditions, or gender. RT asynchrony in mESCs strongly correlated with changes in Hi-C compartments between alleles but not as strongly with SNP density, gene expression, imprinting, or chromatin accessibility. We then tracked mESC RT asynchronous regions during development by analyzing differentiated cell types, including extraembryonic endoderm stem (XEN) cells, four male and female primary mouse embryonic fibroblasts (MEFs), and neural precursor cells (NPCs) differentiated in vitro from mESCs with opposite parental configurations. We found that RT asynchrony and allelic discordance in Hi-C compartments seen in mESCs were largely lost in all differentiated cell types, accompanied by novel sites of allelic asynchrony at a considerably smaller proportion of the genome, suggesting that genome organization of homologs converges to similar folding patterns during cell fate commitment.

SUMMARYThe temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, butcis-elements regulating spatio-temporal control of replication have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discretecis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.Highlightscis-elements (ERCEs) regulate large scale chromosome structure and functionMultiple ERCEs cooperatively control domain-wide replicationERCEs harbor prominent active chromatin features and form CTCF-independent loopsERCEs enable genetic dissection of large-scale chromosome structure-function.

// Brittany S. Sexton 1,2 , Brooke R. Druliner 1,3 , Daniel L. Vera 1,4 , Denis Avey 1 , Fanxiu Zhu 1 , and Jonathan H. Dennis 1,4,5 1 Department of Biological Science, The Florida State University, Tallahassee, FL, United States of America 2 Current address: Massachusetts General Hospital Cancer Center and Department of Medicine, Harvard Medical School, Charlestown, MA, United States of America 3 Current address: Division of Gastroenterology and Hepatology, Mayo Clinic, Rochester, MN, United States of America 4 The Center for Genomics and Personalized Medicine The Florida State University, Tallahassee, FL, United States of America 5 Institute of Molecular Biophysics, The Florida State University, Tallahassee, FL, United States of America Correspondence to: Jonathan H. Dennis, email: // Keywords : nucleosome, chromatin, iSLK.219, KSHV, next generation sequencing, Chromosome Section Received : December 22, 2015 Accepted : December 28, 2015 Published : January 07, 2016 Abstract Nucleosome occupancy is critically important in regulating access to the eukaryotic genome. Few studies in human cells have measured genome-wide nucleosome distributions at high temporal resolution during a response to a common stimulus. We measured nucleosome distributions at high temporal resolution following Kaposi’s-sarcoma-associated herpesvirus (KSHV) reactivation using our newly developed mTSS-seq technology, which maps nucleosome distribution at the transcription start sites (TSS) of all human genes. Nucleosomes underwent widespread changes in organization 24 hours after KSHV reactivation and returned to their basal nucleosomal architecture 48 hours after KSHV reactivation. The widespread changes consisted of an indiscriminate remodeling event resulting in the loss of nucleosome rotational phasing signals. Additionally, one in six TSSs in the human genome possessed nucleosomes that are translationally remodeled. 72% of the loci with translationally remodeled nucleosomes have nucleosomes that moved to positions encoded by the underlying DNA sequence. Finally we demonstrated that these widespread alterations in nucleosomal architecture potentiated regulatory factor binding. These descriptions of nucleosomal architecture changes provide a new framework for understanding the role of chromatin in the genomic response, and have allowed us to propose a hierarchical model for chromatin-based regulation of genome response.

Mammalian genomes are partitioned into domains that replicate in a defined temporal order. These domains can replicate at similar times in all cell types (constitutive) or at cell type-specific times (developmental). Genome-wide chromatin conformation capture (Hi-C) has revealed sub-megabase topologically associating domains (TADs), which are the structural counterparts of replication domains. Hi-C also segregates inter-TAD contacts into defined 3D spatial compartments that align precisely to genome-wide replication timing profiles. Determinants of the replication-timing program are re-established during early G1 phase of each cell cycle and lost in G2 phase, but it is not known when TAD structure and inter-TAD contacts are re-established after their elimination during mitosis. Here, we use multiplexed 4C-seq to study dynamic changes in chromatin organization during early G1. We find that both establishment of TADs and their compartmentalization occur during early G1, within the same time frame as establishment of the replication-timing program. Once established, this 3D organization is preserved either after withdrawal into quiescence or for the remainder of interphase including G2 phase, implying 3D structure is not sufficient to maintain replication timing. Finally, we find that developmental domains are less well compartmentalized than constitutive domains and display chromatin properties that distinguish them from early and late constitutive domains. Overall, this study uncovers a strong connection between chromatin re-organization during G1, establishment of replication timing, and its developmental control.

The temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, but cis-elements regulating spatio-temporal control of replicatio have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discrete cis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.

Robustness in biology is the stability of phenotype under diverse genetic and/or environmental perturbations. The circadian clock has remarkable stability of period and phase that-unlike other biological oscillators-is maintained over a wide range of conditions. Here, we show that the high fidelity of the circadian system stems from robust degradation of the clock protein PERIOD. We show that PERIOD degradation is regulated by a balance between ubiquitination and deubiquitination, and that disruption of this balance can destabilize the clock. In mice with a loss-of-function mutation of the E3 ligase gene β-Trcp2, the balance of PERIOD degradation is perturbed and the clock becomes dramatically unstable, presenting a unique behavioral phenotype unlike other circadian mutant animal models. We believe that our data provide a molecular explanation for how circadian phases, such as wake-sleep onset times, can become unstable in humans, and we present a unique mouse model to study human circadian disorders with unstable circadian rhythm phases.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Background Identification of functional elements of a genome often requires dividing a sequence of measurements along a genome into segments where adjacent segments have different properties, such as different mean values. Despite dozens of algorithms developed to address this problem in genomics research, methods with improved accuracy and speed are still needed to effectively tackle both existing and emerging genomic and epigenomic segmentation problems. Results We designed an efficient algorithm, called iSeg, for segmentation of genomic and epigenomic profiles. iSeg first utilizes dynamic programming to identify candidate segments and test for significance. It then uses a novel data structure based on two coupled balanced binary trees to detect overlapping significant segments and update them simultaneously during searching and refinement stages. Refinement and merging of significant segments are performed at the end to generate the final set of segments. By using an objective function based on the p-values of the segments, the algorithm can serve as a general computational framework to be combined with different assumptions on the distributions of the data. As a general segmentation method, it can segment different types of genomic and epigenomic data, such as DNA copy number variation, nucleosome occupancy, nuclease sensitivity, and differential nuclease sensitivity data. Using simple t-tests to compute p-values across multiple datasets of different types, we evaluate iSeg using both simulated and experimental datasets and show that it performs satisfactorily when compared with some other popular methods, which often employ more sophisticated statistical models. Implemented in C++, iSeg is also very computationally efficient, well suited for large numbers of input profiles and data with very long sequences. Conclusions We have developed an efficient general-purpose segmentation tool and showed that it had comparable or more accurate results than many of the most popular segment-calling algorithms used in contemporary genomic data analysis. iSeg is capable of analyzing datasets that have both positive and negative values. Tunable parameters allow users to readily adjust the statistical stringency to best match the biological nature of individual datasets, including widely or sparsely mapped genomic datasets or those with non-normal distributions. Electronic supplementary material The online version of this article (10.1186/s12859-018-2140-3) contains supplementary material, which is available to authorized users.

MotivationSequencing-based methods to examine fundamental features of the genome, such as gene expression and chromatin structure, rely on inferences from the abundance and distribution of reads derived from Illumina sequencing. Drawing sound inferences from such experiments relies on appropriate mathematical methods to model the distribution of reads along the genome, which has been challenging due to the scale and nature of these data.ResultsWe propose a new framework (SRSFseq) based on Square Root Slope Functions shape analysis to analyse Illumina sequencing data. In the new approach the basic unit of information is the density of mapped reads over region of interest located on the known reference genome. The densities are interpreted as shapes and a new shape analysis model is proposed. An equivalent of a Fisher test is used to quantify the significance of shape differences in read distribution patterns between groups of density functions in different experimental conditions. We evaluated the performance of this new framework to analyze RNA-seq data at the exon level, which enabled the detection of variation in read distributions and abundances between experimental conditions not detected by other methods. Thus, the method is a suitable supplement to the state of the are count based techniques. The variety of density representations and flexibility of mathematical design allow the model to be easily adapted to other data types or problems in which the distribution of reads is to be tested. The functional interpretation and SRSF phase-amplitude separation technique gives an efficient noise reduction procedure improving the sensitivity and specificity of the method.

The nucleosome is a fundamental structural and functional chromatin unit that affects nearly all DNA-templated events in eukaryotic genomes. It is also a biochemical substrate for higher order, cis-acting gene expression codes and the monomeric structural unit for chromatin packaging at multiple scales. To predict the nucleosome landscape of a model plant genome, we used a support vector machine computational algorithm trained on human chromatin to predict the nucleosome occupancy likelihood (NOL) across the maize (Zea mays) genome. Experimentally validated NOL plots provide a novel genomic annotation that highlights gene structures, repetitive elements, and chromosome-scale domains likely to reflect regional gene density. We established a new genome browser (http://www.genomaize.org) for viewing support vector machine-based NOL scores. This annotation provides sequence-based comprehensive coverage across the entire genome, including repetitive genomic regions typically excluded from experimental genomics data. We find that transposable elements often displayed family-specific NOL profiles that included distinct regions, especially near their termini, predicted to have strong affinities for nucleosomes. We examined transcription start site consensus NOL plots for maize gene sets and discovered that most maize genes display a typical +1 nucleosome positioning signal just downstream of the start site but not upstream. This overall lack of a –1 nucleosome positioning signal was also predicted by our method for Arabidopsis (Arabidopsis thaliana) genes and verified by additional analysis of previously published Arabidopsis MNase-Seq data, revealing a general feature of plant promoters. Our study advances plant chromatin research by defining the potential contribution of the DNA sequence to observed nucleosome positioning and provides an invariant baseline annotation against which other genomic data can be compared.

Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early and late replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and sub-nuclear position. Moreover, RT is regulated during development and is altered in disease. Exploring mechanisms linking RT to other cellular processes in normal and diseased cells will be facilitated by rapid and robust methods with which to measure RT genome wide. Here, we describe a rapid, robust and relatively inexpensive protocol to analyze genome-wide RT by next-generation sequencing (NGS). This protocol yields highly reproducible results across laboratories and platforms. We also provide computational pipelines for analysis, parsing phased genomes using single nucleotide polymorphisms (SNP) for analyzing RT allelic asynchrony, and for direct comparison to Repli-chip data obtained by analyzing nascent DNA by microarrays.

Cellular processes mediated through nuclear DNA must contend with chromatin. Chromatin structural assays can efficiently integrate information across diverse regulatory elements, revealing the functional noncoding genome. In this study, we use a differential nuclease sensitivity assay based on micrococcal nuclease (MNase) digestion to discover open chromatin regions in the maize genome. We find that maize MNase-hypersensitive (MNase HS) regions localize around active genes and within recombination hotspots, focusing biased gene conversion at their flanks. Although MNase HS regions map to less than 1% of the genome, they consistently explain a remarkably large amount (∼40%) of heritable phenotypic variance in diverse complex traits. MNase HS regions are therefore on par with coding sequences as annotations that demarcate the functional parts of the maize genome. These results imply that less than 3% of the maize genome (coding and MNase HS regions) may give rise to the overwhelming majority of phenotypic variation, greatly narrowing the scope of the functional genome.

Belle (Bel), the Drosophila homolog of the yeast DEAD-box RNA helicase DED1 and human DDX3, has been shown to be required for oogenesis and female fertility. Here we report a novel role of Bel in regulating the expression of transgenes. Abrogation of Bel by mutations or RNAi induces silencing of a variety of P-element-derived transgenes. This silencing effect depends on downregulation of their RNA levels. Our genetic studies have revealed that the RNA helicase Spindle-E (Spn-E), a nuage RNA helicase that plays a crucial role in regulating RNA processing and PIWI-interacting RNA (piRNA) biogenesis in germline cells, is required for loss-of-bel-induced transgene silencing. Conversely, Bel abrogation alleviates the nuage-protein mislocalization phenotype in spn-E mutants, suggesting a competitive relationship between these two RNA helicases. Additionally, disruption of the chromatin remodeling factor Mod(mdg4) or the microRNA biogenesis enzyme Dcr-1 also rescued the transgene-silencing phenotypes in bel mutants, suggesting the involvement of chromatin remodeling and microRNA biogenesis in loss-of-bel-induced transgene silencing. Finally we showed that genetic inhibition of Bel function led to the de novo generation of piRNAs from the transgene region inserted in the genome, suggesting a potential piRNA-dependent mechanism that may mediate transgene silencing as Bel function is inhibited. Our findings have demonstrated a novel involvement of Bel in regulating transgene expression and its loss triggers a transgene silencing mechanism mediated by protein regulators implicated in RNA processing, piRNA biogenesis, chromatin remodeling and the microRNA pathway.

// Brooke R. Druliner 1,5 , Daniel Vera 1,6 , Ruth Johnson 2 , Xiaoyang Ruan 3 , Lynn M. Apone 4 , Eileen T. Dimalanta 4 , Fiona J. Stewart 4 , Lisa Boardman 5 and Jonathan H. Dennis 1,6,7 1 Department of Biological Science, Florida State University, Tallahassee, Florida, United States of America 2 Department of Laboratory Medicine and Experimental Pathology, Mayo Clinic, Rochester, Minnesota, United States of America 3 Biomedical Statistics and Informatics, Department of Health Sciences Research, Mayo Clinic College of Medicine, Rochester, Minnesota, United States of America 4 New England Biolabs Inc., Ipswich, Massachusetts, United States of America 5 Division of Gastroenterology and Hepatology, Department of Internal Medicine, Mayo Clinic, Rochester, Minnesota, United States of America 6 The Center for Genomics and Personalized Medicine, The Florida State University, Tallahassee, Florida, United States of America 7 Institute of Molecular Biophysics, The Florida State University, Tallahassee, Florida, United States of America Correspondence to: Jonathan H. Dennis, email: // Keywords : cancer, chromatin, nucleosome, MNase, whole genome, Chromosome Section Received : December 18, 2015 Accepted : December 21, 2015 Published : December 31, 2015 Abstract Altered chromatin structure is a hallmark of cancer, and inappropriate regulation of chromatin structure may represent the origin of transformation. Important studies have mapped human nucleosome distributions genome wide, but the role of chromatin structure in cancer progression has not been addressed. We developed a MNase-Transcription Start Site Sequence Capture method (mTSS-seq) to map the nucleosome distribution at human transcription start sites genome-wide in primary human lung and colon adenocarcinoma tissue. Here, we confirm that nucleosome redistribution is an early, widespread event in lung (LAC) and colon (CRC) adenocarcinoma. These altered nucleosome architectures are consistent between LAC and CRC patient samples indicating that they may serve as important early adenocarcinoma markers. We demonstrate that the nucleosome alterations are driven by the underlying DNA sequence and potentiate transcription factor binding. We conclude that DNA-directed nucleosome redistributions are widespread early in cancer progression. We have proposed an entirely new hierarchical model for chromatin-mediated genome regulation.

Genome-wide DNA replication timing (RT) profiles reflect the global 3D chromosome architecture of cells. They also provide a comprehensive and unique megabase-scale picture of the cellular epigenetic state. Thus normal differentiation involves reproducible changes in RT and transformation generally perturbs these, although the potential effects of altered RT on the properties of transformed cells remain largely unknown. A major challenge to interrogating these issues in human acute lymphoid leukemia (ALL) is the low proliferative activity of most of the cells, which may be further reduced in cryopreserved samples and difficult to overcome in vitro. In contrast, the ability of many human ALL cell populations to expand when transplanted in highly immunodeficient mice is well documented. To examine the stability of DNA RT profiles of serially passaged xenografts of primary human B- and T-ALL cells, we first devised a method that circumvents the need for BrdU incorporation to distinguish early versus late S-phase cells. Using this and more standard protocols, we found consistent strong retention in xenografts of the original patient-specific RT features, for all 8 primary human ALL cases surveyed (7 B-ALLs and one T-ALL). Moreover, in a case where genomic analyses indicated changing subclonal dynamics in serial passages, the RT profiles tracked concordantly. These results show that DNA RT is a relatively stable feature of human ALLs propagated in immunodeficient mice. In addition, they suggest the power of this approach for future interrogation of the origin and consequences of altered DNA RT in these diseases.

The eukaryotic genome is organized into nucleosomes, the fundamental units of chromatin. The positions of nucleosomes on DNA regulate protein-DNA interactions and in turn influence DNA-templated events. Despite the increasing number of genome-wide maps of nucleosome position, how global changes in gene expression relate to changes in nucleosome position is poorly understood. We show that in nucleosome occupancy mapping experiments in maize (Zea mays), particular genomic regions are highly susceptible to variation introduced by differences in the extent to which chromatin is digested with micrococcal nuclease (MNase). We exploited this digestion-linked variation to identify protein footprints that are hypersensitive to MNase digestion, an approach we term differential nuclease sensitivity profiling (DNS-chip). Hypersensitive footprints were enriched at the 5′ and 3′ ends of genes, associated with gene expression levels, and significantly overlapped with conserved noncoding sequences and the binding sites of the transcription factor KNOTTED1. We also found that the tissue-specific regulation of gene expression was linked to tissue-specific hypersensitive footprints. These results reveal biochemical features of nucleosome organization that correlate with gene expression levels and colocalize with functional DNA elements. This approach to chromatin profiling should be broadly applicable to other species and should shed light on the relationships among chromatin organization, protein-DNA interactions, and genome regulation.

Nucleosome occupancy plays a key role in regulating access to eukaryotic genomes. Although various chromatin regulatory complexes are known to regulate nucleosome occupancy, the role of DNA sequence in this regulation remains unclear, particularly in mammals. To address this problem, we measured nucleosome distribution at high temporal resolution in human cells at hundreds of genes during the reactivation of Kaposi's sarcoma–associated herpesvirus (KSHV). We show that nucleosome redistribution peaks at 24 h post-KSHV reactivation and that the nucleosomal redistributions are widespread and transient. To clarify the role of DNA sequence in these nucleosomal redistributions, we compared the genes with altered nucleosome distribution to a sequence-based computer model and in vitro–assembled nucleosomes. We demonstrate that both the predicted model and the assembled nucleosome distributions are concordant with the majority of nucleosome redistributions at 24 h post-KSHV reactivation. We suggest a model in which loci are held in an unfavorable chromatin architecture and “spring” to a transient intermediate state directed by DNA sequence information. We propose that DNA sequence plays a more considerable role in the regulation of nucleosome positions than was previously appreciated. The surprising findings that nucleosome redistributions are widespread, transient, and DNA-directed shift the current perspective regarding regulation of nucleosome distribution in humans.

The development and progression of lung adenocarcinoma, one of the most common cancers, is driven by the interplay of genetic and epigenetic changes and the role of chromatin structure in malignant transformation remains poorly understood. We used systematic nucleosome distribution and chromatin accessibility microarray mapping platforms to analyze the genome-wide chromatin structure from normal tissues and from primary lung adenocarcinoma of different grades and stages. We identified chromatin-based patterns across different patients with lung adenocarcinoma of different cancer grade and stage. Low-grade cancers had nucleosome distributions very different compared with the corresponding normal tissue but had nearly identical chromatin accessibility. Conversely, nucleosome distributions of high-grade cancers showed few differences. Substantial disruptions in chromosomal accessibility were seen in a patient with a high-grade and high-stage tumor. These data imply that chromatin structure changes during the progression of lung adenocarcinoma. We have therefore developed a model in which low-grade lung adenocarcinomas are linked to changes in nucleosome distributions, whereas higher-grade tumors are linked to large-scale chromosomal changes. These results provide a foundation for the development of a comprehensive framework linking the general and locus-specific roles of chromatin structure to lung cancer progression. We propose that this strategy has the potential to identify a new class of chromatin-based diagnostic, prognostic and therapeutic markers in cancer progression.

The visual system of a particular species is highly adapted to convey detailed ecological and behavioral information essential for survival. The consequences of structural mutations of opsins upon spectral sensitivity and environmental adaptation have been studied in great detail, but lacking is knowledge of the potential influence of alterations in gene regulatory networks upon the diversity of cone subtypes and the variation in the ratio of rods and cones observed in numerous diurnal and nocturnal species. Exploiting photoreceptor patterning in cone-dominated zebrafish, we uncovered two independent mechanisms by which the sine oculis homeobox homolog 7 (six7) regulates photoreceptor development. In a genetic screen, we isolated the lots-of-rods-junior (ljrp23ahub) mutation that resulted in an increased number and uniform distribution of rods in otherwise normal appearing larvae. Sequence analysis, genome editing using TALENs and knockdown strategies confirm ljrp23ahub as a hypomorphic allele of six7, a teleost orthologue of six3, with known roles in forebrain patterning and expression of opsins. Based on the lack of predicted protein-coding changes and a deletion of a conserved element upstream of the transcription start site, a cis-regulatory mutation is proposed as the basis of the reduced expression of six7 in ljrp23ahub. Comparison of the phenotypes of the hypomorphic and knock-out alleles provides evidence of two independent roles in photoreceptor development. EdU and PH3 labeling show that the increase in rod number is associated with extended mitosis of photoreceptor progenitors, and TUNEL suggests that the lack of green-sensitive cones is the result of cell death of the cone precursor. These data add six7 to the small but growing list of essential genes for specification and patterning of photoreceptors in non-mammalian vertebrates, and highlight alterations in transcriptional regulation as a potential source of photoreceptor variation across species., Author Summary Vision begins when an image is focused on the neural retina where rod and cone photoreceptors convert light into the electrical signals of the brain. The 4 cone subtypes in retinas of the majority of fishes, lizards and birds, provide rich color vision. In contrast, retinas of most mammals are better adapted for dim light conditions with rods vastly outnumbering the sparse and less diverse cone subtypes. However, our understanding of photoreceptor development largely based on findings from mammalian models fails to explain the tremendous diversity of cone subtypes and variation of rod and cone numbers across the majority of vertebrate species. Taking advantage of the cone-rich zebrafish retina, we identified a nuclear factor that suppresses the number of rods and is essential for the development of a cone subtype not present in mammals. Combined with prior studies, the findings provide insight into adaptive mechanisms underlying maintenance of a cone-dominated retina.

Telomere length is a quantitative trait important for many cellular functions. Failure to regulate telomere length contributes to genomic instability, cellular senescence, cancer, and apoptosis in humans, but the functional significance of telomere regulation in plants is much less well understood. To gain a better understanding of telomere biology in plants, we used quantitative trait locus (QTL) mapping to identify genetic elements that control telomere length variation in maize (Zea mays L.). For this purpose, we measured the median and mean telomere lengths from 178 recombinant inbred lines of the IBM mapping population and found multiple regions that collectively accounted for 33–38% of the variation in telomere length. Two-way analysis of variance revealed interaction between the quantitative trait loci at genetic bin positions 2.09 and 5.04. Candidate genes within these and other significant QTL intervals, along with select genes known a priori to regulate telomere length, were tested for correlations between expression levels and telomere length in the IBM population and diverse inbred lines by quantitative real-time PCR. A slight but significant positive correlation between expression levels and telomere length was observed for many of the candidate genes, but Ibp2 was a notable exception, showing instead a negative correlation. A rad51-like protein (TEL-MD_5.04) was strongly supported as a candidate gene by several lines of evidence. Our results highlight the value of QTL mapping plus candidate gene expression analysis in a genetically diverse model system for telomere research.