In vivo antiviral host transcriptional response to SARS-CoV-2 by viral load, sex, and age.

Despite limited genomic diversity, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has shown a wide range of clinical manifestations in different patient populations. The mechanisms behind these host differences are still unclear. Here, we examined host response gene expression across infection status, viral load, age, and sex among shotgun RNA sequencing profiles of nasopharyngeal (NP) swabs from 430 individuals with PCR-confirmed SARS-CoV-2 and 54 negative controls. SARS-CoV-2 induced a strong antiviral response with up-regulation of antiviral factors such as OAS1-3 and IFIT1-3 and T helper type 1 (Th1) chemokines CXCL9/10/11, as well as a reduction in transcription of ribosomal proteins. SARS-CoV-2 culture in human airway epithelial (HAE) cultures replicated the in vivo antiviral host response 7 days post infection, with no induction of interferon-stimulated genes after 3 days. Patient-matched longitudinal specimens (mean elapsed time = 6.3 days) demonstrated reduction in interferon-induced transcription, recovery of transcription of ribosomal proteins, and initiation of wound healing and humoral immune responses. Expression of interferon-responsive genes, including ACE2, increased as a function of viral load, while transcripts for B cell-specific proteins and neutrophil chemokines were elevated in patients with lower viral load. Older individuals had reduced expression of the Th1 chemokines CXCL9/10/11 and their cognate receptor CXCR3, as well as CD8A and granzyme B, suggesting deficiencies in trafficking and/or function of cytotoxic T cells and natural killer (NK) cells. Relative to females, males had reduced B cell-specific and NK cell-specific transcripts and an increase in inhibitors of nuclear factor kappa-B (NF-κB) signaling, possibly inappropriately throttling antiviral responses. Collectively, our data demonstrate that host responses to SARS-CoV-2 are dependent on viral load and infection time course, with observed differences due to age and sex that may contribute to disease severity.

Sensitive Identification of Bacterial DNA in Clinical Specimens by Broad-Range 16S rRNA Gene Enrichment.

The broad-range detection and identification of bacterial DNA from clinical specimens are a foundational approach in the practice of molecular microbiology. However, there are circumstances under which conventional testing may yield false-negative or otherwise uninterpretable results, including the presence of multiple bacterial templates or degraded nucleic acids. Here, we describe an alternative, next-generation sequencing approach for the broad range detection of bacterial DNA using broad-range 16S rRNA gene hybrid capture ("16S Capture"). The method is able to deconvolute multiple bacterial species present in a specimen, is compatible with highly fragmented templates, and can be readily implemented when the overwhelming majority of nucleic acids in a specimen derive from the human host. We find that this approach is sensitive to detecting as few as 17 Staphylococcus aureus genomes from a background of 100 ng of human DNA, providing 19- to 189-fold greater sensitivity for identifying bacterial sequences than standard shotgun metagenomic sequencing, and is able to successfully recover organisms from across the eubacterial tree of life. Application of 16S Capture to a proof-of-principle case series demonstrated its ability to identify bacterial species that were consistent with histological evidence of infection, even when diagnosis could not be established using conventional broad range bacterial detection assays. 16S Capture provides a novel means for the efficient and sensitive detection of bacteria embedded in human tissues and for specimens containing highly fragmented template DNA.

RNA Sequencing of the In Vivo Human Herpesvirus 6B Transcriptome To Identify Targets for Clinical Assays Distinguishing between Latent and Active Infections.

Human herpesvirus 6B (HHV-6B) DNA is frequently detected in human samples. Diagnostic assays distinguishing HHV-6B reactivation from latency are limited. This has impaired strategies to diagnose and treat HHV-6B-associated diseases. We used RNA sequencing to characterize and compare the HHV-6B transcriptome in multiple sample types, including (i) whole blood from hematopoietic cell transplant (HCT) recipients with and without HHV-6B plasma viremia, (ii) tumor tissue samples from subjects with large B cell lymphoma infected with HHV-6B, (iii) lymphoblastoid cell lines (LCLs) from subjects with inherited chromosomally integrated HHV-6B or latent infection with HHV-6B, and (iv) HHV-6B Z29 infected SupT1 CD4+ T cells. We demonstrated substantial overlap in the HHV-6B transcriptome observed in in vivo and in vitro samples, although there was variability in the breadth and quantity of gene expression across samples. The HHV-6B viral polymerase gene U38 was the only HHV-6B transcript detected in all next-generation RNA sequencing (RNA-seq) data sets and was one of the most highly expressed genes. We developed a novel reverse transcription-PCR assay targeting HHV-6B U38, which identified U38 mRNA in all tested whole-blood samples from patients with concurrent HHV-6B viremia. No HHV-6B U38 transcripts were detected by RNA-seq or reverse transcription-real-time quantitative PCR (RT-qPCR) in whole-blood samples from subjects without HHV-6B plasma detection or from latently infected LCLs. A RT-qPCR assay for HHV-6B U38 may be useful to identify lytic HHV-6B infection in nonplasma samples and samples from individuals with inherited chromosomally integrated HHV-6B. This study also demonstrates the feasibility of transcriptomic analyses for HCT recipients.IMPORTANCE Human herpesvirus 6B (HHV-6B) is a DNA virus that infects most children within the first few years of life. After primary infection, HHV-6B persists as a chronic, latent infection in many cell types. Additionally, HHV-6B can integrate into germ line chromosomes, resulting in individuals with viral DNA in every nucleated cell. Given that PCR to detect viral DNA is the mainstay for diagnosing HHV-6B infection, the characteristics of HHV-6B infection complicate efforts to distinguish between latent and active viral infection, particularly in immunocompromised patients who have frequent HHV-6B reactivation. In this study, we used RNA sequencing to characterize the HHV-6B gene expression profile in multiple sample types, and our findings identified evidence-based targets for diagnostic tests that distinguish between latent and active viral infection.

Inherited Chromosomally Integrated Human Herpesvirus 6 Demonstrates Tissue-Specific RNA Expression In Vivo That Correlates with an Increased Antibody Immune Response.

Human herpesviruses 6A and 6B (HHV-6A and HHV-6B) are human viruses capable of chromosomal integration. Approximately 1% of the human population carries one copy of HHV-6A/B integrated into every cell in their body, referred to as inherited chromosomally integrated human herpesvirus 6A/B (iciHHV-6A/B). Whether iciHHV-6A/B is transcriptionally active in vivo and how it shapes the immunological response are still unclear. In this study, we screened DNA sequencing (DNA-seq) and transcriptome sequencing (RNA-seq) data for 650 individuals available through the Genotype-Tissue Expression (GTEx) project and identified 2 iciHHV-6A- and 4 iciHHV-6B-positive candidates. When corresponding tissue-specific gene expression signatures were analyzed, low levels HHV-6A/B gene expression was found across multiple tissues, with the highest levels of gene expression in the brain (specifically for HHV-6A), testis, esophagus, and adrenal gland. U90 and U100 were the most highly expressed HHV-6 genes in both iciHHV-6A- and iciHHV-6B-positive individuals. To assess whether tissue-specific gene expression from iciHHV-6A/B influences the immune response, a cohort of 15,498 subjects was screened and 85 iciHHV-6A/B+ subjects were identified. Plasma samples from iciHHV-6A/B+ and age- and sex-matched controls were analyzed for antibodies to control antigens (cytomegalovirus [CMV], Epstein-Barr virus [EBV], and influenza virus [FLU]) or HHV-6A/B antigens. Our results indicate that iciHHV-6A/B+ subjects have significantly more antibodies against the U90 gene product (IE1) than do non-iciHHV-6-positive individuals. Antibody responses against EBV and FLU antigens or HHV-6A/B gene products either not expressed or expressed at low levels, such as U47, U57, and U72, were identical between controls and iciHHV-6A/B+ subjects. CMV-seropositive individuals with iciHHV-6A/B+ have more antibodies against CMV pp150 than do CMV-seropositive controls. These results argue that spontaneous gene expression from integrated HHV-6A/B leads to an increase in antigenic burden that translates into a more robust HHV-6A/B-specific antibody response.IMPORTANCE HHV-6A and -6B are human herpesviruses that have the unique property of being able to integrate into the telomeric regions of human chromosomes. Approximately 1% of the world's population carries integrated HHV-6A/B genome in every cell of their body. Whether viral genes are transcriptionally active in these individuals is unclear. By taking advantage of a unique tissue-specific gene expression data set, we showed that the majority of tissues from iciHHV-6 individuals do not show HHV-6 gene expression. Brain and testes showed the highest tissue-specific expression of HHV-6 genes in two separate data sets. Two HHV-6 genes, U90 (immediate early 1 protein) and U100 (glycoproteins Q1 and Q2), were found to be selectively and consistently expressed across several human tissues. Expression of U90 translates into an increase in antigen-specific antibody response in iciHHV-6A/B+ subjects relative to controls. Future studies will be needed to determine the mechanism of gene expression, the effects of these genes on human gene transcription networks, and the pathophysiological impact of having increased viral protein expression in tissue in conjunction with increased antigen-specific antibody production.

Metagenomic Analysis Reveals Clinical SARS-CoV-2 Infection and Bacterial or Viral Superinfection and Colonization.

BACKGROUND: More than 2 months separated the initial description of SARS-CoV-2 and discovery of its widespread dissemination in the United States. Despite this lengthy interval, implementation of specific quantitative reverse transcription (qRT)-PCR-based SARS-CoV-2 tests in the US has been slow, and testing is still not widely available. Metagenomic sequencing offers the promise of unbiased detection of emerging pathogens, without requiring prior knowledge of the identity of the responsible agent or its genomic sequence. METHODS: To evaluate metagenomic approaches in the context of the current SARS-CoV-2 epidemic, laboratory-confirmed positive and negative samples from Seattle, WA were evaluated by metagenomic sequencing, with comparison to a 2019 reference genomic database created before the emergence of SARS-CoV-2. RESULTS: Within 36 h our results showed clear identification of a novel human Betacoronavirus, closely related to known Betacoronaviruses of bats, in laboratory-proven cases of SARS-CoV-2. A subset of samples also showed superinfection or colonization with human parainfluenza virus 3 or Moraxella species, highlighting the need to test directly for SARS-CoV-2 as opposed to ruling out an infection using a viral respiratory panel. Samples negative for SARS-CoV-2 by RT-PCR were also negative by metagenomic analysis, and positive for Rhinovirus A and C. Unlike targeted SARS-CoV-2 qRT-PCR testing, metagenomic analysis of these SARS-CoV-2 negative samples identified candidate etiological agents for the patients' respiratory symptoms. CONCLUSION: Taken together, these results demonstrate the value of metagenomic analysis in the monitoring and response to this and future viral pandemics.

Comparative genomic, transcriptomic, and proteomic reannotation of human herpesvirus 6.

BACKGROUND: Human herpesvirus-6A and -6B (HHV-6) are betaherpesviruses that reach > 90% seroprevalence in the adult population. Unique among human herpesviruses, HHV-6 can integrate into the subtelomeric regions of human chromosomes; when this occurs in germ line cells it causes a condition called inherited chromosomally integrated HHV-6 (iciHHV-6). Only two complete genomes are available for replicating HHV-6B, leading to numerous conflicting annotations and little known about the global genomic diversity of this ubiquitous virus. RESULTS: Using a custom capture panel for HHV-6B, we report complete genomes from 61 isolates of HHV-6B from active infections (20 from Japan, 35 from New York state, and 6 from Uganda), and 64 strains of iciHHV-6B (mostly from North America). HHV-6B sequence clustered by geography and illustrated extensive recombination. Multiple iciHHV-6B sequences from unrelated individuals across the United States were found to be completely identical, consistent with a founder effect. Several iciHHV-6B strains clustered with strains from recent active pediatric infection. Combining our genomic analysis with the first RNA-Seq and shotgun proteomics studies of HHV-6B, we completely reannotated the HHV-6B genome, altering annotations for more than 10% of existing genes, with multiple instances of novel splicing and genes that hitherto had gone unannotated. CONCLUSION: Our results are consistent with a model of intermittent de novo integration of HHV-6B into host germline cells during active infection with a large contribution of founder effect in iciHHV-6B. Our data provide a significant advance in the genomic annotation of HHV-6B, which will contribute to the detection, diversity, and control of this virus.

Transposable element RNA dysregulation in mutant KRAS(G12C) 3D lung cancer spheroids.

Mutant KRAS regulates transposable element (TE) RNA and interferon-stimulated gene (ISG) expression, but it remains unclear whether diverse mutations in KRAS affect different TE RNAs throughout the genome. We analyzed the transcriptomes of 3D human lung cancer spheroids that harbor KRAS(G12C) mutations to determine the landscape of TE RNAs regulated by mutant KRAS(G12C). We found that KRAS(G12C) signaling is required for the expression of LINE- and LTR-derived TE RNAs that are distinct from TE RNAs previously shown to be regulated by mutant KRAS(G12D) or KRAS(G12V). Moreover, KRAS(G12C) inhibition specifically upregulates SINE-derived TE RNAs from the youngest Alu subfamily AluY. Our results reveal that TE RNA dysregulation in KRAS-driven lung cancer cells is mutation-dependent, while also highlighting a subset of young, Alu-derived TE RNAs that are coordinately activated with innate immunity genes upon KRAS(G12C) inhibition.

Profiling of repetitive RNA sequences in the blood plasma of patients with cancer.

Liquid biopsies provide a means for the profiling of cell-free RNAs secreted by cells throughout the body. Although well-annotated coding and non-coding transcripts in blood are readily detectable and can serve as biomarkers of disease, the overall diagnostic utility of the cell-free transcriptome remains unclear. Here we show that RNAs derived from transposable elements and other repeat elements are enriched in the cell-free transcriptome of patients with cancer, and that they serve as signatures for the accurate classification of the disease. We used repeat-element-aware liquid-biopsy technology and single-molecule nanopore sequencing to profile the cell-free transcriptome in plasma from patients with cancer and to examine millions of genomic features comprising all annotated genes and repeat elements throughout the genome. By aggregating individual repeat elements to the subfamily level, we found that samples with pancreatic cancer are enriched with specific Alu subfamilies, whereas other cancers have their own characteristic cell-free RNA profile. Our findings show that repetitive RNA sequences are abundant in blood and can be used as disease-specific diagnostic biomarkers.

Mutant KRAS regulates transposable element RNA and innate immunity via KRAB zinc-finger genes.

RAS genes are the most frequently mutated oncogenes in cancer, yet the effects of oncogenic RAS signaling on the noncoding transcriptome remain unclear. We analyzed the transcriptomes of human airway and bronchial epithelial cells transformed with mutant KRAS to define the landscape of KRAS-regulated noncoding RNAs. We find that oncogenic KRAS signaling upregulates noncoding transcripts throughout the genome, many of which arise from transposable elements (TEs). These TE RNAs exhibit differential expression, are preferentially released in extracellular vesicles, and are regulated by KRAB zinc-finger (KZNF) genes, which are broadly downregulated in mutant KRAS cells and lung adenocarcinomas in vivo. Moreover, mutant KRAS induces an intrinsic IFN-stimulated gene (ISG) signature that is often seen across many different cancers. Our results indicate that mutant KRAS remodels the repetitive noncoding transcriptome, demonstrating the broad scope of intracellular and extracellular RNAs regulated by this oncogenic signaling pathway.

Host-pathogen dynamics in longitudinal clinical specimens from patients with COVID-19.

Rapid dissemination of SARS-CoV-2 sequencing data to public repositories has enabled widespread study of viral genomes, but studies of longitudinal specimens from infected persons are relatively limited. Analysis of longitudinal specimens enables understanding of how host immune pressures drive viral evolution in vivo. Here we performed sequencing of 49 longitudinal SARS-CoV-2-positive samples from 20 patients in Washington State collected between March and September of 2020. Viral loads declined over time with an average increase in RT-QPCR cycle threshold of 0.87 per day. We found that there was negligible change in SARS-CoV-2 consensus sequences over time, but identified a number of nonsynonymous variants at low frequencies across the genome. We observed enrichment for a relatively small number of these variants, all of which are now seen in consensus genomes across the globe at low prevalence. In one patient, we saw rapid emergence of various low-level deletion variants at the N-terminal domain of the spike glycoprotein, some of which have previously been shown to be associated with reduced neutralization potency from sera. In a subset of samples that were sequenced using metagenomic methods, differential gene expression analysis showed a downregulation of cytoskeletal genes that was consistent with a loss of ciliated epithelium during infection and recovery. We also identified co-occurrence of bacterial species in samples from multiple hospitalized individuals. These results demonstrate that the intrahost genetic composition of SARS-CoV-2 is dynamic during the course of COVID-19, and highlight the need for continued surveillance and deep sequencing of minor variants.

Molecular Features of the Measles Virus Viral Fusion Complex That Favor Infection and Spread in the Brain.

Measles virus (MeV) bearing a single amino acid change in the fusion protein (F)-L454W-was isolated from two patients who died of MeV central nervous system (CNS) infection. This mutation in F confers an advantage over wild-type virus in the CNS, contributing to disease in these patients. Using murine ex vivo organotypic brain cultures and human induced pluripotent stem cell-derived brain organoids, we show that CNS adaptive mutations in F enhance the spread of virus ex vivo. The spread of virus in human brain organoids is blocked by an inhibitory peptide that targets F, confirming that dissemination in the brain tissue is attributable to F. A single mutation in MeV F thus alters the fusion complex to render MeV more neuropathogenic. IMPORTANCE Measles virus (MeV) infection can cause serious complications in immunocompromised individuals, including measles inclusion body encephalitis (MIBE). In some cases, MeV persistence and subacute sclerosing panencephalitis (SSPE), another severe central nervous system (CNS) complication, develop even in the face of a systemic immune response. Both MIBE and SSPE are relatively rare but lethal. It is unclear how MeV causes CNS infection. We introduced specific mutations that are found in MIBE or SSPE cases into the MeV fusion protein to test the hypothesis that dysregulation of the viral fusion complex-comprising F and the receptor binding protein, H-allows virus to spread in the CNS. Using metagenomic, structural, and biochemical approaches, we demonstrate that altered fusion properties of the MeV H-F fusion complex permit MeV to spread in brain tissue.

Cell free DNA from respiratory pathogens is detectable in the blood plasma of Cystic Fibrosis patients.

Diagnostically informative microbial cell-free DNA (cfDNA) can be detected from blood plasma during fulminant infections such as sepsis. However, the potential for DNA from airway pathogens to enter the circulation of cystic fibrosis (CF) patients during chronic infective states has not yet been evaluated. We assessed whether patient blood contained measurable quantities of cfDNA from CF respiratory microorganisms by sequencing plasma from 21 individuals with CF recruited from outpatient clinics and 12 healthy controls. To account for possible contamination with exogenous microbial nucleic acids, statistical significance of microbe-derived read counts from CF patients was determined relative to the healthy control population. In aggregate, relative abundance of microbial cfDNA was nearly an order of magnitude higher in CF patients than in healthy subjects (p = 8.0×10-3). 15 of 21 (71%) CF patients demonstrated cfDNA from one or more relevant organisms. In contrast, none of the healthy subjects evidenced significant microbial cfDNA for any of the organisms examined. Concordance of cfDNA with standard microbiological culture of contemporaneously collected patient sputum was variable. Our findings provide evidence that cfDNA from respiratory pathogens are present in the bloodstream of most CF patients, which could potentially be exploited for the purposes of noninvasive clinical diagnosis.

In vivo antiviral host response to SARS-CoV-2 by viral load, sex, and age.

Despite limited genomic diversity, SARS-CoV-2 has shown a wide range of clinical manifestations in different patient populations. The mechanisms behind these host differences are still unclear. Here, we examined host response gene expression across infection status, viral load, age, and sex among shotgun RNA-sequencing profiles of nasopharyngeal swabs from 430 individuals with PCR-confirmed SARS-CoV-2 and 54 negative controls. SARS-CoV-2 induced a strong antiviral response with upregulation of antiviral factors such as OAS1-3 and IFIT1-3 , and Th1 chemokines CXCL9/10/11 , as well as a reduction in transcription of ribosomal proteins. SARS-CoV-2 culture in human airway epithelial cultures replicated the in vivo antiviral host response. Patient-matched longitudinal specimens (mean elapsed time = 6.3 days) demonstrated reduction in interferon-induced transcription, recovery of transcription of ribosomal proteins, and initiation of wound healing and humoral immune responses. Expression of interferon-responsive genes, including ACE2 , increased as a function of viral load, while transcripts for B cell-specific proteins and neutrophil chemokines were elevated in patients with lower viral load. Older individuals had reduced expression of Th1 chemokines CXCL9/10/11 and their cognate receptor, CXCR3 , as well as CD8A and granzyme B, suggesting deficiencies in trafficking and/or function of cytotoxic T cells and natural killer (NK) cells. Relative to females, males had reduced B and NK cell-specific transcripts and an increase in inhibitors of NF-κB signaling, possibly inappropriately throttling antiviral responses. Collectively, our data demonstrate that host responses to SARS-CoV-2 are dependent on viral load and infection time course, with observed differences due to age and sex that may contribute to disease severity.

Inhibition of Coronavirus Entry In Vitro and Ex Vivo by a Lipid-Conjugated Peptide Derived from the SARS-CoV-2 Spike Glycoprotein HRC Domain.

The emergence of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), the etiological agent of the 2019 coronavirus disease (COVID-19), has erupted into a global pandemic that has led to tens of millions of infections and hundreds of thousands of deaths worldwide. The development of therapeutics to treat infection or as prophylactics to halt viral transmission and spread is urgently needed. SARS-CoV-2 relies on structural rearrangements within a spike (S) glycoprotein to mediate fusion of the viral and host cell membranes. Here, we describe the development of a lipopeptide that is derived from the C-terminal heptad repeat (HRC) domain of SARS-CoV-2 S that potently inhibits infection by SARS-CoV-2. The lipopeptide inhibits cell-cell fusion mediated by SARS-CoV-2 S and blocks infection by live SARS-CoV-2 in Vero E6 cell monolayers more effectively than previously described lipopeptides. The SARS-CoV-2 lipopeptide exhibits broad-spectrum activity by inhibiting cell-cell fusion mediated by SARS-CoV-1 and Middle East respiratory syndrome coronavirus (MERS-CoV) and blocking infection by live MERS-CoV in cell monolayers. We also show that the SARS-CoV-2 HRC-derived lipopeptide potently blocks the spread of SARS-CoV-2 in human airway epithelial (HAE) cultures, an ex vivo model designed to mimic respiratory viral propagation in humans. While viral spread of SARS-CoV-2 infection was widespread in untreated airways, those treated with SARS-CoV-2 HRC lipopeptide showed no detectable evidence of viral spread. These data provide a framework for the development of peptide therapeutics for the treatment of or prophylaxis against SARS-CoV-2 as well as other coronaviruses.IMPORTANCE SARS-CoV-2, the causative agent of COVID-19, continues to spread globally, placing strain on health care systems and resulting in rapidly increasing numbers of cases and mortalities. Despite the growing need for medical intervention, no FDA-approved vaccines are yet available, and treatment has been limited to supportive therapy for the alleviation of symptoms. Entry inhibitors could fill the important role of preventing initial infection and preventing spread. Here, we describe the design, synthesis, and evaluation of a lipopeptide that is derived from the HRC domain of the SARS-CoV-2 S glycoprotein that potently inhibits fusion mediated by SARS-CoV-2 S glycoprotein and blocks infection by live SARS-CoV-2 in both cell monolayers (in vitro) and human airway tissues (ex vivo). Our results highlight the SARS-CoV-2 HRC-derived lipopeptide as a promising therapeutic candidate for SARS-CoV-2 infections.

High-resolution profiling of human cytomegalovirus cell-free DNA in human plasma highlights its exceptionally fragmented nature.

Human cytomegalovirus (CMV) infections comprise a leading cause of newborn impairments worldwide and are pervasive concerns among the immunocompromised. Quantification of CMV viral loads is increasingly used to guide definitions of CMV disease but standardization of CMV quantitation remains problematic, mostly due to differences in qPCR amplicon sizes between clinical laboratories. Here, we used plasma cfDNA sequencing data from 2,208 samples sent for non-invasive prenatal aneuploidy screening to detect CMV and precisely measure the length of CMV fragments in human plasma. CMV reads were identified in 120 (5.4%) samples. Median cfDNA fragment size derived from CMV was significantly shorter than cfDNA derived from human chromosomes (103 vs 172 bp, p < 0.0001), corresponding to the 3rd percentile of human cfDNA. Sequencing of cfDNA from seven plasma samples from transplant patients positive for CMV confirmed the extraordinarily short nature of CMV cfDNA fragment size with a median length of 149 bp. We further show that these high-resolution measurements of CMV DNA fragment size accurately predict measured discrepancies in serum viral load measurements by different qPCR assays. These results highlight the exceptionally fragmented nature of CMV cfDNA and illustrate the promise of plasma cfDNA sequencing for quantitating viral loads through detection of fragments that would be unrecoverable by qPCR.

NCBI's Virus Discovery Hackathon: Engaging Research Communities to Identify Cloud Infrastructure Requirements.

A wealth of viral data sits untapped in publicly available metagenomic data sets when it might be extracted to create a usable index for the virological research community. We hypothesized that work of this complexity and scale could be done in a hackathon setting. Ten teams comprised of over 40 participants from six countries, assembled to create a crowd-sourced set of analysis and processing pipelines for a complex biological data set in a three-day event on the San Diego State University campus starting 9 January 2019. Prior to the hackathon, 141,676 metagenomic data sets from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) were pre-assembled into contiguous assemblies (contigs) by NCBI staff. During the hackathon, a subset consisting of 2953 SRA data sets (approximately 55 million contigs) was selected, which were further filtered for a minimal length of 1 kb. This resulted in 4.2 million (Mio) contigs, which were aligned using BLAST against all known virus genomes, phylogenetically clustered and assigned metadata. Out of the 4.2 Mio contigs, 360,000 contigs were labeled with domains and an additional subset containing 4400 contigs was screened for virus or virus-like genes. The work yielded valuable insights into both SRA data and the cloud infrastructure required to support such efforts, revealing analysis bottlenecks and possible workarounds thereof. Mainly: (i) Conservative assemblies of SRA data improves initial analysis steps; (ii) existing bioinformatic software with weak multithreading/multicore support can be elevated by wrapper scripts to use all cores within a computing node; (iii) redesigning existing bioinformatic algorithms for a cloud infrastructure to facilitate its use for a wider audience; and (iv) a cloud infrastructure allows a diverse group of researchers to collaborate effectively. The scientific findings will be extended during a follow-up event. Here, we present the applied workflows, initial results, and lessons learned from the hackathon.