Metagenomics for Pathogen Detection During a Mass Mortality Event in Songbirds.

Mass mortality events in wildlife can be indications of an emerging infectious disease. During the spring and summer of 2021, hundreds of dead passerines were reported across the eastern US. Birds exhibited a range of clinical signs including swollen conjunctiva, ocular discharge, ataxia, and nystagmus. As part of the diagnostic investigation, high-throughput metagenomic next-generation sequencing was performed across three molecular laboratories on samples from affected birds. Many potentially pathogenic microbes were detected, with bacteria forming the largest proportion; however, no singular agent was consistently identified, with many of the detected microbes also found in unaffected (control) birds and thus considered to be subclinical infections. Congruent results across laboratories have helped drive further investigation into alternative causes, including environmental contaminants and nutritional deficiencies. This work highlights the utility of metagenomic approaches in investigations of emerging diseases and provides a framework for future wildlife mortality events.

Enabling Data Discovery with the Astrobiology Resource Metadata Standard.

As scientific investigations increasingly adopt Open Science practices, reuse of data becomes paramount. However, despite decades of progress in internet search tools, finding relevant astrobiology datasets for an envisioned investigation remains challenging due to the precise and atypical needs of the astrobiology researcher. In response, we have developed the Astrobiology Resource Metadata Standard (ARMS), a metadata standard designed to uniformly describe astrobiology "resources," that is, virtually any product of astrobiology research. Those resources include datasets, physical samples, software (modeling codes and scripts), publications, websites, images, videos, presentations, and so on. ARMS has been formulated to describe astrobiology resources generated by individual scientists or smaller scientific teams, rather than larger mission teams who may be required to use more complex archival metadata schemes. In the following, we discuss the participatory development process, give an overview of the metadata standard, describe its current use in practice, and close with a discussion of additional possible uses and extensions.

Genetically corrected RAG2-SCID human hematopoietic stem cells restore V(D)J-recombinase and rescue lymphoid deficiency.

ABSTRACT: Recombination-activating genes (RAG1 and RAG2) are critical for lymphoid cell development and function by initiating the variable (V), diversity (D), and joining (J) (V(D)J)-recombination process to generate polyclonal lymphocytes with broad antigen specificity. The clinical manifestations of defective RAG1/2 genes range from immune dysregulation to severe combined immunodeficiencies (SCIDs), causing life-threatening infections and death early in life without hematopoietic cell transplantation (HCT). Despite improvements, haploidentical HCT without myeloablative conditioning carries a high risk of graft failure and incomplete immune reconstitution. The RAG complex is only expressed during the G0-G1 phase of the cell cycle in the early stages of T- and B-cell development, underscoring that a direct gene correction might capture the precise temporal expression of the endogenous gene. Here, we report a feasibility study using the CRISPR/Cas9-based "universal gene-correction" approach for the RAG2 locus in human hematopoietic stem/progenitor cells (HSPCs) from healthy donors and RAG2-SCID patient. V(D)J-recombinase activity was restored after gene correction of RAG2-SCID-derived HSPCs, resulting in the development of T-cell receptor (TCR) αß and γδ CD3+ cells and single-positive CD4+ and CD8+ lymphocytes. TCR repertoire analysis indicated a normal distribution of CDR3 length and preserved usage of the distal TRAV genes. We confirmed the in vivo rescue of B-cell development with normal immunoglobulin M surface expression and a significant decrease in CD56bright natural killer cells. Together, we provide specificity, toxicity, and efficacy data supporting the development of a gene-correction therapy to benefit RAG2-deficient patients.

In host evolution of beta lactam resistance during active treatment for Pseudomonas aeruginosa bacteremia.

Multidrug-resistant (MDR) Pseudomonas aeruginosa has been declared a serious threat by the United States Centers for Disease Control and Prevention. Here, we used whole genome sequencing (WGS) to investigate recurrent P. aeruginosa bloodstream infections in a severely immunocompromised patient. The infections demonstrated unusual, progressive increases in resistance to beta lactam antibiotics in the setting of active treatment with appropriate, guideline-directed agents. WGS followed by comparative genomic analysis of isolates collected over 44 days demonstrated in host evolution of a single P. aeruginosa isolate characterized by stepwise acquisition of two de-novo genetic resistance mechanisms over the course of treatment. We found a novel deletion affecting the ampC repressor ampD and neighboring gene ampE, which associated with initial cefepime treatment failure. This was followed by acquisition of a porin nonsense mutation, OprD, associated with resistance to carbapenems. This study highlights the potential for in-host evolution of P. aeruginosa during bloodstream infections in severely immunocompromised patients despite appropriate antimicrobial therapy. In addition, it demonstrates the utility of WGS for understanding unusual resistance patterns in the clinical context.

Integrated host/microbe metagenomics enables accurate lower respiratory tract infection diagnosis in critically ill children.

BACKGROUNDLower respiratory tract infection (LRTI) is a leading cause of death in children worldwide. LRTI diagnosis is challenging because noninfectious respiratory illnesses appear clinically similar and because existing microbiologic tests are often falsely negative or detect incidentally carried microbes, resulting in antimicrobial overuse and adverse outcomes. Lower airway metagenomics has the potential to detect host and microbial signatures of LRTI. Whether it can be applied at scale and in a pediatric population to enable improved diagnosis and treatment remains unclear.METHODSWe used tracheal aspirate RNA-Seq to profile host gene expression and respiratory microbiota in 261 children with acute respiratory failure. We developed a gene expression classifier for LRTI by training on patients with an established diagnosis of LRTI (n = 117) or of noninfectious respiratory failure (n = 50). We then developed a classifier that integrates the host LRTI probability, abundance of respiratory viruses, and dominance in the lung microbiome of bacteria/fungi considered pathogenic by a rules-based algorithm.RESULTSThe host classifier achieved a median AUC of 0.967 by cross-validation, driven by activation markers of T cells, alveolar macrophages, and the interferon response. The integrated classifier achieved a median AUC of 0.986 and increased the confidence of patient classifications. When applied to patients with an uncertain diagnosis (n = 94), the integrated classifier indicated LRTI in 52% of cases and nominated likely causal pathogens in 98% of those.CONCLUSIONLower airway metagenomics enables accurate LRTI diagnosis and pathogen identification in a heterogeneous cohort of critically ill children through integration of host, pathogen, and microbiome features.FUNDINGSupport for this study was provided by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and the National Heart, Lung, and Blood Institute (UG1HD083171, 1R01HL124103, UG1HD049983, UG01HD049934, UG1HD083170, UG1HD050096, UG1HD63108, UG1HD083116, UG1HD083166, UG1HD049981, K23HL138461, and 5R01HL155418) as well as by the Chan Zuckerberg Biohub.

Prolonged silent carriage, genomic virulence potential and transmission between staff and patients characterize a neonatal intensive care unit (NICU) outbreak of methicillin-resistant Staphylococcus aureus (MRSA).

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) is an important pathogen in neonatal intensive care units (NICU) that confers significant morbidity and mortality. OBJECTIVE: Improving our understanding of MRSA transmission dynamics, especially among high-risk patients, is an infection prevention priority. METHODS: We investigated a cluster of clinical MRSA cases in the NICU using a combination of epidemiologic review and whole-genome sequencing (WGS) of isolates from clinical and surveillance cultures obtained from patients and healthcare personnel (HCP). RESULTS: Phylogenetic analysis identified 2 genetically distinct phylogenetic clades and revealed multiple silent-transmission events between HCP and infants. The predominant outbreak strain harbored multiple virulence factors. Epidemiologic investigation and genomic analysis identified a HCP colonized with the dominant MRSA outbreak strain who cared for most NICU patients who were infected or colonized with the same strain, including 1 NICU patient with severe infection 7 months before the described outbreak. These results guided implementation of infection prevention interventions that prevented further transmission events. CONCLUSIONS: Silent transmission of MRSA between HCP and NICU patients likely contributed to a NICU outbreak involving a virulent MRSA strain. WGS enabled data-driven decision making to inform implementation of infection control policies that mitigated the outbreak. Prospective WGS coupled with epidemiologic analysis can be used to detect transmission events and prompt early implementation of control strategies.

A 2-Gene Host Signature for Improved Accuracy of COVID-19 Diagnosis Agnostic to Viral Variants.

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false-negative viral PCR test results. Such tests are also susceptible to false-positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses, and nonviral conditions (n = 318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to diagnose COVID-19. We find that optimal classifiers include an interferon-stimulated gene that is strongly induced in COVID-19 compared with nonviral conditions, such as IFI6, and a second immune-response gene that is more strongly induced in other viral infections, such as GBP5. The IFI6+GBP5 classifier achieves an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n = 553). We further provide proof-of-concept demonstration that the classifier can be implemented in a clinically relevant RT-qPCR assay. Finally, we show that its performance is robust across common SARS-CoV-2 variants and is unaffected by cross-contamination, demonstrating its utility for improved accuracy of COVID-19 diagnostics. IMPORTANCE In this work, we study upper respiratory tract gene expression to develop and validate a 2-gene host-based COVID-19 diagnostic classifier and then demonstrate its implementation in a clinically practical qPCR assay. We find that the host classifier has utility for mitigating false-negative results, for example due to SARS-CoV-2 variants harboring mutations at primer target sites, and for mitigating false-positive viral PCR results due to laboratory cross-contamination. Both types of error carry serious consequences of either unrecognized viral transmission or unnecessary isolation and contact tracing. This work is directly relevant to the ongoing COVID-19 pandemic given the continued emergence of viral variants and the continued challenges of false-positive PCR assays. It also suggests the feasibility of pan-respiratory virus host-based diagnostics that would have value in congregate settings, such as hospitals and nursing homes, where unrecognized respiratory viral transmission is of particular concern.

Upper airway gene expression shows a more robust adaptive immune response to SARS-CoV-2 in children.

Unlike other respiratory viruses, SARS-CoV-2 disproportionately causes severe disease in older adults whereas disease burden in children is lower. To investigate whether differences in the upper airway immune response may contribute to this disparity, we compare nasopharyngeal gene expression in 83 children (<19-years-old; 38 with SARS-CoV-2, 11 with other respiratory viruses, 34 with no virus) and 154 older adults (>40-years-old; 45 with SARS-CoV-2, 28 with other respiratory viruses, 81 with no virus). Expression of interferon-stimulated genes is robustly activated in both children and adults with SARS-CoV-2 infection compared to the respective non-viral groups, with only subtle distinctions. Children, however, demonstrate markedly greater upregulation of pathways related to B cell and T cell activation and proinflammatory cytokine signaling, including response to TNF and production of IFNγ, IL-2 and IL-4. Cell type deconvolution confirms greater recruitment of B cells, and to a lesser degree macrophages, to the upper airway of children. Only children exhibit a decrease in proportions of ciliated cells, among the primary targets of SARS-CoV-2, upon infection. These findings demonstrate that children elicit a more robust innate and especially adaptive immune response to SARS-CoV-2 in the upper airway that likely contributes to their protection from severe disease in the lower airway.

Discovery of Chlorophyll d: Isolation and Characterization of a Far-Red Cyanobacterium from the Original Site of Manning and Strain (1943) at Moss Beach, California.

We have isolated a chlorophyll-d-containing cyanobacterium from the intertidal field site at Moss Beach, on the coast of Central California, USA, where Manning and Strain (1943) originally discovered this far-red chlorophyll. Here, we present the cyanobacterium's environmental description, culturing procedure, pigment composition, ultrastructure, and full genome sequence. Among cultures of far-red cyanobacteria obtained from red algae from the same site, this strain was an epiphyte on a brown macroalgae. Its Qyin vivo absorbance peak is centered at 704-705 nm, the shortest wavelength observed thus far among the various known Acaryochloris strains. Its Chl a/Chl d ratio was 0.01, with Chl d accounting for 99% of the total Chl d and Chl a mass. TEM imagery indicates the absence of phycobilisomes, corroborated by both pigment spectra and genome analysis. The Moss Beach strain codes for only a single set of genes for producing allophycocyanin. Genomic sequencing yielded a 7.25 Mbp circular chromosome and 10 circular plasmids ranging from 16 kbp to 394 kbp. We have determined that this strain shares high similarity with strain S15, an epiphyte of red algae, while its distinct gene complement and ecological niche suggest that this strain could be the closest known relative to the original Chl d source of Manning and Strain (1943). The Moss Beach strain is designated Acaryochloris sp. (marina) strain Moss Beach.

DiMeLo-seq: a long-read, single-molecule method for mapping protein-DNA interactions genome wide.

Studies of genome regulation routinely use high-throughput DNA sequencing approaches to determine where specific proteins interact with DNA, and they rely on DNA amplification and short-read sequencing, limiting their quantitative application in complex genomic regions. To address these limitations, we developed directed methylation with long-read sequencing (DiMeLo-seq), which uses antibody-tethered enzymes to methylate DNA near a target protein's binding sites in situ. These exogenous methylation marks are then detected simultaneously with endogenous CpG methylation on unamplified DNA using long-read, single-molecule sequencing technologies. We optimized and benchmarked DiMeLo-seq by mapping chromatin-binding proteins and histone modifications across the human genome. Furthermore, we identified where centromere protein A localizes within highly repetitive regions that were unmappable with short sequencing reads, and we estimated the density of centromere protein A molecules along single chromatin fibers. DiMeLo-seq is a versatile method that provides multimodal, genome-wide information for investigating protein-DNA interactions.

Using genomic epidemiology of SARS-CoV-2 to support contact tracing and public health surveillance in rural Humboldt County, California.

BACKGROUND: During the COVID-19 pandemic within the United States, much of the responsibility for diagnostic testing and epidemiologic response has relied on the action of county-level departments of public health. Here we describe the integration of genomic surveillance into epidemiologic response within Humboldt County, a rural county in northwest California. METHODS: Through a collaborative effort, 853 whole SARS-CoV-2 genomes were generated, representing ~58% of the 1,449 SARS-CoV-2-positive cases detected in Humboldt County as of March 12, 2021. Phylogenetic analysis of these data was used to develop a comprehensive understanding of SARS-CoV-2 introductions to the county and to support contact tracing and epidemiologic investigations of all large outbreaks in the county. RESULTS: In the case of an outbreak on a commercial farm, viral genomic data were used to validate reported epidemiologic links and link additional cases within the community who did not report a farm exposure to the outbreak. During a separate outbreak within a skilled nursing facility, genomic surveillance data were used to rule out the putative index case, detect the emergence of an independent Spike:N501Y substitution, and verify that the outbreak had been brought under control. CONCLUSIONS: These use cases demonstrate how developing genomic surveillance capacity within local public health departments can support timely and responsive deployment of genomic epidemiology for surveillance and outbreak response based on local needs and priorities.

Discovering disease-causing pathogens in resource-scarce Southeast Asia using a global metagenomic pathogen monitoring system.

SignificanceMetagenomic pathogen sequencing offers an unbiased approach to characterizing febrile illness. In resource-scarce settings with high biodiversity, it is critical to identify disease-causing pathogens in order to understand burden and to prioritize efforts for control. Here, metagenomic next-generation sequencing (mNGS) characterization of the pathogen landscape in Cambodia revealed diverse vector-borne and zoonotic pathogens irrespective of age and gender as risk factors. Identification of key pathogens led to changes in national program surveillance. This study is a "real world" example of the use of mNGS surveillance of febrile individuals, executed in-country, to identify outbreaks of vector-borne, zoonotic, and other emerging pathogens in a resource-scarce setting.

Early prediction of preeclampsia in pregnancy with cell-free RNA.

Liquid biopsies that measure circulating cell-free RNA (cfRNA) offer an opportunity to study the development of pregnancy-related complications in a non-invasive manner and to bridge gaps in clinical care1-4. Here we used 404 blood samples from 199 pregnant mothers to identify and validate cfRNA transcriptomic changes that are associated with preeclampsia, a multi-organ syndrome that is the second largest cause of maternal death globally5. We find that changes in cfRNA gene expression between normotensive and preeclamptic mothers are marked and stable early in gestation, well before the onset of symptoms. These changes are enriched for genes specific to neuromuscular, endothelial and immune cell types and tissues that reflect key aspects of preeclampsia physiology6-9, suggest new hypotheses for disease progression and correlate with maternal organ health. This enabled the identification and independent validation of a panel of 18 genes that when measured between 5 and 16 weeks of gestation can form the basis of a liquid biopsy test that would identify mothers at risk of preeclampsia long before clinical symptoms manifest themselves. Tests based on these observations could help predict and manage who is at risk for preeclampsia-an important objective for obstetric care10,11.

Upper airway gene expression reveals a more robust innate and adaptive immune response to SARS-CoV-2 in children compared with older adults.

Unlike other respiratory viruses, SARS-CoV-2 disproportionately causes severe disease in older adults and only rarely in children. To investigate whether differences in the upper airway immune response could contribute to this disparity, we compared nasopharyngeal gene expression in 83 children (<19-years-old; 38 with SARS-CoV-2, 11 with other respiratory viruses, 34 with no virus) and 154 adults (>40-years-old; 45 with SARS-CoV-2, 28 with other respiratory viruses, 81 with no virus). Expression of interferon-stimulated genes (ISGs) was robustly activated in both children and adults with SARS-CoV-2 compared to the respective non-viral groups, with only relatively subtle distinctions. Children, however, demonstrated markedly greater upregulation of pathways related to B cell and T cell activation and proinflammatory cytokine signaling, including TNF, IFNγ, IL-2 and IL-4 production. Cell type deconvolution confirmed greater recruitment of B cells, and to a lesser degree macrophages, to the upper airway of children. Only children exhibited a decrease in proportions of ciliated cells, the primary target of SARS-CoV-2, upon infection with the virus. These findings demonstrate that children elicit a more robust innate and adaptive immune response to SARS-CoV-2 infection in the upper airway that likely contributes to their protection from severe disease in the lower airway.

Tracheal aspirate RNA sequencing identifies distinct immunological features of COVID-19 ARDS.

The immunological features that distinguish COVID-19-associated acute respiratory distress syndrome (ARDS) from other causes of ARDS are incompletely understood. Here, we report the results of comparative lower respiratory tract transcriptional profiling of tracheal aspirate from 52 critically ill patients with ARDS from COVID-19 or from other etiologies, as well as controls without ARDS. In contrast to a "cytokine storm," we observe reduced proinflammatory gene expression in COVID-19 ARDS when compared to ARDS due to other causes. COVID-19 ARDS is characterized by a dysregulated host response with increased PTEN signaling and elevated expression of genes with non-canonical roles in inflammation and immunity. In silico analysis of gene expression identifies several candidate drugs that may modulate gene expression in COVID-19 ARDS, including dexamethasone and granulocyte colony stimulating factor. Compared to ARDS due to other types of viral pneumonia, COVID-19 is characterized by impaired interferon-stimulated gene (ISG) expression. The relationship between SARS-CoV-2 viral load and expression of ISGs is decoupled in patients with COVID-19 ARDS when compared to patients with mild COVID-19. In summary, assessment of host gene expression in the lower airways of patients reveals distinct immunological features of COVID-19 ARDS.

Impaired immune signaling and changes in the lung microbiome precede secondary bacterial pneumonia in COVID-19.

Secondary bacterial infections, including ventilator-associated pneumonia (VAP), lead to worse clinical outcomes and increased mortality following viral respiratory infections including in patients with coronavirus disease 2019 (COVID-19). Using a combination of tracheal aspirate bulk and single-cell RNA sequencing we assessed lower respiratory tract immune responses and microbiome dynamics in 23 COVID-19 patients, 10 of whom developed VAP, and eight critically ill uninfected controls. At a median of three days (range: 2-4 days) before VAP onset we observed a transcriptional signature of bacterial infection. At a median of 15 days prior to VAP onset (range: 8-38 days), we observed a striking impairment in immune signaling in COVID-19 patients who developed VAP. Longitudinal metatranscriptomic analysis revealed disruption of lung microbiome community composition in patients with VAP, providing a connection between dysregulated immune signaling and outgrowth of opportunistic pathogens. These findings suggest that COVID-19 patients who develop VAP have impaired antibacterial immune defense detectable weeks before secondary infection onset.

Lettuce (Lactuca sativa) productivity influenced by microbial inocula under nitrogen-limited conditions in aquaponics.

The demand for food will outpace productivity of conventional agriculture due to projected growth of the human population, concomitant with shrinkage of arable land, increasing scarcity of freshwater, and a rapidly changing climate. While aquaponics has potential to sustainably supplement food production with minimal environmental impact, there is a need to better characterize the complex interplay between the various components (fish, plant, microbiome) of these systems to optimize scale up and productivity. Here, we investigated how the commonly-implemented practice of continued microbial community transfer from pre-existing systems might promote or impede productivity of aquaponics. Specifically, we monitored plant growth phenotypes, water chemistry, and microbiome composition of rhizospheres, biofilters, and fish feces over 61-days of lettuce (Lactuca sativa var. crispa) growth in nitrogen-limited aquaponic systems inoculated with bacteria that were either commercially sourced or originating from a pre-existing aquaponic system. Lettuce above- and below-ground growth were significantly reduced across replicates treated with a pre-existing aquaponic system inoculum when compared to replicates treated with a commercial inoculum. Reduced productivity was associated with enrichment in specific bacterial genera in plant roots, including Pseudomonas, following inoculum transfer from pre-existing systems. Increased productivity was associated with enrichment of nitrogen-fixing Rahnella in roots of plants treated with the commercial inoculum. Thus, we show that inoculation from a pre-existing system, rather than from a commercial inoculum, is associated with lower yields. Further work will be necessary to test the putative mechanisms involved.

Upper airway gene expression reveals suppressed immune responses to SARS-CoV-2 compared with other respiratory viruses.

SARS-CoV-2 infection is characterized by peak viral load in the upper airway prior to or at the time of symptom onset, an unusual feature that has enabled widespread transmission of the virus and precipitated a global pandemic. How SARS-CoV-2 is able to achieve high titer in the absence of symptoms remains unclear. Here, we examine the upper airway host transcriptional response in patients with COVID-19 (n = 93), other viral (n = 41) or non-viral (n = 100) acute respiratory illnesses (ARIs). Compared with other viral ARIs, COVID-19 is characterized by a pronounced interferon response but attenuated activation of other innate immune pathways, including toll-like receptor, interleukin and chemokine signaling. The IL-1 and NLRP3 inflammasome pathways are markedly less responsive to SARS-CoV-2, commensurate with a signature of diminished neutrophil and macrophage recruitment. This pattern resembles previously described distinctions between symptomatic and asymptomatic viral infections and may partly explain the propensity for pre-symptomatic transmission in COVID-19. We further use machine learning to build 27-, 10- and 3-gene classifiers that differentiate COVID-19 from other ARIs with AUROCs of 0.981, 0.954 and 0.885, respectively. Classifier performance is stable across a wide range of viral load, suggesting utility in mitigating false positive or false negative results of direct SARS-CoV-2 tests.