Single-Cell RNA Sequencing of the T Helper Cell Response to House Dust Mites Defines a Distinct Gene Expression Signature in Airway Th2 Cells.

Naive CD4+ T cells differentiate into functionally diverse T helper (Th) cell subsets. Th2 cells play a pathogenic role in asthma, yet a clear picture of their transcriptional profile is lacking. We performed single-cell RNA sequencing (scRNA-seq) of T helper cells from lymph node, lung, and airways in the house dust mite (HDM) model of allergic airway disease. scRNA-seq resolved transcriptional profiles of naive CD4+ T, Th1, Th2, regulatory T (Treg) cells, and a CD4+ T cell population responsive to type I interferons. Th2 cells in the airways were enriched for transcription of many genes, including Cd200r1, Il6, Plac8, and Igfbp7, and their mRNA profile was supported by analysis of chromatin accessibility and flow cytometry. Pathways associated with lipid metabolism were enriched in Th2 cells, and experiments with inhibitors of key metabolic pathways supported roles for glucose and lipid metabolism. These findings provide insight into the differentiation of pathogenic Th2 cells in the context of allergy.

Stellate Cells, Hepatocytes, and Endothelial Cells Imprint the Kupffer Cell Identity on Monocytes Colonizing the Liver Macrophage Niche.

Macrophages are strongly adapted to their tissue of residence. Yet, little is known about the cell-cell interactions that imprint the tissue-specific identities of macrophages in their respective niches. Using conditional depletion of liver Kupffer cells, we traced the developmental stages of monocytes differentiating into Kupffer cells and mapped the cellular interactions imprinting the Kupffer cell identity. Kupffer cell loss induced tumor necrosis factor (TNF)- and interleukin-1 (IL-1) receptor-dependent activation of stellate cells and endothelial cells, resulting in the transient production of chemokines and adhesion molecules orchestrating monocyte engraftment. Engrafted circulating monocytes transmigrated into the perisinusoidal space and acquired the liver-associated transcription factors inhibitor of DNA 3 (ID3) and liver X receptor-α (LXR-α). Coordinated interactions with hepatocytes induced ID3 expression, whereas endothelial cells and stellate cells induced LXR-α via a synergistic NOTCH-BMP pathway. This study shows that the Kupffer cell niche is composed of stellate cells, hepatocytes, and endothelial cells that together imprint the liver-specific macrophage identity.

The Transcription Factor ZEB2 Is Required to Maintain the Tissue-Specific Identities of Macrophages.

Heterogeneity between different macrophage populations has become a defining feature of this lineage. However, the conserved factors defining macrophages remain largely unknown. The transcription factor ZEB2 is best described for its role in epithelial to mesenchymal transition; however, its role within the immune system is only now being elucidated. We show here that Zeb2 expression is a conserved feature of macrophages. Using Clec4f-cre, Itgax-cre, and Fcgr1-cre mice to target five different macrophage populations, we found that loss of ZEB2 resulted in macrophage disappearance from the tissues, coupled with their subsequent replenishment from bone-marrow precursors in open niches. Mechanistically, we found that ZEB2 functioned to maintain the tissue-specific identities of macrophages. In Kupffer cells, ZEB2 achieved this by regulating expression of the transcription factor LXRα, removal of which recapitulated the loss of Kupffer cell identity and disappearance. Thus, ZEB2 expression is required in macrophages to preserve their tissue-specific identities.

ConTra v3: a tool to identify transcription factor binding sites across species, update 2017.

Transcription factors are important gene regulators with distinctive roles in development, cell signaling and cell cycling, and they have been associated with many diseases. The ConTra v3 web server allows easy visualization and exploration of predicted transcription factor binding sites (TFBSs) in any genomic region surrounding coding or non-coding genes. In this updated version, with a completely re-implemented user interface using latest web technologies, users can choose from nine reference organisms ranging from human to yeast. ConTra v3 can analyze promoter regions, 5΄-UTRs, 3΄-UTRs and introns or any other genomic region of interest. Thousands of position weight matrices are available to choose from for detecting specific binding sites. Besides this visualization option, additional new exploration functionality is added to the tool that will automatically detect TFBSs having at the same time the highest regulatory potential, the highest conservation scores of the genomic regions covered by the predicted TFBSs and strongest co-localizations with genomic regions exhibiting regulatory activity. The ConTra v3 web server is freely available at http://bioit2.irc.ugent.be/contra/v3.

GATA1-Deficient Dendritic Cells Display Impaired CCL21-Dependent Migration toward Lymph Nodes Due to Reduced Levels of Polysialic Acid.

Dendritic cells (DCs) play a pivotal role in the regulation of the immune response. DC development and activation is finely orchestrated through transcriptional programs. GATA1 transcription factor is required for murine DC development, and data suggest that it might be involved in the fine-tuning of the life span and function of activated DCs. We generated DC-specific Gata1 knockout mice (Gata1-KODC), which presented a 20% reduction of splenic DCs, partially explained by enhanced apoptosis. RNA sequencing analysis revealed a number of deregulated genes involved in cell survival, migration, and function. DC migration toward peripheral lymph nodes was impaired in Gata1-KODC mice. Migration assays performed in vitro showed that this defect was selective for CCL21, but not CCL19. Interestingly, we show that Gata1-KODC DCs have reduced polysialic acid levels on their surface, which is a known determinant for the proper migration of DCs toward CCL21.

Computationally designed liver-specific transcriptional modules and hyperactive factor IX improve hepatic gene therapy.

The development of the next-generation gene therapy vectors for hemophilia requires using lower and thus potentially safer vector doses and augmenting their therapeutic efficacy. We have identified hepatocyte-specific transcriptional cis-regulatory modules (CRMs) by using a computational strategy that increased factor IX (FIX) levels 11- to 15-fold. Vector efficacy could be enhanced by combining these hepatocyte-specific CRMs with a synthetic codon-optimized hyperfunctional FIX-R338L Padua transgene. This Padua mutation boosted FIX activity up to sevenfold, with no apparent increase in thrombotic risk. We then validated this combination approach using self-complementary adenoassociated virus serotype 9 (scAAV9) vectors in hemophilia B mice. This resulted in sustained supraphysiologic FIX activity (400%), correction of the bleeding diathesis at clinically relevant, low vector doses (5 × 10(10) vector genomes [vg]/kg) that are considered safe in patients undergoing gene therapy. Moreover, immune tolerance could be induced that precluded induction of inhibitory antibodies to FIX upon immunization with recombinant FIX protein.

Genome-wide computational analysis reveals cardiomyocyte-specific transcriptional Cis-regulatory motifs that enable efficient cardiac gene therapy.

Gene therapy is a promising emerging therapeutic modality for the treatment of cardiovascular diseases and hereditary diseases that afflict the heart. Hence, there is a need to develop robust cardiac-specific expression modules that allow for stable expression of the gene of interest in cardiomyocytes. We therefore explored a new approach based on a genome-wide bioinformatics strategy that revealed novel cardiac-specific cis-acting regulatory modules (CS-CRMs). These transcriptional modules contained evolutionary-conserved clusters of putative transcription factor binding sites that correspond to a "molecular signature" associated with robust gene expression in the heart. We then validated these CS-CRMs in vivo using an adeno-associated viral vector serotype 9 that drives a reporter gene from a quintessential cardiac-specific α-myosin heavy chain promoter. Most de novo designed CS-CRMs resulted in a >10-fold increase in cardiac gene expression. The most robust CRMs enhanced cardiac-specific transcription 70- to 100-fold. Expression was sustained and restricted to cardiomyocytes. We then combined the most potent CS-CRM4 with a synthetic heart and muscle-specific promoter (SPc5-12) and obtained a significant 20-fold increase in cardiac gene expression compared to the cytomegalovirus promoter. This study underscores the potential of rational vector design to improve the robustness of cardiac gene therapy.

Liver-specific transcriptional modules identified by genome-wide in silico analysis enable efficient gene therapy in mice and non-human primates.

The robustness and safety of liver-directed gene therapy can be substantially improved by enhancing expression of the therapeutic transgene in the liver. To achieve this, we developed a new approach of rational in silico vector design. This approach relies on a genome-wide bio-informatics strategy to identify cis-acting regulatory modules (CRMs) containing evolutionary conserved clusters of transcription factor binding site motifs that determine high tissue-specific gene expression. Incorporation of these CRMs into adeno-associated viral (AAV) and non-viral vectors enhanced gene expression in mice liver 10 to 100-fold, depending on the promoter used. Furthermore, these CRMs resulted in robust and sustained liver-specific expression of coagulation factor IX (FIX), validating their immediate therapeutic and translational relevance. Subsequent translational studies indicated that therapeutic FIX expression levels could be attained reaching 20-35% of normal levels after AAV-based liver-directed gene therapy in cynomolgus macaques. This study underscores the potential of rational vector design using computational approaches to improve their robustness and therefore allows for the use of lower and thus safer vector doses for gene therapy, while maximizing therapeutic efficacy.

PhysBinder: Improving the prediction of transcription factor binding sites by flexible inclusion of biophysical properties.

The most important mechanism in the regulation of transcription is the binding of a transcription factor (TF) to a DNA sequence called the TF binding site (TFBS). Most binding sites are short and degenerate, which makes predictions based on their primary sequence alone somewhat unreliable. We present a new web tool that implements a flexible and extensible algorithm for predicting TFBS. The algorithm makes use of both direct (the sequence) and several indirect readout features of protein-DNA complexes (biophysical properties such as bendability or the solvent-excluded surface of the DNA). This algorithm significantly outperforms state-of-the-art approaches for in silico identification of TFBS. Users can submit FASTA sequences for analysis in the PhysBinder integrative algorithm and choose from >60 different TF-binding models. The results of this analysis can be used to plan and steer wet-lab experiments. The PhysBinder web tool is freely available at http://bioit.dmbr.ugent.be/physbinder/index.php.

A flexible integrative approach based on random forest improves prediction of transcription factor binding sites.

Transcription factor binding sites (TFBSs) are DNA sequences of 6-15 base pairs. Interaction of these TFBSs with transcription factors (TFs) is largely responsible for most spatiotemporal gene expression patterns. Here, we evaluate to what extent sequence-based prediction of TFBSs can be improved by taking into account the positional dependencies of nucleotides (NPDs) and the nucleotide sequence-dependent structure of DNA. We make use of the random forest algorithm to flexibly exploit both types of information. Results in this study show that both the structural method and the NPD method can be valuable for the prediction of TFBSs. Moreover, their predictive values seem to be complementary, even to the widely used position weight matrix (PWM) method. This led us to combine all three methods. Results obtained for five eukaryotic TFs with different DNA-binding domains show that our method improves classification accuracy for all five eukaryotic TFs compared with other approaches. Additionally, we contrast the results of seven smaller prokaryotic sets with high-quality data and show that with the use of high-quality data we can significantly improve prediction performance. Models developed in this study can be of great use for gaining insight into the mechanisms of TF binding.

ConTra v2: a tool to identify transcription factor binding sites across species, update 2011.

Transcription factors are important gene regulators with distinctive roles in development, cell signaling and cell cycling, and they have been associated with many diseases. The ConTra v2 web server allows easy visualization and exploration of predicted transcription factor binding sites in any genomic region surrounding coding or non-coding genes. In this new version, users can choose from nine reference organisms ranging from human to yeast. ConTra v2 can analyze promoter regions, 5'-UTRs, 3'-UTRs and introns or any other genomic region of interest. Hundreds of position weight matrices are available to choose from, but the user can also upload any other matrices for detecting specific binding sites. A typical analysis is run in four simple steps of choosing the gene, the transcript, the region of interest and then selecting one or more transcription factor binding sites. The ConTra v2 web server is freely available at http://bioit.dmbr.ugent.be/contrav2/index.php.

Nanobodies® specific for respiratory syncytial virus fusion protein protect against infection by inhibition of fusion.

Despite the medical importance of respiratory syncytial virus (RSV) infections, there is no vaccine or therapeutic agent available. Prophylactic administration of palivizumab, a humanized monoclonal RSV fusion (F) protein-specific antibody, can protect high-risk children. Previously, we have demonstrated that RSV can be neutralized by picomolar concentrations of a camelid immunoglobulin single-variable domain that binds the RSV protein F (F-VHHb nanobodies). Here, we investigated the mechanism by which these nanobodies neutralize RSV and tested their antiviral activity in vivo. We demonstrate that bivalent RSV F-specific nanobodies neutralize RSV infection by inhibiting fusion without affecting viral attachment. The ability of RSV F-specific nanobodies to protect against RSV infection was investigated in vivo. Intranasal administration of bivalent RSV F-specific nanobodies protected BALB/c mice from RSV infection, and associated pulmonary inflammation. Moreover, therapeutic treatment with these nanobodies after RSV infection could reduce viral replication and reduced pulmonary inflammation. Thus, nanobodies are promising therapeutic molecules for treatment of RSV.

Novel expression cassettes for increasing apolipoprotein AI transgene expression in vascular endothelial cells.

Transduction of endothelial cells (EC) with a vector that expresses apolipoprotein A-I (APOAI) reduces atherosclerosis in arteries of fat-fed rabbits. However, the effects on atherosclerosis are partial and might be enhanced if APOAI expression could be increased. With a goal of developing an expression cassette that generates higher levels of APOAI mRNA in EC, we tested 4 strategies, largely in vitro: addition of 2 types of enhancers, addition of computationally identified EC-specific cis-regulatory modules (CRM), and insertion of the rabbit APOAI gene at the transcription start site (TSS) of sequences cloned from genes that are highly expressed in cultured EC. Addition of a shear stress-responsive enhancer did not increase APOAI expression. Addition of 2 copies of a Mef2c enhancer increased APOAI expression from a moderately active promoter/enhancer but decreased APOAI expression from a highly active promoter/enhancer. Of the 11 CRMs, 3 increased APOAI expression from a moderately active promoter (2-7-fold; P < 0.05); none increased expression from a highly active promoter/enhancer. Insertion of the APOAI gene into the TSS of highly expressed EC genes did not increase expression above levels obtained with a moderately active promoter/enhancer. New strategies are needed to further increase APOAI transgene expression in EC.

CARD14 Signalling Ensures Cell Survival and Cancer Associated Gene Expression in Prostate Cancer Cells.

Prostate cancer (PCa) is one of the most common cancer types in men and represents an increasing global problem due to the modern Western lifestyle. The signalling adapter protein CARD14 is specifically expressed in epithelial cells, where it has been shown to mediate NF-κB signalling, but a role for CARD14 in carcinoma has not yet been described. By analysing existing cancer databases, we found that CARD14 overexpression strongly correlates with aggressive PCa in human patients. Moreover, we showed that CARD14 is overexpressed in the LNCaP PCa cell line and that knockdown of CARD14 severely reduces LNCaP cell survival. Similarly, knockdown of BCL10 and MALT1, which are known to form a signalling complex with CARD14, also induced LNCaP cell death. MALT1 is a paracaspase that mediates downstream signalling by acting as a scaffold, as well as a protease. Recent studies have already indicated a role for the scaffold function of MALT1 in PCa cell growth. Here, we also demonstrated constitutive MALT1 proteolytic activity in several PCa cell lines, leading to cleavage of A20 and CYLD. Inhibition of MALT1 protease activity did not affect PCa cell survival nor activation of NF-κB and JNK signalling, but reduced expression of cancer-associated genes, including the cytokine IL-6. Taken together, our results revealed a novel role for CARD14-induced signalling in regulating PCa cell survival and gene expression. The epithelial cell type-specific expression of CARD14 may offer novel opportunities for more specific therapeutic targeting approaches in PCa.

GTF3A mutations predispose to herpes simplex encephalitis by disrupting biogenesis of the host-derived RIG-I ligand RNA5SP141.

Herpes simplex virus 1 (HSV-1) infects several billion people worldwide and can cause life-threatening herpes simplex encephalitis (HSE) in some patients. Monogenic defects in components of the type I interferon system have been identified in patients with HSE, emphasizing the role of inborn errors of immunity underlying HSE pathogenesis. Here, we identify compound heterozygous loss-of-function mutations in the gene GTF3A encoding for transcription factor IIIA (TFIIIA), a component of the RNA polymerase III complex, in a patient with common variable immunodeficiency and HSE. Patient fibroblasts and GTF3A gene-edited cells displayed impaired HSV-1-induced innate immune responses and enhanced HSV-1 replication. Chromatin immunoprecipitation sequencing analysis identified the 5S ribosomal RNA pseudogene 141 (RNA5SP141), an endogenous ligand of the RNA sensor RIG-I, as a transcriptional target of TFIIIA. GTF3A mutant cells exhibited diminished RNA5SP141 expression and abrogated RIG-I activation upon HSV-1 infection. Our work unveils a crucial role for TFIIIA in transcriptional regulation of a cellular RIG-I agonist and shows that GTF3A genetic defects lead to impaired cell-intrinsic anti-HSV-1 responses and can predispose to HSE.

ConTra: a promoter alignment analysis tool for identification of transcription factor binding sites across species.

Transcription factors (TFs) are key components in signaling pathways, and the presence of their binding sites in the promoter regions of DNA is essential for their regulation of the expression of the corresponding genes. Orthologous promoter sequences are commonly used to increase the specificity with which potentially functional transcription factor binding sites (TFBSs) are recognized and to detect possibly important similarities or differences between the different species. The ConTra (conserved TFBSs) web server provides the biologist at the bench with a user-friendly tool to interactively visualize TFBSs predicted using either TransFac (1) or JASPAR (2) position weight matrix libraries, on a promoter alignment of choice. The visualization can be preceded by a simple scoring analysis to explore which TFs are the most likely to bind to the promoter of interest. The ConTra web server is available at http://bioit.dmbr.ugent.be/ConTra/index.php.

ORegAnno: an open-access community-driven resource for regulatory annotation.

ORegAnno is an open-source, open-access database and literature curation system for community-based annotation of experimentally identified DNA regulatory regions, transcription factor binding sites and regulatory variants. The current release comprises 30 145 records curated from 922 publications and describing regulatory sequences for over 3853 genes and 465 transcription factors from 19 species. A new feature called the 'publication queue' allows users to input relevant papers from scientific literature as targets for annotation. The queue contains 4438 gene regulation papers entered by experts and another 54 351 identified by text-mining methods. Users can enter or 'check out' papers from the queue for manual curation using a series of user-friendly annotation pages. A typical record entry consists of species, sequence type, sequence, target gene, binding factor, experimental outcome and one or more lines of experimental evidence. An evidence ontology was developed to describe and categorize these experiments. Records are cross-referenced to Ensembl or Entrez gene identifiers, PubMed and dbSNP and can be visualized in the Ensembl or UCSC genome browsers. All data are freely available through search pages, XML data dumps or web services at: http://www.oreganno.org.

Comparison of the PU.1 transcriptional regulome and interactome in human and mouse inflammatory dendritic cells.

Dendritic cells (DCs) are key immune modulators and are able to mount immune responses or tolerance. DC differentiation and activation imply a plethora of molecular and cellular responses, including transcriptional changes. PU.1 is a highly expressed transcription factor in DCs and coordinates relevant aspects of DC biology. Due to their role as immune regulators, DCs pose as a promising immunotherapy tool. However, some of their functional features, such as survival, activation, or migration, are compromised due to the limitations to simulate in vitro the physiologic DC differentiation process. A better knowledge of transcriptional programs would allow the identification of potential targets for manipulation with the aim of obtaining "qualified" DCs for immunotherapy purposes. Most of the current knowledge regarding DC biology derives from studies using mouse models, which not always find a parallel in human. In the present study, we dissect the PU.1 transcriptional regulome and interactome in mouse and human DCs, in the steady state or LPS activated. The PU.1 transcriptional regulome was identified by performing PU.1 chromatin immunoprecipitation followed by high-throughput sequencing and pairing these data with RNAsequencing data. The PU.1 interactome was identified by performing PU.1 immunoprecipitation followed by mass spectrometry analysis. Our results portray PU.1 as a pivotal factor that plays an important role in the regulation of genes required for proper DC activation and function, and assures the repression of nonlineage genes. The interspecies differences between human and mouse DCs are surprisingly substantial, highlighting the need to study the biology of human DCs.

Human metallothionein expression under normal and pathological conditions: mechanisms of gene regulation based on in silico promoter analysis.

Metallothioneins (MTs) are ubiquitous metal-binding proteins that have been highly conserved throughout evolution. Although their physiological function is not completely understood, they are involved in diverse processes including metal homeostasis and detoxification, the oxidative stress response, inflammation, and cell proliferation. Te human MT gene family consists of at least 18 isoforms, containing pseudogenes as well as genes encoding functional proteins. Most of the MT isoforms can be induced by a wide variety of substances, such as metals, cytokines, and hormones. Different cell types express discrete MT isoforms, which reflects the specifically adapted functions of MTs and a divergence in their regulation. Te aberrant expression of MTs has been described in a number of diseases, including Crohn's disease, cancer, Alzheimer's disease, amyotrophic lateral sclerosis, Menkes disease, and Wilson's disease. Therefore, a thorough understanding of MT gene regulation is imperative. To date, the transcriptional regulation of MTs has primarily been studied in mice. While only four murine MT isoforms exist, the homology between murine and human MTs allows for the evaluation of the regulatory regions in their respective promoters. Here, we review the aberrant expression of MTs in human diseases and the mechanisms that regulate MT1 expression based on an in silico evaluation of transcription factor binding sites.

Low nucleosome occupancy is encoded around functional human transcription factor binding sites.

BACKGROUND: Transcriptional regulation of genes in eukaryotes is achieved by the interactions of multiple transcription factors with arrays of transcription factor binding sites (TFBSs) on DNA and with each other. Identification of these TFBSs is an essential step in our understanding of gene regulatory networks, but computational prediction of TFBSs with either consensus or commonly used stochastic models such as Position-Specific Scoring Matrices (PSSMs) results in an unacceptably high number of hits consisting of a few true functional binding sites and numerous false non-functional binding sites. This is due to the inability of the models to incorporate higher order properties of sequences including sequences surrounding TFBSs and influencing the positioning of nucleosomes and/or the interactions that might occur between transcription factors. RESULTS: Significant improvement can be expected through the development of a new framework for the modeling and prediction of TFBSs that considers explicitly these higher order sequence properties. It would be particularly interesting to include in the new modeling framework the information present in the nucleosome positioning sequences (NPSs) surrounding TFBSs, as it can be hypothesized that genomes use this information to encode the formation of stable nucleosomes over non-functional sites, while functional sites have a more open chromatin configuration. In this report we evaluate the usefulness of the latter feature by comparing the nucleosome occupancy probabilities around experimentally verified human TFBSs with the nucleosome occupancy probabilities around false positive TFBSs and in random sequences. CONCLUSION: We present evidence that nucleosome occupancy is remarkably lower around true functional human TFBSs as compared to non-functional human TFBSs, which supports the use of this feature to improve current TFBS prediction approaches in higher eukaryotes.