Search

Genetic variants and de novo mutations in regulatory regions of the genome are typically discovered by whole-genome sequencing (WGS), however WGS is expensive and most WGS reads come from non-regulatory regions. The Assay for Transposase-Accessible Chromatin (ATAC-seq) generates reads from regulatory sequences and could potentially be used as a low-cost 'capture' method for regulatory variant discovery, but its use for this purpose has not been systematically evaluated. Here we apply seven variant callers to bulk and single-cell ATAC-seq data and evaluate their ability to identify single nucleotide variants (SNVs) and insertions/deletions (indels). In addition, we develop an ensemble classifier, VarCA, which combines features from individual variant callers to predict variants. The Genome Analysis Toolkit (GATK) is the best-performing individual caller with precision/recall on a bulk ATAC test dataset of 0.92/0.97 for SNVs and 0.87/0.82 for indels within ATAC-seq peak regions with at least 10 reads. On bulk ATAC-seq reads, VarCA achieves superior performance with precision/recall of 0.99/0.95 for SNVs and 0.93/0.80 for indels. On single-cell ATAC-seq reads, VarCA attains precision/recall of 0.98/0.94 for SNVs and 0.82/0.82 for indels. In summary, ATAC-seq reads can be used to accurately discover non-coding regulatory variants in the absence of whole-genome sequencing data and our ensemble method, VarCA, has the best overall performance., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Heterochromatin has essential functions in maintaining chromosome structure, in protecting genome integrity and in stabilizing gene expression programs. Heterochromatin is often nucleated by underlying DNA repeat sequences, such as major satellite repeats (MSR) and long interspersed nuclear elements (LINE). In order to establish heterochromatin, MSR and LINE elements need to be transcriptionally competent and generate non-coding repeat RNA that remain chromatin associated. We explored whether these heterochromatic RNA, similar to DNA and histones, may be methylated, particularly for 5-methylcytosine (5mC) or methyl-6-adenosine (m6A). Our analysis in mouse ES cells identifies only background level of 5mC but significant enrichment for m6A on heterochromatic RNA. Moreover, MSR transcripts are a novel target for m6A RNA modification, and their m6A RNA enrichment is decreased in ES cells that are mutant for Mettl3 or Mettl14, which encode components of a central RNA methyltransferase complex. Importantly, MSR transcripts that are partially deficient in m6A RNA methylation display impaired chromatin association and have a reduced potential to form RNA:DNA hybrids. We propose that m6A modification of MSR RNA will enhance the functions of MSR repeat transcripts to stabilize mouse heterochromatin., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)