maxATAC: Genome-scale transcription-factor binding prediction from ATAC-seq with deep neural networks.

Transcription factors read the genome, fundamentally connecting DNA sequence to gene expression across diverse cell types. Determining how, where, and when TFs bind chromatin will advance our understanding of gene regulatory networks and cellular behavior. The 2017 ENCODE-DREAM in vivo Transcription-Factor Binding Site (TFBS) Prediction Challenge highlighted the value of chromatin accessibility data to TFBS prediction, establishing state-of-the-art methods for TFBS prediction from DNase-seq. However, the more recent Assay-for-Transposase-Accessible-Chromatin (ATAC)-seq has surpassed DNase-seq as the most widely-used chromatin accessibility profiling method. Furthermore, ATAC-seq is the only such technique available at single-cell resolution from standard commercial platforms. While ATAC-seq datasets grow exponentially, suboptimal motif scanning is unfortunately the most common method for TFBS prediction from ATAC-seq. To enable community access to state-of-the-art TFBS prediction from ATAC-seq, we (1) curated an extensive benchmark dataset (127 TFs) for ATAC-seq model training and (2) built "maxATAC", a suite of user-friendly, deep neural network models for genome-wide TFBS prediction from ATAC-seq in any cell type. With models available for 127 human TFs, maxATAC is the largest collection of high-performance TFBS prediction models for ATAC-seq. maxATAC performance extends to primary cells and single-cell ATAC-seq, enabling improved TFBS prediction in vivo. We demonstrate maxATAC's capabilities by identifying TFBS associated with allele-dependent chromatin accessibility at atopic dermatitis genetic risk loci.

An atlas of gene regulatory networks for memory CD4 + T cells in youth and old age.

Aging profoundly affects immune-system function, promoting susceptibility to pathogens, cancers and chronic inflammation. We previously identified a population of IL-10-producing, T follicular helper-like cells (" Tfh10 "), linked to suppressed vaccine responses in aged mice. Here, we integrate single-cell ( sc )RNA-seq, scATAC-seq and genome-scale modeling to characterize Tfh10 - and the full CD4 + memory T cell ( CD4 + TM ) compartment - in young and old mice. We identified 13 CD4 + TM populations, which we validated through cross-comparison to prior scRNA-seq studies. We built gene regulatory networks ( GRNs ) that predict transcription-factor control of gene expression in each T-cell population and how these circuits change with age. Through integration with pan-cell aging atlases, we identified intercellular-signaling networks driving age-dependent changes in CD4 + TM. Our atlas of finely resolved CD4 + TM subsets, GRNs and cell-cell communication networks is a comprehensive resource of predicted regulatory mechanisms operative in memory T cells, presenting new opportunities to improve immune responses in the elderly.