Deconvoluting complex correlates of COVID-19 severity with a multi-omic pandemic tracking strategy.

The SARS-CoV-2 pandemic has differentially impacted populations across race and ethnicity. A multi-omic approach represents a powerful tool to examine risk across multi-ancestry genomes. We leverage a pandemic tracking strategy in which we sequence viral and host genomes and transcriptomes from nasopharyngeal swabs of 1049 individuals (736 SARS-CoV-2 positive and 313 SARS-CoV-2 negative) and integrate them with digital phenotypes from electronic health records from a diverse catchment area in Northern California. Genome-wide association disaggregated by admixture mapping reveals novel COVID-19-severity-associated regions containing previously reported markers of neurologic, pulmonary and viral disease susceptibility. Phylodynamic tracking of consensus viral genomes reveals no association with disease severity or inferred ancestry. Summary data from multiomic investigation reveals metagenomic and HLA associations with severe COVID-19. The wealth of data available from residual nasopharyngeal swabs in combination with clinical data abstracted automatically at scale highlights a powerful strategy for pandemic tracking, and reveals distinct epidemiologic, genetic, and biological associations for those at the highest risk.

Immune imprinting, breadth of variant recognition, and germinal center response in human SARS-CoV-2 infection and vaccination.

During the SARS-CoV-2 pandemic, novel and traditional vaccine strategies have been deployed globally. We investigated whether antibodies stimulated by mRNA vaccination (BNT162b2), including third-dose boosting, differ from those generated by infection or adenoviral (ChAdOx1-S and Gam-COVID-Vac) or inactivated viral (BBIBP-CorV) vaccines. We analyzed human lymph nodes after infection or mRNA vaccination for correlates of serological differences. Antibody breadth against viral variants is lower after infection compared with all vaccines evaluated but improves over several months. Viral variant infection elicits variant-specific antibodies, but prior mRNA vaccination imprints serological responses toward Wuhan-Hu-1 rather than variant antigens. In contrast to disrupted germinal centers (GCs) in lymph nodes during infection, mRNA vaccination stimulates robust GCs containing vaccine mRNA and spike antigen up to 8 weeks postvaccination in some cases. SARS-CoV-2 antibody specificity, breadth, and maturation are affected by imprinting from exposure history and distinct histological and antigenic contexts in infection compared with vaccination.

mRNA vaccination compared to infection elicits an IgG-predominant response with greater SARS-CoV-2 specificity and similar decrease in variant spike recognition.

During the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, new vaccine strategies including lipid nanoparticle delivery of antigen encoding RNA have been deployed globally. The BioNTech/Pfizer mRNA vaccine BNT162b2 encoding SARS-CoV-2 spike protein shows 95% efficacy in preventing disease, but it is unclear how the antibody responses to vaccination differ from those generated by infection. Here we compare the magnitude and breadth of antibodies targeting SARS-CoV-2, SARS-CoV-2 variants of concern, and endemic coronaviruses, in vaccinees and infected patients. We find that vaccination differs from infection in the dominance of IgG over IgM and IgA responses, with IgG reaching levels similar to those of severely ill COVID-19 patients and shows decreased breadth of the antibody response targeting endemic coronaviruses. Viral variants of concern from B.1.1.7 to P.1 to B.1.351 form a remarkably consistent hierarchy of progressively decreasing antibody recognition by both vaccinees and infected patients exposed to Wuhan-Hu-1 antigens.

High-throughput SARS-CoV-2 and host genome sequencing from single nasopharyngeal swabs.

During COVID19 and other viral pandemics, rapid generation of host and pathogen genomic data is critical to tracking infection and informing therapies. There is an urgent need for efficient approaches to this data generation at scale. We have developed a scalable, high throughput approach to generate high fidelity low pass whole genome and HLA sequencing, viral genomes, and representation of human transcriptome from single nasopharyngeal swabs of COVID19 patients.

An Extensive Meta-Metagenomic Search Identifies SARS-CoV-2-Homologous Sequences in Pangolin Lung Viromes.

In numerous instances, tracking the biological significance of a nucleic acid sequence can be augmented through the identification of environmental niches in which the sequence of interest is present. Many metagenomic data sets are now available, with deep sequencing of samples from diverse biological niches. While any individual metagenomic data set can be readily queried using web-based tools, meta-searches through all such data sets are less accessible. In this brief communication, we demonstrate such a meta-metagenomic approach, examining close matches to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in all high-throughput sequencing data sets in the NCBI Sequence Read Archive accessible with the "virome" keyword. In addition to the homology to bat coronaviruses observed in descriptions of the SARS-CoV-2 sequence (F. Wu, S. Zhao, B. Yu, Y. M. Chen, et al., Nature 579:265-269, 2020, https://doi.org/10.1038/s41586-020-2008-3; P. Zhou, X. L. Yang, X. G. Wang, B. Hu, et al., Nature 579:270-273, 2020, https://doi.org/10.1038/s41586-020-2012-7), we note a strong homology to numerous sequence reads in metavirome data sets generated from the lungs of deceased pangolins reported by Liu et al. (P. Liu, W. Chen, and J. P. Chen, Viruses 11:979, 2019, https://doi.org/10.3390/v11110979). While analysis of these reads indicates the presence of a similar viral sequence in pangolin lung, the similarity is not sufficient to either confirm or rule out a role for pangolins as an intermediate host in the recent emergence of SARS-CoV-2. In addition to the implications for SARS-CoV-2 emergence, this study illustrates the utility and limitations of meta-metagenomic search tools in effective and rapid characterization of potentially significant nucleic acid sequences.IMPORTANCE Meta-metagenomic searches allow for high-speed, low-cost identification of potentially significant biological niches for sequences of interest.

Loop-closure kinetics reveal a stable, right-handed DNA intermediate in Cre recombination.

In Cre site-specific recombination, the synaptic intermediate is a recombinase homotetramer containing a pair of loxP DNA target sites. The enzyme system's strand-exchange mechanism proceeds via a Holliday-junction (HJ) intermediate; however, the geometry of DNA segments in the synapse has remained highly controversial. In particular, all crystallographic structures are consistent with an achiral, planar Holliday-junction (HJ) structure, whereas topological assays based on Cre-mediated knotting of plasmid DNAs are consistent with a right-handed chiral junction. We use the kinetics of loop closure involving closely spaced (131-151 bp) loxP sites to investigate the in-aqueo ensemble of conformations for the longest-lived looped DNA intermediate. Fitting the experimental site-spacing dependence of the loop-closure probability, J, to a statistical-mechanical theory of DNA looping provides evidence for substantial out-of-plane HJ distortion, which unequivocally stands in contrast to the square-planar intermediate geometry from Cre-loxP crystal structures and those of other int-superfamily recombinases. J measurements for an HJ-isomerization-deficient Cre mutant suggest that the apparent geometry of the wild-type complex is consistent with temporal averaging of right-handed and achiral structures. Our approach connects the static pictures provided by crystal structures and the natural dynamics of macromolecules in solution, thus advancing a more comprehensive dynamic analysis of large nucleoprotein structures and their mechanisms.

Deconvolution of nucleic-acid length distributions: a gel electrophoresis analysis tool and applications.

Next-generation DNA-sequencing (NGS) technologies, which are designed to streamline the acquisition of massive amounts of sequencing data, are nonetheless dependent on various preparative steps to generate DNA fragments of required concentration, purity and average size (molecular weight). Current automated electrophoresis systems for DNA- and RNA-sample quality control, such as Agilent's Bioanalyzer® and TapeStation® products, are costly to acquire and use; they also provide limited information for samples having broad size distributions. Here, we describe a software tool that helps determine the size distribution of DNA fragments in an NGS library, or other DNA sample, based on gel-electrophoretic line profiles. The software, developed as an ImageJ plug-in, allows for straightforward processing of gel images, including lane selection and fitting of univariate functions to intensity distributions. The user selects the option of fitting either discrete profiles in cases where discrete gel bands are visible or continuous profiles, having multiple bands buried under a single broad peak. The method requires only modest imaging capabilities and is a cost-effective, rigorous alternative characterization method to augment existing techniques for library quality control.

Recompleting the Caenorhabditis elegans genome.

Caenorhabditis elegans was the first multicellular eukaryotic genome sequenced to apparent completion. Although this assembly employed a standard C. elegans strain (N2), it used sequence data from several laboratories, with DNA propagated in bacteria and yeast. Thus, the N2 assembly has many differences from any C. elegans available today. To provide a more accurate C. elegans genome, we performed long-read assembly of VC2010, a modern strain derived from N2. Our VC2010 assembly has 99.98% identity to N2 but with an additional 1.8 Mb including tandem repeat expansions and genome duplications. For 116 structural discrepancies between N2 and VC2010, 97 structures matching VC2010 (84%) were also found in two outgroup strains, implying deficiencies in N2. Over 98% of N2 genes encoded unchanged products in VC2010; moreover, we predicted ≥53 new genes in VC2010. The recompleted genome of C. elegans should be a valuable resource for genetics, genomics, and systems biology.

Intricate and Cell Type-Specific Populations of Endogenous Circular DNA (eccDNA) in Caenorhabditis elegans and Homo sapiens.

Investigations aimed at defining the 3D configuration of eukaryotic chromosomes have consistently encountered an endogenous population of chromosome-derived circular genomic DNA, referred to as extrachromosomal circular DNA (eccDNA). While the production, distribution, and activities of eccDNAs remain understudied, eccDNA formation from specific regions of the linear genome has profound consequences on the regulatory and coding capabilities for these regions. Here, we define eccDNA distributions in Caenorhabditis elegans and in three human cell types, utilizing a set of DNA topology-dependent approaches for enrichment and characterization. The use of parallel biophysical, enzymatic, and informatic approaches provides a comprehensive profiling of eccDNA robust to isolation and analysis methodology. Results in human and nematode systems provide quantitative analysis of the eccDNA loci at both unique and repetitive regions. Our studies converge on and support a consistent picture, in which endogenous genomic DNA circles are present in normal physiological states, and in which the circles come from both coding and noncoding genomic regions. Prominent among the coding regions generating DNA circles are several genes known to produce a diversity of protein isoforms, with mucin proteins and titin as specific examples.

Contribution of fluorophore dynamics and solvation to resonant energy transfer in protein-DNA complexes: a molecular-dynamics study.

In Förster resonance energy transfer (FRET) experiments, extracting accurate structural information about macromolecules depends on knowing the positions and orientations of donor and acceptor fluorophores. Several approaches have been employed to reduce uncertainties in quantitative FRET distance measurements. Fluorophore-position distributions can be estimated by surface accessibility (SA) calculations, which compute the region of space explored by the fluorophore within a static macromolecular structure. However, SA models generally do not take fluorophore shape, dye transition-moment orientation, or dye-specific chemical interactions into account. We present a detailed molecular-dynamics (MD) treatment of fluorophore dynamics for an ATTO donor/acceptor dye pair and specifically consider as case studies dye-labeled protein-DNA intermediates in Cre site-specific recombination. We carried out MD simulations in both an aqueous solution and glycerol/water mixtures to assess the effects of experimental solvent systems on dye dynamics. Our results unequivocally show that MD simulations capture solvent effects and dye-dye interactions that can dramatically affect energy transfer efficiency. We also show that results from SA models and MD simulations strongly diverge in cases where donor and acceptor fluorophores are in close proximity. Although atomistic simulations are computationally more expensive than SA models, explicit MD studies are likely to give more realistic results in both homogeneous and mixed solvents. Our study underscores the model-dependent nature of FRET analyses, but also provides a starting point to develop more realistic in silico approaches for obtaining experimental ensemble and single-molecule FRET data.

The thermodynamics of DNA loop formation, from J to Z.

The formation of DNA loops is a ubiquitous theme in biological processes, including DNA replication, recombination and repair, and gene regulation. These loops are mediated by proteins bound at specific sites along the contour of a single DNA molecule, in some cases many thousands of base pairs apart. Loop formation incurs a thermodynamic cost that is a sensitive function of the length of looped DNA as well as the geometry and elastic properties of the DNA-bound protein. The free energy of DNA looping is logarithmically related to a generalization of the Jacobson-Stockmayer factor for DNA cyclization, termed the J factor. In the present article, we review the thermodynamic origins of this quantity, discuss how it is measured experimentally and connect the macroscopic interpretation of the J factor with a statistical-mechanical description of DNA looping and cyclization.

Measurements of DNA-loop formation via Cre-mediated recombination.

The Cre-recombination system has become an important tool for genetic manipulation of higher organisms and a model for site-specific DNA-recombination mechanisms employed by the λ-Int superfamily of recombinases. We report a novel quantitative approach for characterizing the probability of DNA-loop formation in solution using time-dependent ensemble Förster resonance energy transfer measurements of intra- and inter-molecular Cre-recombination kinetics. Our method uses an innovative technique for incorporating multiple covalent modifications at specific sites in covalently closed DNA. Because the mechanism of Cre recombinase does not conform to a simple kinetic scheme, we employ numerical methods to extract rate constants for fundamental steps that pertain to Cre-mediated loop closure. Cre recombination does not require accessory proteins, DNA supercoiling or particular metal-ion cofactors and is thus a highly flexible system for quantitatively analyzing DNA-loop formation in vitro and in vivo.