External Periodic Driving of Large Systems of Globally Coupled Phase Oscillators

Large systems of coupled oscillators subjected to a periodic external drive occur in many situations in physics and biology. Here the simple, paradigmatic case of equal-strength, all-to-all sine-coupling of phase oscillators subject to a sinusoidal external drive is considered. The stationary states and their stability are determined. Using the stability information and numerical experiments, parameter space phase diagrams showing when different types of system behavior apply are constructed, and the bifurcations marking transitions between different types of behavior are delineated. The analysis is supported by results of direct numerical simulation of an ensemble of oscillators

Identification of an emphysema- associated genetic variant near TGFB2 with regulatory effects in lung fibroblasts

Murine studies have linked TGF-β signaling to emphysema, and human genome-wide association studies (GWAS) studies of lung function and COPD have identified associated regions near genes in the TGF-β superfamily. However, the functional regulatory mechanisms at these loci have not been identified. We performed the largest GWAS of emphysema patterns to date, identifying 10 GWAS loci including an association peak spanning a 200 kb region downstream from TGFB2. Integrative analysis of publicly available eQTL, DNaseI, and chromatin conformation data identified a putative functional variant, rs1690789, that may regulate TGFB2 expression in human fibroblasts. Using chromatin conformation capture, we confirmed that the region containing rs1690789 contacts the TGFB2 promoter in fibroblasts, and CRISPR/Cas-9 targeted deletion of a ~ 100 bp region containing rs1690789 resulted in decreased TGFB2 expression in primary human lung fibroblasts. These data provide novel mechanistic evidence linking genetic variation affecting the TGF-β pathway to emphysema in humans

Transcription-driven genome organization::A model for chromosome structure and the regulation of gene expression tested through simulations

Current models for the folding of the human genome see a hierarchy stretching down from chromosome territories, through A/B compartments and topologically-associating domains (TADs), to contact domains stabilized by cohesin and CTCF. However, molecular mechanisms underlying this folding, and the way folding affects transcriptional activity, remain obscure. Here we review physical principles driving proteins bound to long polymers into clusters surrounded by loops, and present a parsimonious yet comprehensive model for the way the organization determines function. We argue that clusters of active RNA polymerases and their transcription factors are major architectural features; then, contact domains, TADs and compartments just reflect one or more loops and clusters. We suggest tethering a gene close to a cluster containing appropriate factors—a transcription factory—increases the firing frequency, and offer solutions to many current puzzles concerning the actions of enhancers, super-enhancers, boundaries and eQTLs (expression quantitative trait loci). As a result, the activity of any gene is directly influenced by the activity of other transcription units around it in 3D space, and this is supported by Brownian-dynamics simulations of transcription factors binding to cognate sites on long polymers

The Network Zoo: a multilingual package for the inference and analysis of gene regulatory networks

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods

Regulatory network changes between cell lines and their tissues of origin

Background: Cell lines are an indispensable tool in biomedical research and often used as surrogates for tissues. Although there are recognized important cellular and transcriptomic differences between cell lines and tissues, a systematic overview of the differences between the regulatory processes of a cell line and those of its tissue of origin has not been conducted. The RNA-Seq data generated by the GTEx project is the first available data resource in which it is possible to perform a large-scale transcriptional and regulatory network analysis comparing cell lines with their tissues of origin. Results: We compared 127 paired Epstein-Barr virus transformed lymphoblastoid cell lines (LCLs) and whole blood samples, and 244 paired primary fibroblast cell lines and skin samples. While gene expression analysis confirms that these cell lines carry the expression signatures of their primary tissues, albeit at reduced levels, network analysis indicates that expression changes are the cumulative result of many previously unreported alterations in transcription factor (TF) regulation. More specifically, cell cycle genes are over-expressed in cell lines compared to primary tissues, and this alteration in expression is a result of less repressive TF targeting. We confirmed these regulatory changes for four TFs, including SMAD5, using independent ChIP-seq data from ENCODE. Conclusions: Our results provide novel insights into the regulatory mechanisms controlling the expression differences between cell lines and tissues. The strong changes in TF regulation that we observe suggest that network changes, in addition to transcriptional levels, should be considered when using cell lines as models for tissues. Electronic supplementary material The online version of this article (10.1186/s12864-017-4111-x) contains supplementary material, which is available to authorized users

Ensemble genomic analysis in human lung tissue identifies novel genes for chronic obstructive pulmonary disease

Abstract Background Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with chronic obstructive pulmonary disease (COPD). However, many genetic variants show suggestive evidence for association but do not meet the strict threshold for genome-wide significance. Integrative analysis of multiple omics datasets has the potential to identify novel genes involved in disease pathogenesis by leveraging these variants in a functional, regulatory context. Results We performed expression quantitative trait locus (eQTL) analysis using genome-wide SNP genotyping and gene expression profiling of lung tissue samples from 86 COPD cases and 31 controls, testing for SNPs associated with gene expression levels. These results were integrated with a prior COPD GWAS using an ensemble statistical and network methods approach to identify relevant genes and observe them in the context of overall genetic control of gene expression to highlight co-regulated genes and disease pathways. We identified 250,312 unique SNPs and 4997 genes in the cis(local)-eQTL analysis (5% false discovery rate). The top gene from the integrative analysis was MAPT, a gene recently identified in an independent GWAS of lung function. The genes HNRNPAB and PCBP2 with RNA binding activity and the gene ACVR1B were identified in network communities with validated disease relevance. Conclusions The integration of lung tissue gene expression with genome-wide SNP genotyping and subsequent intersection with prior GWAS and omics studies highlighted candidate genes within COPD loci and in communities harboring known COPD genes. This integration also identified novel disease genes in sub-threshold regions that would otherwise have been missed through GWAS

Multisystem Analysis of Mycobacterium tuberculosis Reveals Kinase-Dependent Remodeling of the Pathogen-Environment Interface

ABSTRACT Tuberculosis is the leading killer among infectious diseases worldwide. Increasing multidrug resistance has prompted new approaches for tuberculosis drug development, including targeted inhibition of virulence determinants and of signaling cascades that control many downstream pathways. We used a multisystem approach to determine the effects of a potent small-molecule inhibitor of the essential Mycobacterium tuberculosis Ser/Thr protein kinases PknA and PknB. We observed differential levels of phosphorylation of many proteins and extensive changes in levels of gene expression, protein abundance, cell wall lipids, and intracellular metabolites. The patterns of these changes indicate regulation by PknA and PknB of several pathways required for cell growth, including ATP synthesis, DNA synthesis, and translation. These data also highlight effects on pathways for remodeling of the mycobacterial cell envelope via control of peptidoglycan turnover, lipid content, a SigE-mediated envelope stress response, transmembrane transport systems, and protein secretion systems. Integrated analysis of phosphoproteins, transcripts, proteins, and lipids identified an unexpected pathway whereby threonine phosphorylation of the essential response regulator MtrA decreases its DNA binding activity. Inhibition of this phosphorylation is linked to decreased expression of genes for peptidoglycan turnover, and of genes for mycolyl transferases, with concomitant changes in mycolates and glycolipids in the cell envelope. These findings reveal novel roles for PknA and PknB in regulating multiple essential cell functions and confirm that these kinases are potentially valuable targets for new antituberculosis drugs. In addition, the data from these linked multisystems provide a valuable resource for future targeted investigations into the pathways regulated by these kinases in the M. tuberculosis cell