Using symptom-based case predictions to identify host genetic factors that contribute to COVID-19 susceptibility.

Epidemiological and genetic studies on COVID-19 are currently hindered by inconsistent and limited testing policies to confirm SARS-CoV-2 infection. Recently, it was shown that it is possible to predict COVID-19 cases using cross-sectional self-reported disease-related symptoms. Here, we demonstrate that this COVID-19 prediction model has reasonable and consistent performance across multiple independent cohorts and that our attempt to improve upon this model did not result in improved predictions. Using the existing COVID-19 prediction model, we then conducted a GWAS on the predicted phenotype using a total of 1,865 predicted cases and 29,174 controls. While we did not find any common, large-effect variants that reached genome-wide significance, we do observe suggestive genetic associations at two SNPs (rs11844522, p = 1.9x10-7; rs5798227, p = 2.2x10-7). Explorative analyses furthermore suggest that genetic variants associated with other viral infectious diseases do not overlap with COVID-19 susceptibility and that severity of COVID-19 may have a different genetic architecture compared to COVID-19 susceptibility. This study represents a first effort that uses a symptom-based predicted phenotype as a proxy for COVID-19 in our pursuit of understanding the genetic susceptibility of the disease. We conclude that the inclusion of symptom-based predicted cases could be a useful strategy in a scenario of limited testing, either during the current COVID-19 pandemic or any future viral outbreak.

Decoding diversity in a coral reef fish species complex with restricted range using metagenomic sequencing of gut contents.

AIM: Identification of the processes that generate and maintain species diversity within the same region can provide insight into biogeographic patterns at broader spatiotemporal scales. Hawkfishes in the genus Paracirrhites are a unique taxon to explore with respect to niche differentiation, exhibiting diagnostic differences in coloration, and an apparent center of distribution outside of the Indo-Malay-Philippine (IMP) biodiversity hotspot for coral reef fishes. Our aim is to use next-generation sequencing methods to leverage samples of a taxon at their center of maximum diversity to explore phylogenetic relationships and a possible mechanism of coexistence. LOCATION: Flint Island, Southern Line Islands, Republic of Kiribati. METHODS: A comprehensive review of museum records, the primary literature, and unpublished field survey records was undertaken to determine ranges for four "arc-eye" hawkfish species in the Paracirrhites species complex and a potential hybrid. Fish from four Paracirrhites species were collected from Flint Island in the Southern Line Islands, Republic of Kiribati. Hindgut contents were sequenced, and subsequent metagenomic analyses were used to assess the phylogenetic relatedness of the host fish, the microbiome community structure, and prey remains for each species. RESULTS: Phylogenetic analyses conducted with recovered mitochondrial genomes revealed clustering of P. bicolor with P. arcatus and P. xanthus with P. nisus, which were unexpected on the basis of previous morphological work in this species complex. Differences in taxonomic composition of gut microbial communities and presumed prey remains indicate likely separation of foraging niches. MAIN CONCLUSIONS: Our findings point toward previously unidentified relationships in this cryptic species complex at its proposed center of distribution. The three species endemic to the Polynesian province (P. nisus, P. xanthus, and P. bicolor) cluster separately from the more broadly distributed P. arcatus on the basis of relative abundance of metazoan sequences in the gut (presumed prey remains). Discordance between gut microbial communities and phylogeny of the host fish further reinforce the hypothesis of niche separation.

Molecular and Microbial Microenvironments in Chronically Diseased Lungs Associated with Cystic Fibrosis.

To visualize the personalized distributions of pathogens and chemical environments, including microbial metabolites, pharmaceuticals, and their metabolic products, within and between human lungs afflicted with cystic fibrosis (CF), we generated three-dimensional (3D) microbiome and metabolome maps of six explanted lungs from three cystic fibrosis patients. These 3D spatial maps revealed that the chemical environments differ between patients and within the lungs of each patient. Although the microbial ecosystems of the patients were defined by the dominant pathogen, their chemical diversity was not. Additionally, the chemical diversity between locales in the lungs of the same individual sometimes exceeded interindividual variation. Thus, the chemistry and microbiome of the explanted lungs appear to be not only personalized but also regiospecific. Previously undescribed analogs of microbial quinolones and antibiotic metabolites were also detected. Furthermore, mapping the chemical and microbial distributions allowed visualization of microbial community interactions, such as increased production of quorum sensing quinolones in locations where Pseudomonas was in contact with Staphylococcus and Granulicatella, consistent with in vitro observations of bacteria isolated from these patients. Visualization of microbe-metabolite associations within a host organ in early-stage CF disease in animal models will help elucidate the complex interplay between the presence of a given microbial structure, antibiotics, metabolism of antibiotics, microbial virulence factors, and host responses.IMPORTANCE Microbial infections are now recognized to be polymicrobial and personalized in nature. Comprehensive analysis and understanding of the factors underlying the polymicrobial and personalized nature of infections remain limited, especially in the context of the host. By visualizing microbiomes and metabolomes of diseased human lungs, we reveal how different the chemical environments are between hosts that are dominated by the same pathogen and how community interactions shape the chemical environment or vice versa. We highlight that three-dimensional organ mapping methods represent hypothesis-building tools that allow us to design mechanistic studies aimed at addressing microbial responses to other microbes, the host, and pharmaceutical drugs.

Diel population and functional synchrony of microbial communities on coral reefs.

On coral reefs, microorganisms are essential for recycling nutrients to primary producers through the remineralization of benthic-derived organic matter. Diel investigations of reef processes are required to holistically understand the functional roles of microbial players in these ecosystems. Here we report a metagenomic analysis characterizing microbial communities in the water column overlying 16 remote forereef sites over a diel cycle. Our results show that microbial community composition is more dissimilar between day and night samples collected from the same site than between day or night samples collected across geographically distant reefs. Diel community differentiation is largely driven by the flux of Psychrobacter sp., which is two-orders of magnitude more abundant during the day. Nighttime communities are enriched with species of Roseobacter, Halomonas, and Alteromonas encoding a greater variety of pathways for carbohydrate catabolism, further illustrating temporal patterns of energetic provisioning between different marine microbes. Dynamic diel fluctuations of microbial populations could also support the efficient trophic transfer of energy posited in coral reef food webs.

Cystic Fibrosis Rapid Response: Translating Multi-omics Data into Clinically Relevant Information.

Pulmonary exacerbations are the leading cause of death in cystic fibrosis (CF) patients. To track microbial dynamics during acute exacerbations, a CF rapid response (CFRR) strategy was developed. The CFRR relies on viromics, metagenomics, metatranscriptomics, and metabolomics data to rapidly monitor active members of the viral and microbial community during acute CF exacerbations. To highlight CFRR, a case study of a CF patient is presented, in which an abrupt decline in lung function characterized a fatal exacerbation. The microbial community in the patient's lungs was closely monitored through the multi-omics strategy, which led to the identification of pathogenic shigatoxigenic Escherichia coli (STEC) expressing Shiga toxin. This case study illustrates the potential for the CFRR to deconstruct complicated disease dynamics and provide clinicians with alternative treatments to improve the outcomes of pulmonary exacerbations and expand the life spans of individuals with CF.IMPORTANCE Proper management of polymicrobial infections in patients with cystic fibrosis (CF) has extended their life span. Information about the composition and dynamics of each patient's microbial community aids in the selection of appropriate treatment of pulmonary exacerbations. We propose the cystic fibrosis rapid response (CFRR) as a fast approach to determine viral and microbial community composition and activity during CF pulmonary exacerbations. The CFRR potential is illustrated with a case study in which a cystic fibrosis fatal exacerbation was characterized by the presence of shigatoxigenic Escherichia coli The incorporation of the CFRR within the CF clinic could increase the life span and quality of life of CF patients.

Three-Dimensional Microbiome and Metabolome Cartography of a Diseased Human Lung.

Our understanding of the spatial variation in the chemical and microbial makeup of an entire human organ remains limited, in part due to the size and heterogeneity of human organs and the complexity of the associated metabolome and microbiome. To address this challenge, we developed a workflow to enable the cartography of metabolomic and microbiome data onto a three-dimensional (3D) organ reconstruction built off radiological images. This enabled the direct visualization of the microbial and chemical makeup of a human lung from a cystic fibrosis patient. We detected host-derived molecules, microbial metabolites, medications, and region-specific metabolism of medications and placed it in the context of microbial distributions in the lung. Our tool further created browsable maps of a 3D microbiome/metabolome reconstruction map on a radiological image of a human lung and forms an interactive resource for the scientific community.

Spatial Molecular Architecture of the Microbial Community of a Peltigera Lichen.

Microbes are commonly studied as individual species, but they exist as mixed assemblages in nature. At present, we know very little about the spatial organization of the molecules, including natural products that are produced within these microbial networks. Lichens represent a particularly specialized type of symbiotic microbial assemblage in which the component microorganisms exist together. These composite microbial assemblages are typically comprised of several types of microorganisms representing phylogenetically diverse life forms, including fungi, photosymbionts, bacteria, and other microbes. Here, we employed matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) imaging mass spectrometry to characterize the distributions of small molecules within a Peltigera lichen. In order to probe how small molecules are organized and localized within the microbial consortium, analytes were annotated and assigned to their respective producer microorganisms using mass spectrometry-based molecular networking and metagenome sequencing. The spatial analysis of the molecules not only reveals an ordered layering of molecules within the lichen but also supports the compartmentalization of unique functions attributed to various layers. These functions include chemical defense (e.g., antibiotics), light-harvesting functions associated with the cyanobacterial outer layer (e.g., chlorophyll), energy transfer (e.g., sugars) surrounding the sun-exposed cyanobacterial layer, and carbohydrates that may serve a structural or storage function and are observed with higher intensities in the non-sun-exposed areas (e.g., complex carbohydrates). IMPORTANCE Microbial communities have evolved over centuries to live symbiotically. The direct visualization of such communities at the chemical and functional level presents a challenge. Overcoming this challenge may allow one to visualize the spatial distributions of specific molecules involved in symbiosis and to define their functional roles in shaping the community structure. In this study, we examined the diversity of microbial genes and taxa and the presence of biosynthetic gene clusters by metagenomic sequencing and the compartmentalization of organic chemical components within a lichen using mass spectrometry. This approach allowed the identification of chemically distinct sections within this composite organism. Using our multipronged approach, various fungal natural products, not previously reported from lichens, were identified and two different fungal layers were visualized at the chemical level.

Using viromes to predict novel immune proteins in non-model organisms.

Immunity is mostly studied in a few model organisms, leaving the majority of immune systems on the planet unexplored. To characterize the immune systems of non-model organisms alternative approaches are required. Viruses manipulate host cell biology through the expression of proteins that modulate the immune response. We hypothesized that metagenomic sequencing of viral communities would be useful to identify both known and unknown host immune proteins. To test this hypothesis, a mock human virome was generated and compared to the human proteome using tBLASTn, resulting in 36 proteins known to be involved in immunity. This same pipeline was then applied to reef-building coral, a non-model organism that currently lacks traditional molecular tools like transgenic animals, gene-editing capabilities, and in vitro cell cultures. Viromes isolated from corals and compared with the predicted coral proteome resulted in 2503 coral proteins, including many proteins involved with pathogen sensing and apoptosis. There were also 159 coral proteins predicted to be involved with coral immunity but currently lacking any functional annotation. The pipeline described here provides a novel method to rapidly predict host immune components that can be applied to virtually any system with the potential to discover novel immune proteins.

Metabolomics of pulmonary exacerbations reveals the personalized nature of cystic fibrosis disease.

Background. Cystic fibrosis (CF) is a genetic disease that results in chronic infections of the lungs. CF patients experience intermittent pulmonary exacerbations (CFPE) that are associated with poor clinical outcomes. CFPE involves an increase in disease symptoms requiring more aggressive therapy. Methods. Longitudinal sputum samples were collected from 11 patients (n = 44 samples) to assess the effect of exacerbations on the sputum metabolome using liquid chromatography-tandem mass spectrometry (LC-MS/MS). The data was analyzed with MS/MS molecular networking and multivariate statistics. Results. The individual patient source had a larger influence on the metabolome of sputum than the clinical state (exacerbation, treatment, post-treatment, or stable). Of the 4,369 metabolites detected, 12% were unique to CFPE samples; however, the only known metabolites significantly elevated at exacerbation across the dataset were platelet activating factor (PAF) and a related monacylglycerophosphocholine lipid. Due to the personalized nature of the sputum metabolome, a single patient was followed for 4.2 years (capturing four separate exacerbation events) as a case study for the detection of personalized biomarkers with metabolomics. PAF and related lipids were significantly elevated during CFPEs of this patient and ceramide was elevated during CFPE treatment. Correlating the abundance of bacterial 16S rRNA gene amplicons to metabolomics data from the same samples during a CFPE demonstrated that antibiotics were positively correlated to Stenotrophomonas and Pseudomonas, while ceramides and other lipids were correlated with Streptococcus, Rothia, and anaerobes. Conclusions. This study identified PAF and other inflammatory lipids as potential biomarkers of CFPE, but overall, the metabolome of CF sputum was patient specific, supporting a personalized approach to molecular detection of CFPE onset.

Global microbialization of coral reefs.

Microbialization refers to the observed shift in ecosystem trophic structure towards higher microbial biomass and energy use. On coral reefs, the proximal causes of microbialization are overfishing and eutrophication, both of which facilitate enhanced growth of fleshy algae, conferring a competitive advantage over calcifying corals and coralline algae. The proposed mechanism for this competitive advantage is the DDAM positive feedback loop (dissolved organic carbon (DOC), disease, algae, microorganism), where DOC released by ungrazed fleshy algae supports copiotrophic, potentially pathogenic bacterial communities, ultimately harming corals and maintaining algal competitive dominance. Using an unprecedented data set of >400 samples from 60 coral reef sites, we show that the central DDAM predictions are consistent across three ocean basins. Reef algal cover is positively correlated with lower concentrations of DOC and higher microbial abundances. On turf and fleshy macroalgal-rich reefs, higher relative abundances of copiotrophic microbial taxa were identified. These microbial communities shift their metabolic potential for carbohydrate degradation from the more energy efficient Embden-Meyerhof-Parnas pathway on coral-dominated reefs to the less efficient Entner-Doudoroff and pentose phosphate pathways on algal-dominated reefs. This 'yield-to-power' switch by microorganism directly threatens reefs via increased hypoxia and greater CO2 release from the microbial respiration of DOC.

Metabolomics of reef benthic interactions reveals a bioactive lipid involved in coral defence.

Holobionts are assemblages of microbial symbionts and their macrobial host. As extant representatives of some of the oldest macro-organisms, corals and algae are important for understanding how holobionts develop and interact with one another. Using untargeted metabolomics, we show that non-self interactions altered the coral metabolome more than self-interactions (i.e. different or same genus, respectively). Platelet activating factor (PAF) and Lyso-PAF, central inflammatory modulators in mammals, were major lipid components of the coral holobionts. When corals were damaged during competitive interactions with algae, PAF increased along with expression of the gene encoding Lyso-PAF acetyltransferase; the protein responsible for converting Lyso-PAF to PAF. This shows that self and non-self recognition among some of the oldest extant holobionts involve bioactive lipids identical to those in highly derived taxa like humans. This further strengthens the hypothesis that major players of the immune response evolved during the pre-Cambrian.

Diversity, Productivity, and Stability of an Industrial Microbial Ecosystem.

Managing ecosystems to maintain biodiversity may be one approach to ensuring their dynamic stability, productivity, and delivery of vital services. The applicability of this approach to industrial ecosystems that harness the metabolic activities of microbes has been proposed but has never been tested at relevant scales. We used a tag-sequencing approach with bacterial small subunit rRNA (16S) genes and eukaryotic internal transcribed spacer 2 (ITS2) to measuring the taxonomic composition and diversity of bacteria and eukaryotes in an open pond managed for bioenergy production by microalgae over a year. Periods of high eukaryotic diversity were associated with high and more-stable biomass productivity. In addition, bacterial diversity and eukaryotic diversity were inversely correlated over time, possibly due to their opposite responses to temperature. The results indicate that maintaining diverse communities may be essential to engineering stable and productive bioenergy ecosystems using microorganisms.

Ecological networking of cystic fibrosis lung infections.

In the context of a polymicrobial infection, treating a specific pathogen poses challenges because of unknown consequences on other members of the community. The presence of ecological interactions between microbes can change their physiology and response to treatment. For example, in the cystic fibrosis lung polymicrobial infection, antimicrobial susceptibility testing on clinical isolates is often not predictive of antibiotic efficacy. Novel approaches are needed to identify the interrelationships within the microbial community to better predict treatment outcomes. Here we used an ecological networking approach on the cystic fibrosis lung microbiome characterized using 16S rRNA gene sequencing and metagenomics. This analysis showed that the community is separated into three interaction groups: Gram-positive anaerobes, Pseudomonas aeruginosa, and Staphylococcus aureus. The P. aeruginosa and S. aureus groups both anti-correlate with the anaerobic group, indicating a functional antagonism. When patients are clinically stable, these major groupings were also stable, however, during exacerbation, these communities fragment. Co-occurrence networking of functional modules annotated from metagenomics data supports that the underlying taxonomic structure is driven by differences in the core metabolism of the groups. Topological analysis of the functional network identified the non-mevalonate pathway of isoprenoid biosynthesis as a keystone for the microbial community, which can be targeted with the antibiotic fosmidomycin. This study uses ecological theory to identify novel treatment approaches against a polymicrobial disease with more predictable outcomes.

Microbial, host and xenobiotic diversity in the cystic fibrosis sputum metabolome.

Cystic fibrosis (CF) lungs are filled with thick mucus that obstructs airways and facilitates chronic infections. Pseudomonas aeruginosa is a significant pathogen of this disease that produces a variety of toxic small molecules. We used molecular networking-based metabolomics to investigate the chemistry of CF sputa and assess how the microbial molecules detected reflect the microbiome and clinical culture history of the patients. Metabolites detected included xenobiotics, P. aeruginosa specialized metabolites and host sphingolipids. The clinical culture and microbiome profiles did not correspond to the detection of P. aeruginosa metabolites in the same samples. The P. aeruginosa molecules that were detected in sputum did not match those from laboratory cultures. The pseudomonas quinolone signal (PQS) was readily detectable from cultured strains, but absent from sputum, even when its precursor molecules were present. The lack of PQS production in vivo is potentially due to the chemical nature of the CF lung environment, indicating that culture-based studies of this pathogen may not explain its behavior in the lung. The most differentially abundant molecules between CF and non-CF sputum were sphingolipids, including sphingomyelins, ceramides and lactosylceramide. As these highly abundant molecules contain the inflammatory mediator ceramide, they may have a significant role in CF hyperinflammation. This study demonstrates that the chemical makeup of CF sputum is a complex milieu of microbial, host and xenobiotic molecules. Detection of a bacterium by clinical culturing and 16S rRNA gene profiling do not necessarily reflect the active production of metabolites from that bacterium in a sputum sample.

Mass spectral similarity for untargeted metabolomics data analysis of complex mixtures.

While in nucleotide sequencing, the analysis of DNA from complex mixtures of organisms is common, this is not yet true for mass spectrometric data analysis of complex mixtures. The comparative analyses of mass spectrometry data of microbial communities at the molecular level is difficult to perform, especially in the context of a host. The challenge does not lie in generating the mass spectrometry data, rather much of the difficulty falls in the realm of how to derive relevant information from this data. The informatics based techniques to visualize and organize datasets are well established for metagenome sequencing; however, due to the scarcity of informatics strategies in mass spectrometry, it is currently difficult to cross correlate two very different mass spectrometry data sets from microbial communities and their hosts. We highlight that molecular networking can be used as an organizational tool of tandem mass spectrometry data, automated database search for rapid identification of metabolites, and as a workflow to manage and compare mass spectrometry data from complex mixtures of organisms. To demonstrate this platform, we show data analysis from hard corals and a human lung associated with cystic fibrosis.

A Winogradsky-based culture system shows an association between microbial fermentation and cystic fibrosis exacerbation.

There is a poor understanding of how the physiology of polymicrobial communities in cystic fibrosis (CF) lungs contributes to pulmonary exacerbations and lung function decline. In this study, a microbial culture system based on the principles of the Winogradsky column (WinCF system) was developed to study the physiology of CF microbes. The system used glass capillary tubes filled with artificial sputum medium to mimic a clogged airway bronchiole. Chemical indicators were added to observe microbial physiology within the tubes. Characterization of sputum samples from seven patients showed variation in pH, respiration, biofilm formation and gas production, indicating that the physiology of CF microbial communities varied among patients. Incubation of homogenized tissues from an explant CF lung mirrored responses of a Pseudomonas aeruginosa pure culture, supporting evidence that end-stage lungs are dominated by this pathogen. Longitudinal sputum samples taken through two exacerbation events in a single patient showed that a two-unit drop in pH and a 30% increase in gas production occurred in the tubes prior to exacerbation, which was reversed with antibiotic treatment. Microbial community profiles obtained through amplification and sequencing of the 16S rRNA gene showed that fermentative anaerobes became more abundant during exacerbation and were then reduced during treatment where P. aeruginosa became the dominant bacterium. Results from the WinCF experiments support the model where two functionally different CF microbial communities exist, the persistent Climax Community and the acute Attack Community. Fermentative anaerobes are hypothesized to be the core members of the Attack Community and production of acidic and gaseous products from fermentation may drive developing exacerbations. Treatment targeting the Attack Community may better resolve exacerbations and resulting lung damage.

Purifying the impure: sequencing metagenomes and metatranscriptomes from complex animal-associated samples.

The accessibility of high-throughput sequencing has revolutionized many fields of biology. In order to better understand host-associated viral and microbial communities, a comprehensive workflow for DNA and RNA extraction was developed. The workflow concurrently generates viral and microbial metagenomes, as well as metatranscriptomes, from a single sample for next-generation sequencing. The coupling of these approaches provides an overview of both the taxonomical characteristics and the community encoded functions. The presented methods use Cystic Fibrosis (CF) sputum, a problematic sample type, because it is exceptionally viscous and contains high amount of mucins, free neutrophil DNA, and other unknown contaminants. The protocols described here target these problems and successfully recover viral and microbial DNA with minimal human DNA contamination. To complement the metagenomics studies, a metatranscriptomics protocol was optimized to recover both microbial and host mRNA that contains relatively few ribosomal RNA (rRNA) sequences. An overview of the data characteristics is presented to serve as a reference for assessing the success of the methods. Additional CF sputum samples were also collected to (i) evaluate the consistency of the microbiome profiles across seven consecutive days within a single patient, and (ii) compare the consistency of metagenomic approach to a 16S ribosomal RNA gene-based sequencing. The results showed that daily fluctuation of microbial profiles without antibiotic perturbation was minimal and the taxonomy profiles of the common CF-associated bacteria were highly similar between the 16S rDNA libraries and metagenomes generated from the hypotonic lysis (HL)-derived DNA. However, the differences between 16S rDNA taxonomical profiles generated from total DNA and HL-derived DNA suggest that hypotonic lysis and the washing steps benefit in not only removing the human-derived DNA, but also microbial-derived extracellular DNA that may misrepresent the actual microbial profiles.

Unraveling the unseen players in the ocean - a field guide to water chemistry and marine microbiology.

Here we introduce a series of thoroughly tested and well standardized research protocols adapted for use in remote marine environments. The sampling protocols include the assessment of resources available to the microbial community (dissolved organic carbon, particulate organic matter, inorganic nutrients), and a comprehensive description of the viral and bacterial communities (via direct viral and microbial counts, enumeration of autofluorescent microbes, and construction of viral and microbial metagenomes). We use a combination of methods, which represent a dispersed field of scientific disciplines comprising already established protocols and some of the most recent techniques developed. Especially metagenomic sequencing techniques used for viral and bacterial community characterization, have been established only in recent years, and are thus still subjected to constant improvement. This has led to a variety of sampling and sample processing procedures currently in use. The set of methods presented here provides an up to date approach to collect and process environmental samples. Parameters addressed with these protocols yield the minimum on information essential to characterize and understand the underlying mechanisms of viral and microbial community dynamics. It gives easy to follow guidelines to conduct comprehensive surveys and discusses critical steps and potential caveats pertinent to each technique.