Microbial community organization designates distinct pulmonary exacerbation types and predicts treatment outcome in cystic fibrosis.

Polymicrobial infection of the airways is a hallmark of obstructive lung diseases such as cystic fibrosis (CF), non-CF bronchiectasis, and chronic obstructive pulmonary disease. Pulmonary exacerbations (PEx) in these conditions are associated with accelerated lung function decline and higher mortality rates. Understanding PEx ecology is challenged by high inter-patient variability in airway microbial community profiles. We analyze bacterial communities in 880 CF sputum samples collected during an observational prospective cohort study and develop microbiome descriptors to model community reorganization prior to and during 18 PEx. We identify two microbial dysbiosis regimes with opposing ecology and dynamics. Pathogen-governed PEx show hierarchical community reorganization and reduced diversity, whereas anaerobic bloom PEx display stochasticity and increased diversity. A simulation of antimicrobial treatment predicts better efficacy for hierarchically organized communities. This link between PEx, microbiome organization, and treatment success advances the development of personalized clinical management in CF and, potentially, other obstructive lung diseases.

Microbial community organization designates distinct pulmonary exacerbation types and predicts treatment outcome in cystic fibrosis.

Polymicrobial infection of the airways is a hallmark of obstructive lung diseases such as cystic fibrosis (CF), non-CF bronchiectasis, and chronic obstructive pulmonary disease. Pulmonary exacerbations (PEx) in these conditions are associated with accelerated lung function decline and higher mortality rates. An understanding of the microbial underpinnings of PEx is challenged by high inter-patient variability in airway microbial community profiles. We analyzed bacterial communities in 880 CF sputum samples and developed microbiome descriptors to model community reorganization prior to and during 18 PEx. We identified two microbial dysbiosis regimes with opposing ecology and dynamics. Pathogen-governed PEx showed hierarchical community reorganization and reduced diversity, whereas anaerobic bloom PEx displayed stochasticity and increased diversity. A simulation of antimicrobial treatment predicted better efficacy for hierarchically organized communities. This link between PEx type, microbiome organization, and treatment success advances the development of personalized clinical management in CF and, potentially, other obstructive lung diseases.

Forecasting the dynamics of a complex microbial community using integrated meta-omics.

Predicting the behaviour of complex microbial communities is challenging. However, this is essential for complex biotechnological processes such as those in biological wastewater treatment plants (BWWTPs), which require sustainable operation. Here we summarize 14 months of longitudinal meta-omics data from a BWWTP anaerobic tank into 17 temporal signals, explaining 91.1% of the temporal variance, and link those signals to ecological events within the community. We forecast the signals over the subsequent five years and use 21 extra samples collected at defined time intervals for testing and validation. Our forecasts are correct for six signals and hint on phenomena such as predation cycles. Using all the 17 forecasts and the environmental variables, we predict gene abundance and expression, with a coefficient of determination ≥0.87 for the subsequent three years. Our study demonstrates the ability to forecast the dynamics of open microbial ecosystems using interactions between community cycles and environmental parameters.

Longitudinal microbial and molecular dynamics in the cystic fibrosis lung after Elexacaftor-Tezacaftor-Ivacaftor therapy.

BACKGROUND: Cystic fibrosis (CF) is a genetic disorder causing poor mucociliary clearance in the airways and subsequent respiratory infection. The recently approved triple therapy Elexacaftor-Tezacaftor-Ivacaftor (ETI) has significantly improved lung function and decreased airway infection in persons with CF (pwCF). This improvement has been shown to occur rapidly, within the first few weeks of treatment. The effects of longer term ETI therapy on lung infection dynamics, however, remain mostly unknown. RESULTS: Here, we applied 16S rRNA gene amplicon sequencing, untargeted metabolomics, and neutral models to high-resolution, longitudinally collected sputum samples from pwCF on ETI therapy (162 samples, 7 patients) and compared to similarly collected data set from pwCF not taking ETI (630 samples, 9 patients). Because ETI reduces sputum production, samples were collected in freezers provided in the subject's homes at least 3 months after first taking ETI, with those on ETI collecting a sample approximately weekly. The lung function (%ppFEV1) of those in our longitudinal cohort significantly improved after ETI (6.91, SD = 7.74), indicating our study cohort was responsive to ETI. The daily variation of alpha- and beta-diversity of both the microbiome and metabolome was higher for those on ETI, reflecting a more dynamic microbial community and chemical environment during treatment. Four of the seven subjects on ETI were persistently infected with Pseudomonas or Burkholderia in their sputum throughout the sampling period while the total bacterial load significantly decreased with time (R = - 0.42, p = 0.01) in only one subject. The microbiome and metabolome dynamics on ETI were personalized, where some subjects had a progressive change with time on therapy, whereas others had no association with time on treatment. To further classify the augmented variance of the CF microbiome under therapy, we fit the microbiome data to a Hubbell neutral dynamics model in a patient-stratified manner and found that the subjects on ETI had better fit to a neutral model. CONCLUSION: This study shows that the longitudinal microbiology and chemistry in airway secretions from subjects on ETI has become more dynamic and neutral and that after the initial improvement in lung function, many are still persistently infected with CF pathogens.

Longitudinal Microbial and Molecular Dynamics in the Cystic Fibrosis Lung after Elexacaftor-Tezacaftor-Ivacaftor therapy.

Background: Cystic fibrosis (CF) is a genetic disorder causing poor mucociliary clearance in the airways and subsequent respiratory infection. The recently approved triple therapy Elexacaftor-Tezacaftor-Ivacaftor (ETI) has significantly improved the lung function and decreased airway infection of persons with CF (pwCF). This improvement has been shown to occur rapidly, within the first few weeks of treatment. The effects of longer term ETI therapy on lung infection dynamics, however, remains mostly unknown. Results: Here, we applied 16S rRNA gene amplicon sequencing, untargeted metabolomics, and neutral models to high-resolution, longitudinally collected sputum samples from pwCF on ETI therapy (162 samples, 7 patients) and compared to similarly collected data set of CF subjects not taking ETI (630 samples, 9 patients). Because ETI reduces sputum production, samples were collected in freezers provided in the subject's homes at least 3 months after first taking ETI, with those on ETI collecting a sample approximately weekly. The lung function (%ppFEV1) of those in our longitudinal cohort significantly improved after ETI (6.91, SD = 7.74), indicating our study cohort was responsive to ETI. The daily variation of alpha- and beta-diversity of both the microbiome and metabolome was higher for those on ETI, reflecting a more dynamic microbial community and chemical environment during treatment. Four of the seven subjects on ETI were persistently infected with Pseudomonas or Burkholderia in their sputum throughout the sampling period. The microbiome and metabolome dynamics on ETI were personalized, where some subjects had a progressive change with time on therapy, whereas others had no association with time on treatment. To further classify the augmented variance of the CF microbiome under therapy, we fit the microbiome data to a Hubbell neutral dynamics model in a patient-stratified manner and found that the subjects on ETI had better fit to a neutral model. Conclusion: This study shows that the longitudinal microbiology and chemistry in airway secretions from subjects on ETI has become more dynamic and neutral, and that after the initial improvement in lung function, many are still persistently infected with CF pathogens.

Microbial community organization designates distinct pulmonary exacerbation types and predicts treatment outcome in cystic fibrosis.

Polymicrobial infection of the airways is a hallmark of obstructive lung diseases such as cystic fibrosis (CF), non-CF bronchiectasis, and chronic obstructive pulmonary disease (COPD). Intermittent pulmonary exacerbations (PEx) in these conditions are associated with lung function decline and higher mortality rates. An understanding of the microbial underpinnings of PEx is challenged by high inter-patient variability in airway microbial community profiles. We analyzed 880 near-daily CF sputum samples and developed non-standard microbiome descriptors to model community reorganization prior and during 18 PEx. We identified two communal microbial regimes with opposing ecology and dynamics. Whereas pathogen-governed dysbiosis showed hierarchical community organization and reduced diversity, anaerobic bloom dysbiosis displayed stochasticity and increased diversity. Microbiome organization modulated the relevance of pathogens and a simulation of antimicrobial treatment predicted better efficacy for hierarchically organized microbiota. This causal link between PEx, microbiome organization, and treatment success advances the development of personalized dysbiosis management in CF and, potentially, other obstructive lung diseases.

Metagenomic sequencing reveals time, host, and body compartment-specific viral dynamics after lung transplantation.

BACKGROUND: The virome of lung transplant recipients (LTRs) under immunosuppressive therapy is dominated by non-pathogenic Anelloviridae and further includes several pathogenic viruses such as Herpesviruses or respiratory viruses. It is unclear whether the donor-derived virome in the transplanted lung influences recipient virome dynamics in other body compartments and if so, to which degree. Likewise, it is unknown whether dependencies exist among virus populations that mutually shape viral loads and kinetics. RESULTS: To address these questions, we characterized viral communities in airways and plasma of 49 LTRs and analyzed their abundance patterns in a data modeling approach. We found distinct viral clusters that were specific for body compartments and displayed independent dynamics. These clusters robustly gathered specific viral species across the patient cohort. In the lung, viral cluster abundance associated with time after transplantation and we detected mutual exclusion of viral species within the same human host. In plasma, viral cluster dynamics were associated with the indication for transplantation lacking significant short-time changes. Interestingly, pathogenic viruses in the plasma co-occurred specifically with Alpha torque virus genogroup 4 and Gamma torque virus strains suggesting shared functional or ecological requirements. CONCLUSIONS: In summary, the detailed analysis of virome dynamics after lung transplantation revealed host, body compartment, and time-specific dependency patterns among viruses. Furthermore, our results suggested genetic adaptation to the host microenvironment at the level of the virome and support the hypothesis of functional complementarity between Anellovirus groups and other persistent viruses. Video abstract.

Association of bacterial community types, functional microbial processes and lung disease in cystic fibrosis airways.

Bacterial infection and inflammation of the airways are the leading causes of morbidity and mortality in persons with cystic fibrosis (CF). The ecology of the bacterial communities inhabiting CF airways is poorly understood, especially with respect to how community structure, dynamics, and microbial metabolic activity relate to clinical outcomes. In this study, the bacterial communities in 818 sputum samples from 109 persons with CF were analyzed by sequencing bacterial 16S rRNA gene amplicons. We identified eight alternative community types (pulmotypes) by using a Dirichlet multinomial mixture model and studied their temporal dynamics in the cohort. Across patients, the pulmotypes displayed chronological patterns in the transition among each other. Furthermore, significant correlations between pulmotypes and patient clinical status were detected by using multinomial mixed effects models, principal components regression, and statistical testing. Constructing pulmotype-specific metabolic activity profiles, we found that pulmotype microbiota drive distinct community functions including mucus degradation or increased acid production. These results indicate that pulmotypes are the result of ordered, underlying drivers such as predominant metabolism, ecological competition, and niche construction and can form the basis for quantitative, predictive models supporting clinical treatment decisions.

Multi-omics profiling predicts allograft function after lung transplantation.

RATIONALE: Lung transplantation is the ultimate treatment option for patients with end-stage respiratory diseases but bears the highest mortality rate among all solid organ transplantations due to chronic lung allograft dysfunction (CLAD). The mechanisms leading to CLAD remain elusive due to an insufficient understanding of the complex post-transplant adaptation processes. OBJECTIVES: To better understand these lung adaptation processes after transplantation and to investigate their association with future changes in allograft function. METHODS: We performed an exploratory cohort study of bronchoalveolar lavage samples from 78 lung recipients and donors. We analysed the alveolar microbiome using 16S rRNA sequencing, the cellular composition using flow cytometry, as well as metabolome and lipidome profiling. MEASUREMENTS AND MAIN RESULTS: We established distinct temporal dynamics for each of the analysed data sets. Comparing matched donor and recipient samples, we revealed that recipient-specific as well as environmental factors, rather than the donor microbiome, shape the long-term lung microbiome. We further discovered that the abundance of certain bacterial strains correlated with underlying lung diseases even after transplantation. A decline in forced expiratory volume during the first second (FEV1) is a major characteristic of lung allograft dysfunction in transplant recipients. By using a machine learning approach, we could accurately predict future changes in FEV1 from our multi-omics data, whereby microbial profiles showed a particularly high predictive power. CONCLUSION: Bronchoalveolar microbiome, cellular composition, metabolome and lipidome show specific temporal dynamics after lung transplantation. The lung microbiome can predict future changes in lung function with high precision.

Signatures of ecological processes in microbial community time series.

BACKGROUND: Growth rates, interactions between community members, stochasticity, and immigration are important drivers of microbial community dynamics. In sequencing data analysis, such as network construction and community model parameterization, we make implicit assumptions about the nature of these drivers and thereby restrict model outcome. Despite apparent risk of methodological bias, the validity of the assumptions is rarely tested, as comprehensive procedures are lacking. Here, we propose a classification scheme to determine the processes that gave rise to the observed time series and to enable better model selection. RESULTS: We implemented a three-step classification scheme in R that first determines whether dependence between successive time steps (temporal structure) is present in the time series and then assesses with a recently developed neutrality test whether interactions between species are required for the dynamics. If the first and second tests confirm the presence of temporal structure and interactions, then parameters for interaction models are estimated. To quantify the importance of temporal structure, we compute the noise-type profile of the community, which ranges from black in case of strong dependency to white in the absence of any dependency. We applied this scheme to simulated time series generated with the Dirichlet-multinomial (DM) distribution, Hubbell's neutral model, the generalized Lotka-Volterra model and its discrete variant (the Ricker model), and a self-organized instability model, as well as to human stool microbiota time series. The noise-type profiles for all but DM data clearly indicated distinctive structures. The neutrality test correctly classified all but DM and neutral time series as non-neutral. The procedure reliably identified time series for which interaction inference was suitable. Both tests were required, as we demonstrated that all structured time series, including those generated with the neutral model, achieved a moderate to high goodness of fit to the Ricker model. CONCLUSIONS: We present a fast and robust scheme to classify community structure and to assess the prevalence of interactions directly from microbial time series data. The procedure not only serves to determine ecological drivers of microbial dynamics, but also to guide selection of appropriate community models for prediction and follow-up analysis.

The origin and evolution of cell types.

Cell types are the basic building blocks of multicellular organisms and are extensively diversified in animals. Despite recent advances in characterizing cell types, classification schemes remain ambiguous. We propose an evolutionary definition of a cell type that allows cell types to be delineated and compared within and between species. Key to cell type identity are evolutionary changes in the 'core regulatory complex' (CoRC) of transcription factors, that make emergent sister cell types distinct, enable their independent evolution and regulate cell type-specific traits termed apomeres. We discuss the distinction between developmental and evolutionary lineages, and present a roadmap for future research.

Challenges in microbial ecology: building predictive understanding of community function and dynamics.

The importance of microbial communities (MCs) cannot be overstated. MCs underpin the biogeochemical cycles of the earth's soil, oceans and the atmosphere, and perform ecosystem functions that impact plants, animals and humans. Yet our ability to predict and manage the function of these highly complex, dynamically changing communities is limited. Building predictive models that link MC composition to function is a key emerging challenge in microbial ecology. Here, we argue that addressing this challenge requires close coordination of experimental data collection and method development with mathematical model building. We discuss specific examples where model-experiment integration has already resulted in important insights into MC function and structure. We also highlight key research questions that still demand better integration of experiments and models. We argue that such integration is needed to achieve significant progress in our understanding of MC dynamics and function, and we make specific practical suggestions as to how this could be achieved.

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.

Biogeographic variation in the microbiome of the ecologically important sponge, Carteriospongia foliascens.

Sponges are well known for hosting dense and diverse microbial communities, but how these associations vary with biogeography and environment is less clear. Here we compared the microbiome of an ecologically important sponge species, Carteriospongia foliascens, over a large geographic area and identified environmental factors likely responsible for driving microbial community differences between inshore and offshore locations using co-occurrence networks (NWs). The microbiome of C. foliascens exhibited exceptionally high microbial richness, with more than 9,000 OTUs identified at 97% sequence similarity. A large biogeographic signal was evident at the OTU level despite similar phyla level diversity being observed across all geographic locations. The C. foliascens bacterial community was primarily comprised of Gammaproteobacteria (34.2% ± 3.4%) and Cyanobacteria (32.2% ± 3.5%), with lower abundances of Alphaproteobacteria, Bacteroidetes, unidentified Proteobacteria, Actinobacteria, Acidobacteria and Deltaproteobacteria. Co-occurrence NWs revealed a consistent increase in the proportion of Cyanobacteria over Bacteroidetes between turbid inshore and oligotrophic offshore locations, suggesting that the specialist microbiome of C. foliascens is driven by environmental factors.

Wiring for independence: positive feedback motifs facilitate individuation of traits in development and evolution.

Independent selection response of a trait is contingent on the availability of genetic variation that is not entangled with other traits. Mechanistically, such variational individuation in spite of shared genome results from gene regulation. Changes that increase individuation of traits are likely caused by gene regulatory changes. Yet the effect of regulatory evolution on population variation is understudied. Trait individuation also occurs during development. Developmental differentiation involves two stages-induction of differentiation and the maintenance of differentiated fate. The corresponding gene regulatory transition involves the feed-forward and the regulated feedback motifs. Here we consider analogous transition pattern at the evolutionary scale, establishing an autonomous regulatory sub-network involved in the independent trait variation. A population genetic simulation of regulated feedback loop dynamics under small perturbations shows a decoupling of variation in gene expression between the upstream gene and the responding downstream gene. We furthermore observe that the ranges of dynamics that can be generated by feedback and feed-forward networks overlap. Such phenotypic overlap enables genetic accessibility of network-specific expression dynamics. We suggest that feedback topology may eventually confer selective advantage leading from a gradual process to threshold individuation, i.e., the emergence of a novel trait.

Fluvial network organization imprints on microbial co-occurrence networks.

Recent studies highlight linkages among the architecture of ecological networks, their persistence facing environmental disturbance, and the related patterns of biodiversity. A hitherto unresolved question is whether the structure of the landscape inhabited by organisms leaves an imprint on their ecological networks. We analyzed, based on pyrosequencing profiling of the biofilm communities in 114 streams, how features inherent to fluvial networks affect the co-occurrence networks that the microorganisms form in these biofilms. Our findings suggest that hydrology and metacommunity dynamics, both changing predictably across fluvial networks, affect the fragmentation of the microbial co-occurrence networks throughout the fluvial network. The loss of taxa from co-occurrence networks demonstrates that the removal of gatekeepers disproportionately contributed to network fragmentation, which has potential implications for the functions biofilms fulfill in stream ecosystems. Our findings are critical because of increased anthropogenic pressures deteriorating stream ecosystem integrity and biodiversity.

Deciphering microbial interactions and detecting keystone species with co-occurrence networks.

Co-occurrence networks produced from microbial survey sequencing data are frequently used to identify interactions between community members. While this approach has potential to reveal ecological processes, it has been insufficiently validated due to the technical limitations inherent in studying complex microbial ecosystems. Here, we simulate multi-species microbial communities with known interaction patterns using generalized Lotka-Volterra dynamics. We then construct co-occurrence networks and evaluate how well networks reveal the underlying interactions and how experimental and ecological parameters can affect network inference and interpretation. We find that co-occurrence networks can recapitulate interaction networks under certain conditions, but that they lose interpretability when the effects of habitat filtering become significant. We demonstrate that networks suffer from local hot spots of spurious correlation in the neighborhood of hub species that engage in many interactions. We also identify topological features associated with keystone species in co-occurrence networks. This study provides a substantiated framework to guide environmental microbiologists in the construction and interpretation of co-occurrence networks from microbial survey datasets.

Metagenomic analysis reveals presence of Treponema denticola in a tissue biopsy of the Iceman.

Ancient hominoid genome studies can be regarded by definition as metagenomic analyses since they represent a mixture of both hominoid and microbial sequences in an environment. Here, we report the molecular detection of the oral spirochete Treponema denticola in ancient human tissue biopsies of the Iceman, a 5,300-year-old Copper Age natural ice mummy. Initially, the metagenomic data of the Iceman's genomic survey was screened for bacterial ribosomal RNA (rRNA) specific reads. Through ranking the reads by abundance a relatively high number of rRNA reads most similar to T. denticola was detected. Mapping of the metagenome sequences against the T. denticola genome revealed additional reads most similar to this opportunistic pathogen. The DNA damage pattern of specifically mapped reads suggests an ancient origin of these sequences. The haematogenous spread of bacteria of the oral microbiome often reported in the recent literature could already explain the presence of metagenomic reads specific for T. denticola in the Iceman's bone biopsy. We extended, however, our survey to an Iceman gingival tissue sample and a mouth swab sample and could thereby detect T. denticola and Porphyrimonas gingivalis, another important member of the human commensal oral microflora. Taken together, this study clearly underlines the opportunity to detect disease-associated microorganisms when applying metagenomics-enabled approaches on datasets of ancient human remains.