Allele specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs.

Transcription Factors (TFs) regulate gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Since TF occupancy is driven in part by recognition of DNA sequence, genetic variation can influence TF-DNA associations and gene regulation. To identify variants that impact TF binding in human brain tissues, we assessed allele specific binding (ASB) at heterozygous variants for 94 TFs in 9 brain regions from two donors. Leveraging graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signals between alleles at heterozygous variants within each brain region and identified thousands of variants exhibiting ASB for at least one TF. ASB reproducibility was measured by comparisons between independent experiments both within and between donors. We found that rarer alleles in the general population more frequently led to reduced TF binding, whereas common variation had an equal likelihood of increasing or decreasing binding. Motif analysis revealed TF-specific effects, with ASB variants for certain TFs displaying a greater incidence of motif alterations, as well as enrichments for variants under purifying selection. Notably, neuron-specific cis-regulatory elements (cCREs) showed depletion for ASB variants. We identified 2,670 ASB variants with prior evidence of allele-specific gene expression in the brain from GTEx data and observed increasing eQTL effect direction concordance as ASB significance increases. These results provide a valuable and unique resource for mechanistic analysis of cis-regulatory variation in human brain tissue.

Characterization of human transcription factor function and patterns of gene regulation in HepG2 cells.

Transcription factors (TFs) are trans-acting proteins that bind cis-regulatory elements (CREs) in DNA to control gene expression. Here, we analyzed the genomic localization profiles of 529 sequence-specific TFs and 151 cofactors and chromatin regulators in the human cancer cell line HepG2, for a total of 680 broadly termed DNA-associated proteins (DAPs). We used this deep collection to model each TF's impact on gene expression, and identified a cohort of 26 candidate transcriptional repressors. We examine high occupancy target (HOT) sites in the context of three-dimensional genome organization and show biased motif placement in distal-promoter connections involving HOT sites. We also found a substantial number of closed chromatin regions with multiple DAPs bound, and explored their properties, finding that a MAFF/MAFK TF pair correlates with transcriptional repression. Altogether, these analyses provide novel insights into the regulatory logic of the human cell line HepG2 genome and show the usefulness of large genomic analyses for elucidation of individual TF functions.

Multiomic profiling of transcription factor binding and function in human brain.

Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.

Contributions of rare and common variation to early-onset and atypical dementia risk.

We collected and analyzed genomic sequencing data from individuals with clinician- diagnosed early-onset or atypical dementia. Thirty-two patients were previously described, with sixty-eight newly described in this report. Of those sixty-eight, sixty-two patients reported Caucasian, non-Hispanic ethnicity and six reported as African American, non-Hispanic. Fifty-three percent of patients had a returnable variant. Five patients harbored a pathogenic variant as defined by the American College of Medical Genetics criteria for pathogenicity. A polygenic risk score was calculated for Alzheimer's patients in the total cohort and compared to the scores of a late-onset Alzheimer's cohort and a control set. Patients with early-onset Alzheimer's had higher non- APOE polygenic risk scores than patients with late onset Alzheimer's, supporting the conclusion that both rare and common genetic variation associate with early-onset neurodegenerative disease risk.

Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs.

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

Contributions of rare and common variation to early-onset and atypical dementia risk.

We collected and analyzed genomic sequencing data from individuals with clinician-diagnosed early-onset or atypical dementia. Thirty-two patients were previously described, with 68 newly described in this report. Of those 68, 62 patients self-reported white, non-Hispanic ethnicity and 6 reported as African-American, non-Hispanic. Fifty-three percent of patients had a returnable variant. Five patients harbored a pathogenic variant as defined by the American College of Medical Genetics criteria for pathogenicity. A polygenic risk score (PRS) was calculated for Alzheimer's patients in the total cohort and compared to the scores of a late-onset Alzheimer's cohort and a control set. Patients with early-onset Alzheimer's had higher non-APOE PRSs than patients with late-onset Alzheimer's, supporting the conclusion that both rare and common genetic variation associate with early-onset neurodegenerative disease risk.

A subset of chicken statoacoustic ganglion neurites are repelled by Slit1 and Slit2.

Mechanosensory hair cells in the chicken inner ear are innervated by bipolar afferent neurons of the statoacoustic ganglion (SAG). During development, individual SAG neurons project their peripheral process to only one of eight distinct sensory organs. These neuronal subtypes may respond differently to guidance cues as they explore the periphery in search of their target. Previous gene expression data suggested that Slit repellants might channel SAG neurites into the sensory primordia, based on the presence of robo transcripts in the neurons and the confinement of slit transcripts to the flanks of the prosensory domains. This led to the prediction that excess Slit proteins would impede the outgrowth of SAG neurites. As predicted, axonal projections to the primordium of the anterior crista were reduced 2-3 days after electroporation of either slit1 or slit2 expression plasmids into the anterior pole of the otocyst on embryonic day 3 (E3). The posterior crista afferents, which normally grow through and adjacent to slit expression domains as they are navigating towards the posterior pole of the otocyst, did not show Slit responsiveness when similarly challenged by ectopic delivery of slit to their targets. The sensitivity to ectopic Slits shown by the anterior crista afferents was more the exception than the rule: responsiveness to Slits was not observed when the entire E4 SAG was challenged with Slits for 40 h in vitro. The corona of neurites emanating from SAG explants was unaffected by the presence of purified human Slit1 and Slit2 in the culture medium. Reduced axon outgrowth from E8 olfactory bulbs cultured under similar conditions for 24 h confirmed bioactivity of purified human Slits on chicken neurons. In summary, differential sensitivity to Slit repellents may influence the directional outgrowth of otic axons toward either the anterior or posterior otocyst.