N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2.

The SARS-CoV-2 spike (S) glycoprotein contains an immunodominant receptor-binding domain (RBD) targeted by most neutralizing antibodies (Abs) in COVID-19 patient plasma. Little is known about neutralizing Abs binding to epitopes outside the RBD and their contribution to protection. Here, we describe 41 human monoclonal Abs (mAbs) derived from memory B cells, which recognize the SARS-CoV-2 S N-terminal domain (NTD) and show that a subset of them neutralize SARS-CoV-2 ultrapotently. We define an antigenic map of the SARS-CoV-2 NTD and identify a supersite (designated site i) recognized by all known NTD-specific neutralizing mAbs. These mAbs inhibit cell-to-cell fusion, activate effector functions, and protect Syrian hamsters from SARS-CoV-2 challenge, albeit selecting escape mutants in some animals. Indeed, several SARS-CoV-2 variants, including the B.1.1.7, B.1.351, and P.1 lineages, harbor frequent mutations within the NTD supersite, suggesting ongoing selective pressure and the importance of NTD-specific neutralizing mAbs for protective immunity and vaccine design.

Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity.

SARS-CoV-2 can mutate and evade immunity, with consequences for efficacy of emerging vaccines and antibody therapeutics. Here, we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is a highly variable region of S and provide epidemiological, clinical, and molecular characterization of a prevalent, sentinel RBM mutation, N439K. We demonstrate N439K S protein has enhanced binding affinity to the hACE2 receptor, and N439K viruses have similar in vitro replication fitness and cause infections with similar clinical outcomes as compared to wild type. We show the N439K mutation confers resistance against several neutralizing monoclonal antibodies, including one authorized for emergency use by the US Food and Drug Administration (FDA), and reduces the activity of some polyclonal sera from persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

Native elongating transcript sequencing reveals human transcriptional activity at nucleotide resolution.

Major features of transcription by human RNA polymerase II (Pol II) remain poorly defined due to a lack of quantitative approaches for visualizing Pol II progress at nucleotide resolution. We developed a simple and powerful approach for performing native elongating transcript sequencing (NET-seq) in human cells that globally maps strand-specific Pol II density at nucleotide resolution. NET-seq exposes a mode of antisense transcription that originates downstream and converges on transcription from the canonical promoter. Convergent transcription is associated with a distinctive chromatin configuration and is characteristic of lower-expressed genes. Integration of NET-seq with genomic footprinting data reveals stereotypic Pol II pausing coincident with transcription factor occupancy. Finally, exons retained in mature transcripts display Pol II pausing signatures that differ markedly from skipped exons, indicating an intrinsic capacity for Pol II to recognize exons with different processing fates. Together, human NET-seq exposes the topography and regulatory complexity of human gene expression.

A pan-influenza antibody inhibiting neuraminidase via receptor mimicry.

Rapidly evolving influenza A viruses (IAVs) and influenza B viruses (IBVs) are major causes of recurrent lower respiratory tract infections. Current influenza vaccines elicit antibodies predominantly to the highly variable head region of haemagglutinin and their effectiveness is limited by viral drift1 and suboptimal immune responses2. Here we describe a neuraminidase-targeting monoclonal antibody, FNI9, that potently inhibits the enzymatic activity of all group 1 and group 2 IAVs, as well as Victoria/2/87-like, Yamagata/16/88-like and ancestral IBVs. FNI9 broadly neutralizes seasonal IAVs and IBVs, including the immune-evading H3N2 strains bearing an N-glycan at position 245, and shows synergistic activity when combined with anti-haemagglutinin stem-directed antibodies. Structural analysis reveals that D107 in the FNI9 heavy chain complementarity-determinant region 3 mimics the interaction of the sialic acid carboxyl group with the three highly conserved arginine residues (R118, R292 and R371) of the neuraminidase catalytic site. FNI9 demonstrates potent prophylactic activity against lethal IAV and IBV infections in mice. The unprecedented breadth and potency of the FNI9 monoclonal antibody supports its development for the prevention of influenza illness by seasonal and pandemic viruses.

Broadly neutralizing antibodies overcome SARS-CoV-2 Omicron antigenic shift.

The recently emerged SARS-CoV-2 Omicron variant encodes 37 amino acid substitutions in the spike protein, 15 of which are in the receptor-binding domain (RBD), thereby raising concerns about the effectiveness of available vaccines and antibody-based therapeutics. Here we show that the Omicron RBD binds to human ACE2 with enhanced affinity, relative to the Wuhan-Hu-1 RBD, and binds to mouse ACE2. Marked reductions in neutralizing activity were observed against Omicron compared to the ancestral pseudovirus in plasma from convalescent individuals and from individuals who had been vaccinated against SARS-CoV-2, but this loss was less pronounced after a third dose of vaccine. Most monoclonal antibodies that are directed against the receptor-binding motif lost in vitro neutralizing activity against Omicron, with only 3 out of 29 monoclonal antibodies retaining unaltered potency, including the ACE2-mimicking S2K146 antibody1. Furthermore, a fraction of broadly neutralizing sarbecovirus monoclonal antibodies neutralized Omicron through recognition of antigenic sites outside the receptor-binding motif, including sotrovimab2, S2X2593 and S2H974. The magnitude of Omicron-mediated immune evasion marks a major antigenic shift in SARS-CoV-2. Broadly neutralizing monoclonal antibodies that recognize RBD epitopes that are conserved among SARS-CoV-2 variants and other sarbecoviruses may prove key to controlling the ongoing pandemic and future zoonotic spillovers.

Broad sarbecovirus neutralization by a human monoclonal antibody.

The recent emergence of SARS-CoV-2 variants of concern1-10 and the recurrent spillovers of coronaviruses11,12 into the human population highlight the need for broadly neutralizing antibodies that are not affected by the ongoing antigenic drift and that can prevent or treat future zoonotic infections. Here we describe a human monoclonal antibody designated S2X259, which recognizes a highly conserved cryptic epitope of the receptor-binding domain and cross-reacts with spikes from all clades of sarbecovirus. S2X259 broadly neutralizes spike-mediated cell entry of SARS-CoV-2, including variants of concern (B.1.1.7, B.1.351, P.1, and B.1.427/B.1.429), as well as a wide spectrum of human and potentially zoonotic sarbecoviruses through inhibition of angiotensin-converting enzyme 2 (ACE2) binding to the receptor-binding domain. Furthermore, deep-mutational scanning and in vitro escape selection experiments demonstrate that S2X259 possesses an escape profile that is limited to a single substitution, G504D. We show that prophylactic and therapeutic administration of S2X259 protects Syrian hamsters (Mesocricetus auratus) against challenge with the prototypic SARS-CoV-2 and the B.1.351 variant of concern, which suggests that this monoclonal antibody is a promising candidate for the prevention and treatment of emergent variants and zoonotic infections. Our data reveal a key antigenic site that is targeted by broadly neutralizing antibodies and will guide the design of vaccines that are effective against all sarbecoviruses.

SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape.

An ideal therapeutic anti-SARS-CoV-2 antibody would resist viral escape1-3, have activity against diverse sarbecoviruses4-7, and be highly protective through viral neutralization8-11 and effector functions12,13. Understanding how these properties relate to each other and vary across epitopes would aid the development of therapeutic antibodies and guide vaccine design. Here we comprehensively characterize escape, breadth and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD). Despite a trade-off between in vitro neutralization potency and breadth of sarbecovirus binding, we identify neutralizing antibodies with exceptional sarbecovirus breadth and a corresponding resistance to SARS-CoV-2 escape. One of these antibodies, S2H97, binds with high affinity across all sarbecovirus clades to a cryptic epitope and prophylactically protects hamsters from viral challenge. Antibodies that target the angiotensin-converting enzyme 2 (ACE2) receptor-binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency. Nevertheless, we also characterize a potent RBM antibody (S2E128) with breadth across sarbecoviruses related to SARS-CoV-2 and a high barrier to viral escape. These data highlight principles underlying variation in escape, breadth and potency among antibodies that target the RBD, and identify epitopes and features to prioritize for therapeutic development against the current and potential future pandemics.

Transcription Factor-Directed Re-wiring of Chromatin Architecture for Somatic Cell Nuclear Reprogramming toward trans-Differentiation.

MYOD-directed fibroblast trans-differentiation into skeletal muscle provides a unique model to investigate how one transcription factor (TF) reconfigures the three-dimensional chromatin architecture to control gene expression, which is otherwise achieved by the combinatorial activities of multiple TFs. Integrative analysis of genome-wide high-resolution chromatin interactions, MYOD and CTCF DNA-binding profile, and gene expression, revealed that MYOD directs extensive re-wiring of interactions involving cis-regulatory and structural genomic elements, including promoters, enhancers, and insulated neighborhoods (INs). Re-configured INs were hot-spots of differential interactions, whereby MYOD binding to highly constrained sequences at IN boundaries and/or inside INs led to alterations of promoter-enhancer interactions to repress cell-of-origin genes and to activate muscle-specific genes. Functional evidence shows that MYOD-directed re-configuration of chromatin interactions temporally preceded the effect on gene expression and was mediated by direct MYOD-DNA binding. These data illustrate a model whereby a single TF alters multi-loop hubs to drive somatic cell trans-differentiation.

BET Bromodomain Proteins Function as Master Transcription Elongation Factors Independent of CDK9 Recruitment.

Processive elongation of RNA Polymerase II from a proximal promoter paused state is a rate-limiting event in human gene control. A small number of regulatory factors influence transcription elongation on a global scale. Prior research using small-molecule BET bromodomain inhibitors, such as JQ1, linked BRD4 to context-specific elongation at a limited number of genes associated with massive enhancer regions. Here, the mechanistic characterization of an optimized chemical degrader of BET bromodomain proteins, dBET6, led to the unexpected identification of BET proteins as master regulators of global transcription elongation. In contrast to the selective effect of bromodomain inhibition on transcription, BET degradation prompts a collapse of global elongation that phenocopies CDK9 inhibition. Notably, BRD4 loss does not directly affect CDK9 localization. These studies, performed in translational models of T cell leukemia, establish a mechanism-based rationale for the development of BET bromodomain degradation as cancer therapy.

Human gene essentiality.

A gene can be defined as essential when loss of its function compromises viability of the individual (for example, embryonic lethality) or results in profound loss of fitness. At the population level, identification of essential genes is accomplished by observing intolerance to loss-of-function variants. Several computational methods are available to score gene essentiality, and recent progress has been made in defining essentiality in the non-coding genome. Haploinsufficiency is emerging as a critical aspect of gene essentiality: approximately 3,000 human genes cannot tolerate loss of one of the two alleles. Genes identified as essential in human cell lines or knockout mice may be distinct from those in living humans. Reconciling these discrepancies in how we evaluate gene essentiality has applications in clinical genetics and may offer insights for drug development.

Transfer transcriptomic signatures for infectious diseases.

The modulation of the transcriptome is among the earliest responses to infection. However, defining the transcriptomic signatures of disease is challenging because logistic, technical, and cost factors limit the size and representativeness of samples in clinical studies. These limitations lead to a poor performance of signatures when applied to new datasets. Although the study focuses on infection, the central hypothesis of the work is the generalization of sets of signatures across diseases. We use a machine learning approach to identify common elements in datasets and then test empirically whether they are informative about a second dataset from a disease or process distinct from the original dataset. We identify sets of genes, which we name transfer signatures, that are predictive across diverse datasets and/or species (e.g., rhesus to humans). We demonstrate the usefulness of transfer signatures in two use cases: the progression of latent to active tuberculosis and the severity of COVID-19 and influenza A H1N1 infection. This indicates that transfer signatures can be deployed in settings that lack disease-specific biomarkers. The broad significance of our work lies in the concept that a small set of archetypal human immunophenotypes, captured by transfer signatures, can explain a larger set of responses to diverse diseases.

Functional characterization of 3D protein structures informed by human genetic diversity.

Sequence variation data of the human proteome can be used to analyze 3D protein structures to derive functional insights. We used genetic variant data from nearly 140,000 individuals to analyze 3D positional conservation in 4,715 proteins and 3,951 homology models using 860,292 missense and 465,886 synonymous variants. Sixty percent of protein structures harbor at least one intolerant 3D site as defined by significant depletion of observed over expected missense variation. Structural intolerance data correlated with deep mutational scanning functional readouts for PPARG, MAPK1/ERK2, UBE2I, SUMO1, PTEN, CALM1, CALM2, and TPK1 and with shallow mutagenesis data for 1,026 proteins. The 3D structural intolerance analysis revealed different features for ligand binding pockets and orthosteric and allosteric sites. Large-scale data on human genetic variation support a definition of functional 3D sites proteome-wide.

The yeast exoribonuclease Xrn1 and associated factors modulate RNA polymerase II processivity in 5' and 3' gene regions.

mRNA levels are determined by the balance between mRNA synthesis and decay. Protein factors that mediate both processes, including the 5'-3' exonuclease Xrn1, are responsible for a cross-talk between the two processes that buffers steady-state mRNA levels. However, the roles of these proteins in transcription remain elusive and controversial. Applying native elongating transcript sequencing (NET-seq) to yeast cells, we show that Xrn1 functions mainly as a transcriptional activator and that its disruption manifests as a reduction of RNA polymerase II (Pol II) occupancy downstream of transcription start sites. By combining our sequencing data and mathematical modeling of transcription, we found that Xrn1 modulates transcription initiation and elongation of its target genes. Furthermore, Pol II occupancy markedly increased near cleavage and polyadenylation sites in xrn1Δ cells, whereas its activity decreased, a characteristic feature of backtracked Pol II. We also provide indirect evidence that Xrn1 is involved in transcription termination downstream of polyadenylation sites. We noted that two additional decay factors, Dhh1 and Lsm1, seem to function similarly to Xrn1 in transcription, perhaps as a complex, and that the decay factors Ccr4 and Rpb4 also perturb transcription in other ways. Interestingly, the decay factors could differentiate between SAGA- and TFIID-dominated promoters. These two classes of genes responded differently to XRN1 deletion in mRNA synthesis and were differentially regulated by mRNA decay pathways, raising the possibility that one distinction between these two gene classes lies in the mechanisms that balance mRNA synthesis with mRNA decay.

Regulatory genome variants in human susceptibility to infection.

Genome studies have accelerated the discovery of common and rare genetic variants associated with susceptibility to infection and with disease severity. Genome-wide association studies identified many common genetic variants associated with modest risk for infection. Over 80% of these common variants map to the non-coding genome and are thought to modulate the regulatory networks. Exome sequencing has rapidly expanded the number of recognized primary immunodeficiencies through the identification of rare coding variants. In contrast, less than 29 primary immunodeficiencies have causative rare variation mapped outside protein-coding regions. In the future, whole genome sequencing will accelerate the identification of rare variants of substantial phenotypic impact that disrupt essential regulatory elements and the three-dimensional structure of chromatin.

Deep sequencing of 10,000 human genomes.

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

Dynamics of HIV latency and reactivation in a primary CD4+ T cell model.

HIV latency is a major obstacle to curing infection. Current strategies to eradicate HIV aim at increasing transcription of the latent provirus. In the present study we observed that latently infected CD4+ T cells from HIV-infected individuals failed to produce viral particles upon ex vivo exposure to SAHA (vorinostat), despite effective inhibition of histone deacetylases. To identify steps that were not susceptible to the action of SAHA or other latency reverting agents, we used a primary CD4+ T cell model, joint host and viral RNA sequencing, and a viral-encoded reporter. This model served to investigate the characteristics of latently infected cells, the dynamics of HIV latency, and the process of reactivation induced by various stimuli. During latency, we observed persistence of viral transcripts but only limited viral translation. Similarly, the reactivating agents SAHA and disulfiram successfully increased viral transcription, but failed to effectively enhance viral translation, mirroring the ex vivo data. This study highlights the importance of post-transcriptional blocks as one mechanism leading to HIV latency that needs to be relieved in order to purge the viral reservoir.

Bioinformatics and HIV latency.

Despite effective treatment, HIV is not completely eliminated from the infected organism because of the existence of viral reservoirs. A major reservoir consists of infected resting CD4+ T cells, mostly of memory type, that persist over time due to the stable proviral insertion and a long cellular lifespan. Resting cells do not produce viral particles and are protected from viral-induced cytotoxicity or immune killing. However, these latently infected cells can be reactivated by stochastic events or by external stimuli. The present review focuses on novel genome-wide technologies applied to the study of integration, transcriptome, and proteome characteristics and their recent contribution to the understanding of HIV latency.

Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows.

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.