RNA polymerases reshape chromatin architecture and couple transcription on individual fibers.

RNA polymerases must initiate and pause within a complex chromatin environment, surrounded by nucleosomes and other transcriptional machinery. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address this, we employed long-read chromatin fiber sequencing (Fiber-seq) in Drosophila to visualize RNA polymerase (Pol) within its native chromatin context with single-molecule precision along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of individual Pol II, nucleosome, and transcription factor footprints, revealing Pol II pausing-driven destabilization of downstream nucleosomes. Furthermore, we demonstrate pervasive direct distance-dependent transcriptional coupling between nearby Pol II genes, Pol III genes, and transcribed enhancers, modulated by local chromatin architecture. Overall, transcription initiation reshapes surrounding nucleosome architecture and couples nearby transcriptional machinery along individual chromatin fibers.

Developmental fate and cellular maturity encoded in human regulatory DNA landscapes.

Cellular-state information between generations of developing cells may be propagated via regulatory regions. We report consistent patterns of gain and loss of DNase I-hypersensitive sites (DHSs) as cells progress from embryonic stem cells (ESCs) to terminal fates. DHS patterns alone convey rich information about cell fate and lineage relationships distinct from information conveyed by gene expression. Developing cells share a proportion of their DHS landscapes with ESCs; that proportion decreases continuously in each cell type as differentiation progresses, providing a quantitative benchmark of developmental maturity. Developmentally stable DHSs densely encode binding sites for transcription factors involved in autoregulatory feedback circuits. In contrast to normal cells, cancer cells extensively reactivate silenced ESC DHSs and those from developmental programs external to the cell lineage from which the malignancy derives. Our results point to changes in regulatory DNA landscapes as quantitative indicators of cell-fate transitions, lineage relationships, and dysfunction.

DNA-m6A calling and integrated long-read epigenetic and genetic analysis with fibertools.

Long-read DNA sequencing has recently emerged as a powerful tool for studying both genetic and epigenetic architectures at single-molecule and single-nucleotide resolution. Long-read epigenetic studies encompass both the direct identification of native cytosine methylation and the identification of exogenously placed DNA N 6 -methyladenine (DNA-m6A). However, detecting DNA-m6A modifications using single-molecule sequencing, as well as coprocessing single-molecule genetic and epigenetic architectures, is limited by computational demands and a lack of supporting tools. Here, we introduce fibertools, a state-of-the-art toolkit that features a semisupervised convolutional neural network for fast and accurate identification of m6A-marked bases using Pacific Biosciences (PacBio) single-molecule long-read sequencing, as well as the coprocessing of long-read genetic and epigenetic data produced using either the PacBio or Oxford Nanopore Technologies (ONT) sequencing platforms. We demonstrate accurate DNA-m6A identification (>90% precision and recall) along >20 kb long DNA molecules with an ∼1000-fold improvement in speed. In addition, we demonstrate that fibertools can readily integrate genetic and epigenetic data at single-molecule resolution, including the seamless conversion between molecular and reference coordinate systems, allowing for accurate genetic and epigenetic analyses of long-read data within structurally and somatically variable genomic regions.

Circuitry and dynamics of human transcription factor regulatory networks.

The combinatorial cross-regulation of hundreds of sequence-specific transcription factors (TFs) defines a regulatory network that underlies cellular identity and function. Here we use genome-wide maps of in vivo DNaseI footprints to assemble an extensive core human regulatory network comprising connections among 475 sequence-specific TFs and to analyze the dynamics of these connections across 41 diverse cell and tissue types. We find that human TF networks are highly cell selective and are driven by cohorts of factors that include regulators with previously unrecognized roles in control of cellular identity. Moreover, we identify many widely expressed factors that impact transcriptional regulatory networks in a cell-selective manner. Strikingly, in spite of their inherent diversity, all cell-type regulatory networks independently converge on a common architecture that closely resembles the topology of living neuronal networks. Together, our results provide an extensive description of the circuitry, dynamics, and organizing principles of the human TF regulatory network.

Evaluation of N 6-methyldeoxyadenosine antibody-based genomic profiling in eukaryotes.

Low-level DNA N 6-methyldeoxyadenosine (DNA-m6A) has recently been reported across various eukaryotes. Although anti-m6A antibody-based approaches are commonly used to measure DNA-m6A levels, this approach is known to be confounded by DNA secondary structures, RNA contamination, and bacterial contamination. To evaluate for these confounding features, we introduce an approach for systematically validating the selectivity of antibody-based DNA-m6A methods and use a highly selective anti-DNA-m6A antibody to reexamine patterns of DNA-m6A in C. reinhardtii, A. thaliana, and D. melanogaster Our findings raise caution about the use of antibody-based methods for endogenous m6A quantification and mapping in eukaryotes.

Functional categorization of gene regulatory variants that cause Mendelian conditions.

Much of our current understanding of rare human diseases is driven by coding genetic variants. However, non-coding genetic variants play a pivotal role in numerous rare human diseases, resulting in diverse functional impacts ranging from altered gene regulation, splicing, and/or transcript stability. With the increasing use of genome sequencing in clinical practice, it is paramount to have a clear framework for understanding how non-coding genetic variants cause disease. To this end, we have synthesized the literature on hundreds of non-coding genetic variants that cause rare Mendelian conditions via the disruption of gene regulatory patterns and propose a functional classification system. Specifically, we have adapted the functional classification framework used for coding variants (i.e., loss-of-function, gain-of-function, and dominant-negative) to account for features unique to non-coding gene regulatory variants. We identify that non-coding gene regulatory variants can be split into three distinct categories by functional impact: (1) non-modular loss-of-expression (LOE) variants; (2) modular loss-of-expression (mLOE) variants; and (3) gain-of-ectopic-expression (GOE) variants. Whereas LOE variants have a direct corollary with coding loss-of-function variants, mLOE and GOE variants represent disease mechanisms that are largely unique to non-coding variants. These functional classifications aim to provide a unified terminology for categorizing the functional impact of non-coding variants that disrupt gene regulatory patterns in Mendelian conditions.

Implementing evidence-based assertions of clinical actionability in the context of secondary findings: Updates from the ClinGen Actionability Working Group.

PURPOSE: The ClinGen Actionability Working Group (AWG) developed an evidence-based framework to generate actionability reports and scores of gene-condition pairs in the context of secondary findings from genome sequencing. Here we describe the expansion of the framework to include actionability assertions. METHODS: Initial development of the actionability rubric was based on previously scored adult gene-condition pairs and individual expert evaluation. Rubric refinement was iterative and based on evaluation, feedback, and discussion. The final rubric was pragmatically evaluated via integration into actionability assessments for 27 gene-condition pairs. RESULTS: The resulting rubric has a 4-point scale (limited, moderate, strong, and definitive) and uses the highest-scoring outcome-intervention pair of each gene-condition pair to generate a preliminary assertion. During AWG discussions, predefined criteria and factors guide discussion to produce a consensus assertion for a gene-condition pair, which may differ from the preliminary assertion. The AWG has retrospectively generated assertions for all previously scored gene-condition pairs and are prospectively asserting on gene-condition pairs under assessment, having completed over 170 adult and 188 pediatric gene-condition pairs. CONCLUSION: The AWG expanded its framework to provide actionability assertions to enhance the clinical value of their resources and increase their utility as decision aids regarding return of secondary findings.

Recurrent SLC1A2 variants cause epilepsy via a dominant negative mechanism.

SLC1A2 is a trimeric transporter essential for clearing glutamate from neuronal synapses. Recurrent de novo SLC1A2 missense variants cause a severe, early onset developmental and epileptic encephalopathy via an unclear mechanism. We demonstrate that all 3 variants implicated in this condition localize to the trimerization domain of SLC1A2, and that the Leu85Pro variant acts via a dominant negative mechanism to reduce, but not eliminate, wild-type SLC1A2 protein localization and function. Finally, we demonstrate that treatment of a 20-month-old SLC1A2-related epilepsy patient with the SLC1A2-modulating agent ceftriaxone did not result in a significant change in daily spasm count. ANN NEUROL 2019;85:921-926.

A retrospective study of adult patients with noncirrhotic hyperammonemia.

Adult-onset noncirrhotic hyperammonemia (NCH) is poorly understood and has a high morbidity and mortality. To elucidate the etiology and management of NCH, we performed a retrospective analysis of 23 adults (median age 51) with NCH treated between 2014 and 2020 at two academic medical centers. Hyperammonemia was diagnosed in all cases during the evaluation of altered mental status, with 22% presenting with seizures. Peak ammonia levels were >200 µmol/L in 70% of cases. Defects in ammonia metabolism were assessed using urea cycle biochemical testing, germline genetic testing, and testing for urease-producing infectious agents. Ammonia metabolism defects in these cases appear attributable to four major sources: (a) infection with urease-producing organism (n = 5); (b) previously undiagnosed inborn errors of metabolism (IEMs) (n = 4); (c) clinical exposures causing acquired urea cycle dysfunction (n = 6); and (d) unexplained acquired urea cycle dysfunction (uaUCD) (n = 8), as evidenced by biochemical signatures of urea cycle dysfunction without a genetic or clinical exposure. Severe protein malnutrition appeared to be a reversible risk factor for uaUCD. Overall, 13% of our cohort died prior to resolution of hyperammonemia, 26% died after hyperammonemia resolution, 57% survived after having reversible neurological changes, and 4% survived with irreversible neurological changes. Renal replacement therapy for ammonia clearance was often utilized for patients with an ammonia level above 250 µmol/L and patients were frequently empirically treated with antibiotics targeting urea-splitting organisms. Our study demonstrates that acquired urea cycle dysfunction, IEMs and urease-producing infections are major sources of adult-onset NCH and highlights successful management strategies for adult-onset NCH.

Conservation of trans-acting circuitry during mammalian regulatory evolution.

The basic body plan and major physiological axes have been highly conserved during mammalian evolution, yet only a small fraction of the human genome sequence appears to be subject to evolutionary constraint. To quantify cis- versus trans-acting contributions to mammalian regulatory evolution, we performed genomic DNase I footprinting of the mouse genome across 25 cell and tissue types, collectively defining ∼8.6 million transcription factor (TF) occupancy sites at nucleotide resolution. Here we show that mouse TF footprints conjointly encode a regulatory lexicon that is ∼95% similar with that derived from human TF footprints. However, only ∼20% of mouse TF footprints have human orthologues. Despite substantial turnover of the cis-regulatory landscape, nearly half of all pairwise regulatory interactions connecting mouse TF genes have been maintained in orthologous human cell types through evolutionary innovation of TF recognition sequences. Furthermore, the higher-level organization of mouse TF-to-TF connections into cellular network architectures is nearly identical with human. Our results indicate that evolutionary selection on mammalian gene regulation is targeted chiefly at the level of trans-regulatory circuitry, enabling and potentiating cis-regulatory plasticity.

Integrative Genomics Analysis Identifies ACVR1B as a Candidate Causal Gene of Emphysema Distribution.

Genome-wide association studies (GWAS) have identified multiple associations with emphysema apicobasal distribution (EABD), but the biological functions of these variants are unknown. To characterize the functions of EABD-associated variants, we integrated GWAS results with 1) expression quantitative trait loci (eQTL) from the Genotype Tissue Expression (GTEx) project and subjects in the COPDGene (Genetic Epidemiology of COPD) study and 2) cell type epigenomic marks from the Roadmap Epigenomics project. On the basis of these analyses, we selected a variant near ACVR1B (activin A receptor type 1B) for functional validation. SNPs from 168 loci with P values less than 5 × 10-5 in the largest GWAS meta-analysis of EABD were analyzed. Eighty-four loci overlapped eQTL, with 12 of these loci showing greater than 80% likelihood of harboring a single, shared GWAS and eQTL causal variant. Seventeen cell types were enriched for overlap between EABD loci and Roadmap Epigenomics marks (permutation P < 0.05), with the strongest enrichment observed in CD4+, CD8+, and regulatory T cells. We selected a putative causal variant, rs7962469, associated with ACVR1B expression in lung tissue for additional functional investigation, and reporter assays confirmed allele-specific regulatory activity for this variant in human bronchial epithelial and Jurkat immune cell lines. ACVR1B expression levels exhibit a nominally significant association with emphysema distribution. EABD-associated loci are preferentially enriched in regulatory elements of multiple cell types, most notably T-cell subsets. Multiple EABD loci colocalize to regulatory elements that are active across multiple tissues and cell types, and functional analyses confirm the presence of an EABD-associated functional variant that regulates ACVR1B expression, indicating that transforming growth factor-ß signaling plays a role in the EABD phenotype. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Noncirrhotic hyperammonemia after deceased donor kidney transplantation: A case report.

A 72-year-old woman with end-stage kidney disease due to recurrent urinary tract infections and obstructive uropathy of a solitary kidney presented to our hospital for renal transplantation. She underwent successful transplantation of a deceased donor allograft, but developed acute mental status deterioration on the fifth postoperative day. Her serum ammonia was found to be markedly elevated to 447 µmol/L in the setting of normal hepatic function. She was treated with emergent dialysis and empiric antibiotics targeting urea-splitting organisms, and ultimately made a full neurologic recovery with stable renal allograft function. Noncirrhotic hyperammonemia (NCH) is an exceedingly rare clinical entity but seems to have a predilection for patients who have undergone solid organ transplantation. This report emphasizes the importance of rapid diagnosis and initiation of treatment for NCH, which is associated with a high rate of mortality and irreversible neurological morbidity. We outline the successful workup and management approach for this patient.

A single nucleotide polymorphism of cyclin-dependent kinase inhibitor 1B (p27Kip1) associated with human vein graft failure affects growth of human venous adventitial cells but not smooth muscle cells.

BACKGROUND: Cyclin-dependent kinase inhibitor 1B (p27Kip1) is a cell-cycle inhibitor whose -838C>A single nucleotide polymorphism (rs36228499; hereafter called p27 SNP) has been associated with the clinical failure of peripheral vein grafts, but the functional effects of this SNP have not been demonstrated. METHODS: Human saphenous vein adventitial cells and intimal/medial smooth muscle cells (SMCs) were derived from explants obtained at the time of lower extremity bypass operations. We determined the following in adventitial cells and SMCs as a function of the p27 SNP genotype: (1) p27 promoter activity, (2) p27 messenger (m)RNA and protein levels, and (3) growth and collagen gel contraction. Deoxyribonuclease I footprinting was also performed in adventitial cells and SMCs. RESULTS: p27 promoter activity, deoxyribonuclease I footprinting, p27 mRNA levels, and p27 protein levels demonstrated that the p27 SNP is functional in adventitial cells and SMCs. Both cell types with the graft failure protective AA genotype had more p27 mRNA and protein. As predicted because of higher levels of p27 protein, adventitial cells with the AA genotype grew slower than those of the CC genotype. Unexpectedly, SMCs did not show this genotype-dependent growth response. CONCLUSIONS: These results support the functionality of the p27 SNP in venous SMCs and adventitial cells, but an effect of the SNP on cell proliferation is limited to only adventitial cells. These data point to a potential role for adventitial cells in human vein graft failure and also suggest that SMCs express factors that interfere with the activity of p27.

The accessible chromatin landscape of the human genome.

DNase I hypersensitive sites (DHSs) are markers of regulatory DNA and have underpinned the discovery of all classes of cis-regulatory elements including enhancers, promoters, insulators, silencers and locus control regions. Here we present the first extensive map of human DHSs identified through genome-wide profiling in 125 diverse cell and tissue types. We identify ∼2.9 million DHSs that encompass virtually all known experimentally validated cis-regulatory sequences and expose a vast trove of novel elements, most with highly cell-selective regulation. Annotating these elements using ENCODE data reveals novel relationships between chromatin accessibility, transcription, DNA methylation and regulatory factor occupancy patterns. We connect ∼580,000 distal DHSs with their target promoters, revealing systematic pairing of different classes of distal DHSs and specific promoter types. Patterning of chromatin accessibility at many regulatory regions is organized with dozens to hundreds of co-activated elements, and the transcellular DNase I sensitivity pattern at a given region can predict cell-type-specific functional behaviours. The DHS landscape shows signatures of recent functional evolutionary constraint. However, the DHS compartment in pluripotent and immortalized cells exhibits higher mutation rates than that in highly differentiated cells, exposing an unexpected link between chromatin accessibility, proliferative potential and patterns of human variation.

An expansive human regulatory lexicon encoded in transcription factor footprints.

Regulatory factor binding to genomic DNA protects the underlying sequence from cleavage by DNase I, leaving nucleotide-resolution footprints. Using genomic DNase I footprinting across 41 diverse cell and tissue types, we detected 45 million transcription factor occupancy events within regulatory regions, representing differential binding to 8.4 million distinct short sequence elements. Here we show that this small genomic sequence compartment, roughly twice the size of the exome, encodes an expansive repertoire of conserved recognition sequences for DNA-binding proteins that nearly doubles the size of the human cis-regulatory lexicon. We find that genetic variants affecting allelic chromatin states are concentrated in footprints, and that these elements are preferentially sheltered from DNA methylation. High-resolution DNase I cleavage patterns mirror nucleotide-level evolutionary conservation and track the crystallographic topography of protein-DNA interfaces, indicating that transcription factor structure has been evolutionarily imprinted on the human genome sequence. We identify a stereotyped 50-base-pair footprint that precisely defines the site of transcript origination within thousands of human promoters. Finally, we describe a large collection of novel regulatory factor recognition motifs that are highly conserved in both sequence and function, and exhibit cell-selective occupancy patterns that closely parallel major regulators of development, differentiation and pluripotency.

Using Data Independent Acquisition (DIA) to Model High-responding Peptides for Targeted Proteomics Experiments.

Targeted mass spectrometry is an essential tool for detecting quantitative changes in low abundant proteins throughout the proteome. Although selected reaction monitoring (SRM) is the preferred method for quantifying peptides in complex samples, the process of designing SRM assays is laborious. Peptides have widely varying signal responses dictated by sequence-specific physiochemical properties; one major challenge is in selecting representative peptides to target as a proxy for protein abundance. Here we present PREGO, a software tool that predicts high-responding peptides for SRM experiments. PREGO predicts peptide responses with an artificial neural network trained using 11 minimally redundant, maximally relevant properties. Crucial to its success, PREGO is trained using fragment ion intensities of equimolar synthetic peptides extracted from data independent acquisition experiments. Because of similarities in instrumentation and the nature of data collection, relative peptide responses from data independent acquisition experiments are a suitable substitute for SRM experiments because they both make quantitative measurements from integrated fragment ion chromatograms. Using an SRM experiment containing 12,973 peptides from 724 synthetic proteins, PREGO exhibits a 40-85% improvement over previously published approaches at selecting high-responding peptides. These results also represent a dramatic improvement over the rules-based peptide selection approaches commonly used in the literature.

Personal and population genomics of human regulatory variation.

The characteristics and evolutionary forces acting on regulatory variation in humans remains elusive because of the difficulty in defining functionally important noncoding DNA. Here, we combine genome-scale maps of regulatory DNA marked by DNase I hypersensitive sites (DHSs) from 138 cell and tissue types with whole-genome sequences of 53 geographically diverse individuals in order to better delimit the patterns of regulatory variation in humans. We estimate that individuals likely harbor many more functionally important variants in regulatory DNA compared with protein-coding regions, although they are likely to have, on average, smaller effect sizes. Moreover, we demonstrate that there is significant heterogeneity in the level of functional constraint in regulatory DNA among different cell types. We also find marked variability in functional constraint among transcription factor motifs in regulatory DNA, with sequence motifs for major developmental regulators, such as HOX proteins, exhibiting levels of constraint comparable to protein-coding regions. Finally, we perform a genome-wide scan of recent positive selection and identify hundreds of novel substrates of adaptive regulatory evolution that are enriched for biologically interesting pathways such as melanogenesis and adipocytokine signaling. These data and results provide new insights into patterns of regulatory variation in individuals and populations and demonstrate that a large proportion of functionally important variation lies beyond the exome.

Panorama: a targeted proteomics knowledge base.

Panorama is a web application for storing, sharing, analyzing, and reusing targeted assays created and refined with Skyline,1 an increasingly popular Windows client software tool for targeted proteomics experiments. Panorama allows laboratories to store and organize curated results contained in Skyline documents with fine-grained permissions, which facilitates distributed collaboration and secure sharing of published and unpublished data via a web-browser interface. It is fully integrated with the Skyline workflow and supports publishing a document directly to a Panorama server from the Skyline user interface. Panorama captures the complete Skyline document information content in a relational database schema. Curated results published to Panorama can be aggregated and exported as chromatogram libraries. These libraries can be used in Skyline to pick optimal targets in new experiments and to validate peak identification of target peptides. Panorama is open-source and freely available. It is distributed as part of LabKey Server,2 an open source biomedical research data management system. Laboratories and organizations can set up Panorama locally by downloading and installing the software on their own servers. They can also request freely hosted projects on https://panoramaweb.org , a Panorama server maintained by the Department of Genome Sciences at the University of Washington.

Rapid empirical discovery of optimal peptides for targeted proteomics.

We report a method for high-throughput, cost-efficient empirical discovery of optimal proteotypic peptides and fragment ions for targeted proteomics applications using in vitro-synthesized proteins. We demonstrate the approach using human transcription factors, which are typically difficult, low-abundance targets and empirically derived proteotypic peptides for 98% of the target proteins. We show that targeted proteomic assays developed using our approach facilitate robust in vivo quantification of human transcription factors.

Single-nucleoid architecture reveals heterogeneous packaging of mitochondrial DNA.

Cellular metabolism relies on the regulation and maintenance of mitochondrial DNA (mtDNA). Hundreds to thousands of copies of mtDNA exist in each cell, yet because mitochondria lack histones or other machinery important for nuclear genome compaction, it remains unresolved how mtDNA is packaged into individual nucleoids. In this study, we used long-read single-molecule accessibility mapping to measure the compaction of individual full-length mtDNA molecules at near single-nucleotide resolution. We found that, unlike the nuclear genome, human mtDNA largely undergoes all-or-none global compaction, with most nucleoids existing in an inaccessible, inactive state. Highly accessible mitochondrial nucleoids are co-occupied by transcription and replication components and selectively form a triple-stranded displacement loop structure. In addition, we showed that the primary nucleoid-associated protein TFAM directly modulates the fraction of inaccessible nucleoids both in vivo and in vitro, acting consistently with a nucleation-and-spreading mechanism to coat and compact mitochondrial nucleoids. Together, these findings reveal the primary architecture of mtDNA packaging and regulation in human cells.