Abstract:In the human genome, it has been estimated that considerably more sequence is under natural selection in non-coding regions [such as transcription-factor binding sites (TF-binding sites) and non-coding RNAs (ncRNAs)] compared to protein-coding ones. However, less attention has been paid to them. To study selective pressure on non-coding elements, we use next-generation sequencing data from the recently completed pilot phase of the 1000 Genomes Project, which, compared to traditional methods, allows for the characterization of a full spectrum of genomic variations, including single-nucleotide polymorphisms (SNPs), short insertions and deletions (indels) and structural variations (SVs). We develop a framework for combining these variation data with non-coding elements, calculating various population-based metrics to compare classes and subclasses of elements, and developing element-aware aggregation procedures to probe the internal structure of an element. Overall, we find that TF-binding sites and ncRNAs are less selectively constrained for SNPs than coding sequences (CDSs), but more constrained than a neutral reference. We also determine that the relative amounts of constraint for the three types of variations are, in general, correlated, but there are some differences: counter-intuitively, TF-binding sites and ncRNAs are more selectively constrained for indels than for SNPs, compared to CDSs. After inspecting the overall properties of a class of elements, we analyze selective pressure on subclasses within an element class, and show that the extent of selection is associated with the genomic properties of each subclass. We find, for instance, that ncRNAs with higher expression levels tend to be under stronger purifying selection, and the actual regions of TF-binding motifs are under stronger selective pressure than the corresponding peak regions. Further, we develop element-aware aggregation plots to analyze selective pressure across the linear structure of an element, with the confidence intervals evaluated using both simple bootstrapping and block bootstrapping techniques. We find, for example, that both micro-RNAs (particularly the seed regions) and their binding targets are under stronger selective pressure for SNPs than their immediate genomic surroundings. In addition, we demonstrate that substitutions in TF-binding motifs inversely correlate with site conservation, and SNPs unfavorable for motifs are under more selective constraints than favorable SNPs. Finally, to further investigate intra-element differences, we show that SVs have the tendency to use distinctive modes and mechanisms when they interact with genomic elements, such as enveloping whole gene(s) rather than disrupting them partially, as well as duplicating TF motifs in tandem.

Identifying deleterious noncoding variation through gain and loss of CTCF binding activity

Functional Assessment of CTCF Sites at Cytokine-Sensing Mammary Enhancers Using CRISPR/Cas9 Gene Editing in Mice.

Systematic Analysis of Binding of Transcription Factors to Noncoding Variants.

Systematic Identification and Annotation of Multiple-Variant Compound Effects at Transcription Factor Binding Sites in Human Genome.

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

Analysis of genomic variation in non-coding elements using population-scale sequencing data from the 1000 Genomes Project.

A complex CTCF binding code defines TAD boundary structure and function

Exploring the Changing Landscape of Cell-to-cell Variation after CTCF Knockdown Via Single Cell RNA-seq

The Functional Consequences of Variation in Transcription Factor Binding

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

DeepCBS: shedding light on the impact of mutations occurring at CTCF binding sites

Abstract 75: CTCF and Rad21 Mediated Regulation of Gene Expression in Human-Induced Pluripotent Stem Cells-derived Endothelial Cells

Contextualising transcription factor binding during embryogenesis using natural sequence variation

DeFine: Deep Convolutional Neural Networks Accurately Quantify Intensities of Transcription Factor-Dna Binding and Facilitate Evaluation of Functional Non-Coding Variants

A deep catalogue of protein-coding variation in 983,578 individuals

Brain and cancer associated binding domain mutations provide insight into CTCF′s relationship with chromatin and its contribution to gene regulation

A method for scoring the cell type-specific impacts of noncoding variants in personal genomes

Regsnps-Asb: A Computational Framework for Identifying Allele-Specific Transcription Factor Binding from ATAC-seq Data

Distinct properties and functions of CTCF revealed by a rapidly inducible degron system

Clustered CTCF binding is an evolutionary mechanism to maintain topologically associating domains

CADD v1.7: using protein language models, regulatory CNNs and other nucleotide-level scores to improve genome-wide variant predictions