Novel endolithic bacteria of phylum Chloroflexota reveal a myriad of potential survival strategies in the Antarctic desert.

The ice-free McMurdo Dry Valleys of Antarctica are dominated by nutrient-poor mineral soil and rocky outcrops. The principal habitat for microorganisms is within rocks (endolithic). In this environment, microorganisms are provided with protection against sub-zero temperatures, rapid thermal fluctuations, extreme dryness, and ultraviolet and solar radiation. Endolithic communities include lichen, algae, fungi, and a diverse array of bacteria. Chloroflexota is among the most abundant bacterial phyla present in these communities. Among the Chloroflexota are four novel classes of bacteria, here named Candidatus Spiritibacteria class. nov. (=UBA5177), Candidatus Martimicrobia class. nov. (=UBA4733), Candidatus Tarhunnaeia class. nov. (=UBA6077), and Candidatus Uliximicrobia class. nov. (=UBA2235). We retrieved 17 high-quality metagenome-assembled genomes (MAGs) that represent these four classes. Based on genome predictions, all these bacteria are inferred to be aerobic heterotrophs that encode enzymes for the catabolism of diverse sugars. These and other organic substrates are likely derived from lichen, algae, and fungi, as metabolites (including photosynthate), cell wall components, and extracellular matrix components. The majority of MAGs encode the capacity for trace gas oxidation using high-affinity uptake hydrogenases, which could provide energy and metabolic water required for survival and persistence. Furthermore, some MAGs encode the capacity to couple the energy generated from H2 and CO oxidation to support carbon fixation (atmospheric chemosynthesis). All encode mechanisms for the detoxification and efflux of heavy metals. Certain MAGs encode features that indicate possible interactions with other organisms, such as Tc-type toxin complexes, hemolysins, and macroglobulins.IMPORTANCEThe ice-free McMurdo Dry Valleys of Antarctica are the coldest and most hyperarid desert on Earth. It is, therefore, the closest analog to the surface of the planet Mars. Bacteria and other microorganisms survive by inhabiting airspaces within rocks (endolithic). We identify four novel classes of phylum Chloroflexota, and, based on interrogation of 17 metagenome-assembled genomes, we predict specific metabolic and physiological adaptations that facilitate the survival of these bacteria in this harsh environment-including oxidation of trace gases and the utilization of nutrients (including sugars) derived from lichen, algae, and fungi. We propose that such adaptations allow these endolithic bacteria to eke out an existence in this cold and extremely dry habitat.

Distinct and Temporally Stable Assembly Mechanisms Shape Bacterial and Fungal Communities in Vineyard Soils.

Microbial communities in agricultural soils are fundamental for plant growth and in vineyard ecosystems contribute to defining regional wine quality. Managing soil microbes towards beneficial outcomes requires knowledge of how community assembly processes vary across taxonomic groups, spatial scales, and through time. However, our understanding of microbial assembly remains limited. To quantify the contributions of stochastic and deterministic processes to bacterial and fungal assembly across spatial scales and through time, we used 16 s rRNA gene and ITS sequencing in the soil of an emblematic wine-growing region of Italy.Combining null- and neutral-modelling, we found that assembly processes were consistent through time, but bacteria and fungi were governed by different processes. At the within-vineyard scale, deterministic selection and homogenising dispersal dominated bacterial assembly, while neither selection nor dispersal had clear influence over fungal assembly. At the among-vineyard scale, the influence of dispersal limitation increased for both taxonomic groups, but its contribution was much larger for fungal communities. These null-model-based inferences were supported by neutral modelling, which estimated a dispersal rate almost two orders-of-magnitude lower for fungi than bacteria.This indicates that while stochastic processes are important for fungal assembly, bacteria were more influenced by deterministic selection imposed by the biotic and/or abiotic environment. Managing microbes in vineyard soils could thus benefit from strategies that account for dispersal limitation of fungi and the importance of environmental conditions for bacteria. Our results are consistent with theoretical expectations whereby larger individual size and smaller populations can lead to higher levels of stochasticity.

Tissue age, orchard location and disease management influence the composition of fungal and bacterial communities present on the bark of apple trees.

Plants host microbial communities that can be affected by environmental conditions and agronomic practices. Despite the role of bark as a reservoir of plant pathogens and beneficial microorganisms, no information is available on the effects of disease management on the taxonomic composition of the bark-associated communities of apple trees. We assessed the impact of disease management strategies on fungal and bacterial communities on the bark of a scab-resistant apple cultivar in two orchard locations and for two consecutive seasons. The amplicon sequencing revealed that bark age and orchard location strongly affected fungal and bacterial diversity. Microbiota dissimilarity between orchards evolved during the growing season and showed specific temporal series for fungal and bacterial populations in old and young bark. Disease management did not induce global changes in the microbial populations across locations and seasons, but specifically affected the abundance of some taxa according to bark age, orchard location and sampling time. Therefore, the disease management applied to scab-resistant cultivars, which is based on a limited use of fungicides, partially changed the taxonomic composition of bark-associated fungal and bacterial communities, suggesting the need for a more accurate risk assessment regarding possible pathogen outbreaks.

MetaMLST: multi-locus strain-level bacterial typing from metagenomic samples.

Metagenomic characterization of microbial communities has the potential to become a tool to identify pathogens in human samples. However, software tools able to extract strain-level typing information from metagenomic data are needed. Low-throughput molecular typing schema such as Multilocus Sequence Typing (MLST) are still widely used and provide a wealth of strain-level information that is currently not exploited by metagenomic methods. We introduce MetaMLST, a software tool that reconstructs the MLST loci of microorganisms present in microbial communities from metagenomic data. Tested on synthetic and spiked-in real metagenomes, the pipeline was able to reconstruct the MLST sequences with >98.5% accuracy at coverages as low as 1×. On real samples, the pipeline showed higher sensitivity than assembly-based approaches and it proved successful in identifying strains in epidemic outbreaks as well as in intestinal, skin and gastrointestinal microbiome samples.

Diversity and Cyclical Seasonal Transitions in the Bacterial Community in a Large and Deep Perialpine Lake.

High-throughput sequencing (HTS) was used to analyze the seasonal variations in the bacterioplankton community composition (BCC) in the euphotic layer of a large and deep lake south of the Alps (Lake Garda). The BCC was analyzed throughout two annual cycles by monthly samplings using the amplification and sequencing of the V3-V4 hypervariable region of the 16S rRNA gene by the MiSeq Illumina platform. The dominant and most diverse bacterioplankton phyla were among the more frequently reported in freshwater ecosystems, including the Proteobacteria, Cyanobacteria, Bacteroidetes, Verrucomicrobia, Actinobacteria, and Planctomycetes. As a distinctive feature, the development of the BCC showed a cyclical temporal pattern in the two analyzed years and throughout the euphotic layer. The recurring temporal development was controlled by the strong seasonality in water temperature and thermal stratification, and by cyclical temporal changes in nutrients and, possibly, by the remarkable annual cyclical development of cyanobacteria and eukaryotic phytoplankton hosting bacterioplankton that characterizes Lake Garda. Further downstream analyses of operational taxonomic units associated to cyanobacteria allowed confirming the presence of the most abundant taxa previously identified by microscopy and/or phylogenetic analyses, as well as the presence of other small Synechococcales/Chroococcales and rare Nostocales never identified so far in the deep lakes south of the Alps. The implications of the high diversity and strong seasonality are relevant, opening perspectives for the definition of common and discriminating patterns characterizing the temporal and spatial distribution in the BCC, and for the application of the new sequencing technologies in the monitoring of water quality in large and deep lakes.

Intestinal Candida parapsilosis isolates from Rett syndrome subjects bear potential virulent traits and capacity to persist within the host.

BACKGROUND: Rett syndrome (RTT) is a neurological disorder mainly caused by mutations in MeCP2 gene. It has been shown that MeCP2 impairments can lead to cytokine dysregulation due to MeCP2 regulatory role in T-helper and T-reg mediated responses, thus contributing to the pro-inflammatory status associated with RTT. Furthermore, RTT subjects suffer from an intestinal dysbiosis characterized by an abnormal expansion of the Candida population, a known factor responsible for the hyper-activation of pro-inflammatory immune responses. Therefore, we asked whether the intestinal fungal population of RTT subjects might contribute the sub-inflammatory status triggered by MeCP2 deficiency. METHODS: We evaluated the cultivable gut mycobiota from a cohort of 50 RTT patients and 29 healthy controls characterizing the faecal fungal isolates for their virulence-related traits, antifungal resistance and immune reactivity in order to elucidate the role of fungi in RTT's intestinal dysbiosis and gastrointestinal physiology. RESULTS: Candida parapsilosis, the most abundant yeast species in RTT subjects, showed distinct genotypic profiles if compared to healthy controls' isolates as measured by hierarchical clustering analysis from RAPD genotyping. Their phenotypical analysis revealed that RTT's isolates produced more biofilm and were significantly more resistant to azole antifungals compared to the isolates from the healthy controls. In addition, the high levels of IL-1ß and IL-10 produced by peripheral blood mononuclear cells and the mixed Th1/Th17 cells population induced by RTT C. parapsilosis isolates suggest the capacity of these intestinal fungi to persist within the host, being potentially involved in chronic, pro-inflammatory responses. CONCLUSIONS: Here we demonstrated that intestinal C. parapsilosis isolates from RTT subjects hold phenotypic traits that might favour the previously observed low-grade intestinal inflammatory status associated with RTT. Therefore, the presence of putative virulent, pro-inflammatory C. parapsilosis strains in RTT could represent an additional factor in RTT's gastrointestinal pathophysiology, whose mechanisms are not yet clearly understood.

Temperature drives the assembly of endophytic communities' seasonal succession.

Endophytic microorganisms asymptomatically colonise plant tissues. Exploring the assembly dynamics of bacterial endophytic communities is essential to understand the functioning of the plant holobiont and to optimise their possible use as biopesticides or plant biostimulants. The variation in endophytic communities in above and below-ground organs in Vitis vinifera in the field were studied. To understand the specific effect of temperature on endophytic communities, a separate experiment was set up where grapevine cuttings were grown under controlled conditions at three different temperatures. The findings revealed the succession of endophytic communities over the year. Endophytic communities of roots and stems differ in terms of composition and dynamic response to temperature. Noticeably, compositional differences during the seasons affected bacterial taxa more in stems than in roots, suggesting that roots offer a more stable and less easily perturbed environment. Correlation abundance networks showed that the presence of several taxa (including Bradyrhizobium, Burkholderia, Dyella, Mesorhizobium, Propionibacterium and Ralstonia) is linked in both the field and the greenhouse.

Genome sequencing of disease and carriage isolates of nontypeable Haemophilus influenzae identifies discrete population structure.

One of the main hurdles for the development of an effective and broadly protective vaccine against nonencapsulated isolates of Haemophilus influenzae (NTHi) lies in the genetic diversity of the species, which renders extremely difficult the identification of cross-protective candidate antigens. To assess whether a population structure of NTHi could be defined, we performed genome sequencing of a collection of diverse clinical isolates representative of both carriage and disease and of the diversity of the natural population. Analysis of the distribution of polymorphic sites in the core genome and of the composition of the accessory genome defined distinct evolutionary clades and supported a predominantly clonal evolution of NTHi, with the majority of genetic information transmitted vertically within lineages. A correlation between the population structure and the presence of selected surface-associated proteins and lipooligosaccharide structure, known to contribute to virulence, was found. This high-resolution, genome-based population structure of NTHi provides the foundation to obtain a better understanding, of NTHi adaptation to the host as well as its commensal and virulence behavior, that could facilitate intervention strategies against disease caused by this important human pathogen.

Explaining diversity in metagenomic datasets by phylogenetic-based feature weighting.

Metagenomics is revolutionizing our understanding of microbial communities, showing that their structure and composition have profound effects on the ecosystem and in a variety of health and disease conditions. Despite the flourishing of new analysis methods, current approaches based on statistical comparisons between high-level taxonomic classes often fail to identify the microbial taxa that are differentially distributed between sets of samples, since in many cases the taxonomic schema do not allow an adequate description of the structure of the microbiota. This constitutes a severe limitation to the use of metagenomic data in therapeutic and diagnostic applications. To provide a more robust statistical framework, we introduce a class of feature-weighting algorithms that discriminate the taxa responsible for the classification of metagenomic samples. The method unambiguously groups the relevant taxa into clades without relying on pre-defined taxonomic categories, thus including in the analysis also those sequences for which a taxonomic classification is difficult. The phylogenetic clades are weighted and ranked according to their abundance measuring their contribution to the differentiation of the classes of samples, and a criterion is provided to define a reduced set of most relevant clades. Applying the method to public datasets, we show that the data-driven definition of relevant phylogenetic clades accomplished by our ranking strategy identifies features in the samples that are lost if phylogenetic relationships are not considered, improving our ability to mine metagenomic datasets. Comparison with supervised classification methods currently used in metagenomic data analysis highlights the advantages of using phylogenetic information.

Gut microbiota variations in wild yellow baboons (Papio cynocephalus) are associated with sex and habitat disturbance.

Although male and female mammals differ in biological traits and functional needs, the contribution of this sexual dimorphism to variations in gut bacteria and fungi (gut microbiota) in relation to habitat type has not been fully examined. To understand whether the combination of sex and habitat affects gut microbiota variation, we analyzed 40 fecal samples of wild yellow baboons (Papio cynocephalus) living in contrasting habitat types (intact, well-protected vs. fragmented, less protected forests) in the Udzungwa Mountains of Tanzania. Sex determination was performed using the marker genes SRY (Sex-determining Region Y) and DDX3X-DDX3Y (DEAD-Box Helicase 3). Samples were attributed to 34 individuals (19 females and 15 males) belonging to five social groups. Combining the results of sex determination with two amplicon sequencing datasets on bacterial (V1-V3 region of the 16S rRNA gene) and fungal (ITS2) gut communities, we found that overall, baboon females had a significantly higher gut bacterial richness compared to males. Beta diversity estimates indicated that bacterial composition was significantly different between males and females, and this was true for individuals from both well- and less protected forests. Our results highlight the combined role of sex and habitat type in shaping variation in gut microbial communities in wild non-human primates.

Metagenomics untangles potential adaptations of Antarctic endolithic bacteria at the fringe of habitability.

Survival and growth strategies of Antarctic endolithic microbes residing in Earth's driest and coldest desert remain virtually unknown. From 109 endolithic microbiomes, 4539 metagenome-assembled genomes were generated, 49.3 % of which were novel candidate bacterial species. We present evidence that trace gas oxidation and atmospheric chemosynthesis may be the prevalent strategies supporting metabolic activity and persistence of these ecosystems at the fringe of life and the limits of habitability.

A novel computational method identifies intra- and inter-species recombination events in Staphylococcus aureus and Streptococcus pneumoniae.

Advances in high-throughput DNA sequencing technologies have determined an explosion in the number of sequenced bacterial genomes. Comparative sequence analysis frequently reveals evidences of homologous recombination occurring with different mechanisms and rates in different species, but the large-scale use of computational methods to identify recombination events is hampered by their high computational costs. Here, we propose a new method to identify recombination events in large datasets of whole genome sequences. Using a filtering procedure of the gene conservation profiles of a test genome against a panel of strains, this algorithm identifies sets of contiguous genes acquired by homologous recombination. The locations of the recombination breakpoints are determined using a statistical test that is able to account for the differences in the natural rate of evolution between different genes. The algorithm was tested on a dataset of 75 genomes of Staphylococcus aureus and 50 genomes comprising different streptococcal species, and was able to detect intra-species recombination events in S. aureus and in Streptococcus pneumoniae. Furthermore, we found evidences of an inter-species exchange of genetic material between S. pneumoniae and Streptococcus mitis, a closely related commensal species that colonizes the same ecological niche. The method has been implemented in an R package, Reco, which is freely available from supplementary material, and provides a rapid screening tool to investigate recombination on a genome-wide scale from sequence data.

Draft genome sequence of the anatoxin-a producing cyanobacterium Tychonema bourrellyi B0820 isolated from the epilimnion of the deep Alpine Lake Garda.

We report the draft genome sequence of strain B0820 of the cyanobacterium Tychonema bourrellyi isolated from the epilimnion of Lake Garda and assembled from a metagenome of a non-axenic culture. The strain analyzed was shown to produce anatoxin-a, a potent neurotoxin that can cause fatal intoxication in exposed organisms.

Highly diverse and unknown viruses may enhance Antarctic endoliths' adaptability.

BACKGROUND: Rock-dwelling microorganisms are key players in ecosystem functioning of Antarctic ice free-areas. Yet, little is known about their diversity and ecology, and further still, viruses in these communities have been largely unexplored despite important roles related to host metabolism and nutrient cycling. To begin to address this, we present a large-scale viral catalog from Antarctic rock microbial communities. RESULTS: We performed metagenomic analyses on rocks from across Antarctica representing a broad range of environmental and spatial conditions, and which resulted in a predicted viral catalog comprising > 75,000 viral operational taxonomic units (vOTUS). We found largely undescribed, highly diverse and spatially structured virus communities which had predicted auxiliary metabolic genes (AMGs) with functions indicating that they may be potentially influencing bacterial adaptation and biogeochemistry. CONCLUSION: This catalog lays the foundation for expanding knowledge of virosphere diversity, function, spatial ecology, and dynamics in extreme environments. This work serves as a step towards exploring adaptability of microbial communities in the face of a changing climate. Video Abstract.

Intergenerational Transfer of Persistent Bacterial Communities in Female Nile Tilapia.

Resident microbial communities that can support various host functions play a key role in their development and health. In fishes, microbial symbionts are vertically transferred from the parents to their progeny. Such transfer of microbes in mouthbrooder fish species has not been reported yet. Here, we employed Nile tilapia (Oreochromis niloticus) to investigate the vertical transmission of microbes across generations using a 16S rRNA amplicon sequencing approach, based on the presence of bacteria in different generations. Our analysis revealed that the core microbiome in the buccal cavity and posterior intestine of parents shapes the gut microbiome of the progeny across generations. We speculate that the route of this transmission is via the buccal cavity. The identified core microbiome bacteria, namely Nocardioides, Propionibacterium, and Sphingomonas have been reported to play an essential role in the health and development of offspring. These core microbiome members could have specific functions in fish, similar to mammals.

Genome sequencing provides new insights on the distribution of Erwinia amylovora lineages in northern Italy.

Erwinia amylovora is a Gram-negative bacterium that colonizes a wide variety of plant species causing recurrent local outbreaks of fire blight in crops of the Rosaceae family. Recent genomic surveys have documented the limited genomic diversity of this species, possibly related to a recent evolutionary bottleneck and a strong correlation between geography and phylogenetic structure of the species. Despite its economic importance, little is known about the genetic variability of co-circulating strains during local outbreaks. Here, we report the genome sequences of 82 isolates of E. amylovora, collected from different host plants in a period of 16 years in Trentino, a small region in the Northeastern Italian Alps that has been characterized by recurrent outbreaks of fire blight in apple orchards. While the genome isolated before 2018 are closely related to other strains already present in Europe, we found a novel subclade composed only by isolates that were sampled starting from 2018 and demonstrate that the endemic population of this pathogen can be composed by mixture of strains.

Apyrase-mediated amplification of secretory IgA promotes intestinal homeostasis.

Secretory immunoglobulin A (SIgA) interaction with commensal bacteria conditions microbiota composition and function. However, mechanisms regulating reciprocal control of microbiota and SIgA are not defined. Bacteria-derived adenosine triphosphate (ATP) limits T follicular helper (Tfh) cells in the Peyer's patches (PPs) via P2X7 receptor (P2X7R) and thereby SIgA generation. Here we show that hydrolysis of extracellular ATP (eATP) by apyrase results in amplification of the SIgA repertoire. The enhanced breadth of SIgA in mice colonized with apyrase-releasing Escherichia coli influences topographical distribution of bacteria and expression of genes involved in metabolic versus immune functions in the intestinal epithelium. SIgA-mediated conditioning of bacteria and enterocyte function is reflected by differences in nutrient absorption in mice colonized with apyrase-expressing bacteria. Apyrase-induced SIgA improves intestinal homeostasis and attenuates barrier impairment and susceptibility to infection by enteric pathogens in antibiotic-induced dysbiosis. Therefore, amplification of SIgA by apyrase can be leveraged to restore intestinal fitness in dysbiotic conditions.

On the Origin and Propagation of the COVID-19 Outbreak in the Italian Province of Trento, a Tourist Region of Northern Italy.

BACKGROUND: Trentino is an Italian province with a tourism-based economy, bordering the regions of Lombardy and Veneto, where the two earliest and largest outbreaks of COVID-19 occurred in Italy. The earliest cases in Trentino were reported in the first week of March 2020, with most of the cases occurring in the winter sport areas in the Dolomites mountain range. The number of reported cases decreased over the summer months and was followed by a second wave in the autumn and winter of 2020. METHODS: we performed high-coverage Oxford Nanopore sequencing of 253 positive SARS-CoV-2 swabs collected in Trentino between March and December 2020. RESULTS: in this work, we analyzed genome sequences to trace the routes through which the virus entered the area, and assessed whether the autumnal resurgence could be attributed to lineages persisting undetected during summer, or as a consequence of new introductions. CONCLUSIONS: Comparing the draft genomes analyzed with a large selection of European sequences retrieved from GISAID we found that multiple introductions of the virus occurred at the early stage of the epidemics; the two epidemic waves were unrelated; the second wave was due to reintroductions of the virus in summer when traveling restrictions were uplifted.

DNA sequence and taxonomic gap analyses to quantify the coverage of aquatic cyanobacteria and eukaryotic microalgae in reference databases: Results of a survey in the Alpine region.

The taxonomic identification of organisms based on the amplification of specific genetic markers (metabarcoding) implicitly requires adequate discriminatory information and taxonomic coverage of environmental DNA sequences in taxonomic databases. These requirements were quantitatively examined by comparing the determination of cyanobacteria and microalgae obtained by metabarcoding and light microscopy. We used planktic and biofilm samples collected in 37 lakes and 22 rivers across the Alpine region. We focused on two of the most used and best represented genetic markers in the reference databases, namely the 16S rRNA and 18S rRNA genes. A sequence gap analysis using blastn showed that, in the identity range of 99-100%, approximately 30% (plankton) and 60% (biofilm) of the sequences did not find any close counterpart in the reference databases (NCBI GenBank). Similarly, a taxonomic gap analysis showed that approximately 50% of the cyanobacterial and eukaryotic microalgal species identified by light microscopy were not represented in the reference databases. In both cases, the magnitude of the gaps differed between the major taxonomic groups. Even considering the species determined under the microscope and represented in the reference databases, 22% and 26% were still not included in the results obtained by the blastn at percentage levels of identity ≥95% and ≥97%, respectively. The main causes were the absence of matching sequences due to amplification and/or sequencing failure and potential misidentification in the microscopy step. Our results quantitatively demonstrated that in metabarcoding the main obstacles in the classification of 16S rRNA and 18S rRNA sequences and interpretation of high-throughput sequencing biomonitoring data were due to the existence of important gaps in the taxonomic completeness of the reference databases and the short length of reads. The study focused on the Alpine region, but the extent of the gaps could be much greater in other less investigated geographic areas.

Population genetics and evolution of the pan-genome of Streptococcus pneumoniae.

The genetic variability in bacterial species is much larger than in other kingdoms of life. The gene content between pairs of isolates can diverge by as much as 30% in species like Escherichia coli or Streptococcus pneumoniae. This unexpected finding led to the introduction of the concept of the pan-genome, the set of genes that can be found in a given bacterial species. The genome of any isolate is thus composed by a "core genome" shared by all strains and characteristic of the species, and a "dispensable genome" that accounts for many of the phenotypic differences between strains. The pan-genome is usually much larger than the genome of any single isolate and, given the ability of many bacteria to exchange genetic material with the environment, constitutes a reservoir that could enhance their ability to survive in a mutating environment. To understand the evolution of the pan-genome of an important pathogen and its interactions with the commensal microbial flora, we have analyzed the genomes of 44 strains of Streptococcus pneumoniae, one of the most important causes of microbial diseases in humans. Despite evidence of extensive homologous recombination, the S. pneumoniae phylogenetic tree reconstructed from polymorphisms in the core genome identified major groups of genetically related strains. With the exception of serotype 1, the tree correlated poorly with capsular serotype, geographical site of isolation and disease outcome. The distribution of dispensable genes was consistent with phylogeny, although horizontal gene transfer events attenuated this correlation in the case of ancient lineages. Homologous recombination, involving short stretches of DNA, was the dominant evolutionary process of the core genome of S. pneumoniae. Genetic exchange with related species sharing the same ecological niche was the main mechanism of evolution of S. pneumonia; and S. mitis was the main reservoir of genetic diversity of S. pneumoniae. The pan-genome of S. pneumoniae increased logarithmically with the number of strains and linearly with the variability of the sample, suggesting that acquired genes accumulate proportionately to the age of clones.