A recurrent chromosomal inversion suffices for driving escape from oncogene-induced senescence via subTAD reorganization.

Oncogene-induced senescence (OIS) is an inherent and important tumor suppressor mechanism. However, if not removed timely via immune surveillance, senescent cells also have detrimental effects. Although this has mostly been attributed to the senescence-associated secretory phenotype (SASP) of these cells, we recently proposed that "escape" from the senescent state is another unfavorable outcome. The mechanism underlying this phenomenon remains elusive. Here, we exploit genomic and functional data from a prototypical human epithelial cell model carrying an inducible CDC6 oncogene to identify an early-acquired recurrent chromosomal inversion that harbors a locus encoding the circadian transcription factor BHLHE40. This inversion alone suffices for BHLHE40 activation upon CDC6 induction and driving cell cycle re-entry of senescent cells, and malignant transformation. Ectopic overexpression of BHLHE40 prevented induction of CDC6-triggered senescence. We provide strong evidence in support of replication stress-induced genomic instability being a causative factor underlying "escape" from oncogene-induced senescence.

Human Endogenous Retrovirus-K HML-2 integration within RASGRF2 is associated with intravenous drug abuse and modulates transcription in a cell-line model.

HERV-K HML-2 (HK2) has been proliferating in the germ line of humans at least as recently as 250,000 years ago, with some integrations that remain polymorphic in the modern human population. One of the solitary HK2 LTR polymorphic integrations lies between exons 17 and 18 of RASGRF2, a gene that affects dopaminergic activity and is thus related to addiction. Here we show that this antisense HK2 integration (namely RASGRF2-int) is found more frequently in persons who inject drugs compared with the general population. In a Greek HIV-1-positive population (n = 202), we found RASGRF2-int 2.5 times (14 versus 6%) more frequently in patients infected through i.v. drug use compared with other transmission route controls (P = 0.03). Independently, in a United Kingdom-based hepatitis C virus-positive population (n = 184), we found RASGRF2-int 3.6 times (34 versus 9.5%) more frequently in patients infected during chronic drug abuse compared with controls (P < 0.001). We then tested whether RASGRF2-int could be mechanistically responsible for this association by modulating transcription of RASGRF2 We show that the CRISPR/Cas9-mediated insertion of HK2 in HEK293 cells in the exact RASGRF2 intronic position found in the population resulted in significant transcriptional and phenotypic changes. We also explored mechanistic features of other intronic HK2 integrations and show that HK2 LTRs can be responsible for generation of cis-natural antisense transcripts, which could interfere with the transcription of nearby genes. Our findings suggest that RASGRF2-int is a strong candidate for dopaminergic manipulation, and emphasize the importance of accurate mapping of neglected HERV polymorphisms in human genomic studies.

Nanopore sequencing and full genome de novo assembly of human cytomegalovirus TB40/E reveals clonal diversity and structural variations.

BACKGROUND: Human cytomegalovirus (HCMV) has a double-stranded DNA genome of approximately 235 Kbp that is structurally complex including extended GC-rich repeated regions. Genomic recombination events are frequent in HCMV cultures but have also been observed in vivo. Thus, the assembly of HCMV whole genomes from technologies producing shorter than 500 bp sequences is technically challenging. Here we improved the reconstruction of HCMV full genomes by means of a hybrid, de novo genome-assembly bioinformatics pipeline upon data generated from the recently released MinION MkI B sequencer from Oxford Nanopore Technologies. RESULTS: The MinION run of the HCMV (strain TB40/E) library resulted in ~ 47,000 reads from a single R9 flowcell and in ~ 100× average read depth across the virus genome. We developed a novel, self-correcting bioinformatics algorithm to assemble the pooled HCMV genomes in three stages. In the first stage of the bioinformatics algorithm, long contigs (N50 = 21,892) of lower accuracy were reconstructed. In the second stage, short contigs (N50 = 5686) of higher accuracy were assembled, while in the final stage the high quality contigs served as template for the correction of the longer contigs resulting in a high-accuracy, full genome assembly (N50 = 41,056). We were able to reconstruct a single representative haplotype without employing any scaffolding steps. The majority (98.8%) of the genomic features from the reference strain were accurately annotated on this full genome construct. Our method also allowed the detection of multiple alternative sub-genomic fragments and non-canonical structures suggesting rearrangement events between the unique (UL /US) and the repeated (T/IRL/S) genomic regions. CONCLUSIONS: Third generation high-throughput sequencing technologies can accurately reconstruct full-length HCMV genomes including their low-complexity and highly repetitive regions. Full-length HCMV genomes could prove crucial in understanding the genetic determinants and viral evolution underpinning drug resistance, virulence and pathogenesis.

An Innovative Study Design to Assess the Community Effect of Interventions to Mitigate HIV Epidemics Using Transmission-Chain Phylodynamics.

Given globalization and other social phenomena, controlling the spread of infectious diseases has become an imperative public health priority. A plethora of interventions that in theory can mitigate the spread of pathogens have been proposed and applied. Evaluating the effectiveness of such interventions is costly and in many circumstances unrealistic. Most important, the community effect (i.e., the ability of the intervention to minimize the spread of the pathogen from people who received the intervention to other community members) can rarely be evaluated. Here we propose a study design that can build and evaluate evidence in support of the community effect of an intervention. The approach exploits molecular evolutionary dynamics of pathogens in order to track new infections as having arisen from either a control or an intervention group. It enables us to evaluate whether an intervention reduces the number and length of new transmission chains in comparison with a control condition, and thus lets us estimate the relative decrease in new infections in the community due to the intervention. We provide as an example one working scenario of a way the approach can be applied with a simulation study and associated power calculations.

A novel method for the multiplexed target enrichment of MinION next generation sequencing libraries using PCR-generated baits.

The enrichment of targeted regions within complex next generation sequencing libraries commonly uses biotinylated baits to capture the desired sequences. This method results in high read coverage over the targets and their flanking regions. Oxford Nanopore Technologies recently released an USB3.0-interfaced sequencer, the MinION. To date no particular method for enriching MinION libraries has been standardized. Here, using biotinylated PCR-generated baits in a novel approach, we describe a simple and efficient way for multiplexed enrichment of MinION libraries, overcoming technical limitations related with the chemistry of the sequencing-adapters and the length of the DNA fragments. Using Phage Lambda and Escherichia coli as models we selectively enrich for specific targets, significantly increasing the corresponding read-coverage, eliminating unwanted regions. We show that by capturing genomic fragments, which contain the target sequences, we recover reads extending targeted regions and thus can be used for the determination of potentially unknown flanking sequences. By pooling enriched libraries derived from two distinct E. coli strains and analyzing them in parallel, we demonstrate the efficiency of this method in multiplexed format. Crucially we evaluated the optimal bait size for large fragment libraries and we describe for the first time a standardized method for target enrichment in MinION platform.

The interferon receptor-1 promoter polymorphisms affect the outcome of Caucasians with HBeAg-negative chronic HBV infection.

BACKGROUND & AIMS: The outcome of HBeAg-negative chronic hepatitis B virus (HBV) patients who may remain in the inactive carrier state (IC) or progress to HBeAg-negative chronic hepatitis B may be affected by the host genetic profile. Genetic polymorphisms within not only the promoter but also the coding sequence of the interferon receptor 1 (INFAR1) gene have been associated with susceptibility to chronic HBV infection, but their role on the outcomes of HBeAg-negative patients has not been evaluated. We examined the association of INFAR1 promoter polymorphisms with the phase of chronic HBV infection in a demographically characterized Caucasian cohort of 183 consecutive HBeAg-negative chronic HBV patients. METHODS: Using a combination of conventional and allele-specific polymerase chain reactions, bidirectional sequencing and DNA-fragment analysis, we performed typing of three Single Nucleotide Polymorphisms (SNPs -568G/C, -408C/T, -3C/T) and one Variable Number Tandem Repeat [VNTR -77(GT)n] within the INFR1 promoter sequence. RESULTS: The genetic polymorphisms examined were found to be associated with the phase of HBeAg-negative chronic HBV patients. Using a multiple logistic regression model adjusting for age, gender and origin of the individuals, we found that patients with linked genotypes -408CT_-3CT were more likely to be ICs (OR = 2.42 vs. CC, P = 0.036). Also, given the partial linkage between SNP -568G/C and VNTR -77(GT)n, we found that linked genotypes -77(GT)n ≤ 8/≤8_-568GC and -77(GT)n ≤ 8/≤8_-568CC were detected more frequently among ICs (OR = 11.69, P = 0.005 and OR = 7.56, P = 0.001 vs. -77(GT)n >8/>8_-568GG respectively). CONCLUSIONS: These findings suggest that these genetic variations represent important factors associated with the clinical phase of HBeAg-negative chronic HBV infection.

Dual RNA-seq Analysis of Patients' Cells and Viral Genome After Measles Infection.

During the infection of a host cell by an infectious agent, a series of gene expression changes occurs as a consequence of host-pathogen interactions. Unraveling this complex interplay is the key for understanding of microbial virulence and host response pathways, thus providing the basis for new molecular insights into the mechanisms of pathogenesis and the corresponding immune response. Dual RNA sequencing (dual RNA-seq) has been developed to simultaneously determine pathogen and host transcriptomes enabling both differential and coexpression analyses between the two partners as well as genome characterization in the case of RNA viruses. Here, we provide a detailed laboratory protocol and bioinformatics analysis guidelines for dual RNA-seq experiments focusing on - but not restricted to - measles virus (MeV) as a pathogen of interest. The application of dual RNA-seq technologies in MeV-infected patients can potentially provide valuable information on the structure of the viral RNA genome and on cellular innate immune responses and drive the discovery of new targets for antiviral therapy.

Ultrasensitive amplification refractory mutation system real-time PCR (ARMS RT-PCR) assay for detection of minority hepatitis B virus-resistant strains in the era of personalized medicine.

Resistance to antiviral treatment for chronic hepatitis B virus (HBV) has been associated with mutations in the HBV polymerase region. This study aimed at developing an ultrasensitive method for quantifying viral populations with all major HBV resistance-associated mutations, combining the amplification refractory mutation system real-time PCR (ARMS RT-PCR) with a molecular beacon using a LightCycler. The discriminatory ability of this method, the ARMS RT-PCR with molecular beacon assay, was 0.01 to 0.25% for the different HBV resistance-associated mutations, as determined by laboratory-synthesized wild-type (WT) and mutant (Mut) target sequences. The assay showed 100% sensitivity for the detection of mutant variants A181V, T184A, and N236T in samples from 41 chronically HBV-infected patients under antiviral therapy who had developed resistance-associated mutations detected by direct PCR Sanger sequencing. The ratio of mutant to wild-type viral populations (the Mut/WT ratio) was >1% in 38 (63.3%) of 60 samples from chronically HBV-infected nucleos(t)ide analogue-naive patients; combinations of mutations were also detected in half of these samples. The ARMS RT-PCR with molecular beacon assay achieved high sensitivity and discriminatory ability compared to the gold standard of direct PCR Sanger sequencing in identifying resistant viral populations in chronically HBV-infected patients receiving antiviral therapy. Apart from the dominant clones, other quasispecies were also quantified. In samples from chronically HBV-infected nucleos(t)ide analogue-naive patients, the assay proved to be a useful tool in detecting minor variant populations before the initiation of the treatment. These observations need further evaluation with prospective studies before they can be implemented in daily practice.

Nanopore-Sequencing Metabarcoding for Identification of Phytopathogenic and Endophytic Fungi in Olive (Olea europaea) Twigs.

Metabarcoding approaches for the identification of plant disease pathogens and characterization of plant microbial populations constitute a rapidly evolving research field. Fungal plant diseases are of major phytopathological concern; thus, the development of metabarcoding approaches for the detection of phytopathogenic fungi is becoming increasingly imperative in the context of plant disease prognosis. We developed a multiplex metabarcoding method for the identification of fungal phytopathogens and endophytes in olive young shoots, using the MinION sequencing platform (Oxford Nanopore Technologies). Selected fungal-specific primers were used to amplify three different genomic DNA loci (ITS, beta-tubulin, and 28S LSU) originating from olive twigs. A multiplex metabarcoding approach was initially evaluated using healthy olive twigs, and further assessed with naturally infected olive twig samples. Bioinformatic analysis of basecalled reads was carried out using MinKNOW, BLAST+ and R programming, and results were also evaluated using the BugSeq cloud platform. Data analysis highlighted the approaches based on ITS and their combination with beta-tubulin as the most informative ones according to diversity estimations. Subsequent implementation of the method on symptomatic samples identified major olive pathogens and endophytes including genera such as Cladosporium, Didymosphaeria, Paraconiothyrium, Penicillium, Phoma, Verticillium, and others.

A miR-124-mediated post-transcriptional mechanism controlling the cell fate switch of astrocytes to induced neurons.

The microRNA (miRNA) miR-124 has been employed supplementary to neurogenic transcription factors (TFs) and other miRNAs to enhance direct neurogenic conversion. The aim of this study was to investigate whether miR-124 is sufficient to drive direct reprogramming of astrocytes to induced neurons (iNs) on its own and elucidate its independent mechanism of reprogramming action. Our data show that miR-124 is a potent driver of the reprogramming switch of astrocytes toward an immature neuronal fate by directly targeting the RNA-binding protein Zfp36L1 implicated in ARE-mediated mRNA decay and subsequently derepressing Zfp36L1 neurogenic interactome. To this end, miR-124 contribution in iNs' production largely recapitulates endogenous neurogenesis pathways, being further enhanced upon addition of the neurogenic compound ISX9, which greatly improves iNs' differentiation and functional maturation. Importantly, miR-124 is potent in guiding direct conversion of reactive astrocytes to immature iNs in vivo following cortical trauma, while ISX9 supplementation confers a survival advantage to newly produced iNs.

Hepatitis B Virus Integration into Transcriptionally Active Loci and HBV-Associated Hepatocellular Carcinoma.

Hepatitis B Virus (HBV) DNA integrations into the human genome are considered major causative factors to HBV-associated hepatocellular carcinoma development. In the present study, we investigated whether HBV preferentially integrates parts of its genome in specific genes and evaluated the contribution of the integrations in HCC development per gene. We applied dedicated in-house developed pipelines on all of the available HBV DNA integration data and performed a statistical analysis to identify genes that could be characterized as hotspots of integrations, along with the evaluation of their association with HBV-HCC. Our results suggest that 15 genes are recurrently affected by HBV integrations and they are significantly associated with HBV-HCC. Further studies that focus on HBV integrations disrupting these genes are mandatory in order to understand the role of HBV integrations in clonal advantage gain and oncogenesis promotion, as well as to determine whether inhibition of the HBV-disrupted genes can provide a therapy strategy for HBV-HCC.

Molecular and Clinical Prognostic Biomarkers of COVID-19 Severity and Persistence.

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), poses several challenges to clinicians, due to its unpredictable clinical course. The identification of laboratory biomarkers, specific cellular, and molecular mediators of immune response could contribute to the prognosis and management of COVID-19 patients. Of utmost importance is also the detection of differentially expressed genes, which can serve as transcriptomic signatures, providing information valuable to stratify patients into groups, based on the severity of the disease. The role of biomarkers such as IL-6, procalcitonin, neutrophil-lymphocyte ratio, white blood cell counts, etc. has already been highlighted in recently published studies; however, there is a notable amount of new evidence that has not been summarized yet, especially regarding transcriptomic signatures. Hence, in this review, we assess the latest cellular and molecular data and determine the significance of abnormalities in potential biomarkers for COVID-19 severity and persistence. Furthermore, we applied Gene Ontology (GO) enrichment analysis using the genes reported as differentially expressed in the literature in order to investigate which biological pathways are significantly enriched. The analysis revealed a number of processes, such as inflammatory response, and monocyte and neutrophil chemotaxis, which occur as part of the complex immune response to SARS-CoV-2.

Early Signs of Molecular Defects in iPSC-Derived Neural Stems Cells from Patients with Familial Parkinson's Disease.

Parkinson's disease (PD) is the second most common neurodegenerative disorder, classically associated with extensive loss of dopaminergic neurons of the substantia nigra pars compacta. The hallmark of the disease is the accumulation of pathogenic conformations of the presynaptic protein, α-synuclein (αSyn), and the formation of intraneuronal protein aggregate inclusions. Neurodegeneration of dopamine neurons leads to a prominent dopaminergic deficiency in the basal ganglia, responsible for motor disturbances. However, it is now recognized that the disease involves more widespread neuronal dysfunction, leading to early and late non-motor symptoms. The development of in vitro systems based on the differentiation of human-induced pluripotent stem cells provides us the unique opportunity to monitor alterations at the cellular and molecular level throughout the differentiation procedure and identify perturbations that occur early, even at the neuronal precursor stage. Here we aim to identify whether p.A53T-αSyn induced disturbances at the molecular level are already present in neural precursors. Towards this, we present data from transcriptomics analysis of control and p.A53T-αSyn NPCs showing altered expression in transcripts involved in axon guidance, adhesion, synaptogenesis, ion transport, and metabolism. The comparative analysis with the transcriptomics profile of p.A53T-αSyn neurons shows both distinct and overlapping pathways leading to neurodegeneration while meta-analysis with transcriptomics data from both neurodegenerative and neurodevelopmental disorders reveals that p.A53T-pathology has a significant overlap with the latter category. This is the first study showing that molecular dysregulation initiates early at the p.A53T-αSyn NPC level, suggesting that synucleinopathies may have a neurodevelopmental component.

Targeted Virome Sequencing Enhances Unbiased Detection and Genome Assembly of Known and Emerging Viruses-The Example of SARS-CoV-2.

Targeted virome enrichment and sequencing (VirCapSeq-VERT) utilizes a pool of oligos (baits) to enrich all known­up to 2015­vertebrate-infecting viruses, increasing their detection sensitivity. The hybridisation of the baits to the target sequences can be partial, thus enabling the detection and genomic reconstruction of novel pathogens with <40% genetic diversity compared to the strains used for the baits' design. In this study, we deploy this method in multiplexed mixes of viral extracts, and we assess its performance in the unbiased detection of DNA and RNA viruses after cDNA synthesis. We further assess its efficiency in depleting various background genomic material. Finally, as a proof-of-concept, we explore the potential usage of the method for the characterization of unknown, emerging human viruses, such as SARS-CoV-2, which may not be included in the baits' panel. We mixed positive samples of equimolar DNA/RNA viral extracts from SARS-CoV-2, coronavirus OC43, cytomegalovirus, influenza A virus H3N2, parvovirus B19, respiratory syncytial virus, adenovirus C and coxsackievirus A16. Targeted virome enrichment was performed on a dsDNA mix, followed by sequencing on the NextSeq500 (Illumina) and the portable MinION sequencer, to evaluate its usability as a point-of-care (PoC) application. Genome mapping assembly was performed using viral reference sequences. The untargeted libraries contained less than 1% of total reads mapped on most viral genomes, while RNA viruses remained undetected. In the targeted libraries, the percentage of viral-mapped reads were substantially increased, allowing full genome assembly in most cases. Targeted virome sequencing can enrich a broad range of viruses, potentially enabling the discovery of emerging viruses.

Tracing the First Days of the SARS-CoV-2 Pandemic in Greece and the Role of the First Imported Group of Travelers.

The first SARS-CoV-2 case in Greece was confirmed on February 26, 2020, and since then, multiple strains have circulated the country, leading to regional and country-wide outbreaks. Our aim is to enlighten the events that took place during the first days of the SARS-CoV-2 pandemic in Greece, focusing on the role of the first imported group of travelers. We used whole-genome SARS-CoV-2 sequences obtained from the infected travelers of the group as well as Greece-derived and globally subsampled sequences and applied dedicated phylogenetics and phylodynamics tools as well as in-house-developed bioinformatics pipelines. Our analyses reveal the genetic variants circulating in Greece during the first days of the pandemic and the role of the group's imported strains in the course of the first pandemic wave in Greece. The strain that dominated in Greece throughout the first wave, bearing the D614G mutation, was primarily imported from a certain group of travelers, while molecular and clinical data suggest that the infection of the travelers occurred in Egypt. Founder effects early in the pandemic are important for the success of certain strains, as those arriving early, several times, and to diverse locations lead to the formation of large transmission clusters that can be estimated using molecular epidemiology approaches and can be a useful surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks. IMPORTANCE The strain that dominated in Greece during the first pandemic wave was primarily imported from a group of returning travelers in February 2020, while molecular and clinical data suggest that the origin of the transmission was Egypt. The observed molecular transmission clusters reflect the transmission dynamics of this particular strain bearing the D614G mutation while highlighting the necessity of their use as a surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks.

Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations.

Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori from Europe and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks by migration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.

Hepatitis B Virus DNA Integration, Chronic Infections and Hepatocellular Carcinoma.

Hepatitis B Virus (HBV) is an Old World virus with a high mutation rate, which puts its origins in Africa alongside the origins of Homo sapiens, and is a member of the Hepadnaviridae family that is characterized by a unique viral replication cycle. It targets human hepatocytes and can lead to chronic HBV infection either after acute infection via horizontal transmission usually during infancy or childhood or via maternal-fetal transmission. HBV has been found in ~85% of HBV-related Hepatocellular Carcinomas (HCC), and it can integrate the whole or part of its genome into the host genomic DNA. The molecular mechanisms involved in the HBV DNA integration is not yet clear; thus, multiple models have been described with respect to either the relaxed-circular DNA (rcDNA) or the double-stranded linear DNA (dslDNA) of HBV. Various genes have been found to be affected by HBV DNA integration, including cell-proliferation-related genes, oncogenes and long non-coding RNA genes (lincRNAs). The present review summarizes the advances in the research of HBV DNA integration, focusing on the evolutionary and molecular side of the integration events along with the arising clinical aspects in the light of WHO's commitment to eliminate HBV and viral hepatitis by 2030.

Dual RNA-Seq Enables Full-Genome Assembly of Measles Virus and Characterization of Host-Pathogen Interactions.

Measles virus (MeV) has a negative-sense 15 kb long RNA genome, which is generally conserved. Recent advances in high-throughput sequencing (HTS) and Dual RNA-seq allow the analysis of viral RNA genomes and the discovery of viral infection biomarkers, via the simultaneous characterization of the host transcriptome. However, these host-pathogen interactions remain largely unexplored in MeV infections. We performed untargeted Dual RNA-seq in 6 pharyngeal and 6 peripheral blood mononuclear cell (PBMCs) specimens from patients with MeV infection, as confirmed via routine real-time PCR testing. Following optimised DNase treatment of total nucleic acids, we used the pharyngeal samples to build poly-A-enriched NGS libraries. We reconstructed the viral genomes using the pharyngeal datasets and we further conducted differential expression, gene-ontology and pathways enrichment analysis to compare both the pharyngeal and the peripheral blood transcriptomes of the MeV-infected patients vs. control groups of healthy individuals. We obtained 6 MeV genotype-B3 full-genome sequences. We minutely analyzed the transcriptome of the MeV-infected pharyngeal epithelium, detecting all known viral infection biomarkers, but also revealing a functional cluster of local antiviral and inflammatory immune responses, which differ substantially from those observed in the PBMCs transcriptome. The application of Dual RNA-seq technologies in MeV-infected patients can potentially provide valuable information on the virus genome structure and the cellular innate immune responses and drive the discovery of new targets for antiviral therapy.

Upregulation of Human Endogenous Retroviruses in Bronchoalveolar Lavage Fluid of COVID-19 Patients.

Severe COVID-19 pneumonia has been associated with the development of intense inflammatory responses during the course of infections with SARS-CoV-2. Given that human endogenous retroviruses (HERVs) are known to be activated during and participate in inflammatory processes, we examined whether HERV dysregulation signatures are present in COVID-19 patients. By comparing transcriptomes of bronchoalveolar lavage fluid (BALF) of COVID-19 patients and healthy controls, and peripheral blood monocytes (PBMCs) from patients and controls, we have shown that HERVs are intensely dysregulated in BALF of COVID-19 patients compared to those in BALF of healthy control patients but not in PBMCs. In particular, upregulation in the expression of specific HERV families was detected in BALF samples of COVID-19 patients, with HERV-FRD being the most highly upregulated family among the families analyzed. In addition, we compared the expression of HERVs in human bronchial epithelial cells (HBECs) without and after senescence induction in an oncogene-induced senescence model in order to quantitatively measure changes in the expression of HERVs in bronchial cells during the process of cellular senescence. This apparent difference of HERV dysregulation between PBMCs and BALF warrants further studies in the involvement of HERVs in inflammatory pathogenetic mechanisms as well as exploration of HERVs as potential biomarkers for disease progression. Furthermore, the increase in the expression of HERVs in senescent HBECs in comparison to that in noninduced HBECs provides a potential link for increased COVID-19 severity and mortality in aged populations. IMPORTANCE SARS-CoV-2 emerged in late 2019 in China, causing a global pandemic. Severe COVID-19 is characterized by intensive inflammatory responses, and older age is an important risk factor for unfavorable outcomes. HERVs are remnants of ancient infections whose expression is upregulated in multiple conditions, including cancer and inflammation, and their expression is increased with increasing age. The significance of this work is that we were able to recognize dysregulated expression of endogenous retroviral elements in BALF samples but not in PBMCs of COVID-19 patients. At the same time, we were able to identify upregulated expression of multiple HERV families in senescence-induced HBECs in comparison to that in noninduced HBECs, a fact that could possibly explain the differences in disease severity among age groups. These results indicate that HERV expression might play a pathophysiological role in local inflammatory pathways in lungs afflicted by SARS-CoV-2 and their expression could be a potential therapeutic target.

SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave.

BACKGROUND: The spatiotemporal profiling of molecular transmission clusters (MTCs) using viral genomic data can effectively identify transmission networks in order to inform public health actions targeting SARS-CoV-2 spread. METHODS: We used whole genome SARS-CoV-2 sequences derived from ten European regions belonging to eight countries to perform phylogenetic and phylodynamic analysis. We developed dedicated bioinformatics pipelines to identify regional MTCs and to assess demographic factors potentially associated with their formation. RESULTS: The total number and the scale of MTCs varied from small household clusters identified in all regions, to a super-spreading event found in Uusimaa-FI. Specific age groups were more likely to belong to MTCs in different regions. The clustered sequences referring to the age groups 50-100 years old (y.o.) were increased in all regions two weeks after the establishment of the lockdown, while those referring to the age group 0-19 y.o. decreased only in those regions where schools' closure was combined with a lockdown. CONCLUSIONS: The spatiotemporal profiling of the SARS-CoV-2 MTCs can be a useful tool to monitor the effectiveness of the interventions and to reveal cryptic transmissions that have not been identified through contact tracing.