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 four-point scale (limited, moderate, strong, definitive) and uses the highest-scoring outcome-intervention pair of each gene-condition pair to generate a preliminary assertion. During AWG discussions, pre-defined 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.

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.

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.

RNA polymerases reshape chromatin and coordinate transcription on individual fibers

During eukaryotic transcription, RNA polymerases must initiate and pause within a crowded, complex environment, surrounded by nucleosomes and other transcriptional activity. 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 these limitations, we employed long-read chromatin fiber sequencing (Fiber-seq) to visualize RNA polymerases within their native chromatin context at single-molecule and near single-nucleotide resolution along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of single-molecule RNA Polymerase (Pol) II and III transcription associated footprints, which, in aggregate, mirror bulk short-read sequencing-based measurements of transcription. We show that Pol II pausing destabilizes downstream nucleosomes, with frequently paused genes maintaining a short-term memory of these destabilized nucleosomes. Furthermore, we demonstrate pervasive direct coordination and anti-coordination between nearby Pol II genes, Pol III genes, transcribed enhancers, and insulator elements. This coordination is largely limited to spatially organized elements within 5 kb of each other, implicating short-range chromatin environments as a predominant determinant of coordinated polymerase initiation. Overall, transcription initiation reshapes surrounding nucleosome architecture and coordinates nearby transcriptional machinery along individual chromatin fibers.

Multi-site validation of a functional assay to adjudicate SCN5A Brugada Syndrome-associated variants

Brugada Syndrome (BrS) is an inheritable arrhythmia condition that is associated with rare, loss-of-function variants in the cardiac sodium channel gene, SCN5A. Interpreting the pathogenicity of SCN5A missense variants is challenging and ~79% of SCN5A missense variants in ClinVar are currently classified as Variants of Uncertain Significance (VUS). An in vitro SCN5A-BrS automated patch clamp assay was generated for high-throughput functional studies of NaV1.5. The assay was independently studied at two separate research sites - Vanderbilt University Medical Center and Victor Chang Cardiac Research Institute - revealing strong correlations, including peak INa density (R2=0.86). The assay was calibrated according to ClinGen Sequence Variant Interpretation recommendations using high-confidence variant controls (n=49). Normal and abnormal ranges of function were established based on the distribution of benign variant assay results. The assay accurately distinguished benign controls (24/25) from pathogenic controls (23/24). Odds of Pathogenicity values derived from the experimental results yielded 0.042 for normal function (BS3 criterion) and 24.0 for abnormal function (PS3 criterion), resulting in up to strong evidence for both ACMG criteria. The calibrated assay was then used to study SCN5A VUS observed in four families with BrS and other arrhythmia phenotypes associated with SCN5A loss-of-function. The assay revealed loss-of-function for three of four variants, enabling reclassification to likely pathogenic. This validated APC assay provides clinical-grade functional evidence for the reclassification of current VUS and will aid future SCN5A-BrS variant classification.

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 as well as the identification of exogenously placed DNA N6-methyladenine (DNA-m6A). However, detecting DNA-m6A modifications using single-molecule sequencing, as well as co-processing 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 semi-supervised convolutional neural network for fast and accurate identification of m6A-marked bases using PacBio single-molecule long-read sequencing, as well as the co-processing of long-read genetic and epigenetic data produced using either PacBio or Oxford Nanopore sequencing platforms. We demonstrate accurate DNA-m6A identification (>90% precision and recall) along >20 kilobase long DNA molecules with a ~1,000-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.

Full-length Isoform Sequencing for Resolving the Molecular Basis of Charcot-Marie-Tooth 2A

Objectives

Transcript sequencing of patient-derived samples has been shown to improve the diagnostic yield for solving cases of suspected Mendelian conditions, yet the added benefit of full-length long-read transcript sequencing is largely unexplored.

Methods

We applied short-read and full-length transcript sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis.

Results

We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts the branch point critical for intron 6 splicing. Full-length long-read isoform complementary DNA (cDNA) sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates 5 distinct altered splicing transcripts. All 5 altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A).

Discussion

This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Dominant-negative variant in SLC1A4 causes an autosomal dominant epilepsy syndrome

SLC1A4 is a trimeric neutral amino acid transporter essential for shuttling L-serine from astrocytes into neurons. Individuals with biallelic variants in SLC1A4 are known to have spastic tetraplegia, thin corpus callosum, and progressive microcephaly (SPATCCM) syndrome, but individuals with heterozygous variants are not thought to have disease. We identify an 8-year-old patient with global developmental delay, spasticity, epilepsy, and microcephaly who has a de novo heterozygous three amino acid duplication in SLC1A4 (L86_M88dup). We demonstrate that L86_M88dup causes a dominant-negative N-glycosylation defect of SLC1A4, which in turn reduces the plasma membrane localization of SLC1A4 and the transport rate of SLC1A4 for L-serine.

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.

Full-length isoform sequencing for resolving the molecular basis of Charcot-Marie-Tooth 2A

Objectives

Transcript sequencing of patient derived samples has been shown to improve the diagnostic yield for solving cases of likely Mendelian disorders, yet the added benefit of full-length long-read transcript sequencing is largely unexplored.

Methods

We applied short-read and full-length isoform cDNA sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis.

Results

We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts a branch point critical for intron 6 spicing. Full-length long-read isoform cDNA sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates five distinct altered splicing transcripts. All five altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 protein levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A).

Discussion

This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Rewiring of human neurodevelopmental gene regulatory programs by human accelerated regions

Human accelerated regions (HARs) are the fastest-evolving regions of the human genome, and many are hypothesized to function as regulatory elements that drive human-specific gene regulatory programs. We interrogate the in vitro enhancer activity and in vivo epigenetic landscape of more than 3,100 HARs during human neurodevelopment, demonstrating that many HARs appear to act as neurodevelopmental enhancers and that sequence divergence at HARs has largely augmented their neuronal enhancer activity. Furthermore, we demonstrate PPP1R17 to be a putative HAR-regulated gene that has undergone remarkable rewiring of its cell type and developmental expression patterns between non-primates and primates and between non-human primates and humans. Finally, we show that PPP1R17 slows neural progenitor cell cycle progression, paralleling the cell cycle length increase seen predominantly in primate and especially human neurodevelopment. Our findings establish HARs as key components in rewiring human-specific neurodevelopmental gene regulatory programs and provide an integrated resource to study enhancer activity of specific HARs.

Single cell biology-a Keystone Symposia report

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Clinical utility of brain MRS imaging of patients with adult-onset non-cirrhotic hyperammonemia

Adult-onset non-cirrhotic hyperammonemia (NCH) is a rare, but often fatal condition that can result in both reversible and irreversible neurological defects. Here we present five cases of adult-onset non-cirrhotic hyperammonemia wherein brain magnetic resonance spectroscopy (MRS) scans for cerebral glutamine (Gln) and myo-inositol (mI) levels helped guide clinical management. Specifically, we demonstrate that when combined with traditional brain magnetic resonance imaging (MRI) scans, cerebral Gln and mI MRS can help disentangle the reversible from irreversible neurological defects associated with hyperammonemic crisis. Specifically, we demonstrate that whereas an elevated brain MRS Gln level is associated with reversible neurological defects, markedly low mI levels are associated with a risk for irreversible neurological defects such as central pontine myelinolysis. Overall, our findings indicate the utility of brain MRS in guiding clinical care and prognosis in patients with adult-onset non-cirrhotic hyperammonemia.

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.

Single-molecule regulatory architectures captured by chromatin fiber sequencing

Gene regulation is chiefly determined at the level of individual linear chromatin molecules, yet our current understanding of cis-regulatory architectures derives from fragmented sampling of large numbers of disparate molecules. We developed an approach for precisely stenciling the structure of individual chromatin fibers onto their composite DNA templates using nonspecific DNA N⁶-adenine methyltransferases. Single-molecule long-read sequencing of chromatin stencils enabled nucleotide-resolution readout of the primary architecture of multikilobase chromatin fibers (Fiber-seq). Fiber-seq exposed widespread plasticity in the linear organization of individual chromatin fibers and illuminated principles guiding regulatory DNA actuation, the coordinated actuation of neighboring regulatory elements, single-molecule nucleosome positioning, and single-molecule transcription factor occupancy. Our approach and results open new vistas on the primary architecture of gene regulation.

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.

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.

Expanding the clinical spectrum of biallelic ZNF335 variants

ZNF335 plays an essential role in neurogenesis and biallelic variants in ZNF335 have been identified as the cause of severe primary autosomal recessive microcephaly in 2 unrelated families. We describe, herein, 2 additional affected individuals with biallelic ZNF335 variants, 1 individual with a homozygous c.1399 T > C, p.(Cys467Arg) variant, and a second individual with compound heterozygous c.2171_2173delTCT, p.(Phe724del) and c.3998A > G, p.(Glu1333Gly) variants with the latter variant predicted to affect splicing. Whereas the first case presented with early death and a severe phenotype characterized by anterior agyria with prominent extra-axial spaces, absent basal ganglia, and hypoplasia of the brainstem and cerebellum, the second case had a milder clinical presentation with hypomyelination and otherwise preserved brain structures on MRI. Our findings expand the clinical spectrum of ZNF335-associated microcephaly.

Selecting Optimal Peptides for Targeted Proteomic Experiments in Human Plasma Using In Vitro Synthesized Proteins as Analytical Standards

In targeted proteomics, the development of robust methodologies is dependent upon the selection of a set of optimal peptides for each protein-of-interest. Unfortunately, predicting which peptides and respective product ion transitions provide the greatest signal-to-noise ratio in a particular assay matrix is complicated. Using in vitro synthesized proteins as analytical standards, we report here an empirically driven method for the selection of said peptides in a human plasma assay matrix.

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.

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.

Mouse genomic footprinting reveals conservation of transcription factor (TF) recognition repertoires and trans-regulatory circuitry despite massive turnover of DNA elements that contact TFs in vivo.

Having generated genomic DNase I footprinting data of the mouse genome across 25 cell and tissue types, these authors use these data to quantify cis-versus-trans regulatory contributions to mammalian regulatory evolution. They describe more than 600 motifs that collectively are over 95% similar to that recognized in vivo by human transcription factors (TFs). Despite substantial turnover of the cis-regulatory landscape around each TF gene, 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. Conservation between mouse and human TF regulatory networks is particularly similar at the highest organization level. The work was performed as part of the mouse ENCODE project.

Mapping and dynamics of regulatory DNA and transcription factor networks in A. thaliana

Our understanding of gene regulation in plants is constrained by our limited knowledge of plant cis-regulatory DNA and its dynamics. We mapped DNase I hypersensitive sites (DHSs) in A. thaliana seedlings and used genomic footprinting to delineate ∼ 700,000 sites of in vivo transcription factor (TF) occupancy at nucleotide resolution. We show that variation associated with 72 diverse quantitative phenotypes localizes within DHSs. TF footprints encode an extensive cis-regulatory lexicon subject to recent evolutionary pressures, and widespread TF binding within exons may have shaped codon usage patterns. The architecture of A. thaliana TF regulatory networks is strikingly similar to that of animals in spite of diverged regulatory repertoires. We analyzed regulatory landscape dynamics during heat shock and photomorphogenesis, disclosing thousands of environmentally sensitive elements and enabling mapping of key TF regulatory circuits underlying these fundamental responses. Our results provide an extensive resource for the study of A. thaliana gene regulation and functional biology.

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.

Exonic transcription factor binding directs codon choice and affects protein evolution

Genomes contain both a genetic code specifying amino acids and a regulatory code specifying transcription factor (TF) recognition sequences. We used genomic deoxyribonuclease I footprinting to map nucleotide resolution TF occupancy across the human exome in 81 diverse cell types. We found that ~15% of human codons are dual-use codons ("duons") that simultaneously specify both amino acids and TF recognition sites. Duons are highly conserved and have shaped protein evolution, and TF-imposed constraint appears to be a major driver of codon usage bias. Conversely, the regulatory code has been selectively depleted of TFs that recognize stop codons. More than 17% of single-nucleotide variants within duons directly alter TF binding. Pervasive dual encoding of amino acid and regulatory information appears to be a fundamental feature of genome evolution.

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.

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.

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.

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.

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 with an overall success rate of 98%. We show further that targeted proteomic assays developed using our approach facilitate robust in vivo quantification of human transcription factors.