Search

In 2019 the National Science Foundation in the United States began operation of its Frontera supercomputer at the Texas Advanced Computing Center, debuting at number 5 on the TOP500 list. Unlike more general purpose supercomputers, Frontera was designed for the focused community of users doing computational science at the largest scale. This special issue features articles from Frontera's science users spanning the sciences from climate change to tornadoes and organic solar cells. We also include two computational science articles highlighting the tools and techniques that are needed to enable users running single jobs beyond 8,000 nodes.

Computational tools addressing various components of design–build–test–learn (DBTL) loops for the construction of synthetic genetic networks exist but do not generally cover the entire DBTL loop. This manuscript introduces an end-to-end sequence of tools that together form a DBTL loop called Design Assemble Round Trip (DART). DART provides rational selection and refinement of genetic parts to construct and test a circuit. Computational support for experimental process, metadata management, standardized data collection and reproducible data analysis is provided via the previously published Round Trip (RT) test–learn loop. The primary focus of this work is on the Design Assemble (DA) part of the tool chain, which improves on previous techniques by screening up to thousands of network topologies for robust performance using a novel robustness score derived from dynamical behavior based on circuit topology only. In addition, novel experimental support software is introduced for the assembly of genetic circuits. A complete design-through-analysis sequence is presented using several OR and NOR circuit designs, with and without structural redundancy, that are implemented in budding yeast. The execution of DART tested the predictions of the design tools, specifically with regard to robust and reproducible performance under different experimental conditions. The data analysis depended on a novel application of machine learning techniques to segment bimodal flow cytometry distributions. Evidence is presented that, in some cases, a more complex build may impart more robustness and reproducibility across experimental conditions. Graphical Abstract

We describe an experimental campaign that replicated the performance assessment of logic gates engineered into cells of S. cerevisiae by Gander, et al. Our experimental campaign used a novel high throughput experimentation framework developed under DARPA’s Synergistic Discovery and Design (SD2) program: a remote robotic lab at Strateos executed a parameterized experimental protocol. Using this protocol and robotic execution, we generated two orders of magnitude more flow cytometry data than the original experiments. We discuss our results, which largely, but not completely, agree with the original report, and make some remarks about lessons learned.

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.

Synthetic biology is a complex discipline that involves creating detailed, purpose-built designs from genetic parts. This process is often phrased as a Design-Build-Test-Learn loop, where iterative design improvements can be made, implemented, measured, and analyzed. Automation can potentially improve both the end-to-end duration of the process and the utility of data produced by the process. One of the most important considerations for the development of effective automation and quality data is a rigorous description of implicit knowledge encoded as a formal knowledge representation. The development of knowledge representation for the process poses a number of challenges, including developing effective human-machine interfaces, protecting against and repairing user error, providing flexibility for terminological mismatches, and supporting extensibility to new experimental types. We address these challenges with the DARPA SD2 Round Trip software architecture. The Round Trip is an open architecture that automates many of the key steps in the Test and Learn phases of a Design-Build-Test-Learn loop for high-throughput laboratory science. The primary contribution of the Round Trip is to assist with and otherwise automate metadata creation, curation, standardization, and linkage with experimental data. The Round Trip's focus on metadata supports fast, automated, and replicable analysis of experiments as well as experimental situational awareness and experimental interpretability. We highlight the major software components and data representations that enable the Round Trip to speed up the design and analysis of experiments by 2 orders of magnitude over prior

We describe an experimental campaign that replicated the performance assessment of logic gates engineered into cells of Saccharomyces cerevisiae by Gander et al. Our experimental campaign used a novel high-throughput experimentation framework developed under Defense Advanced Research Projects Agency’s Synergistic Discovery and Design program: a remote robotic lab at Strateos executed a parameterized experimental protocol. Using this protocol and robotic execution, we generated two orders of magnitude more flow cytometry data than the original experiments. We discuss our results, which largely, but not completely, agree with the original report and make some remarks about lessons learned. Graphical Abstract

Summary Improvements in next‐generation sequencing technologies have resulted in dramatically reduced sequencing costs. This has led to an explosion of ‘‐seq’‐based methods, of which RNA sequencing (RNA‐seq) for generating transcriptomic data is the most popular. By analysing global patterns of gene expression in organs/tissues/cells of interest or in response to chemical or environmental perturbations, researchers can better understand an organism's biology. Tools designed to work with large RNA‐seq data sets enable analyses and visualizations to help generate hypotheses about a gene's function. We present here a user‐friendly RNA‐seq data exploration tool, called the ‘eFP‐Seq Browser’, that shows the read map coverage of a gene of interest in each of the samples along with ‘electronic fluorescent pictographic’ (eFP) images that serve as visual representations of expression levels. The tool also summarizes the details of each RNA‐seq experiment, providing links to archival databases and publications. It automatically computes the reads per kilobase per million reads mapped expression‐level summaries and point biserial correlation scores to sort the samples based on a gene's expression level or by how dissimilar the read map profile is from a gene splice variant, to quickly identify samples with the strongest expression level or where alternative splicing might be occurring. Links to the Integrated Genome Browser desktop visualization tool allow researchers to visualize and explore the details of RNA‐seq alignments summarized in eFP‐Seq Browser as coverage graphs. We present four cases of use of the eFP‐Seq Browser for ABI3,SR34,SR45a and U2AF65B, where we examine expression levels and identify alternative splicing. The URL for the browser is https://bar.utoronto.ca/eFP-Seq_Browser/. Open research badges This article has earned an Open Data Badge for making publicly available the digitally‐shareable data necessary to reproduce the reported results. Tool is at http://sps:urlprefix::https; RNA‐seq data at http://sps:urlprefix::https and http://sps:urlprefix::https. Code is available at http://sps:urlprefix::https, Significance Statement We present a tool, the eFP‐Seq Browser, for rapidly identifying RNA sequencing samples with strong expression levels of a given gene, or where the read maps for a given gene/sample best match a particular gene model. Details can be called up with convenient links to the Integrated Genome Browser.

Jetstream2 will be a category I production cloud resource that is part of the National Science Foundation’s Innovative HPC Program. The project’s aim is to accelerate science and engineering by providing “on-demand” programmable infrastructure built around a core system at Indiana University and four regional sites. Jetstream2 is an evolution of the Jetstream platform, which functions primarily as an Infrastructure-as-a-Service cloud. The lessons learned in cloud architecture, distributed storage, and container orchestration have inspired changes in both hardware and software for Jetstream2. These lessons have wide implications as institutions converge HPC and cloud technology while building on prior work when deploying their own cloud environments. Jetstream2’s next-generation hardware, robust open-source software, and enhanced virtualization will provide a significant platform to further cloud adoption within the US research and education communities.

Engineered proteins generally must possess a stable structure in order to achieve their designed function. Stable designs, however, are astronomically rare within the space of all possible amino acid sequences. As a consequence, many designs must be tested computationally and experimentally in order to find stable ones, which is expensive in terms of time and resources. Here we report a neural network model that predicts protein stability based only on sequences of amino acids, and demonstrate its performance by evaluating the stability of almost 200,000 novel proteins. These include a wide range of sequence perturbations, providing a baseline for future work in the field. We also report a second neural network model that is able to generate novel stable proteins. Finally, we show that the predictive model can be used to substantially increase the stability of both expert-designed and model-generated proteins.

Containerization technologies provide a mechanism to encapsulate applications and many of their dependencies, facilitating software portability and reproducibility on HPC systems. However, in order to access many of the architectural features that enable HPC system performance, compatibility between certain components of the container and host is required, resulting in a trade-off between portability and performance. In this work, we discuss our experiences running three state-of-the-art containerization technologies on five leading petascale systems. We present how we build the containers to ensure performance and security and their performance at scale. We ran microbenchmarks at a scale of 6,144 nodes containing 0.35 M MPI processes and baseline the performance of container technologies. We establish the near-native performance and minimal memory overheads by the containerized environments using MILC - a lattice quantum chromodynamics code at 139,968 processes and using VPIC - a 3d electromagnetic relativistic Vector Particle-In-Cell code for modeling kinetic plasmas at 32,768 processes. We demonstrate an on-par performance trend at a large scale on Intel, AMD, and three NVIDIA architectures for both HPC applications.

Beginning with the initial release of the DesignSafe JupyterHub in late 2015, TACC has been building and maintaining custom JupyterHub clusters for research groups across different domains of science and engineering. Today, TACC maintains five production JupyterHub systems utilizing over half a terabyte of memory and hundreds of CPU cores supporting nearly 1,600 unique users combined. In this paper, we describe our approach to utilizing JupyterHub in these cyberinfrastructure projects and our collaborative approach to integrating Jupyter into different research communities. For two such groups, we present an in-depth discussion of the science use cases and technical drivers that informed and evolved the design of our offering, and our outreach and engagement efforts to promote adoption. We discuss the implementation of custom features we have added to better integrate JupyterHub into other components of the cyberinfrastructure platforms, and we conclude with a description of our plans for future architecture and usage of JupyterHub at TACC.

DNA methylation is a chromatin modification that can provide epigenetic regulation of gene and transposon expression. Plants utilize several pathways to establish and maintain DNA methylation in specific sequence contexts. The chromomethylase (CMT) genes maintain CHG (where H = A, C or T) methylation. The RNA-directed DNA methylation (RdDM) pathway is important for CHH methylation. Transcriptome analysis was performed in a collection of Zea mays lines carrying mutant alleles for CMT or RdDM-associated genes. While the majority of the transcriptome was not affected, we identified sets of genes and transposon families sensitive to context-specific decreases in DNA methylation in mutant lines. Many of the genes that are up-regulated in CMT mutant lines have high levels of CHG methylation, while genes that are differentially expressed in RdDM mutants are enriched for having nearby mCHH islands, implicating context-specific DNA methylation in the regulation of expression for a small number of genes. Many genes regulated by CMTs exhibit natural variation for DNA methylation and transcript abundance in a panel of diverse inbred lines. Transposon families with differential expression in the mutant genotypes show few defining features, though several families up-regulated in RdDM mutants show enriched expression in endosperm tissue, highlighting the potential importance for this pathway during reproduction. Taken together, our findings suggest that while the number of genes and transposon families whose expression is reproducibly affected by mild perturbations in context-specific methylation is small, there are distinct patterns for loci impacted by RdDM and CMT mutants.

Eukaryotes use a temporally regulated process, known as the replication timing program, to ensure that their genomes are fully and accurately duplicated during S phase. Replication timing programs are predictive of genomic features and activity and are considered to be functional readouts of chromatin organization. Although replication timing programs have been described for yeast and animal systems, much less is known about the temporal regulation of plant DNA replication or its relationship to genome sequence and chromatin structure. We used the thymidine analog, 5-ethynyl-2'-deoxyuridine, in combination with flow sorting and Repli-Seq to describe, at high-resolution, the genome-wide replication timing program for Arabidopsis (Arabidopsis thaliana) Col-0 suspension cells. We identified genomic regions that replicate predominantly during early, mid, and late S phase, and correlated these regions with genomic features and with data for chromatin state, accessibility, and long-distance interaction. Arabidopsis chromosome arms tend to replicate early while pericentromeric regions replicate late. Early and mid-replicating regions are gene-rich and predominantly euchromatic, while late regions are rich in transposable elements and primarily heterochromatic. However, the distribution of chromatin states across the different times is complex, with each replication time corresponding to a mixture of states. Early and mid-replicating sequences interact with each other and not with late sequences, but early regions are more accessible than mid regions. The replication timing program in Arabidopsis reflects a bipartite genomic organization with early/mid-replicating regions and late regions forming separate, noninteracting compartments. The temporal order of DNA replication within the early/mid compartment may be modulated largely by chromatin accessibility.

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed., Author summary In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

Epigenetic modification of DNA through methylation is known to be involved in multiple biological processes such as gene suppression. However, the exact mechanism of how DNA methylations play their part is yet unclear. In mammals, CpG islands (CGI) have been studied extensively for their involvement in cancer. Whereas in plants, despite the fact that there are not only CpG but also CHG and CHH contexts of methylation, an efficient and easy-to-use pipeline to decipher these phenomena is still to be developed. Both ZED-align and BisuKit are user-friendly apps deployed on CyVerse infrastructure where users can use their bisulfite sequence files to run multiple command line-based packages with minimal intervention. © 2016 by John Wiley & Sons, Inc.

Biomedical data are quickly growing in volume and in variety, providing clinicians an opportunity for better clinical decision support. Here, we demonstrate a robust platform that uses software automation and high performance computing (HPC) resources to achieve real-time analytics of clinical data, specifically magnetic resonance imaging (MRI) data. We used the Agave application programming interface (API) to facilitate communication, data transfer, and job control between an MRI scanner and an off-site HPC resource. In this use case, Agave executed the graphical pipeline tool GRAPE (GRAphical Pipeline Environment) to perform automated, real-time, quantitative analysis of MRI scans. Same-session image processing will open the door for adaptive scanning and real-time quality control, potentially accelerating the discovery of pathologies and minimizing patient callbacks. We envision this platform can be adapted to other medical instruments, HPC resources, and analytics tools.

DNA methylation is a chromatin modification that can provide epigenetic regulation of gene and transposon expression. Plants utilize several pathways to establish and maintain DNA methylation in specific sequence contexts. The chromomethylase (CMT) genes maintain CHG (where H = A, C or T) methylation. The RNA-directed DNA methylation (RdDM) pathway is important for CHH methylation. Transcriptome analysis was performed in a collection ofZea mayslines carrying mutant alleles for CMT or RdDM-associated genes. While the majority of the transcriptome was not affected, we identified sets of genes and transposon families sensitive to context-specific decreases in DNA methylation in mutant lines. Many of the genes that are up-regulated in CMT mutant lines have high levels of CHG methylation, while genes that are differentially expressed in RdDM mutants are enriched for having nearby mCHH islands, providing evidence that context-specific DNA methylation directly regulates expression of a small number of genes. The analysis of a diverse set of inbred lines revealed that many genes regulated by CMTs exhibit natural variation for DNA methylation and gene expression. Transposon families with differential expression in the mutant genotypes show few defining features, though several families up-regulated in RdDM mutants show enriched expression in endosperm, highlighting the importance for this pathway during reproduction. Taken together, our findings suggest that while the number of genes and transposon families whose expression is reproducibly affected by mild perturbations in context-specific methylation is small, there are distinct patterns for loci impacted by RdDM and CMT mutants.

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that include several different cell lineages, and present replication profiles from cells in early, mid, and late S phase of the mitotic cell cycle. This represents the first genomic analysis of replication timing in a crop species, and the first in any plant species to characterize replication at the whole-genome level. Maize root tips have a complex replication timing program, including regions of distinct early, mid and late-S replication that each constitute between 20-24% of the genome, as well as other loci corresponding to ~32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

The Arabidopsis Information Portal (https://www. araport.org) is a new online resource for plant biology research. It houses the Arabidopsis thaliana genome sequence and associated annotation. It was conceived as a framework that allows the research community to develop and release ‘modules’ that integrate, analyze and visualize Arabidopsis data that may reside at remote sites. The current implementation provides an indexed database of core genomic information. These data are made available through feature-rich web applications that provide search, data mining, and genome browser functionality, and also by bulk download and web services. Araport uses software from the InterMine and JBrowse projects to expose curated data from TAIR, GO, BAR, EBI, UniProt, PubMed and EPIC CoGe. The site also hosts ‘science apps,’ developed as prototypes for community modules that use dynamic web pages to present data obtained on-demand from third-party servers via RESTful web services. Designed for sustainability, the Arabidopsis Information Portal strategy exploits existing scientific computing infrastructure, adopts a practical mixture of data integration technologies and encourages collaborative enhancement of the resource by its user community.

In the USA, the national cyberinfrastructure refers to a system of research supercomputer and other IT facilities and the high speed networks that connect them. These resources have been heavily leveraged by scientists in disciplines such as high energy physics, astronomy, and climatology, but until recently they have been little used by biomedical researchers. We suggest that many of the ‘Big Data’ challenges facing the medical informatics community can be efficiently handled using national-scale cyberinfrastructure. Resources such as the Extreme Science and Discovery Environment, the Open Science Grid, and Internet2 provide economical and proven infrastructures for Big Data challenges, but these resources can be difficult to approach. Specialized web portals, support centers, and virtual organizations can be constructed on these resources to meet defined computational challenges, specifically for genomics. We provide examples of how this has been done in basic biology as an illustration for the biomedical informatics community.

Scaffold or matrix attachment regions (S/MARs) are found in all eukaryotes. The pattern of distribution and genomic context of S/MARs is thought to be important for processes such as chromatin organization and modulation of gene expression. Despite the importance of such processes, much is unknown about the large-scale distribution and sequence content of S/MARs in vivo. Here, we report the use of tiling microarrays to map 1358 S/MARs on Arabidopsis thaliana chromosome 4 (chr4). S/MARs occur throughout chr4, spaced much more closely than in the large plant and animal genomes that have been studied to date. Arabidopsis S/MARs can be divided into five clusters based on their association with other genomic features, suggesting a diversity of functions. While some Arabidopsis S/MARs may define structural domains, most occur near the transcription start sites of genes. Genes associated with these S/MARs have an increased probability of expression, which is particularly pronounced in the case of transcription factor genes. Analysis of sequence motifs and 6-mer enrichment patterns show that S/MARs are preferentially enriched in poly(dA:dT) tracts, sequences that resist nucleosome formation, and the majority of S/MARs contain at least one nucleosome-depleted region. This global view of S/MARs provides a framework to begin evaluating genome-scale models for S/MAR function.

In plants, a subset of genes exhibit imprinting in endosperm tissue such that expression is primarily from the maternal or paternal allele. Imprinting may arise as a consequence of mechanisms for silencing of transposons during reproduction, and in some cases imprinted expression of particular genes may provide a selective advantage such that it is conserved across species. Separate mechanisms for the origin of imprinted expression patterns and maintenance of these patterns may result in substantial variation in the targets of imprinting in different species. Here we present deep sequencing of RNAs isolated from reciprocal crosses of four diverse maize genotypes, providing a comprehensive analysis that allows evaluation of imprinting at more than 95% of endosperm-expressed genes. We find that over 500 genes exhibit statistically significant parent-of-origin effects in maize endosperm tissue, but focused our analyses on a subset of these genes that had >90% expression from the maternal allele (69 genes) or from the paternal allele (108 genes) in at least one reciprocal cross. Over 10% of imprinted genes show evidence of allelic variation for imprinting. A comparison of imprinting in maize and rice reveals that 13% of genes with syntenic orthologs in both species exhibit conserved imprinting. Genes that exhibit conserved imprinting between maize and rice have elevated nonsynonymous to synonymous substitution ratios compared with other imprinted genes, suggesting a history of more rapid evolution. Together, these data suggest that imprinting only has functional relevance at a subset of loci that currently exhibit imprinting in maize.

Trimethylation of histone H3 Lys-27 (H3K27me3) plays a critical role in regulating gene expression during plant and animal development. We characterized the genome-wide distribution of H3K27me3 in five developmentally distinct tissues in maize (Zea mays) plants of two genetic backgrounds, B73 and Mo17. There were more substantial differences in the genome-wide profile of H3K27me3 between different tissues than between the two genotypes. The tissue-specific patterns of H3K27me3 were often associated with differences in gene expression among the tissues and most of the imprinted genes that are expressed solely from the paternal allele in endosperm are targets of H3K27me3. A comparison of the H3K27me3 targets in rice (Oryza sativa), maize, and Arabidopsis thaliana provided evidence for conservation of the H3K27me3 targets among plant species. However, there was limited evidence for conserved targeting of H3K27me3 in the two maize subgenomes derived from whole-genome duplication, suggesting the potential for subfunctionalization of chromatin regulation of paralogs. Genomic profiling of H3K27me3 in loss-of-function mutant lines for Maize Enhancer of zeste-like2 (Mez2) and Mez3, two of the three putative H3K27me3 methyltransferases present in the maize genome, suggested partial redundancy of this gene family for maintaining H3K27me3 patterns. Only a portion of the targets of H3K27me3 required Mez2 and/or Mez3, and there was limited evidence for functional consequences of H3K27me3 at these targets.

ThaleMine (https://apps.araport.org/thalemine/) is a comprehensive data warehouse that integrates a wide array of genomic information of the model plant Arabidopsis thaliana. The data collection currently includes the latest structural and functional annotation from the Araport11 update, the Col-0 genome sequence, RNA-seq and array expression, co-expression, protein interactions, homologs, pathways, publications, alleles, germplasm and phenotypes. The data are collected from a wide variety of public resources. Users can browse gene-specific data through Gene Report pages, identify and create gene lists based on experiments or indexed keywords, and run GO enrichment analysis to investigate the biological significance of selected gene sets. Developed by the Arabidopsis Information Portal project (Araport, https://www.araport.org/), ThaleMine uses the InterMine software framework, which builds well-structured data, and provides powerful data query and analysis functionality. The warehoused data can be accessed by users via graphical interfaces, as well as programmatically via web-services. Here we describe recent developments in ThaleMine including new features and extensions, and discuss future improvements. InterMine has been broadly adopted by the model organism research community including nematode, rat, mouse, zebrafish, budding yeast, the modENCODE project, as well as being used for human data. ThaleMine is the first InterMine developed for a plant model. As additional new plant InterMines are developed by the legume and other plant research communities, the potential of cross-organism integrative data analysis will be further enabled.

Jetstream is a first-of-a-kind system for the NSF - a distributed production cloud resource. The NSF awarded funds to create Jetstream in November 2014. Here we review the purpose for creating Jetstream, present the acceptance test results that define Jetstream's key characteristics, describe our experiences in standing up an OpenStack-based cloud environment, and share some of the early scientific results that have been obtained by researchers and students using this system. Jetstream offers unique capability within the XSEDE-supported US national cyberinfrastructure, delivering interactive virtual machines (VMs) via the Atmosphere interface developed by the University of Arizona. As a multi-region deployment that operates as a single integrated system, Jetstream is proving effective in supporting modes and disciplines of research traditionally underrepresented on larger XSEDE-supported clusters and supercomputers. Already, researchers in biology, network science, economics, earth science, and computer science have used Jetstream to perform research -- much of it research in the "long tail of science."

We have entered a new era in agricultural and biomedical science made possible by remarkable advances in DNA sequencing technologies. The complete sequence of an individual's set of chromosomes (collectively, its genome) provides a primary genetic code for what makes that individual unique, just as the contents of every personal computer reflect the unique attributes of its owner. But a second code, composed of "epigenetic" layers of information, affects the accessibility of the stored information and the execution of specific tasks. Nature's second code is enigmatic and must be deciphered if we are to fully understand and optimize the genetic potential of crop plants. The goal of the Epigenomics of Plants International Consortium is to crack this second code, and ultimately master its control, to help catalyze a new green revolution.

Fission yeast is an important model for epigenetic studies due to the ease with which genetic mutants can be isolated. However, it can be difficult to complement epigenetic phenotypes with genomic libraries in order to identify the genes responsible. This is because epigenetic phenotypes are typically unstable, and can prohibit complementation if silencing cannot be reestablished. Here we have resequenced the fission yeast genome following mutagenesis to readily identify a novel mutant involved in heterochromatic silencing. Candidate genes were identified as functional single base changes linked to the mutation, which were then reconstituted in a wild-type strain to recapitulate the mutant phenotype. By this procedure we identified a weak allele of ubc4, which encodes an essential E2 ubiquitin ligase, as responsible for the swi*603 mutant phenotype. In combination with a large collection of mutants and suppressor plasmids, next-generation genomic resequencing promises to dramatically enhance the power of yeast genetics, permitting the isolation of subtle alleles of essential genes, alleles with quantitative effects, and enhancers and suppressors of heterochromatic silencing.

Plant cells grown in culture exhibit genetic and epigenetic instability. Using a combination of chromatin immunoprecipitation and DNA methylation profiling on tiling microarrays, we have mapped the location and abundance of histone and DNA modifications in a continuously proliferating, dedifferentiated cell suspension culture of Arabidopsis. We have found that euchromatin becomes hypermethylated in culture and that a small percentage of the hypermethylated genes become associated with heterochromatic marks. In contrast, the heterochromatin undergoes dramatic and very precise DNA hypomethylation with transcriptional activation of specific transposable elements (TEs) in culture. High throughput sequencing of small interfering RNA (siRNA) revealed that TEs activated in culture have increased levels of 21-nucleotide (nt) siRNA, sometimes at the expense of the 24-nt siRNA class. In contrast, TEs that remain silent, which match the predominant 24-nt siRNA class, do not change significantly in their siRNA profiles. These results implicate RNA interference and chromatin modification in epigenetic restructuring of the genome following the activation of TEs in immortalized cell culture., Author Summary Cultured plant cells are genetically and epigenetically unstable. We investigated the epigenomic consequences of long-term plant cell culture and found that some genes show an increase in DNA methylation, reminiscent of immortalized animal cell lines and cancer cells. By contrast, in the heterochromatic portion of the genome, some transposable elements (TE), undergo dramatic and very precise loss of DNA methylation and transcriptional activation. The reactivated TEs in culture are also accompanied by a production of 21-nucleotide (nt) small RNAs. In contrast, TEs that remain methylated and silent retain the predominant class of 24-nt small RNAs, and do not change significantly in their small RNA profiles. Our results implicate RNA interference in epigenetic restructuring of the genome following the activation of TEs in immortalized cell culture., Immortalized plant cells undergo epigenomic changes, reminiscent of immortalized animal cell lines and cancer cells, accompanied by shifts in small RNA classes, implicating a role for the RNAi machinery in regulating the epigenome.

In most eukaryotes, histone methylation patterns regulate chromatin architecture and function: methylation of histone H3 lysine-9 (H3K9) demarcates heterochromatin while H3K4 methylation demarcates euchromatin. We show here that the S. pombe JmjC-domain protein Lid2 is a trimethyl H3K4 demethylase responsible for H3K4 hypomethylation in heterochromatin. Lid2 interacts with the histone lysine-9 methyltransferase, Clr4, through the Dos1/Clr8-Rik1 complex, which also functions in the RNA interference pathway. Disruption of the JmjC-domain alone results in severe heterochromatin defects and depletion of siRNA, while overexpressing Lid2 enhances heterochromatin silencing. The physical and functional link between H3K4 demethylation and H3K9 methylation suggests that the two reactions act in a coordinated manner. Surprisingly, cross-regulation of H3K4 and H3K9 methylation in euchromatin also requires Lid2. We provide evidence suggesting that Lid2 enzymatic activity in euchromatin is regulated through dynamic interplay with other histone modification enzymes. Our findings provide a mechanistic insight into the coordination of H3K4 and H3K9 methylation.

Histone methylation regulates transcription, chromatin structure, and the epigenetic state of the cell. Recent studies identified the JmjC domain as a catalytic module for histone demethylation. Schizosaccharomyces pombe contains seven JmjC proteins, but it was unclear whether any of them functioned as histone demethylases. In this report, we show that the JmjC protein Jmj2, which is evolutionarily conserved from yeast to human, reversed trimethylated H3-Lys-4 to di- and mono-but not unmethylated products. Overexpression of Jmj2 but not a catalytically inactive mutant reduced H3-Lys-4 trimethylation levels in vivo and suppressed the toxicity caused by overexpression of the H3-Lys-4-me3-binding protein Yng1 in budding yeast. Genome-wide analysis showed that the loss of jmj2 was associated with an increase in the H3-Lys-4-me3 signal, which was enriched near the transcriptional start sites and the coding regions. At the mating-type locus, the loss of jmj2 or substitution of jmj2 with a catalytically inactive form is correlated with increased reporter gene transcription and H3-Lys-4-me3/2 levels, suggesting that Jmj2 and its demethylase activity may play a role in heterochromatin biology. Our findings identified a novel S. pombe histone demethylase with specificity toward di- and trimethylated histone H3-Lys-4 and a possible role in heterochromatin regulation.

International audience; LSD1 represses and activates transcription by demethylating histone H3K4me and H3K9me, respectively. Genetic ablation of the S. pombe homologs, splsd1 and splsd2, resulted in slow growth and lethality, respectively, underscoring their physiological importance. spLsd1 and spLsd2 form a stable protein complex, which exhibits demethylase activity toward methylated H3K9 in vitro. Both proteins were associated with the heterochromatin boundary regions and euchromatic gene promoters. Loss of spLsd1 resulted in increased H3K9 methylation accompanied by reduced euchromatic gene transcription and heterochromatin propagation. Removal of the H3K9 methylase Clr4 partially suppressed the slow growth phenotype of splsd1Delta. Conversely, catalytically inactivating point mutations in the splsd1 and splsd2 genes partially mimicked the growth and heterochromatin propagation phenotypes. Taken together, these findings suggest the importance of both enzymatic and nonenzymatic roles of spLsd1 in regulating heterochromatin propagation and euchromatic transcription and also suggest that misregulation of spLsd1/2 is likely to impact the epigenetic state of the cell.

The maize genome is relatively large (∼ 2.3 Gb) and has a complex organization of interspersed genes and transposable elements, which necessitates frequent boundaries between different types of chromatin. The examination of maize genes and conserved noncoding sequences revealed that many of these are flanked by regions of elevated asymmetric CHH (where H is A, C, or T) methylation (termed mCHH islands). These mCHH islands are quite short (∼ 100 bp), are enriched near active genes, and often occur at the edge of the transposon that is located nearest to genes. The analysis of DNA methylation in other sequence contexts and several chromatin modifications revealed that mCHH islands mark the transition from heterochromatin-associated modifications to euchromatin-associated modifications. The presence of an mCHH island is fairly consistent in several distinct tissues that were surveyed but shows some variation among different haplotypes. The presence of insertion/deletions in promoters often influences the presence and position of an mCHH island. The mCHH islands are dependent upon RNA-directed DNA methylation activities and are lost in mop1 and mop3 mutants, but the nearby genes rarely exhibit altered expression levels. Instead, loss of an mCHH island is often accompanied by additional loss of DNA methylation in CG and CHG contexts associated with heterochromatin in nearby transposons. This suggests that mCHH islands and RNA-directed DNA methylation near maize genes may act to preserve the silencing of transposons from activity of nearby genes.

ARGONAUTE1 (AGO1) of Arabidopsis thaliana mediates the cleavage of microRNA (miRNA)-targeted mRNAs, and it has also been implicated in the posttranscriptional silencing of transgenes and the maintenance of chromatin structure. Mutations in AGO1 severely disrupt plant development, indicating that miRNA function and possibly other aspects of RNA interference are essential for maintaining normal patterns of gene expression. Using microarrays, we found that 1 to 6% of genes display significant expression changes in several alleles of ago1 at multiple developmental stages, with the majority showing higher levels. Several classes of known miRNA targets increased markedly in ago1, whereas others showed little or no change. Cleavage of mRNAs within miRNA-homologous sites was reduced but not abolished in an ago1 -null background, indicating that redundant slicer activity exists in Arabidopsis. Small interfering RNAs and larger 30- to 60-nucleotide RNA fragments corresponding to highly upregulated miRNA target genes accumulated in wild-type plants but not in ago1, the RNA-dependent RNA polymerase mutants rdr2 and rdr6, or the Dicer-like mutants dcl1 and dcl3. Both sense and antisense RNAs corresponding to these miRNA targets accumulated in the ago1 and dcl1 backgrounds. These results indicate that a subset of endogenous mRNA targets of RNA interference may be regulated through a mechanism of second-strand RNA synthesis and degradation initiated by or in addition to miRNA-mediated cleavage.

Arabidopsis is a model system not only for studying numerous aspects of plant biology, but also for understanding mechanisms of the rapid evolutionary process associated with genome duplication and polyploidization. Although in animals interspecific hybrids are often sterile and aneuploids are related to disease syndromes, both Arabidopsis autopolyploids and allopolyploids occur in nature and can be readily formed in the laboratory, providing an attractive system for comparing changes in gene expression and genome structure among relatively 'young' and 'established' or 'ancient' polyploids. Powerful reverse and forward genetics in Arabidopsis offer an exceptional means by which regulatory mechanisms of gene and genome duplication may be revealed. Moreover, the Arabidopsis genome is completely sequenced; both coding and non-coding sequences are available. We have developed spotted oligo-gene and chromosome microarrays using the complete Arabidopsis genome sequence. The oligo-gene microarray consists of ~26 000 70-mer oligonucleotides that are designed from all annotated genes in Arabidopsis, and the chromosome microarray contains 1 kb genomic tiling fragments amplified from a chromosomal region or the complete sequence of chromosome 4. We have demonstrated the utility of microarrays for genome-wide analysis of changes in gene expression, genome organization and chromatin structure in Arabidopsis polyploids and related species.

Assimilate partitioning refers to the systemic distribution of sugars and amino acids from sites of primary assimilation (source tissue) to import-dependent tissues and organs (sinks). One of the defining questions in this area is how plants balance source productivity with sink demand. Recent results from our laboratory showed that sucrose transport activity is directly proportional to the transcription rate of the phloem-specific proton-sucrose symporter BvSUT1 in Beta vulgaris L. Moreover, symporter gene transcription is regulated by sucrose levels in the leaf. Here we show that sucrose-dependent regulation of BvSUT1 transcription is mediated, at least in part, by a protein phosphorylation relay pathway. Protein phosphatase inhibitors decreased sucrose transport activity, symporter protein and mRNA abundance, and the relative transcription rate of the symporter gene. In contrast, protein kinase inhibitors had no effect or increased sucrose transport, protein and mRNA abundance, and transcription. Furthermore, pre-treating leaves with kinase inhibitors before feeding with sucrose blocked the sucrose-dependent decrease in symporter transcription and transport activity. The latter observation provides direct evidence for a protein phosphorylation cascade operating between the sucrose-sensor and the transcriptional regulator that controls BvSUT1 expression and, ultimately, phloem loading.