A lineage-specific requirement for YY1 Polycomb Group protein function in early T cell development.

Yin Yang 1 (YY1) is a ubiquitous transcription factor and mammalian Polycomb Group protein (PcG) with important functions for regulating lymphocyte development and stem cell self-renewal. YY1 mediates stable PcG-dependent transcriptional repression via recruitment of PcG proteins that result in histone modifications. Many questions remain unanswered regarding how cell- and tissue-specificity is achieved by PcG proteins. Here, we demonstrate that a conditional knockout of Yy1 in the hematopoietic system results in an early T cell developmental blockage at the double negative (DN) 1 stage with reduced Notch1 signaling. There is a lineage-specific requirement for YY1 PcG function. YY1 PcG domain is required for T and B cell development but not necessary for myeloid cells. YY1 functions in early T cell development are multicomponent and involve both PcG-dependent and -independent regulations. Although YY1 promotes early T cell survival through its PcG function, its function to promote the DN1-to-DN2 transition and Notch1 expression and signaling is independent of its PcG function. Our results reveal how a ubiquitously expressed PcG protein mediates lineage-specific and context-specific functions to control early T cell development.

Data exploration, quality control and statistical analysis of ChIP-exo/nexus experiments.

ChIP-exo/nexus experiments rely on innovative modifications of the commonly used ChIP-seq protocol for high resolution mapping of transcription factor binding sites. Although many aspects of the ChIP-exo data analysis are similar to those of ChIP-seq, these high throughput experiments pose a number of unique quality control and analysis challenges. We develop a novel statistical quality control pipeline and accompanying R/Bioconductor package, ChIPexoQual, to enable exploration and analysis of ChIP-exo and related experiments. ChIPexoQual evaluates a number of key issues including strand imbalance, library complexity, and signal enrichment of data. Assessment of these features are facilitated through diagnostic plots and summary statistics computed over regions of the genome with varying levels of coverage. We evaluated our QC pipeline with both large collections of public ChIP-exo/nexus data and multiple, new ChIP-exo datasets from Escherichia coli. ChIPexoQual analysis of these datasets resulted in guidelines for using these QC metrics across a wide range of sequencing depths and provided further insights for modelling ChIP-exo data.

Predicting gene expression in massively parallel reporter assays: A comparative study.

In many human diseases, associated genetic changes tend to occur within noncoding regions, whose effect might be related to transcriptional control. A central goal in human genetics is to understand the function of such noncoding regions: given a region that is statistically associated with changes in gene expression (expression quantitative trait locus [eQTL]), does it in fact play a regulatory role? And if so, how is this role "coded" in its sequence? These questions were the subject of the Critical Assessment of Genome Interpretation eQTL challenge. Participants were given a set of sequences that flank eQTLs in humans and were asked to predict whether these are capable of regulating transcription (as evaluated by massively parallel reporter assays), and whether this capability changes between alternative alleles. Here, we report lessons learned from this community effort. By inspecting predictive properties in isolation, and conducting meta-analysis over the competing methods, we find that using chromatin accessibility and transcription factor binding as features in an ensemble of classifiers or regression models leads to the most accurate results. We then characterize the loci that are harder to predict, putting the spotlight on areas of weakness, which we expect to be the subject of future studies.

Epstein-Barr Virus Nuclear Antigen 3 (EBNA3) Proteins Regulate EBNA2 Binding to Distinct RBPJ Genomic Sites.

UNLABELLED: Latent infection of B lymphocytes by Epstein-Barr virus (EBV) in vitro results in their immortalization into lymphoblastoid cell lines (LCLs); this latency program is controlled by the EBNA2 viral transcriptional activator, which targets promoters via RBPJ, a DNA binding protein in the Notch signaling pathway. Three other EBNA3 proteins (EBNA3A, EBNA3B, and EBNA3C) interact with RBPJ to regulate cell gene expression. The mechanism by which EBNAs regulate different genes via RBPJ remains unclear. Our chromatin immunoprecipitation with deep sequencing (ChIP-seq) analysis of the EBNA3 proteins analyzed in concert with prior EBNA2 and RBPJ data demonstrated that EBNA3A, EBNA3B, and EBNA3C bind to distinct, partially overlapping genomic locations. Although RBPJ interaction is critical for EBNA3A and EBNA3C growth effects, only 30 to 40% of EBNA3-bound sites colocalize with RBPJ. Using LCLs conditional for EBNA3A or EBNA3C activity, we demonstrate that EBNA2 binding at sites near EBNA3A- or EBNA3C-regulated genes is specifically regulated by the respective EBNA3. To investigate EBNA3 binding specificity, we identified sequences and transcription factors enriched at EBNA3A-, EBNA3B-, and EBNA3C-bound sites. This confirmed the prior observation that IRF4 is enriched at EBNA3A- and EBNA3C-bound sites and revealed IRF4 enrichment at EBNA3B-bound sites. Using IRF4-negative BJAB cells, we demonstrate that IRF4 is essential for EBNA3C, but not EBNA3A or EBNA3B, binding to specific sites. These results support a model in which EBNA2 and EBNA3s compete for distinct subsets of RBPJ sites to regulate cell genes and where EBNA3 subset specificity is determined by interactions with other cell transcription factors. IMPORTANCE: Epstein-Barr virus (EBV) latent gene products cause human cancers and transform B lymphocytes into immortalized lymphoblastoid cell lines in vitro. EBV nuclear antigens (EBNAs) and membrane proteins constitutively activate pathways important for lymphocyte growth and survival. An important unresolved question is how four different EBNAs (EBNA2, -3A, -3B, and -3C) exert unique effects via a single transcription factor, RBPJ. Here, we report that each EBNA binds to distinct but partially overlapping sets of genomic sites. EBNA3A and EBNA3C specifically regulate EBNA2's access to different RBPJ sites, providing a mechanism by which each EBNA can regulate distinct cell genes. We show that IRF4, an essential regulator of B cell differentiation, is critical for EBNA3C binding specificity; EBNA3A and EBNA3B specificities are likely due to interactions with other cell transcription factors. EBNA3 titration of EBNA2 transcriptional function at distinct sites likely limits cell defenses that would be triggered by unchecked EBNA2 prooncogenic activity.

Perm-seq: Mapping Protein-DNA Interactions in Segmental Duplication and Highly Repetitive Regions of Genomes with Prior-Enhanced Read Mapping.

Segmental duplications and other highly repetitive regions of genomes contribute significantly to cells' regulatory programs. Advancements in next generation sequencing enabled genome-wide profiling of protein-DNA interactions by chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq). However, interactions in highly repetitive regions of genomes have proven difficult to map since short reads of 50-100 base pairs (bps) from these regions map to multiple locations in reference genomes. Standard analytical methods discard such multi-mapping reads and the few that can accommodate them are prone to large false positive and negative rates. We developed Perm-seq, a prior-enhanced read allocation method for ChIP-seq experiments, that can allocate multi-mapping reads in highly repetitive regions of the genomes with high accuracy. We comprehensively evaluated Perm-seq, and found that our prior-enhanced approach significantly improves multi-read allocation accuracy over approaches that do not utilize additional data types. The statistical formalism underlying our approach facilitates supervising of multi-read allocation with a variety of data sources including histone ChIP-seq. We applied Perm-seq to 64 ENCODE ChIP-seq datasets from GM12878 and K562 cells and identified many novel protein-DNA interactions in segmental duplication regions. Our analysis reveals that although the protein-DNA interactions sites are evolutionarily less conserved in repetitive regions, they share the overall sequence characteristics of the protein-DNA interactions in non-repetitive regions.

Genome-wide association study (GWAS) of circulating vitamin D outcomes among individuals of African ancestry.

BACKGROUND: Vitamin D deficiency is more common among African-ancestry individuals and may be associated with adverse health outcomes. Vitamin D binding protein (VDBP) regulates concentrations of biologically active vitamin D. OBJECTIVE: We conducted genome-wide association study (GWAS) of VDBP and 25-hydroxyvitamin D among African-ancestry individuals. METHODS: Data were collected from 2,602 African American adults from the Southern Community Cohort Study (SCCS) and 6,934 African- or Caribbean-ancestry adults from the UK Biobank. Serum VDBP concentrations were available only in the SCCS and were measured by using the Polyclonal Human VDBP ELISA kit. Serum 25-hydroxyvitamin D concentrations for both study samples were measured by using Diasorin Liason, a chemiluminescent immunoassay. Participants were genotyped for single nucleotide polymorphisms (SNPs) with genome-wide coverage by using Illumina or Affymetrix platforms. Fine-mapping analysis was performed by using forward stepwise linear regression models including all variants with P value < 5 × 10-8 and within 250 kbps of a lead SNP. RESULTS: We identified 4 loci notably associated with VDBP concentrations in the SCCS population: rs7041 (per allele ß = 0.61 µg/mL, SE = 0.05, P = 1.4 × 10-48) and rs842998 (per allele ß = 0.39 µg/mL, SE = 0.03, P = 4.0 × 10-31) in GC, rs8427873 (per allele ß = 0.31 µg/mL, SE = 0.04, P = 3.0 × 10-14) near GC and rs11731496 (per allele ß = 0.21 µg/mL, SE = 0.03, P = 3.6 × 10-11) in between GC and NPFFR2. In conditional analyses, which included the above-mentioned SNPs, only rs7041 remained notable (P = 4.1 × 10-21). SNP rs4588 in GC was the only GWAS-identified SNP associated with 25-hydroxyvitamin D concentration. Among UK Biobank participants: per allele ß = -0.11 µg/mL, SE = 0.01, P = 1.5 × 10-13; in the SCCS: per allele ß = -0.12 µg/mL, SE = 0.06, P = 2.8 × 10-02). rs7041 and rs4588 are functional SNPs that influence the binding affinity of VDBP to 25-hydroxyvitamin D. CONCLUSIONS: Our results were in line with previous studies conducted in European-ancestry populations, showing that GC, the gene that directly encodes for VDBP, would be important for VDBP and 25-hydroxyvitamin D concentrations. The current study extends our knowledge of the genetics of vitamin D in diverse populations.

Epstein-Barr Virus Infection Promotes Epithelial Cell Growth by Attenuating Differentiation-Dependent Exit from the Cell Cycle.

Epstein-Barr virus (EBV) is a human herpesvirus that is associated with lymphomas as well as nasopharyngeal and gastric carcinomas. Although carcinomas account for almost 90% of EBV-associated cancers, progress in examining EBV's role in their pathogenesis has been limited by difficulty in establishing latent infection in nontransformed epithelial cells. Recently, EBV infection of human telomerase reverse transcriptase (hTERT)-immortalized normal oral keratinocytes (NOKs) has emerged as a model that recapitulates aspects of EBV infection in vivo, such as differentiation-associated viral replication. Using uninfected NOKs and NOKs infected with the Akata strain of EBV (NOKs-Akata), we examined changes in gene expression due to EBV infection and differentiation. Latent EBV infection produced very few significant gene expression changes in undifferentiated NOKs but significantly reduced the extent of differentiation-induced gene expression changes. Gene set enrichment analysis revealed that differentiation-induced downregulation of the cell cycle and metabolism pathways was markedly attenuated in NOKs-Akata relative to that in uninfected NOKs. We also observed that pathways induced by differentiation were less upregulated in NOKs-Akata. We observed decreased differentiation markers and increased suprabasal MCM7 expression in NOKs-Akata versus NOKs when both were grown in raft cultures, consistent with our transcriptome sequencing (RNA-seq) results. These effects were also observed in NOKs infected with a replication-defective EBV mutant (AkataΔRZ), implicating mechanisms other than lytic-gene-induced host shutoff. Our results help to define the mechanisms by which EBV infection alters keratinocyte differentiation and provide a basis for understanding the role of EBV in epithelial cancers.IMPORTANCE Latent infection by Epstein-Barr virus (EBV) is an early event in the development of EBV-associated carcinomas. In oral epithelial tissues, EBV establishes a lytic infection of differentiated epithelial cells to facilitate the spread of the virus to new hosts. Because of limitations in existing model systems, the effects of latent EBV infection on undifferentiated and differentiating epithelial cells are poorly understood. Here, we characterize latent infection of an hTERT-immortalized oral epithelial cell line (NOKs). We find that although EBV expresses a latency pattern similar to that seen in EBV-associated carcinomas, infection of undifferentiated NOKs results in differential expression of a small number of host genes. In differentiating NOKs, however, EBV has a more substantial effect, reducing the extent of differentiation and delaying the exit from the cell cycle. This effect may synergize with preexisting cellular abnormalities to prevent exit from the cell cycle, representing a critical step in the development of cancer.