C5a-licensed phagocytes drive sterilizing immunity during systemic fungal infection.

Systemic candidiasis is a common, high-mortality, nosocomial fungal infection. Unexpectedly, it has emerged as a complication of anti-complement C5-targeted monoclonal antibody treatment, indicating a critical niche for C5 in antifungal immunity. We identified transcription of complement system genes as the top biological pathway induced in candidemic patients and as predictive of candidemia. Mechanistically, C5a-C5aR1 promoted fungal clearance and host survival in a mouse model of systemic candidiasis by stimulating phagocyte effector function and ERK- and AKT-dependent survival in infected tissues. C5ar1 ablation rewired macrophage metabolism downstream of mTOR, promoting their apoptosis and enhancing mortality through kidney injury. Besides hepatocyte-derived C5, local C5 produced intrinsically by phagocytes provided a key substrate for antifungal protection. Lower serum C5a concentrations or a C5 polymorphism that decreases leukocyte C5 expression correlated independently with poor patient outcomes. Thus, local, phagocyte-derived C5 production licenses phagocyte antimicrobial function and confers innate protection during systemic fungal infection.

Autocrine vitamin D signaling switches off pro-inflammatory programs of TH1 cells.

The molecular mechanisms governing orderly shutdown and retraction of CD4+ type 1 helper T (TH1) cell responses remain poorly understood. Here we show that complement triggers contraction of TH1 responses by inducing intrinsic expression of the vitamin D (VitD) receptor and the VitD-activating enzyme CYP27B1, permitting T cells to both activate and respond to VitD. VitD then initiated the transition from pro-inflammatory interferon-γ+ TH1 cells to suppressive interleukin-10+ cells. This process was primed by dynamic changes in the epigenetic landscape of CD4+ T cells, generating super-enhancers and recruiting several transcription factors, notably c-JUN, STAT3 and BACH2, which together with VitD receptor shaped the transcriptional response to VitD. Accordingly, VitD did not induce interleukin-10 expression in cells with dysfunctional BACH2 or STAT3. Bronchoalveolar lavage fluid CD4+ T cells of patients with COVID-19 were TH1-skewed and showed de-repression of genes downregulated by VitD, from either lack of substrate (VitD deficiency) and/or abnormal regulation of this system.

TCF-1 controls Treg cell functions that regulate inflammation, CD8+ T cell cytotoxicity and severity of colon cancer.

The transcription factor TCF-1 is essential for the development and function of regulatory T (Treg) cells; however, its function is poorly understood. Here, we show that TCF-1 primarily suppresses transcription of genes that are co-bound by Foxp3. Single-cell RNA-sequencing analysis identified effector memory T cells and central memory Treg cells with differential expression of Klf2 and memory and activation markers. TCF-1 deficiency did not change the core Treg cell transcriptional signature, but promoted alternative signaling pathways whereby Treg cells became activated and gained gut-homing properties and characteristics of the TH17 subset of helper T cells. TCF-1-deficient Treg cells strongly suppressed T cell proliferation and cytotoxicity, but were compromised in controlling CD4+ T cell polarization and inflammation. In mice with polyposis, Treg cell-specific TCF-1 deficiency promoted tumor growth. Consistently, tumor-infiltrating Treg cells of patients with colorectal cancer showed lower TCF-1 expression and increased TH17 expression signatures compared to adjacent normal tissue and circulating T cells. Thus, Treg cell-specific TCF-1 expression differentially regulates TH17-mediated inflammation and T cell cytotoxicity, and can determine colorectal cancer outcome.

Loss of CD4+ T cell-intrinsic arginase 1 accelerates Th1 response kinetics and reduces lung pathology during influenza infection.

Arginase 1 (Arg1), the enzyme catalyzing the conversion of arginine to ornithine, is a hallmark of IL-10-producing immunoregulatory M2 macrophages. However, its expression in T cells is disputed. Here, we demonstrate that induction of Arg1 expression is a key feature of lung CD4+ T cells during mouse in vivo influenza infection. Conditional ablation of Arg1 in CD4+ T cells accelerated both virus-specific T helper 1 (Th1) effector responses and its resolution, resulting in efficient viral clearance and reduced lung pathology. Using unbiased transcriptomics and metabolomics, we found that Arg1-deficiency was distinct from Arg2-deficiency and caused altered glutamine metabolism. Rebalancing this perturbed glutamine flux normalized the cellular Th1 response. CD4+ T cells from rare ARG1-deficient patients or CRISPR-Cas9-mediated ARG1-deletion in healthy donor cells phenocopied the murine cellular phenotype. Collectively, CD4+ T cell-intrinsic Arg1 functions as an unexpected rheostat regulating the kinetics of the mammalian Th1 lifecycle with implications for Th1-associated tissue pathologies.

Human retinoic acid-regulated CD161+ regulatory T cells support wound repair in intestinal mucosa.

Repair of tissue damaged during inflammatory processes is key to the return of local homeostasis and restoration of epithelial integrity. Here we describe CD161+ regulatory T (Treg) cells as a distinct, highly suppressive population of Treg cells that mediate wound healing. These Treg cells were enriched in intestinal lamina propria, particularly in Crohn's disease. CD161+ Treg cells had an all-trans retinoic acid (ATRA)-regulated gene signature, and CD161 expression on Treg cells was induced by ATRA, which directly regulated the CD161 gene. CD161 was co-stimulatory, and ligation with the T cell antigen receptor induced cytokines that accelerated the wound healing of intestinal epithelial cells. We identified a transcription-factor network, including BACH2, RORγt, FOSL2, AP-1 and RUNX1, that controlled expression of the wound-healing program, and found a CD161+ Treg cell signature in Crohn's disease mucosa associated with reduced inflammation. These findings identify CD161+ Treg cells as a population involved in controlling the balance between inflammation and epithelial barrier healing in the gut.

BACH2 immunodeficiency illustrates an association between super-enhancers and haploinsufficiency.

The transcriptional programs that guide lymphocyte differentiation depend on the precise expression and timing of transcription factors (TFs). The TF BACH2 is essential for T and B lymphocytes and is associated with an archetypal super-enhancer (SE). Single-nucleotide variants in the BACH2 locus are associated with several autoimmune diseases, but BACH2 mutations that cause Mendelian monogenic primary immunodeficiency have not previously been identified. Here we describe a syndrome of BACH2-related immunodeficiency and autoimmunity (BRIDA) that results from BACH2 haploinsufficiency. Affected subjects had lymphocyte-maturation defects that caused immunoglobulin deficiency and intestinal inflammation. The mutations disrupted protein stability by interfering with homodimerization or by causing aggregation. We observed analogous lymphocyte defects in Bach2-heterozygous mice. More generally, we observed that genes that cause monogenic haploinsufficient diseases were substantially enriched for TFs and SE architecture. These findings reveal a previously unrecognized feature of SE architecture in Mendelian diseases of immunity: heterozygous mutations in SE-regulated genes identified by whole-exome/genome sequencing may have greater significance than previously recognized.

Diapedesis-Induced Integrin Signaling via LFA-1 Facilitates Tissue Immunity by Inducing Intrinsic Complement C3 Expression in Immune Cells.

Intrinsic complement C3 activity is integral to human T helper type 1 (Th1) and cytotoxic T cell responses. Increased or decreased intracellular C3 results in autoimmunity and infections, respectively. The mechanisms regulating intracellular C3 expression remain undefined. We identified complement, including C3, as among the most significantly enriched biological pathway in tissue-occupying cells. We generated C3-reporter mice and confirmed that C3 expression was a defining feature of tissue-immune cells, including T cells and monocytes, occurred during transendothelial diapedesis, and depended on integrin lymphocyte-function-associated antigen 1 (LFA-1) signals. Immune cells from patients with leukocyte adhesion deficiency type 1 (LAD-1) had reduced C3 transcripts and diminished effector activities, which could be rescued proportionally by intracellular C3 provision. Conversely, increased C3 expression by T cells from arthritis patients correlated with disease severity. Our study defines integrins as key controllers of intracellular complement, demonstrates that perturbations in the LFA-1-C3-axis contribute to primary immunodeficiency, and identifies intracellular C3 as biomarker of severity in autoimmunity.

Long-term organic carbon preservation enhanced by iron and manganese.

The balance between degradation and preservation of sedimentary organic carbon (OC) is important for global carbon and oxygen cycles1. The relative importance of different mechanisms and environmental conditions contributing to marine sedimentary OC preservation, however, remains unclear2-8. Simple organic molecules can be geopolymerized into recalcitrant forms by means of the Maillard reaction5, although reaction kinetics at marine sedimentary temperatures are thought to be slow9,10. More recent work in terrestrial systems suggests that the reaction can be catalysed by manganese minerals11-13, but the potential for the promotion of geopolymerized OC formation at marine sedimentary temperatures is uncertain. Here we present incubation experiments and find that iron and manganese ions and minerals abiotically catalyse the Maillard reaction by up to two orders of magnitude at temperatures relevant to continental margins where most preservation occurs4. Furthermore, the chemical signature of the reaction products closely resembles dissolved and total OC found in continental margin sediments globally. With the aid of a pore-water model14, we estimate that iron- and manganese-catalysed transformation of simple organic molecules into complex macromolecules might generate on the order of approximately 4.1 Tg C yr-1 for preservation in marine sediments. In the context of perhaps only about 63 Tg C yr-1 variation in sedimentary organic preservation over the past 300 million years6, we propose that variable iron and manganese inputs to the ocean could exert a substantial but hitherto unexplored impact on global OC preservation over geological time.

Kaposi's sarcoma herpesvirus latency-associated nuclear antigen broadly regulates viral gene expression and is essential for lytic infection.

Kaposi's sarcoma herpesvirus (KSHV) is a leading cause of malignancy in AIDS and current therapies are limited. Like all herpesviruses, KSHV infection can be latent or lytic. KSHV latency-associated nuclear antigen (LANA) is essential for viral genome persistence during latent infection. LANA also maintains latency by antagonizing expression and function of the KSHV lytic switch protein, RTA. Here, we find LANA null KSHV is not capable of lytic replication, indicating a requirement for LANA. While LANA promoted both lytic and latent gene expression in cells partially permissive for lytic infection, it repressed expression in non-permissive cells. Importantly, forced RTA expression in non-permissive cells led to induction of lytic infection and LANA switched to promote, rather than repress, most lytic viral gene expression. When basal viral gene expression levels were high, LANA promoted expression, but repressed expression at low basal levels unless RTA expression was forcibly induced. LANA's effects were broad, but virus gene specific, extending to an engineered, recombinant viral GFP under control of host EF1α promoter, but not to host EF1α. Together, these results demonstrate that, in addition to its essential role in genome maintenance, LANA broadly regulates viral gene expression, and is required for high levels of lytic gene expression during lytic infection. Strategies that target LANA are expected to abolish KSHV infection.

Non-coding RNAs in immunoregulation and autoimmunity: Technological advances and critical limitations.

Immune cell function is critically dependent on precise control over transcriptional output from the genome. In this respect, integration of environmental signals that regulate gene expression, specifically by transcription factors, enhancer DNA elements, genome topography and non-coding RNAs (ncRNAs), are key components. The first three have been extensively investigated. Even though non-coding RNAs represent the vast majority of cellular RNA species, this class of RNA remains historically understudied. This is partly because of a lag in technological and bioinformatic innovations specifically capable of identifying and accurately measuring their expression. Nevertheless, recent progress in this domain has enabled a profusion of publications identifying novel sub-types of ncRNAs and studies directly addressing the function of ncRNAs in human health and disease. Many ncRNAs, including circular and enhancer RNAs, have now been demonstrated to play key functions in the regulation of immune cells and to show associations with immune-mediated diseases. Some ncRNAs may function as biomarkers of disease, aiding in diagnostics and in estimating response to treatment, while others may play a direct role in the pathogenesis of disease. Importantly, some are relatively stable and are amenable to therapeutic targeting, for example through gene therapy. Here, we provide an overview of ncRNAs and review technological advances that enable their study and hold substantial promise for the future. We provide context-specific examples by examining the associations of ncRNAs with four prototypical human autoimmune diseases, specifically rheumatoid arthritis, psoriasis, inflammatory bowel disease and multiple sclerosis. We anticipate that the utility and mechanistic roles of these ncRNAs in autoimmunity will be further elucidated in the near future.

A comprehensive single cell data analysis of lymphoblastoid cells reveals the role of super-enhancers in maintaining EBV latency.

We probed the lifecycle of Epstein-Barr virus (EBV) on a cell-by-cell basis using single cell RNA sequencing (scRNA-seq) data from nine publicly available lymphoblastoid cell lines (LCLs). While the majority of LCLs comprised cells containing EBV in the latent phase, two other clusters of cells were clearly evident and were distinguished by distinct expression of host and viral genes. Notably, both were high expressors of EBV LMP1/BNLF2 and BZLF1 compared to another cluster that expressed neither gene. The two novel clusters differed from each other in their expression of EBV lytic genes, including glycoprotein gene GP350. The first cluster, comprising GP350- LMP1hi cells, expressed high levels of HIF1A and was transcriptionally regulated by HIF1-α. Treatment of LCLs with Pevonedistat, a drug that enhances HIF1-α signaling, markedly induced this cluster. The second cluster, containing GP350+ LMP1hi cells, expressed EBV lytic genes. Host genes that are controlled by super-enhancers (SEs), such as transcription factors MYC and IRF4, had the lowest expression in this cluster. Functionally, the expression of genes regulated by MYC and IRF4 in GP350+ LMP1hi cells were lower compared to other cells. Indeed, induction of EBV lytic reactivation in EBV+ AKATA reduced the expression of these SE-regulated genes. Furthermore, CRISPR-mediated perturbation of the MYC or IRF4 SEs in LCLs induced the lytic EBV gene expression, suggesting that host SEs and/or SE target genes are required for maintenance of EBV latency. Collectively, our study revealed EBV-associated heterogeneity among LCLs that may have functional consequence on host and viral biology.

STAT5 Represses a STAT3-Independent Th17-like Program during Th9 Cell Differentiation.

IL-9-producing Th cells, termed Th9 cells, contribute to immunity against parasites and cancers but have detrimental roles in allergic disease and colitis. Th9 cells differentiate in response to IL-4 and TGF-ß, but these signals are insufficient to drive Th9 differentiation in the absence of IL-2. IL-2-induced STAT5 activation is required for chromatin accessibility within Il9 enhancer and promoter regions and directly transactivates the Il9 locus. STAT5 also suppresses gene expression during Th9 cell development, but these roles are less well defined. In this study, we demonstrate that human allergy-associated Th9 cells exhibited a signature of STAT5-mediated gene repression that is associated with the silencing of a Th17-like transcriptional signature. In murine Th9 cell differentiation, blockade of IL-2/STAT5 signaling induced the expression of IL-17 and the Th17-associated transcription factor Rorγt. However, IL-2-deprived Th9 cells did not exhibit a significant Th17- or STAT3-associated transcriptional signature. Consistent with these observations, differentiation of IL-17-producing cells under these conditions was STAT3-independent but did require Rorγt and BATF. Furthermore, ectopic expression of Rorγt and BATF partially rescued IL-17 production in STAT3-deficient Th17 cells, highlighting the importance of these factors in this process. Although STAT3 was not required for the differentiation of IL-17-producing cells under IL-2-deprived Th9 conditions, their prolonged survival was STAT3-dependent, potentially explaining why STAT3-independent IL-17 production is not commonly observed in vivo. Together, our data suggest that IL-2/STAT5 signaling plays an important role in controlling the balance of a Th9 versus a Th17-like differentiation program in vitro and in allergic disease.

Folic acid-mediated fibrosis is driven by C5a receptor 1-mediated activation of kidney myeloid cells.

We have previously reported that increased expression and activation of kidney cell complement components play an important role in the pathogenesis of renal scarring. Here, we used floxed green fluorescent protein (GFP)-C5a receptor 1 (C5aR1) knockin mice (GFP-C5ar1fl/fl) and the model of folic acid (FA)-induced kidney injury to define the cell types and potential mechanisms by which increased C5aR1 activation leads to fibrosis. Using flow cytometry and confocal microscopy, we identified macrophages as the major interstitial cell type showing increased expression of C5aR1 in FA-treated mice. C5ar1fl/fl.Lyz2Cre+/- mice, in which C5aR1 has been specifically deleted in lysozyme M-expressing myeloid cells, experienced reduced fibrosis compared with control C5ar1fl/fl mice. Examination of C5aR1-expressing macrophage transcriptomes by gene set enrichment analysis demonstrated that these cells were enriched in pathways corresponding to the complement cascade, collagen formation, and the NABA matrisome, strongly pointing to their critical roles in tissue repair/scarring. Since C5aR1 was also detected in a small population of platelet-derived growth factor receptor-ß+ GFP+ cells, we developed C5ar1fl/fl.Foxd1Cre+/- mice, in which C5aR1 is deleted specifically in pericytes, and found reduced FA-induced fibrosis. Primary cell cultures of platelet-derived growth factor receptor-ß+ pericytes isolated from FA-treated C5ar1fl/fl.Foxd1Cre+/- mice showed reduced secretion of several cytokines, including IL-6 and macrophage inflammatory protein-2, compared with pericytes isolated from FA-treated control GFP-C5ar1fl/fl mice. Collectively, these data imply that C5a/C5aR1 axis activation primarily in interstitial cells contributes to the development of renal fibrosis.NEW & NOTEWORTHY This study used novel green fluorescent protein C5a receptor 1 floxed mice and the model of folic acid-mediated kidney fibrosis to demonstrate the pathogenic role of increased expression of this complement receptor on macrophages.

Host-Virus Chimeric Events in SARS-CoV-2-Infected Cells Are Infrequent and Artifactual.

The pathogenic mechanisms underlying severe SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection remain largely unelucidated. High-throughput sequencing technologies that capture genome and transcriptome information are key approaches to gain detailed mechanistic insights from infected cells. These techniques readily detect both pathogen- and host-derived sequences, providing a means of studying host-pathogen interactions. Recent studies have reported the presence of host-virus chimeric (HVC) RNA in transcriptome sequencing (RNA-seq) data from SARS-CoV-2-infected cells and interpreted these findings as evidence of viral integration in the human genome as a potential pathogenic mechanism. Since SARS-CoV-2 is a positive-sense RNA virus that replicates in the cytoplasm, it does not have a nuclear phase in its life cycle. Thus, it is biologically unlikely to be in a location where splicing events could result in genome integration. Therefore, we investigated the biological authenticity of HVC events. In contrast to true biological events like mRNA splicing and genome rearrangement events, which generate reproducible chimeric sequencing fragments across different biological isolates, we found that HVC events across >100 RNA-seq libraries from patients with coronavirus disease 2019 (COVID-19) and infected cell lines were highly irreproducible. RNA-seq library preparation is inherently error prone due to random template switching during reverse transcription of RNA to cDNA. By counting chimeric events observed when constructing an RNA-seq library from human RNA and spiked-in RNA from an unrelated species, such as the fruit fly, we estimated that ∼1% of RNA-seq reads are artifactually chimeric. In SARS-CoV-2 RNA-seq, we found that the frequency of HVC events was, in fact, not greater than this background "noise." Finally, we developed a novel experimental approach to enrich SARS-CoV-2 sequences from bulk RNA of infected cells. This method enriched viral sequences but did not enrich HVC events, suggesting that the majority of HVC events are, in all likelihood, artifacts of library construction. In conclusion, our findings indicate that HVC events observed in RNA-sequencing libraries from SARS-CoV-2-infected cells are extremely rare and are likely artifacts arising from random template switching of reverse transcriptase and/or sequence alignment errors. Therefore, the observed HVC events do not support SARS-CoV-2 fusion to cellular genes and/or integration into human genomes. IMPORTANCE The pathogenic mechanisms underlying SARS-CoV-2, the virus responsible for COVID-19, are not fully understood. In particular, relatively little is known about the reasons some individuals develop life-threatening or persistent COVID-19. Recent studies identified host-virus chimeric (HVC) reads in RNA-sequencing data from SARS-CoV-2-infected cells and suggested that HVC events support potential "human genome invasion" and "integration" by SARS-CoV-2. This suggestion has fueled concerns about the long-term effects of current mRNA vaccines that incorporate elements of the viral genome. SARS-CoV-2 is a positive-sense, single-stranded RNA virus that does not encode a reverse transcriptase and does not include a nuclear phase in its life cycle, so some doubts have rightfully been expressed regarding the authenticity of HVCs and the role played by endogenous retrotransposons in this phenomenon. Thus, it is important to independently authenticate these HVC events. Here, we provide several lines of evidence suggesting that the observed HVC events are likely artifactual.

Epstein-Barr Virus Episome Physically Interacts with Active Regions of the Host Genome in Lymphoblastoid Cells.

The Epstein-Barr virus (EBV) episome is known to interact with the three-dimensional structure of the human genome in infected cells. However, the exact locations of these interactions and their potential functional consequences remain unclear. Recently, high-resolution chromatin conformation capture (Hi-C) assays in lymphoblastoid cells have become available, enabling us to precisely map the contacts between the EBV episome(s) and the human host genome. Using available Hi-C data at a 10-kb resolution, we have identified 15,000 reproducible contacts between EBV episome(s) and the human genome. These contacts are highly enriched in chromatin regions denoted by typical or super enhancers and active markers, including histone H3K27ac and H3K4me1. Additionally, these contacts are highly enriched at loci bound by host transcription factors that regulate B cell growth (e.g., IKZF1 and RUNX3), factors that enhance cell proliferation (e.g., HDGF), or factors that promote viral replication (e.g., NBS1 and NFIC). EBV contacts show nearly 2-fold enrichment in host regions bound by EBV nuclear antigen 2 (EBNA2) and EBNA3 transcription factors. Circular chromosome conformation capture followed by sequencing (4C-seq) using the EBV origin of plasmid replication (oriP) as a "bait" in lymphoblastoid cells further confirmed contacts with active chromatin regions. Collectively, our analysis supports interactions between EBV episome(s) and active regions of the human genome in lymphoblastoid cells.IMPORTANCE EBV is associated with ∼200,000 cancers each year. In vitro, EBV can transform primary human B lymphocytes into immortalized cell lines. EBV-encoded proteins, along with noncoding RNAs and microRNAs, hijack cellular proteins and pathways to control cell growth. EBV nuclear proteins usurp normal transcriptional programs to activate the expression of key oncogenes, including MYC, to provide a proliferation signal. EBV nuclear antigens also repress CDKN2A to suppress senescence. EBV membrane protein activates NF-κB to provide survival signals. EBV genomes are maintained by EBNA1, which tethers EBV episomes to the host chromosomes during mitosis. However, little is known about where EBV episomes are located in interphase cells. In interphase cells, EBV promoters drive the expression of latency genes, while oriP functions as an enhancer for these promoters. In this study, integrative analyses of published lymphoblastoid cell line (LCL) Hi-C data and our 4C-seq experiments position EBV episomes to host genomes with active epigenetic marks. These contact points were significantly enriched for super enhancers. The close proximity of EBV episomes and the super enhancers that are enriched for transcription cofactors or mediators in lymphoblasts may benefit EBV gene expression, suggesting a novel mechanism of transcriptional activation.

Opposing roles of STAT1 and STAT3 in IL-21 function in CD4+ T cells.

IL-21 is a type I cytokine essential for immune cell differentiation and function. Although IL-21 can activate several STAT family transcription factors, previous studies focused mainly on the role of STAT3 in IL-21 signaling. Here, we investigated the role of STAT1 and show that STAT1 and STAT3 have at least partially opposing roles in IL-21 signaling in CD4(+) T cells. IL-21 induced STAT1 phosphorylation, and this was augmented in Stat3-deficient CD4(+) T cells. RNA-Seq analysis of CD4(+) T cells from Stat1- and Stat3-deficient mice revealed that both STAT1 and STAT3 are critical for IL-21-mediated gene regulation. Expression of some genes, including Tbx21 and Ifng, was differentially regulated by STAT1 and STAT3. Moreover, opposing actions of STAT1 and STAT3 on IFN-γ expression in CD4(+) T cells were demonstrated in vivo during chronic lymphocytic choriomeningitis infection. Finally, IL-21-mediated induction of STAT1 phosphorylation, as well as IFNG and TBX21 expression, were higher in CD4(+) T cells from patients with autosomal dominant hyper-IgE syndrome, which is caused by STAT3 deficiency, as well as in cells from STAT1 gain-of-function patients. These data indicate an interplay between STAT1 and STAT3 in fine-tuning IL-21 actions.

Integrating motif, DNA accessibility and gene expression data to build regulatory maps in an organism.

Characterization of cell type specific regulatory networks and elements is a major challenge in genomics, and emerging strategies frequently employ high-throughput genome-wide assays of transcription factor (TF) to DNA binding, histone modifications or chromatin state. However, these experiments remain too difficult/expensive for many laboratories to apply comprehensively to their system of interest. Here, we explore the potential of elucidating regulatory systems in varied cell types using computational techniques that rely on only data of gene expression, low-resolution chromatin accessibility, and TF-DNA binding specificities ('motifs'). We show that static computational motif scans overlaid with chromatin accessibility data reasonably approximate experimentally measured TF-DNA binding. We demonstrate that predicted binding profiles and expression patterns of hundreds of TFs are sufficient to identify major regulators of ∼200 spatiotemporal expression domains in the Drosophila embryo. We are then able to learn reliable statistical models of enhancer activity for over 70 expression domains and apply those models to annotate domain specific enhancers genome-wide. Throughout this work, we apply our motif and accessibility based approach to comprehensively characterize the regulatory network of fruitfly embryonic development and show that the accuracy of our computational method compares favorably to approaches that rely on data from many experimental assays.

Global analysis of Drosophila Cys2-His2 zinc finger proteins reveals a multitude of novel recognition motifs and binding determinants.

Cys2-His2 zinc finger proteins (ZFPs) are the largest group of transcription factors in higher metazoans. A complete characterization of these ZFPs and their associated target sequences is pivotal to fully annotate transcriptional regulatory networks in metazoan genomes. As a first step in this process, we have characterized the DNA-binding specificities of 129 zinc finger sets from Drosophila using a bacterial one-hybrid system. This data set contains the DNA-binding specificities for at least one encoded ZFP from 70 unique genes and 23 alternate splice isoforms representing the largest set of characterized ZFPs from any organism described to date. These recognition motifs can be used to predict genomic binding sites for these factors within the fruit fly genome. Subsets of fingers from these ZFPs were characterized to define their orientation and register on their recognition sequences, thereby allowing us to define the recognition diversity within this finger set. We find that the characterized fingers can specify 47 of the 64 possible DNA triplets. To confirm the utility of our finger recognition models, we employed subsets of Drosophila fingers in combination with an existing archive of artificial zinc finger modules to create ZFPs with novel DNA-binding specificity. These hybrids of natural and artificial fingers can be used to create functional zinc finger nucleases for editing vertebrate genomes.

Possible Human Papillomavirus 38 Contamination of Endometrial Cancer RNA Sequencing Samples in The Cancer Genome Atlas Database.

UNLABELLED: Viruses are causally associated with a number of human malignancies. In this study, we sought to identify new virus-cancer associations by searching RNA sequencing data sets from >2,000 patients, encompassing 21 cancers from The Cancer Genome Atlas (TCGA), for the presence of viral sequences. In agreement with previous studies, we found human papillomavirus 16 (HPV16) and HPV18 in oropharyngeal cancer and hepatitis B and C viruses in liver cancer. Unexpectedly, however, we found HPV38, a cutaneous form of HPV associated with skin cancer, in 32 of 168 samples from endometrial cancer. In 12 of the HPV38-positive (HPV38(+)) samples, we observed at least one paired read that mapped to both human and HPV38 genomes, indicative of viral integration into the host DNA, something not previously demonstrated for HPV38. The expression levels of HPV38 transcripts were relatively low, and all 32 HPV38(+) samples belonged to the same experimental batch of 40 samples, whereas none of the other 128 endometrial carcinoma samples were HPV38(+), raising doubts about the significance of the HPV38 association. Moreover, the HPV38(+) samples contained the same 10 novel single nucleotide variations (SNVs), leading us to hypothesize that one patient was infected with this new isolate of HPV38, which was integrated into his/her genome and may have cross-contaminated other TCGA samples within batch 228. Based on our analysis, we propose guidelines to examine the batch effect, virus expression level, and SNVs as part of next-generation sequencing (NGS) data analysis for evaluating the significance of viral/pathogen sequences in clinical samples. IMPORTANCE: High-throughput RNA sequencing (RNA-Seq), followed by computational analysis, has vastly accelerated the identification of viral and other pathogenic sequences in clinical samples, but cross-contamination during the processing of the samples remain a major problem that can lead to erroneous conclusions. We found HPV38 sequences specifically present in RNA-Seq samples from endometrial cancer patients from TCGA, a virus not previously associated with this type of cancer. However, multiple lines of evidence suggest possible cross-contamination in these samples, which were processed together in the same batch. Despite this potential cross-contamination, our data indicate that we have detected a new isolate of HPV38 that appears to be integrated into the human genome. We also provide general guidelines for computational detection and interpretation of pathogen-disease associations.