Pairs and Pairix: a file format and a tool for efficient storage and retrieval for Hi-C read pairs.

SUMMARY: As the amount of 3D chromosomal interaction data continues to increase, storing and accessing such data efficiently becomes paramount. We introduce Pairs, a block-compressed text file format for storing paired genomic coordinates from Hi-C data, and Pairix, an open-source C application to index and query Pairs files. Pairix (also available in Python and R) extends the functionalities of Tabix to paired coordinates data. We have also developed PairsQC, a collapsible HTML quality control report generator for Pairs files. AVAILABILITY AND IMPLEMENTATION: The format specification and source code are available at https://github.com/4dn-dcic/pairix, https://github.com/4dn-dcic/Rpairix and https://github.com/4dn-dcic/pairsqc.

Tibanna: software for scalable execution of portable pipelines on the cloud.

SUMMARY: We introduce Tibanna, an open-source software tool for automated execution of bioinformatics pipelines on Amazon Web Services (AWS). Tibanna accepts reproducible and portable pipeline standards including Common Workflow Language (CWL), Workflow Description Language (WDL) and Docker. It adopts a strategy of isolation and optimization of individual executions, combined with a serverless scheduling approach. Pipelines are executed and monitored using local commands or the Python Application Programming Interface (API) and cloud configuration is automatically handled. Tibanna is well suited for projects with a range of computational requirements, including those with large and widely fluctuating loads. Notably, it has been used to process terabytes of data for the 4D Nucleome (4DN) Network. AVAILABILITY AND IMPLEMENTATION: Source code is available on GitHub at https://github.com/4dn-dcic/tibanna. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PaSD-qc: quality control for single cell whole-genome sequencing data using power spectral density estimation.

Single cell whole-genome sequencing (scWGS) is providing novel insights into the nature of genetic heterogeneity in normal and diseased cells. However, the whole-genome amplification process required for scWGS introduces biases into the resulting sequencing that can confound downstream analysis. Here, we present a statistical method, with an accompanying package PaSD-qc (Power Spectral Density-qc), that evaluates the properties and quality of single cell libraries. It uses a modified power spectral density to assess amplification uniformity, amplicon size distribution, autocovariance and inter-sample consistency as well as to identify chromosomes with aberrant read-density profiles due either to copy alterations or poor amplification. These metrics provide a standard way to compare the quality of single cell samples as well as yield information necessary to improve variant calling strategies. We demonstrate the usefulness of this tool in comparing the properties of scWGS protocols, identifying potential chromosomal copy number variation, determining chromosomal and subchromosomal regions of poor amplification, and selecting high-quality libraries from low-coverage data for deep sequencing. The software is available free and open-source at https://github.com/parklab/PaSDqc.

The 4D Nucleome Data Portal as a resource for searching and visualizing curated nucleomics data.

The 4D Nucleome (4DN) Network aims to elucidate the complex structure and organization of chromosomes in the nucleus and the impact of their disruption in disease biology. We present the 4DN Data Portal ( https://data.4dnucleome.org/ ), a repository for datasets generated in the 4DN network and relevant external datasets. Datasets were generated with a wide range of experiments, including chromosome conformation capture assays such as Hi-C and other innovative sequencing and microscopy-based assays probing chromosome architecture. All together, the 4DN data portal hosts more than 1800 experiment sets and 36000 files. Results of sequencing-based assays from different laboratories are uniformly processed and quality-controlled. The portal interface allows easy browsing, filtering, and bulk downloads, and the integrated HiGlass genome browser allows interactive visualization and comparison of multiple datasets. The 4DN data portal represents a primary resource for chromosome contact and other nuclear architecture data for the scientific community.

Aging and neurodegeneration are associated with increased mutations in single human neurons.

It has long been hypothesized that aging and neurodegeneration are associated with somatic mutation in neurons; however, methodological hurdles have prevented testing this hypothesis directly. We used single-cell whole-genome sequencing to perform genome-wide somatic single-nucleotide variant (sSNV) identification on DNA from 161 single neurons from the prefrontal cortex and hippocampus of 15 normal individuals (aged 4 months to 82 years), as well as 9 individuals affected by early-onset neurodegeneration due to genetic disorders of DNA repair (Cockayne syndrome and xeroderma pigmentosum). sSNVs increased approximately linearly with age in both areas (with a higher rate in hippocampus) and were more abundant in neurodegenerative disease. The accumulation of somatic mutations with age-which we term genosenium-shows age-related, region-related, and disease-related molecular signatures and may be important in other human age-associated conditions.

Developmental expression of Kv1 voltage-gated potassium channels in the avian hypothalamus.

Specialized hypothalamic neurons integrate the homeostatic balance between food intake and energy expenditure, processes that may become dysregulated during the development of diabetes, obesity, and other metabolic disorders. Shaker family voltage-gated potassium channels (Kv1) contribute to the maintenance of resting membrane potential, action potential characteristics, and neurotransmitter release in many populations of neurons, although hypothalamic Kv1 channel expression has been largely unexplored. Whole-cell patch clamp recordings from avian hypothalamic brain slices demonstrate a developmental shift in the electrophysiological properties of avian arcuate nucleus neurons, identifying an increase in outward ionic current that corresponds with action potential maturation. Additionally, RT-PCR experiments identified the early expression of Kv1.2, Kv1.3, and Kv1.5 mRNA in the embryonic avian hypothalamus, suggesting that these channels may underlie the electrophysiological changes observed in these neurons. Real-time quantitative PCR analysis on intact microdissections of embryonic hypothalamic tissue revealed a concomitant increase in Kv1.2 and Kv1.5 gene expression at key electrophysiological time points during development. This study is the first to demonstrate hypothalamic mRNA expression of Kv1 channels in developing avian embryos and may suggest a role for voltage-gated ion channel regulation in the physiological patterning of embryonic hypothalamic circuits governing energy homeostasis.

Exposure to the lampricide 3-trifluoromethyl-4-nitrophenol results in increased expression of carbohydrate transporters in Saccharomyces cerevisiae.

The lampricide 3-trifluoromethyl-4-nitrophenol (TFM) is used to control sea lamprey (Petromyzon marinus) populations in freshwater lakes. Although TFM can have sublethal and lethal effects, little is known about gene expression changes with TFM exposure. Microarray analysis was used to determine differential gene expression over 4 h of exposure in Saccharomyces cerevisiae. Among the most significantly up-regulated genes were regulators of carbohydrate transport, including HXT1, HXT3, HXT4, IMA5, MIG2, and YKR075C. Environ Toxicol Chem 2016;35:1727-1732. © 2015 SETAC.