Superspreading events shaped the Coronavirus Disease 2019 (COVID-19) pandemic, and their rapid identification and containment are essential for disease control. Here we provide a national-scale analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) superspreading during the first wave of infections in Austria, a country that played a major role in initial virus transmissions in Europe. Capitalizing on Austria’s well-developed epidemiological surveillance system, we identified major SARS-CoV-2 clusters during the first wave of infections and performed deep whole-genome sequencing of more than 500 virus samples. Phylogenetic-epidemiological analysis enabled the reconstruction of superspreading events and charts a map of tourism-related viral spread originating from Austria in spring 2020. Moreover, we exploited epidemiologically well-defined clusters to quantify SARS-CoV-2 mutational dynamics, including the observation of a low-frequency mutation that progressed to fixation within the infection chain. Time-resolved virus sequencing unveiled viral mutation dynamics within individuals with COVID-19, and epidemiologically validated infector-infectee pairs enabled us to determine an average transmission bottleneck size of 103 SARS-CoV-2 particles. In conclusion, this study illustrates the power of combining epidemiological analysis with deep viral genome sequencing to unravel the spread of SARS-CoV-2, and to gain fundamental insights into mutational dynamics and transmission properties.

Fragile X syndrome (FXS), the most common form of inherited intellectual disability, is due to the functional deficiency of the fragile X mental retardation protein (FMRP), an RNA-binding protein involved in translational regulation of many messenger RNAs, playing key roles in synaptic morphology and plasticity. To date, no effective treatment for FXS is available. We searched for FMRP targets by HITS-CLIP during early development of multiple mouse brain regions (hippocampus, cortex and cerebellum) at a time of brain development when FMRP is most highly expressed and synaptogenesis reaches a peak. We identified the largest dataset of mRNA targets of FMRP available in brain and we defined their cellular origin. We confirmed the G-quadruplex containing structure as an enriched motif in FMRP RNA targets. In addition to four less represented motifs, our study points out that, in the brain, CTGKA is the prominent motif bound by FMRP, which recognizes it when not engaged in Watson–Crick pairing. All of these motifs negatively modulated the expression level of a reporter protein. While the repertoire of FMRP RNA targets in cerebellum is quite divergent, the ones of cortex and hippocampus are vastly overlapping. In these two brain regions, the Phosphodiesterase 2a (Pde2a) mRNA is a prominent target of FMRP, which modulates its translation and intracellular transport. This enzyme regulates the homeostasis of cAMP and cGMP and represents a novel and attractive therapeutic target to treat FXS.

CD8+ T cell immunity to SARS-CoV-2 has been implicated in COVID-19 severity and virus control. Here, we identified nonsynonymous mutations in MHC-I-restricted CD8+ T cell epitopes after deep sequencing of 747 SARS-CoV-2 virus isolates. Mutant peptides exhibited diminished or abrogated MHC-I binding in a cell-free in vitro assay. Reduced MHC-I binding of mutant peptides was associated with decreased proliferation, IFN-γ production and cytotoxic activity of CD8+ T cells isolated from HLA-matched COVID-19 patients. Single cell RNA sequencing of ex vivo expanded, tetramer-sorted CD8+ T cells from COVID-19 patients further revealed qualitative differences in the transcriptional response to mutant peptides. Our findings highlight the capacity of SARS-CoV-2 to subvert CD8+ T cell surveillance through point mutations in MHC-I-restricted viral epitopes.

GC-biased gene conversion (gBGC) is a process that tends to increase the GC content of recombining DNA over evolutionary time and is thought to explain the evolution of GC content in mammals and yeasts. Evidence for gBGC outside these two groups is growing but is still limited. Here, we analyzed 36 completely sequenced genomes representing four of the five major groups in eukaryotes (Unikonts, Excavates, Chromalveolates and Plantae). gBGC was investigated by directly comparing GC content and recombination rates in species where recombination data are available, that is, half of them. To study all species of our dataset, we used chromosome size as a proxy for recombination rate and compared it with GC content. Among the 17 species showing a significant relationship between GC content and chromosome size, 15 are consistent with the predictions of the gBGC model. Importantly, the species showing a pattern consistent with gBGC are found in all the four major groups of eukaryotes studied, which suggests that gBGC may be widespread in eukaryotes.

Graphical Abstract Highlights d Chronic viral infection causes long-term transcriptome and proteome changes in the liver d Hepatocyte-intrinsic type I interferon (IFN-I) signaling regulates hepatic metabolism d IFN-I signaling reprograms the urea cycle in hepatocytes and alters serum metabolites d Serum levels of arginine and ornithine modulate T cell responses and pathology

BackgroundMost eukaryotic genes are subject to alternative splicing (AS), which may contribute to the production of protein variants or to the regulation of gene expression via nonsense-mediated messenger RNA (mRNA) decay (NMD). However, a fraction of splice variants might correspond to spurious transcripts and the question of the relative proportion of splicing errors to functional splice variants remains highly debated.ResultsWe propose a test to quantify the fraction of AS events corresponding to errors. This test is based on the fact that the fitness cost of splicing errors increases with the number of introns in a gene and with expression level. We analyzed the transcriptome of the intron-rich eukaryote Paramecium tetraurelia. We show that in both normal and in NMD-deficient cells, AS rates strongly decrease with increasing expression level and with increasing number of introns. This relationship is observed for AS events that are detectable by NMD as well as for those that are not, which invalidates the hypothesis of a link with the regulation of gene expression. Our results show that in genes with a median expression level, 92–98% of observed splice variants correspond to errors. We observed the same patterns in human transcriptomes and we further show that AS rates correlate with the fitness cost of splicing errors.ConclusionsThese observations indicate that genes under weaker selective pressure accumulate more maladaptive substitutions and are more prone to splicing errors. Thus, to a large extent, patterns of gene expression variants simply reflect the balance between selection, mutation, and drift.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1344-6) contains supplementary material, which is available to authorized users.

Cachexia represents a leading cause of morbidity and mortality in various cancers, chronic inflammation and infections. Understanding of the mechanisms that drive cachexia has remained limited, especially for infection-associated cachexia (IAC). Here we describe a model of reversible cachexia in mice with chronic viral infection and identify an essential role for CD8+ T cells in IAC. Cytokines linked to cancer-associated cachexia did not contribute to IAC. Instead, virus-specific CD8+ T cells caused morphological and molecular changes in the adipose tissue, which led to depletion of lipid stores. These changes occurred at a time point that preceded the peak of the CD8+ T cell response and required T cell–intrinsic type 1 interferon signaling and antigen-specific priming. Our results link systemic antiviral immune responses to adipose-tissue remodeling and reveal an underappreciated role of CD8+ T cells in IAC.

Incidence of cutaneous squamous cell carcinomas (cSCCs) constantly increases in the Caucasian population. Developing preferentially on precancerous lesions such as actinic keratoses due to chronic sunlight exposure, cSCCs result from the malignant transformation of keratinocytes. Although a resection of the primary tumor is usually curative, a subset of aggressive cSCCs shows a high risk of recurrence and metastases. The characterization of the molecular dysfunctions involved in cSCC development should help to identify new relevant targets against these aggressive cSCCs. In that context, we have used small RNA sequencing to identify 100 microRNAs (miRNAs) whose expression was altered during chemically induced mouse skin tumorigenesis. The decreased expression of the miR-193b/365a cluster during tumor progression suggests a tumor suppressor role. Ectopic expression of these miRNAs in tumor cells indeed inhibited their proliferation, clonogenic potential and migration, which were stimulated in normal keratinocytes when these miRNAs were blocked with antisense oligonucleotides. A combination of in silico predictions and transcriptome analyses identified several target genes of interest. We validated KRAS and MAX as direct targets of miR-193b and miR-365a. Repression of these targets using siRNAs mimicked the effects of miR-193b and miR-365a, suggesting that these genes might mediate, at least in part, the tumor-suppressive action of these miRNAs.