HERVs characterize normal and leukemia stem cells and represent a source of shared epitopes for cancer immunotherapy.

Human endogenous retroviruses (HERVs) represent 8% of the human genome. The expression of HERVs and their immune impact have not been extensively studied in Acute Myeloid Leukemia (AML). In this study, we used a reference of 14 968 HERV functional units to provide a thorough analysis of HERV expression in normal and AML bone marrow cells. We show that the HERV retrotranscriptome accurately characterizes normal and leukemic cell subpopulations, including leukemia stem cells, in line with different epigenetic profiles. We then show that HERV expression delineates AML subtypes with different prognoses. We finally propose a method to select and prioritize CD8+ T cell epitopes derived from AML-specific HERVs and we show that lymphocytes infiltrating patient bone marrow at diagnosis contain naturally occurring CD8+ T cells against these HERV epitopes. We also provide in vitro data supporting the functionality of HERV-specific CD8+ T-cells against AML cells. These results show that HERVs represent an important source of genetic information that can help enhancing disease stratification or biomarker identification and an important reservoir of alternative tumor-specific T cell epitopes relevant for cancer immunotherapy.

SARS-CoV-2 historical global testing and genomic variability.

The world has been dealing with the SARS-COV-2 pandemic since December 2019 and a lot of effort has focused on tracking the spread of the virus by gathering information regarding testing statistics and generating viral genomic sequences. Unfortunately, there is neither a single comprehensive resource with global historical testing data nor a centralized database with summary statistics of the identified genomic variants. We merged different pre-aggregated historical testing data and complemented them with our manually extracted ones, which consist of 6852 historical test statistics from 76 countries/states unreported in any other dataset, at the date of submission, making our dataset the most comprehensive to date. We also analyzed all publicly deposited SARS-CoV-2 genomic sequences in GISAID and annotated their variants. Both datasets can be accessed through our interactive dashboard which also provides important insights on different outbreak trends across countries and states. The dashboard is available at https://bioinfo.lau.edu.lb/gkhazen/covid19 . A daily updated version of the datasets can be downloaded from github.com/KhazenLab/covid19-data.

Multi-scale structural community organisation of the human genome.

BACKGROUND: Structural interaction frequency matrices between all genome loci are now experimentally achievable thanks to high-throughput chromosome conformation capture technologies. This ensues a new methodological challenge for computational biology which consists in objectively extracting from these data the structural motifs characteristic of genome organisation. RESULTS: We deployed the fast multi-scale community mining algorithm based on spectral graph wavelets to characterise the networks of intra-chromosomal interactions in human cell lines. We observed that there exist structural domains of all sizes up to chromosome length and demonstrated that the set of structural communities forms a hierarchy of chromosome segments. Hence, at all scales, chromosome folding predominantly involves interactions between neighbouring sites rather than the formation of links between distant loci. CONCLUSIONS: Multi-scale structural decomposition of human chromosomes provides an original framework to question structural organisation and its relationship to functional regulation across the scales. By construction the proposed methodology is independent of the precise assembly of the reference genome and is thus directly applicable to genomes whose assembly is not fully determined.

Tumour burden and antigen-specific T cell magnitude represent major parameters for clinical response to cancer vaccine and TCR-engineered T cell therapy.

BACKGROUND: Cancer vaccines and T-cell receptor (TCR) engineered T cells (Tg-T cell) represent two different therapeutic strategies that can target the same tumour epitopes. The first approach requires the induction of a specific immune response in patients, while the second relies on the efficacy of adoptively transferred T cells. Because the ratio of antigen-specific T cells to tumour cells engaged by these strategies may influence the clinical outcome, we evaluated the efficacy of these two therapeutic approaches in solid tumours according to the tumour burden. METHODS: We performed a meta-analysis restricted to the therapeutic vaccine and Tg-T cell trials, presenting annotated individual clinical data. We adapted a previously published mathematical model for tumour immune dynamics to estimate the clinical impact of the number of specific T cells in regard to the tumour burden. RESULTS: A focused analysis of Tg-T cell studies revealed that clinical responses were mostly observed with the highest doses of infused T cells, suggesting that exceeding a threshold of effector T cells may be required for clinical efficacy. In silico modelling of cancer vaccine and Tg-T cell therapies starting at different tumour burdens showed that therapeutic vaccines control low or moderate tumour burdens, whereas increasing the amount of infused Tg-T cells succeeds in controlling high tumour masses. CONCLUSION: We propose that therapeutic vaccines should be considered in the context of low or moderate tumour burden, whereas Tg-T cell strategies may be more adapted for the treatment of advanced metastatic diseases.

Identification of shared tumor epitopes from endogenous retroviruses inducing high-avidity cytotoxic T cells for cancer immunotherapy.

Human endogenous retroviruses (HERVs) represent 8% of the human genome. HERV products may represent tumor antigens relevant for cancer immunotherapy. We developed a bioinformatic approach to identify shared CD8+ T cell epitopes derived from cancer-associated HERVs in solid tumors. Six candidates among the most commonly shared HLA-A2 epitopes with evidence of translation were selected for immunological evaluation. In vitro priming assays confirmed the immunogenicity of these epitopes, which induced high-avidity CD8+ T cell clones. These T cells specifically recognize and kill HLA-A2+ tumor cells presenting HERV epitopes on HLA molecules, as demonstrated by mass spectrometry. Furthermore, epitope-specific CD8+ T cells were identified by dextramer staining among tumor-infiltrating lymphocytes from HLA-A2+ patients with breast cancer. Last, we showed that HERV-specific T cells lyse patient-derived organoids. These shared virus-like epitopes are of major interest for the development of cancer vaccines or T cell-based immunotherapies, especially in tumors with low/intermediate mutational burden.

Cross-species analysis of viral nucleic acid interacting proteins identifies TAOKs as innate immune regulators.

The cell intrinsic antiviral response of multicellular organisms developed over millions of years and critically relies on the ability to sense and eliminate viral nucleic acids. Here we use an affinity proteomics approach in evolutionary distant species (human, mouse and fly) to identify proteins that are conserved in their ability to associate with diverse viral nucleic acids. This approach shows a core of orthologous proteins targeting viral genetic material and species-specific interactions. Functional characterization of the influence of 181 candidates on replication of 6 distinct viruses in human cells and flies identifies 128 nucleic acid binding proteins with an impact on virus growth. We identify the family of TAO kinases (TAOK1, -2 and -3) as dsRNA-interacting antiviral proteins and show their requirement for type-I interferon induction. Depletion of TAO kinases in mammals or flies leads to an impaired response to virus infection characterized by a reduced induction of interferon stimulated genes in mammals and impaired expression of srg1 and diedel in flies. Overall, our study shows a larger set of proteins able to mediate the interaction between viral genetic material and host factors than anticipated so far, attesting to the ancestral roots of innate immunity and to the lineage-specific pressures exerted by viruses.

Structural organization of human replication timing domains.

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.

Large replication skew domains delimit GC-poor gene deserts in human.

Besides their large-scale organization in isochores, mammalian genomes display megabase-sized regions, spanning both genes and intergenes, where the strand nucleotide composition asymmetry decreases linearly, possibly due to replication activity. These so-called skew-N domains cover about a third of the human genome and are bordered by two skew upward jumps that were hypothesized to compose a subset of "master" replication origins active in the germline. Skew-N domains were shown to exhibit a particular gene organization. Genes with CpG-rich promoters likely expressed in the germline are over represented near the master replication origins, with large genes being co-oriented with replication fork progression, which suggests some coordination of replication and transcription. In this study, we describe another skew structure that covers ∼13% of the human genome and that is bordered by putative master replication origins similar to the ones flanking skew-N domains. These skew-split-N domains have a shape reminiscent of a N, but split in half, leaving in the center a region of null skew whose length increases with domain size. These central regions (median size ∼860 kb) have a homogeneous composition, i.e. both a null and constant skew and a constant and low GC content. They correspond to heterochromatin gene deserts found in low-GC isochores with an average gene density of 0.81 promoters/Mb as compared to 7.73 promoters/Mb genome wide. The analysis of epigenetic marks and replication timing data confirms that, in these late replicating heterochomatic regions, the initiation of replication is likely to be random. This contrasts with the transcriptionally active euchromatin state found around the bordering well positioned master replication origins. Altogether skew-N domains and skew-split-N domains cover about 50% of the human genome.