Evaluation of EPISEQ SARS-CoV-2 and a Fully Integrated Application to Identify SARS-CoV-2 Variants from Several Next-Generation Sequencing Approaches.

Whole-genome sequencing has become an essential tool for real-time genomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide. The handling of raw next-generation sequencing (NGS) data is a major challenge for sequencing laboratories. We developed an easy-to-use web-based application (EPISEQ SARS-CoV-2) to analyse SARS-CoV-2 NGS data generated on common sequencing platforms using a variety of commercially available reagents. This application performs in one click a quality check, a reference-based genome assembly, and the analysis of the generated consensus sequence as to coverage of the reference genome, mutation screening and variant identification according to the up-to-date Nextstrain clade and Pango lineage. In this study, we validated the EPISEQ SARS-CoV-2 pipeline against a reference pipeline and compared the performance of NGS data generated by different sequencing protocols using EPISEQ SARS-CoV-2. We showed a strong agreement in SARS-CoV-2 clade and lineage identification (>99%) and in spike mutation detection (>99%) between EPISEQ SARS-CoV-2 and the reference pipeline. The comparison of several sequencing approaches using EPISEQ SARS-CoV-2 revealed 100% concordance in clade and lineage classification. It also uncovered reagent-related sequencing issues with a potential impact on SARS-CoV-2 mutation reporting. Altogether, EPISEQ SARS-CoV-2 allows an easy, rapid and reliable analysis of raw NGS data to support the sequencing efforts of laboratories with limited bioinformatics capacity and those willing to accelerate genomic surveillance of SARS-CoV-2.

TET2 regulates immune tolerance in chronically activated mast cells.

Mutation of the TET2 DNA-hydroxymethylase has been associated with a number of immune pathologies. The disparity in phenotype and clinical presentation among these pathologies leads to questions regarding the role of TET2 mutation in promoting disease evolution in different immune cell types. Here we show that, in primary mast cells, Tet2 expression is induced in response to chronic and acute activation signals. In TET2-deficient mast cells, chronic activation via the oncogenic KITD816V allele associated with mastocytosis, selects for a specific epigenetic signature characterized by hypermethylated DNA regions (HMR) at immune response genes. H3K27ac and transcription factor binding is consistent with priming or more open chromatin at both HMR and non-HMR in proximity to immune genes in these cells, and this signature coincides with increased pathological inflammation signals. HMR are also associated with a subset of immune genes that are direct targets of TET2 and repressed in TET2-deficient cells. Repression of these genes results in immune tolerance to acute stimulation that can be rescued with vitamin C treatment or reiterated with a Tet inhibitor. Overall, our data support a model where TET2 plays a direct role in preventing immune tolerance in chronically activated mast cells, supporting TET2 as a viable target to reprogram the innate immune response for innovative therapies.

High-Risk International Clones of Carbapenem-Nonsusceptible Pseudomonas aeruginosa Endemic to Indonesian Intensive Care Units: Impact of a Multifaceted Infection Control Intervention Analyzed at the Genomic Level.

Infection control effectiveness evaluations require detailed epidemiological and microbiological data. We analyzed the genomic profiles of carbapenem-nonsusceptible Pseudomonas aeruginosa (CNPA) strains collected from two intensive care units (ICUs) in the national referral hospital in Jakarta, Indonesia, where a multifaceted infection control intervention was applied. We used clinical data combined with whole-genome sequencing (WGS) of systematically collected CNPA to infer the transmission dynamics of CNPA strains and to characterize their resistome. We found that the number of CNPA transmissions and acquisitions by patients was highly variable over time but that, overall, the rates were not significantly reduced by the intervention. Environmental sources were involved in these transmissions and acquisitions. Four high-risk international CNPA clones (ST235, ST823, ST357, and ST446) dominated, but the distribution of these clones changed significantly after the intervention was implemented. Using resistome analysis, carbapenem resistance was explained by the presence of various carbapenemase-encoding genes (blaGES-5, blaVIM-2-8, and blaIMP-1-7-43) and by mutations within the porin OprD. Our results reveal for the first time the dynamics of P. aeruginosa antimicrobial resistance (AMR) profiles in Indonesia and additionally show the utility of WGS in combination with clinical data to evaluate the impact of an infection control intervention. (This study has been registered at www.trialregister.nl under registration no. NTR5541).IMPORTANCE In low-to-middle-income countries such as Indonesia, work in intensive care units (ICUs) can be hampered by lack of resources. Conducting large epidemiological studies in such settings using genomic tools is rather challenging. Still, we were able to systematically study the transmissions of carbapenem-nonsusceptible strains of P. aeruginosa (CNPA) within and between ICUs, before and after an infection control intervention. Our data show the importance of the broad dissemination of the internationally recognized CNPA clones, the relevance of environmental reservoirs, and the mixed effects of the implemented intervention; it led to a profound change in the clonal make-up of CNPA, but it did not reduce the patients' risk of CNPA acquisitions. Thus, CNPA epidemiology in Indonesian ICUs is part of a global expansion of multiple CNPA clones that remains difficult to control by infection prevention measures.

Genome-wide characterization of mammalian promoters with distal enhancer functions.

Gene expression in mammals is precisely regulated by the combination of promoters and gene-distal regulatory regions, known as enhancers. Several studies have suggested that some promoters might have enhancer functions. However, the extent of this type of promoters and whether they actually function to regulate the expression of distal genes have remained elusive. Here, by exploiting a high-throughput enhancer reporter assay, we unravel a set of mammalian promoters displaying enhancer activity. These promoters have distinct genomic and epigenomic features and frequently interact with other gene promoters. Extensive CRISPR-Cas9 genomic manipulation demonstrated the involvement of these promoters in the cis regulation of expression of distal genes in their natural loci. Our results have important implications for the understanding of complex gene regulation in normal development and disease.

Dynamic recruitment of Ets1 to both nucleosome-occupied and -depleted enhancer regions mediates a transcriptional program switch during early T-cell differentiation.

Ets1 is a sequence-specific transcription factor that plays an important role during hematopoiesis, and is essential for the transition of CD4(-)/CD8(-) double negative (DN) to CD4(+)/CD8(+) double positive (DP) thymocytes. Using genome-wide and functional approaches, we investigated the binding properties, transcriptional role and chromatin environment of Ets1 during this transition. We found that while Ets1 binding at distal sites was associated with active genes at both DN and DP stages, its enhancer activity was attained at the DP stage, as reflected by levels of the core transcriptional hallmarks H3K4me1/3, RNA Polymerase II and eRNA. This dual, stage-specific ability reflected a switch from non-T hematopoietic toward T-cell specific gene expression programs during the DN-to-DP transition, as indicated by transcriptome analyses of Ets1(-/-) thymic cells. Coincidentally, Ets1 associates more specifically with Runx1 in DN and with TCF1 in DP cells. We also provide evidence that Ets1 predominantly binds distal nucleosome-occupied regions in DN and nucleosome-depleted regions in DP. Finally and importantly, we demonstrate that Ets1 induces chromatin remodeling by displacing H3K4me1-marked nucleosomes. Our results thus provide an original model whereby the ability of a transcription factor to bind nucleosomal DNA changes during differentiation with consequences on its cognate enhancer activity.

High-throughput and quantitative assessment of enhancer activity in mammals by CapStarr-seq.

Cell-type specific regulation of gene expression requires the activation of promoters by distal genomic elements defined as enhancers. The identification and the characterization of enhancers are challenging in mammals due to their genome complexity. Here we develop CapStarr-Seq, a novel high-throughput strategy to quantitatively assess enhancer activity in mammals. This approach couples capture of regions of interest to previously developed Starr-seq technique. Extensive assessment of CapStarr-seq demonstrates accurate quantification of enhancer activity. Furthermore, we find that enhancer strength is associated with binding complexity of tissue-specific transcription factors and super-enhancers, while additive enhancer activity isolates key genes involved in cell identity and function. The CapStarr-Seq thus provides a fast and cost-effective approach to assess the activity of potential enhancers for a given cell type and will be helpful in decrypting transcription regulation mechanisms.

Integrative analysis of public ChIP-seq experiments reveals a complex multi-cell regulatory landscape.

The large collections of ChIP-seq data rapidly accumulating in public data warehouses provide genome-wide binding site maps for hundreds of transcription factors (TFs). However, the extent of the regulatory occupancy space in the human genome has not yet been fully apprehended by integrating public ChIP-seq data sets and combining it with ENCODE TFs map. To enable genome-wide identification of regulatory elements we have collected, analysed and retained 395 available ChIP-seq data sets merged with ENCODE peaks covering a total of 237 TFs. This enhanced repertoire complements and refines current genome-wide occupancy maps by increasing the human genome regulatory search space by 14% compared to ENCODE alone, and also increases the complexity of the regulatory dictionary. As a direct application we used this unified binding repertoire to annotate variant enhancer loci (VELs) from H3K4me1 mark in two cancer cell lines (MCF-7, CRC) and observed enrichments of specific TFs involved in biological key functions to cancer development and proliferation. Those enrichments of TFs within VELs provide a direct annotation of non-coding regions detected in cancer genomes. Finally, full access to this catalogue is available online together with the TFs enrichment analysis tool (http://tagc.univ-mrs.fr/remap/).

Divergent transcription is associated with promoters of transcriptional regulators.

BACKGROUND: Divergent transcription is a wide-spread phenomenon in mammals. For instance, short bidirectional transcripts are a hallmark of active promoters, while longer transcripts can be detected antisense from active genes in conditions where the RNA degradation machinery is inhibited. Moreover, many described long non-coding RNAs (lncRNAs) are transcribed antisense from coding gene promoters. However, the general significance of divergent lncRNA/mRNA gene pair transcription is still poorly understood. Here, we used strand-specific RNA-seq with high sequencing depth to thoroughly identify antisense transcripts from coding gene promoters in primary mouse tissues. RESULTS: We found that a substantial fraction of coding-gene promoters sustain divergent transcription of long non-coding RNA (lncRNA)/mRNA gene pairs. Strikingly, upstream antisense transcription is significantly associated with genes related to transcriptional regulation and development. Their promoters share several characteristics with those of transcriptional developmental genes, including very large CpG islands, high degree of conservation and epigenetic regulation in ES cells. In-depth analysis revealed a unique GC skew profile at these promoter regions, while the associated coding genes were found to have large first exons, two genomic features that might enforce bidirectional transcription. Finally, genes associated with antisense transcription harbor specific H3K79me2 epigenetic marking and RNA polymerase II enrichment profiles linked to an intensified rate of early transcriptional elongation. CONCLUSIONS: We concluded that promoters of a class of transcription regulators are characterized by a specialized transcriptional control mechanism, which is directly coupled to relaxed bidirectional transcription.

Active STAT5 regulates T-bet and eomesodermin expression in CD8 T cells and imprints a T-bet-dependent Tc1 program with repressed IL-6/TGF-ß1 signaling.

In adoptive therapy, CD8 T cells expressing active STAT5 (STAT5CA) transcription factors were found to be superior to unmanipulated counterparts in long-term persistence, capacity to infiltrate autochthonous mouse melanomas, thrive in their microenvironment, and induce their regression. However, the molecular mechanisms sustaining these properties were undefined. In this study, we report that STAT5CA induced sustained expression of genes controlling tissue homing, cytolytic granule composition, type 1 CD8 cytotoxic T cell-associated effector molecules granzyme B(+), IFN-γ(+), TNF-α(+), and CCL3(+), but not IL-2, and transcription factors T-bet and eomesodermin (Eomes). Chromatin immunoprecipitation sequencing analyses identified the genes possessing regulatory regions to which STAT5 bound in long-term in vivo maintained STAT5CA-expressing CD8 T cells. This analysis identified 34% of the genes differentially expressed between STAT5CA-expressing and nonexpressing effector T cells as direct STAT5CA target genes, including those encoding T-bet, Eomes, and granzyme B. Additionally, genes encoding the IL-6R and TGFbRII subunits were stably repressed, resulting in dampened IL-17-producing CD8 T cell polarization in response to IL-6 and TGF-ß1. The absence of T-bet did not affect STAT5CA-driven accumulation of the T cells in tissue or their granzyme B expression but restored IL-2 secretion and IL-6R and TGFbRII expression and signaling, as illustrated by IL-17 induction. Therefore, concerted STAT5/T-bet/Eomes regulation controls homing, long-term maintenance, recall responses, and resistance to polarization towards IL-17-producing CD8 T cells while maintaining expression of an efficient type 1 CD8 cytotoxic T cell program (granzyme B(+), IFN-γ(+)).

A comprehensive analysis of gene expression changes provoked by bacterial and fungal infection in C. elegans.

While Caenorhabditis elegans specifically responds to infection by the up-regulation of certain genes, distinct pathogens trigger the expression of a common set of genes. We applied new methods to conduct a comprehensive and comparative study of the transcriptional response of C. elegans to bacterial and fungal infection. Using tiling arrays and/or RNA-sequencing, we have characterized the genome-wide transcriptional changes that underlie the host's response to infection by three bacterial (Serratia marcescens, Enterococcus faecalis and otorhabdus luminescens) and two fungal pathogens (Drechmeria coniospora and Harposporium sp.). We developed a flexible tool, the WormBase Converter (available at http://wormbasemanager.sourceforge.net/), to allow cross-study comparisons. The new data sets provided more extensive lists of differentially regulated genes than previous studies. Annotation analysis confirmed that genes commonly up-regulated by bacterial infections are related to stress responses. We found substantial overlaps between the genes regulated upon intestinal infection by the bacterial pathogens and Harposporium, and between those regulated by Harposporium and D. coniospora, which infects the epidermis. Among the fungus-regulated genes, there was a significant bias towards genes that are evolving rapidly and potentially encode small proteins. The results obtained using new methods reveal that the response to infection in C. elegans is determined by the nature of the pathogen, the site of infection and the physiological imbalance provoked by infection. They form the basis for future functional dissection of innate immune signaling. Finally, we also propose alternative methods to identify differentially regulated genes that take into account the greater variability in lowly expressed genes.