Genomic surveillance of SARS-CoV-2 evolution by a centralised pipeline and weekly focused sequencing, Austria, January 2021 to March 2023.

BackgroundThe COVID-19 pandemic was largely driven by genetic mutations of SARS-CoV-2, leading in some instances to enhanced infectiousness of the virus or its capacity to evade the host immune system. To closely monitor SARS-CoV-2 evolution and resulting variants at genomic-level, an innovative pipeline termed SARSeq was developed in Austria.AimWe discuss technical aspects of the SARSeq pipeline, describe its performance and present noteworthy results it enabled during the pandemic in Austria.MethodsThe SARSeq pipeline was set up as a collaboration between private and public clinical diagnostic laboratories, a public health agency, and an academic institution. Representative SARS-CoV-2 positive specimens from each of the nine Austrian provinces were obtained from SARS-CoV-2 testing laboratories and processed centrally in an academic setting for S-gene sequencing and analysis.ResultsSARS-CoV-2 sequences from up to 2,880 cases weekly resulted in 222,784 characterised case samples in January 2021-March 2023. Consequently, Austria delivered the fourth densest genomic surveillance worldwide in a very resource-efficient manner. While most SARS-CoV-2 variants during the study showed comparable kinetic behaviour in all of Austria, some, like Beta, had a more focused spread. This highlighted multifaceted aspects of local population-level acquired immunity. The nationwide surveillance system enabled reliable nowcasting. Measured early growth kinetics of variants were predictive of later incidence peaks.ConclusionWith low automation, labour, and cost requirements, SARSeq is adaptable to monitor other pathogens and advantageous even for resource-limited countries. This multiplexed genomic surveillance system has potential as a rapid response tool for future emerging threats.

Loss of cohesin regulator PDS5A reveals repressive role of Polycomb loops.

Polycomb Repressive Complexes 1 and 2 (PRC1, PRC2) are conserved epigenetic regulators that promote transcriptional gene silencing. PRC1 and PRC2 converge on shared targets, catalyzing repressive histone modifications. Additionally, a subset of PRC1/PRC2 targets engage in long-range interactions whose functions in gene silencing are poorly understood. Using a CRISPR screen in mouse embryonic stem cells, we found that the cohesin regulator PDS5A links transcriptional silencing by Polycomb and 3D genome organization. PDS5A deletion impairs cohesin unloading and results in derepression of a subset of endogenous PRC1/PRC2 target genes. Importantly, derepression is not linked to loss of Polycomb chromatin domains. Instead, PDS5A removal causes aberrant cohesin activity leading to ectopic insulation sites, which disrupt the formation of ultra-long Polycomb loops. We show that these loops are important for robust silencing at a subset of PRC1/PRC2 target genes and that maintenance of cohesin-dependent genome architecture is critical for Polycomb regulation.

Identification of CCZ1 as an essential lysosomal trafficking regulator in Marburg and Ebola virus infections.

Marburg and Ebola filoviruses are two of the deadliest infectious agents and several outbreaks have occurred in the last decades. Although several receptors and co-receptors have been reported for Ebola virus, key host factors remain to be elucidated. In this study, using a haploid cell screening platform, we identify the guanine nucleotide exchange factor CCZ1 as a key host factor in the early stage of filovirus replication. The critical role of CCZ1 for filovirus infections is validated in 3D primary human hepatocyte cultures and human blood-vessel organoids, both critical target sites for Ebola and Marburg virus tropism. Mechanistically, CCZ1 controls early to late endosomal trafficking of these viruses. In addition, we report that CCZ1 has a role in the endosomal trafficking of endocytosis-dependent SARS-CoV-2 infections, but not in infections by Lassa virus, which enters endo-lysosomal trafficking at the late endosome stage. Thus, we have identified an essential host pathway for filovirus infections in cell lines and engineered human target tissues. Inhibition of CCZ1 nearly completely abolishes Marburg and Ebola infections. Thus, targeting CCZ1 could potentially serve as a promising drug target for controlling infections caused by various viruses, such as SARS-CoV-2, Marburg, and Ebola.

Ascl1 and Ngn2 convert mouse embryonic stem cells to neurons via functionally distinct paths.

Ascl1 and Ngn2, closely related proneural transcription factors, are able to convert mouse embryonic stem cells into induced neurons. Despite their similarities, these factors elicit only partially overlapping transcriptional programs, and it remains unknown whether cells are converted via distinct mechanisms. Here we show that Ascl1 and Ngn2 induce mutually exclusive side populations by binding and activating distinct lineage drivers. Furthermore, Ascl1 rapidly dismantles the pluripotency network and installs neuronal and trophoblast cell fates, while Ngn2 generates a neural stem cell-like intermediate supported by incomplete shutdown of the pluripotency network. Using CRISPR-Cas9 knockout screening, we find that Ascl1 relies more on factors regulating pluripotency and the cell cycle, such as Tcf7l1. In the absence of Tcf7l1, Ascl1 still represses core pluripotency genes but fails to exit the cell cycle. However, overexpression of Cdkn1c induces cell cycle exit and restores the generation of neurons. These findings highlight that cell type conversion can occur through two distinct mechanistic paths, even when induced by closely related transcription factors.

Single-cell brain organoid screening identifies developmental defects in autism.

The development of the human brain involves unique processes (not observed in many other species) that can contribute to neurodevelopmental disorders1-4. Cerebral organoids enable the study of neurodevelopmental disorders in a human context. We have developed the CRISPR-human organoids-single-cell RNA sequencing (CHOOSE) system, which uses verified pairs of guide RNAs, inducible CRISPR-Cas9-based genetic disruption and single-cell transcriptomics for pooled loss-of-function screening in mosaic organoids. Here we show that perturbation of 36 high-risk autism spectrum disorder genes related to transcriptional regulation uncovers their effects on cell fate determination. We find that dorsal intermediate progenitors, ventral progenitors and upper-layer excitatory neurons are among the most vulnerable cell types. We construct a developmental gene regulatory network of cerebral organoids from single-cell transcriptomes and chromatin modalities and identify autism spectrum disorder-associated and perturbation-enriched regulatory modules. Perturbing members of the BRG1/BRM-associated factor (BAF) chromatin remodelling complex leads to enrichment of ventral telencephalon progenitors. Specifically, mutating the BAF subunit ARID1B affects the fate transition of progenitors to oligodendrocyte and interneuron precursor cells, a phenotype that we confirmed in patient-specific induced pluripotent stem cell-derived organoids. Our study paves the way for high-throughput phenotypic characterization of disease susceptibility genes in organoid models with cell state, molecular pathway and gene regulatory network readouts.

FAM3C/ILEI protein is elevated in psoriatic lesions and triggers psoriasiform hyperproliferation in mice.

FAM3C/ILEI is an important cytokine for tumor progression and metastasis. However, its involvement in inflammation remains elusive. Here, we show that ILEI protein is highly expressed in psoriatic lesions. Inducible keratinocyte-specific ILEI overexpression in mice (K5-ILEIind ) recapitulates many aspects of psoriasis following TPA challenge, primarily manifested by impaired epidermal differentiation and increased neutrophil recruitment. Mechanistically, ILEI triggers Erk and Akt signaling, which then activates STAT3 via Ser727 phosphorylation. Keratinocyte-specific ILEI deletion ameliorates TPA-induced skin inflammation. A transcriptomic ILEI signature obtained from the K5-ILEIind model shows enrichment in several signaling pathways also found in psoriasis and identifies urokinase as a targetable enzyme to counteract ILEI activity. Pharmacological inhibition of urokinase in TPA-induced K5-ILEIind mice results in significant improvement of psoriasiform symptoms by reducing ILEI secretion. The ILEI signature distinguishes psoriasis from healthy skin with uPA ranking among the top "separator" genes. Our study identifies ILEI as a key driver in psoriasis, indicates the relevance of ILEI-regulated genes for disease manifestation, and shows the clinical impact of ILEI and urokinase as novel potential therapeutic targets in psoriasis.

Low-density lipoprotein receptor-related protein 1 (LRP1) as an auxiliary host factor for RNA viruses.

Viruses with an RNA genome are often the cause of zoonotic infections. In order to identify novel pro-viral host cell factors, we screened a haploid insertion-mutagenized mouse embryonic cell library for clones that are resistant to Rift Valley fever virus (RVFV). This screen returned the low-density lipoprotein receptor-related protein 1 (LRP1) as a top hit, a plasma membrane protein involved in a wide variety of cell activities. Inactivation of LRP1 in human cells reduced RVFV RNA levels already at the attachment and entry stages of infection. Moreover, the role of LRP1 in promoting RVFV infection was dependent on physiological levels of cholesterol and on endocytosis. In the human cell line HuH-7, LRP1 also promoted early infection stages of sandfly fever Sicilian virus and La Crosse virus, but had a minor effect on late infection by vesicular stomatitis virus, whereas encephalomyocarditis virus was entirely LRP1-independent. Moreover, siRNA experiments in human Calu-3 cells demonstrated that also SARS-CoV-2 infection benefitted from LRP1. Thus, we identified LRP1 as a host factor that supports infection by a spectrum of RNA viruses.

A whole-genome scan for Artemisinin cytotoxicity reveals a novel therapy for human brain tumors.

The natural compound Artemisinin is the most widely used antimalarial drug worldwide. Based on its cytotoxicity, it is also used for anticancer therapy. Artemisinin and its derivates are endoperoxides that damage proteins in eukaryotic cells; their definite mechanism of action and host cell targets, however, have remained largely elusive. Using yeast and haploid stem cell screening, we demonstrate that a single cellular pathway, namely porphyrin (heme) biosynthesis, is required for the cytotoxicity of Artemisinins. Genetic or pharmacological modulation of porphyrin production is sufficient to alter its cytotoxicity in eukaryotic cells. Using multiple model systems of human brain tumor development, such as cerebral glioblastoma organoids, and patient-derived tumor spheroids, we sensitize cancer cells to dihydroartemisinin using the clinically approved porphyrin enhancer and surgical fluorescence marker 5-aminolevulinic acid, 5-ALA. A combination treatment of Artemisinins and 5-ALA markedly and specifically killed brain tumor cells in all model systems tested, including orthotopic patient-derived xenografts in vivo. These data uncover the critical molecular pathway for Artemisinin cytotoxicity and a sensitization strategy to treat different brain tumors, including drug-resistant human glioblastomas.

ZFP462 safeguards neural lineage specification by targeting G9A/GLP-mediated heterochromatin to silence enhancers.

ZNF462 haploinsufficiency is linked to Weiss-Kruszka syndrome, a genetic disorder characterized by neurodevelopmental defects, including autism. Though conserved in vertebrates and essential for embryonic development, the molecular functions of ZNF462 remain unclear. We identified its murine homologue ZFP462 in a screen for mediators of epigenetic gene silencing. Here we show that ZFP462 safeguards neural lineage specification of mouse embryonic stem cells (ESCs) by targeting the H3K9-specific histone methyltransferase complex G9A/GLP to silence meso-endodermal genes. ZFP462 binds to transposable elements that are potential enhancers harbouring pluripotency and meso-endoderm transcription factor binding sites. Recruiting G9A/GLP, ZFP462 seeds heterochromatin, restricting transcription factor binding. Loss of ZFP462 in ESCs results in increased chromatin accessibility at target sites and ectopic expression of meso-endodermal genes. Taken together, ZFP462 confers lineage and locus specificity to the broadly expressed epigenetic regulator G9A/GLP. Our results suggest that aberrant activation of lineage non-specific genes in the neuronal lineage underlies ZNF462-associated neurodevelopmental pathology.

The HUSH complex controls brain architecture and protocadherin fidelity.

The HUSH (human silencing hub) complex contains the H3K9me3 binding protein M-phase phosphoprotein 8 (MPP8) and recruits the histone methyltransferase SETDB1 as well as Microrchidia CW-type zinc finger protein 2 (MORC2). Functional and mechanistic studies of the HUSH complex have hitherto been centered around SETDB1 while the in vivo functions of MPP8 and MORC2 remain elusive. Here, we show that genetic inactivation of Mphosph8 or Morc2a in the nervous system of mice leads to increased brain size, altered brain architecture, and behavioral changes. Mechanistically, in both mouse brains and human cerebral organoids, MPP8 and MORC2 suppress the repetitive-like protocadherin gene cluster in an H3K9me3-dependent manner. Our data identify MPP8 and MORC2, previously linked to silencing of repetitive elements via the HUSH complex, as key epigenetic regulators of protocadherin expression in the nervous system and thereby brain development and neuronal individuality in mice and humans.

DNA sequence and chromatin modifiers cooperate to confer epigenetic bistability at imprinting control regions.

Genomic imprinting is regulated by parental-specific DNA methylation of imprinting control regions (ICRs). Despite an identical DNA sequence, ICRs can exist in two distinct epigenetic states that are memorized throughout unlimited cell divisions and reset during germline formation. Here, we systematically study the genetic and epigenetic determinants of this epigenetic bistability. By iterative integration of ICRs and related DNA sequences to an ectopic location in the mouse genome, we first identify the DNA sequence features required for maintenance of epigenetic states in embryonic stem cells. The autonomous regulatory properties of ICRs further enabled us to create DNA-methylation-sensitive reporters and to screen for key components involved in regulating their epigenetic memory. Besides DNMT1, UHRF1 and ZFP57, we identify factors that prevent switching from methylated to unmethylated states and show that two of these candidates, ATF7IP and ZMYM2, are important for the stability of DNA and H3K9 methylation at ICRs in embryonic stem cells.

Transient upregulation of IRF1 during exit from naive pluripotency confers viral protection.

Stem cells intrinsically express a subset of genes which are normally associated with interferon stimulation and the innate immune response. However, the expression of these interferon-stimulated genes (ISG) in stem cells is independent from external stimuli such as viral infection. Here, we show that the interferon regulatory factor 1, Irf1, is directly controlled by the murine formative pluripotency gene regulatory network and transiently upregulated during the transition from naive to formative pluripotency. IRF1 binds to regulatory regions of a conserved set of ISGs and is required for their faithful expression upon exit from naive pluripotency. We show that in the absence of IRF1, cells exiting the naive pluripotent stem cell state are more susceptible to viral infection. Irf1 therefore acts as a link between the formative pluripotency network, regulation of innate immunity genes, and defense against viral infections during formative pluripotency.

Redirecting Imipramine against Bluetongue Virus Infection: Insights from a Genome-wide Haploid Screening Study.

Bluetongue virus (BTV), an arbovirus of ruminants, is a causative agent of numerous epidemics around the world. Due to the emergence of novel reassortant BTV strains and new outbreaks, there is an unmet need for efficacious antivirals. In this study, we used an improved haploid screening platform to identify the relevant host factors for BTV infection. Our screening tool identified and validated the host factor Niemann-Pick C1 (NPC1), a lysosomal membrane protein that is involved in lysosomal cholesterol transport, as a critical factor in BTV infection. This finding prompted us to investigate the possibility of testing imipramine, an antidepressant drug known to inhibit NPC1 function by interfering with intracellular cholesterol trafficking. In this study, we evaluated the sensitivity of BTV to imipramine using in vitro assays. Our results demonstrate that imipramine pretreatment inhibited in vitro replication and progeny release of BTV-4, BTV-8, and BTV-16. Collectively, our findings highlight the importance of NPC1 for BTV infection and recommend the reprofiling of imipramine as a potential antiviral drug against BTV.

GFPT2/GFAT2 and AMDHD2 act in tandem to control the hexosamine pathway.

The hexosamine biosynthetic pathway (HBP) produces the essential metabolite UDP-GlcNAc and plays a key role in metabolism, health, and aging. The HBP is controlled by its rate-limiting enzyme glutamine fructose-6-phosphate amidotransferase (GFPT/GFAT) that is directly inhibited by UDP-GlcNAc in a feedback loop. HBP regulation by GFPT is well studied but other HBP regulators have remained obscure. Elevated UDP-GlcNAc levels counteract the glycosylation toxin tunicamycin (TM), and thus we screened for TM resistance in haploid mouse embryonic stem cells (mESCs) using random chemical mutagenesis to determine alternative HBP regulation. We identified the N-acetylglucosamine deacetylase AMDHD2 that catalyzes a reverse reaction in the HBP and its loss strongly elevated UDP-GlcNAc. To better understand AMDHD2, we solved the crystal structure and found that loss-of-function (LOF) is caused by protein destabilization or interference with its catalytic activity. Finally, we show that mESCs express AMDHD2 together with GFPT2 instead of the more common paralog GFPT1. Compared with GFPT1, GFPT2 had a much lower sensitivity to UDP-GlcNAc inhibition, explaining how AMDHD2 LOF resulted in HBP activation. This HBP configuration in which AMDHD2 serves to balance GFPT2 activity was also observed in other mESCs and, consistently, the GFPT2:GFPT1 ratio decreased with differentiation of human embryonic stem cells. Taken together, our data reveal a critical function of AMDHD2 in limiting UDP-GlcNAc production in cells that use GFPT2 for metabolite entry into the HBP.

The BTB-domain transcription factor ZBTB2 recruits chromatin remodelers and a histone chaperone during the exit from pluripotency.

Transcription factors (TFs) harboring broad-complex, tramtrack, and bric-a-brac (BTB) domains play important roles in development and disease. These BTB domains are thought to recruit transcriptional modulators to target DNA regions. However, a systematic molecular understanding of the mechanism of action of this TF family is lacking. Here, we identify the zinc finger BTB-TF Zbtb2 from a genetic screen for regulators of exit from pluripotency and demonstrate that its absence perturbs embryonic stem cell differentiation and the gene expression dynamics underlying peri-implantation development. We show that ZBTB2 binds the chromatin remodeler Ep400 to mediate downstream transcription. Independently, the BTB domain directly interacts with nucleosome remodeling and deacetylase and histone chaperone histone regulator A. Nucleosome remodeling and deacetylase recruitment is a common feature of BTB TFs, and based on phylogenetic analysis, we propose that this is a conserved evolutionary property. Binding to UBN2, in contrast, is specific to ZBTB2 and requires a C-terminal extension of the BTB domain. Taken together, this study identifies a BTB-domain TF that recruits chromatin modifiers and a histone chaperone during a developmental cell state transition and defines unique and shared molecular functions of the BTB-domain TF family.

Multiplexed detection of SARS-CoV-2 and other respiratory infections in high throughput by SARSeq.

The COVID-19 pandemic has demonstrated the need for massively-parallel, cost-effective tests monitoring viral spread. Here we present SARSeq, saliva analysis by RNA sequencing, a method to detect SARS-CoV-2 and other respiratory viruses on tens of thousands of samples in parallel. SARSeq relies on next generation sequencing of multiple amplicons generated in a multiplexed RT-PCR reaction. Two-dimensional, unique dual indexing, using four indices per sample, enables unambiguous and scalable assignment of reads to individual samples. We calibrate SARSeq on SARS-CoV-2 synthetic RNA, virions, and hundreds of human samples of various types. Robustness and sensitivity were virtually identical to quantitative RT-PCR. Double-blinded benchmarking to gold standard quantitative-RT-PCR performed by human diagnostics laboratories confirms this high sensitivity. SARSeq can be used to detect Influenza A and B viruses and human rhinovirus in parallel, and can be expanded for detection of other pathogens. Thus, SARSeq is ideally suited for differential diagnostic of infections during a pandemic.

A human tissue screen identifies a regulator of ER secretion as a brain-size determinant.

Loss-of-function (LOF) screens provide a powerful approach to identify regulators in biological processes. Pioneered in laboratory animals, LOF screens of human genes are currently restricted to two-dimensional cell cultures, which hinders the testing of gene functions requiring tissue context. Here, we present CRISPR-lineage tracing at cellular resolution in heterogeneous tissue (CRISPR-LICHT), which enables parallel LOF studies in human cerebral organoid tissue. We used CRISPR-LICHT to test 173 microcephaly candidate genes, revealing 25 to be involved in known and uncharacterized microcephaly-associated pathways. We characterized IER3IP1, which regulates the endoplasmic reticulum (ER) function and extracellular matrix protein secretion crucial for tissue integrity, the dysregulation of which results in microcephaly. Our human tissue screening technology identifies microcephaly genes and mechanisms involved in brain-size control.

Multilayered VBC score predicts sgRNAs that efficiently generate loss-of-function alleles.

CRISPR-Cas9 screens have emerged as a transformative approach to systematically probe gene functions. The quality and success of these screens depends on the frequencies of loss-of-function alleles, particularly in negative-selection screens widely applied for probing essential genes. Using optimized screening workflows, we performed essentialome screens in cancer cell lines and embryonic stem cells and achieved dropout efficiencies that could not be explained by common frameshift frequencies. We find that these superior effect sizes are mainly determined by the impact of in-frame mutations on protein function, which can be predicted based on amino acid composition and conservation. We integrate protein features into a 'Bioscore' and fuse it with improved predictors of single-guide RNA activity and indel formation to establish a score that captures all relevant processes in CRISPR-Cas9 mutagenesis. This Vienna Bioactivity CRISPR score (www.vbc-score.org) outperforms previous prediction tools and enables the selection of sgRNAs that effectively produce loss-of-function alleles.