MRSA strains with distinct accessory genes predominate at different ages in cystic fibrosis

RATIONALE

Methicillin resistant Staphylococcus aureus (MRSA) is prevalent and consequential in cystic fibrosis (CF). Whole genome sequencing (WGS) could reveal genomic differences in MRSA associated with poorer outcomes or detect MRSA transmission.

OBJECTIVES

To identify MRSA genes associated with low lung function and potential MRSA transmission in CF.

METHODS

We collected 97 MRSA isolates from 74 individuals with CF from 2017 and performed short-read WGS. We determined sequence type (ST) and the phylogenetic relationship between isolates. We aligned accessory genes from 25 reference genomes to genome assemblies, classified isolates by accessory gene content, and correlated the accessory genome to clinical outcomes.

RESULTS

The most prevalent ST were ST5 (N = 55), ST8 (N = 15), and ST105 (N = 14). Closely related MRSA strains were shared by family members with CF, but rarely between unrelated individuals. Three clusters of MRSA were identified by accessory genome content. Cluster A, including ST5 and ST105, was highly prevalent at all ages. Cluster B, including ST8, was more limited to younger patients. Cluster C included 6 distantly related strains. Patients 20 years old and younger infected with Cluster A had lower forced expiratory volume in the first second (FEV₁ ) and higher sputum biomass compared to similar-aged patients with Cluster B.

CONCLUSIONS

In this CF cohort, we identified MRSA subtypes that predominate at different ages and differ by accessory gene content. The most prevalent cluster of MRSA, including ST5 and ST105, was associated with lower FEV₁ . ST8 MRSA was more common in younger patients and thus has the potential to rise in prevalence as these patients age.

Contemporaneous radiations of fungi and plants linked to symbiosis

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Plants and fungi interact widely and in diverse ways, from mutualism to parasitism and decomposition. Here, Lutzoni et al. analyse the timing of plant and fungal evolutionary radiations and identify four major periods in which plant-fungal interactions likely drove lineage diversification.

Cyanidiales diversity in Yellowstone National Park

The Cyanidiales are unicellular red algae that are unique among phototrophs. They thrive in acidic, moderately high-temperature habitats typically associated with geothermally active regions, although much remains to be learned about their distribution and diversity within such extreme environments. We focused on Yellowstone National Park (YNP), using culture-dependent efforts in combination with a park-wide environmental polymerase chain reaction (PCR) survey to examine Cyanidiales diversity and distribution in aqueous (i.e. submerged), soil and endolithic environments. Phylogenetic reconstruction of Cyanidiales biodiversity demonstrated the presence of Cyanidioschyzon and Galdieria lineages exhibiting distinct habitat preferences. Cyanidioschyzon was the only phylotype detected in aqueous environments, but was also prominent in moist soil and endolithic habitats, environments where this genus was thought to be scarce. Galdieria was found in soil and endolithic samples, but absent in aqueous habitats. Interestingly, Cyanidium could not be found in the surveys, suggesting this genus may be absent or rare in YNP. Direct microscopic counts and viable counts from soil samples collected along a moisture gradient were positively correlated with moisture content, providing the first in situ evidence that gravimetric moisture is an important environmental parameter controlling distribution of these algae.

Good to the bone: microbial community thrives within bone cavities of a bison carcass at Yellowstone National Park

The discovery of unanticipated microbial diversity in remote, often hostile environments has led to a greater appreciation of the complexity and richness of the natural world. Yellowstone National Park (YNP) has long been a focus of work on taxa that inhabit extreme environments. Here we report the finding of microbial flora that inhabit an unexpected niche: the cavities of bone remnants from a bison carcass in Norris Geyser Basin in YNP. Although bleached white on the surface, the bone cavities are bright green due to the presence of Stichococcus-like trebouxiophyte green algae. The cavities also harbour different fungi and bacteria. Stichococcus species are common lichen photobionts and the Thelebolales fungi present in the bone cavities have previously been found in association with animal remains. Scanning electron microscope analysis suggests the fungi and algae do not form lichen-like associations in the bone. Rather these taxa and the bacteria appear to be opportunists that have colonized an isolated oasis that provides nutrients and protection from desiccation and UV radiation.

A five-gene phylogeny of Pezizomycotina

Pezizomycotina is the largest subphylum of Ascomycota and includes the vast majority of filamentous, ascoma-producing species. Here we report the results from weighted parsimony, maximum likelihood and Bayesian phylogenetic analyses of five nuclear loci (SSU rDNA, LSU rDNA, RPB1, RPB2 and EF-lalpha) from 191 taxa. Nine of the 10 Pezizomycotina classes currently recognized were represented in the sampling. These data strongly supported the monophyly of Pezizomycotina, Arthoniomycetes, Eurotiomycetes, Orbiliomycetes and Sordariomycetes. Pezizomycetes and Dothideomycetes also were resolved as monophyletic but not strongly supported by the data. Lecanoromycetes was resolved as paraphyletic in parsimony analyses but monophyletic in maximum likelihood and Bayesian analyses. Leotiomycetes was polyphyletic due to exclusion of Geoglossaceae. The two most basal classes of Pezizomycotina were Orbiliomycetes and Pezizomycetes, both of which comprise species that produce apothecial ascomata. The seven remaining classes formed a monophyletic group that corresponds to Leotiomyceta. Within Leotiomyceta, the supraclass clades of Leotiomycetes s.s. plus Sordariomycetes and Arthoniomycetes plus Dothideomycetes were resolved with moderate support.

A higher-level phylogenetic classification of the Fungi

A comprehensive phylogenetic classification of the kingdom Fungi is proposed, with reference to recent molecular phylogenetic analyses, and with input from diverse members of the fungal taxonomic community. The classification includes 195 taxa, down to the level of order, of which 16 are described or validated here: Dikarya subkingdom nov.; Chytridiomycota, Neocallimastigomycota phyla nov.; Monoblepharidomycetes, Neocallimastigomycetes class. nov.; Eurotiomycetidae, Lecanoromycetidae, Mycocaliciomycetidae subclass. nov.; Acarosporales, Corticiales, Baeomycetales, Candelariales, Gloeophyllales, Melanosporales, Trechisporales, Umbilicariales ords. nov. The clade containing Ascomycota and Basidiomycota is classified as subkingdom Dikarya, reflecting the putative synapomorphy of dikaryotic hyphae. The most dramatic shifts in the classification relative to previous works concern the groups that have traditionally been included in the Chytridiomycota and Zygomycota. The Chytridiomycota is retained in a restricted sense, with Blastocladiomycota and Neocallimastigomycota representing segregate phyla of flagellated Fungi. Taxa traditionally placed in Zygomycota are distributed among Glomeromycota and several subphyla incertae sedis, including Mucoromycotina, Entomophthoromycotina, Kickxellomycotina, and Zoopagomycotina. Microsporidia are included in the Fungi, but no further subdivision of the group is proposed. Several genera of 'basal' Fungi of uncertain position are not placed in any higher taxa, including Basidiobolus, Caulochytrium, Olpidium, and Rozella.

Evolution of Pleopsidium (Lichenized Ascomycota) S943 Group I Introns and the Phylogeography of an Intron-Encoded Putative Homing Endonuclease

The sporadic distribution of nuclear group I introns among different fungal lineages can be explained by vertical inheritance of the introns followed by successive losses, or horizontal transfers from one lineage to another through intron homing or reverse splicing. Homing is mediated by an intron-encoded homing endonuclease (HE) and recent studies suggest that the introns and their associated HE gene (HEG) follow a recurrent cyclical model of invasion, degeneration, loss, and reinvasion. The purpose of this study was to compare this model to the evolution of HEGs found in the group I intron at position S943 of the nuclear ribosomal DNA of the lichen-forming fungus Pleopsidium. Forty-eight S943 introns were found in the 64 Pleopsidium samples from a worldwide screen, 22 of which contained a full-length HEG that encodes a putative 256-amino acid HE, and 2 contained HE pseudogenes. The HEGs are divided into two closely related types (as are the introns that encode them) that differ by 22.6% in their nucleotide sequences. The evolution of the Pleopsidium intron-HEG element shows strong evidence for a cyclical model of evolution. The intron was likely acquired twice in the genus and then transmitted via two or three interspecific horizontal transfers. Close geographical proximity plays an important role in intron-HEG horizontal transfer because most of these mobile elements were found in Europe. Once acquired in a lineage, the intron-HEG element was also vertically transmitted, and occasionally degenerated or was lost.

New insights into classification and evolution of the Lecanoromycetes (Pezizomycotina, Ascomycota) from phylogenetic analyses of three ribosomal RNA- and two protein-coding genes

The Lecanoromycetes includes most of the lichen-forming fungal species (> 13500) and is therefore one of the most diverse class of all Fungi in terms of phenotypic complexity. We report phylogenetic relationships within the Lecanoromycetes resulting from Bayesian and maximum likelihood analyses with complementary posterior probabilities and bootstrap support values based on three combined multilocus datasets using a supermatrix approach. Nine of 10 orders and 43 of 64 families currently recognized in Eriksson's classification of the Lecanoromycetes (Outline of Ascomycota--2006 Myconet 12:1-82) were represented in this sampling. Our analyses strongly support the Acarosporomycetidae and Ostropomycetidae as monophyletic, whereas the delimitation of the largest subclass, the Lecanoromycetidae, remains uncertain. Independent of future delimitation of the Lecanoromycetidae, the Rhizocarpaceae and Umbilicariaceae should be elevated to the ordinal level. This study shows that recent classifications include several nonmonophyletic taxa at different ranks that need to be recircumscribed. Our phylogenies confirm that ascus morphology cannot be applied consistently to shape the classification of lichen-forming fungi. The increasing amount of missing data associated with the progressive addition of taxa resulted in some cases in the expected loss of support, but we also observed an improvement in statistical support for many internodes. We conclude that a phylogenetic synthesis for a chosen taxonomic group should include a comprehensive assessment of phylogenetic confidence based on multiple estimates using different methods and on a progressive taxon sampling with an increasing number of taxa, even if it involves an increasing amount of missing data.

Reconstructing the early evolution of Fungi using a six-gene phylogeny

The ancestors of fungi are believed to be simple aquatic forms with flagellated spores, similar to members of the extant phylum Chytridiomycota (chytrids). Current classifications assume that chytrids form an early-diverging clade within the kingdom Fungi and imply a single loss of the spore flagellum, leading to the diversification of terrestrial fungi. Here we develop phylogenetic hypotheses for Fungi using data from six gene regions and nearly 200 species. Our results indicate that there may have been at least four independent losses of the flagellum in the kingdom Fungi. These losses of swimming spores coincided with the evolution of new mechanisms of spore dispersal, such as aerial dispersal in mycelial groups and polar tube eversion in the microsporidia (unicellular forms that lack mitochondria). The enigmatic microsporidia seem to be derived from an endoparasitic chytrid ancestor similar to Rozella allomycis, on the earliest diverging branch of the fungal phylogenetic tree.

Phylogenetic analyses suggest reverse splicing spread of group I introns in fungal ribosomal DNA

BACKGROUND

Group I introns have spread into over 90 different sites in nuclear ribosomal DNA (rDNA) with greater than 1700 introns reported in these genes. These ribozymes generally spread through endonuclease-mediated intron homing. Another putative pathway is reverse splicing whereby a free group I intron inserts into a homologous or heterologous RNA through complementary base-pairing between the intron and exon RNA. Reverse-transcription of the RNA followed by general recombination results in intron spread. Here we used phylogenetics to test for reverse splicing spread in a taxonomically broadly sampled data set of fungal group I introns including 9 putatively ancient group I introns in the rDNA of the yeast-like symbiont Symbiotaphrina buchneri.

RESULTS

Our analyses reveal a complex evolutionary history of the fungal introns with many cases of vertical inheritance (putatively for the 9 introns in S. buchneri) and intron lateral transfer. There are several examples in which introns, many of which are still present in S. buchneri, may have spread through reverse splicing into heterologous rDNA sites. If the S. buchneri introns are ancient as we postulate, then group I intron loss was widespread in fungal rDNA evolution.

CONCLUSION

On the basis of these results, we suggest that the extensive distribution of fungal group I introns is at least partially explained by the reverse splicing movement of existing introns into ectopic rDNA sites.

Assembling the fungal tree of life: progress, classification, and evolution of subcellular traits

Based on an overview of progress in molecular systematics of the true fungi (Fungi/Eumycota) since 1990, little overlap was found among single-locus data matrices, which explains why no large-scale multilocus phylogenetic analysis had been undertaken to reveal deep relationships among fungi. As part of the project "Assembling the Fungal Tree of Life" (AFTOL), results of four Bayesian analyses are reported with complementary bootstrap assessment of phylogenetic confidence based on (1) a combined two-locus data set (nucSSU and nucLSU rDNA) with 558 species representing all traditionally recognized fungal phyla (Ascomycota, Basidiomycota, Chytridiomycota, Zygomycota) and the Glomeromycota, (2) a combined three-locus data set (nucSSU, nucLSU, and mitSSU rDNA) with 236 species, (3) a combined three-locus data set (nucSSU, nucLSU rDNA, and RPB2) with 157 species, and (4) a combined four-locus data set (nucSSU, nucLSU, mitSSU rDNA, and RPB2) with 103 species. Because of the lack of complementarity among single-locus data sets, the last three analyses included only members of the Ascomycota and Basidiomycota. The four-locus analysis resolved multiple deep relationships within the Ascomycota and Basidiomycota that were not revealed previously or that received only weak support in previous studies. The impact of this newly discovered phylogenetic structure on supraordinal classifications is discussed. Based on these results and reanalysis of subcellular data, current knowledge of the evolution of septal features of fungal hyphae is synthesized, and a preliminary reassessment of ascomal evolution is presented. Based on previously unpublished data and sequences from GenBank, this study provides a phylogenetic synthesis for the Fungi and a framework for future phylogenetic studies on fungi.

Contribution of RPB2 to multilocus phylogenetic studies of the euascomycetes (Pezizomycotina, Fungi) with special emphasis on the lichen-forming Acarosporaceae and evolution of polyspory

Despite the recent progress in molecular phylogenetics, many of the deepest relationships among the main lineages of the largest fungal phylum, Ascomycota, remain unresolved. To increase both resolution and support on a large-scale phylogeny of lichenized and non-lichenized ascomycetes, we combined the protein coding-gene RPB2 with the traditionally used nuclear ribosomal genes SSU and LSU. Our analyses resulted in the naming of the new subclasses Acarosporomycetidae and Ostropomycetidae, and the new class Lichinomycetes, as well as the establishment of the phylogenetic placement and novel circumscription of the lichen-forming fungi family Acarosporaceae. The delimitation of this family has been problematic over the past century, because its main diagnostic feature, true polyspory (numerous spores issued from multiple post-meiosis mitoses) with over 100 spores per ascus, is probably not restricted to the Acarosporaceae. This observation was confirmed by our reconstruction of the origin and evolution of this form of true polyspory using maximum likelihood as the optimality criterion. The various phylogenetic analyses carried out on our data sets allowed us to conclude that: (1) the inclusion of phylogenetic signal from ambiguously aligned regions into the maximum parsimony analyses proved advantageous in reconstructing phylogeny; however, when more data become available, Bayesian analysis using different models of evolution is likely to be more efficient; (2) neighbor-joining bootstrap proportions seem to be more appropriate in detecting topological conflict between data partitions of large-scale phylogenies than posterior probabilities; and (3) Bayesian bootstrap proportion provides a compromise between posterior probability outcomes (i.e., higher accuracy, but with a higher number of significantly supported wrong internodes) vs. maximum likelihood bootstrap proportion outcomes (i.e., lower accuracy, with a lower number of significantly supported wrong internodes).

The evolution of homing endonuclease genes and group I introns in nuclear rDNA

Group I introns are autonomous genetic elements that can catalyze their own excision from pre-RNA. Understanding how group I introns move in nuclear ribosomal (r)DNA remains an important question in evolutionary biology. Two models are invoked to explain group I intron movement. The first is termed homing and results from the action of an intron-encoded homing endonuclease that recognizes and cleaves an intronless allele at or near the intron insertion site. Alternatively, introns can be inserted into RNA through reverse splicing. Here, we present the sequences of two large group I introns from fungal nuclear rDNA, which both encode putative full-length homing endonuclease genes (HEGs). Five remnant HEGs in different fungal species are also reported. This brings the total number of known nuclear HEGs from 15 to 22. We determined the phylogeny of all known nuclear HEGs and their associated introns. We found evidence for intron-independent HEG invasion into both homologous and heterologous introns in often distantly related lineages, as well as the "switching" of HEGs between different intron peripheral loops and between sense and antisense strands of intron DNA. These results suggest that nuclear HEGs are frequently mobilized. HEG invasion appears, however, to be limited to existing introns in the same or neighboring sites. To study the intron-HEG relationship in more detail, the S943 group I intron in fungal small-subunit rDNA was used as a model system. The S943 HEG is shown to be widely distributed as functional, inactivated, or remnant ORFs in S943 introns.

Major fungal lineages are derived from lichen symbiotic ancestors

About one-fifth of all known extant fungal species form obligate symbiotic associations with green algae, cyanobacteria or with both photobionts. These symbioses, known as lichens, are one way for fungi to meet their requirement for carbohydrates. Lichens are widely believed to have arisen independently on several occasions, accounting for the high diversity and mixed occurrence of lichenized and non-lichenized (42 and 58%, respectively) fungal species within the Ascomycota. Depending on the taxonomic classification chosen, 15-18 orders of the Ascomycota include lichen-forming taxa, and 8-11 of these orders (representing about 60% of the Ascomycota species) contain both lichenized and non-lichenized species. Here we report a phylogenetic comparative analysis of the Ascomycota, a phylum that includes greater than 98% of known lichenized fungal species. Using a Bayesian phylogenetic tree sampling methodology combined with a statistical model of trait evolution, we take into account uncertainty about the phylogenetic tree and ancestral state reconstructions. Our results show that lichens evolved earlier than believed, and that gains of lichenization have been infrequent during Ascomycota evolution, but have been followed by multiple independent losses of the lichen symbiosis. As a consequence, major Ascomycota lineages of exclusively non-lichen-forming species are derived from lichen-forming ancestors. These species include taxa with important benefits and detriments to humans, such as Penicillium and Aspergillus.

Widespread occurrence of spliceosomal introns in the rDNA genes of ascomycetes

Spliceosomal (pre-mRNA) introns have previously been found in eukaryotic protein-coding genes, in the small nuclear RNAs of some fungi, and in the small- and large-subunit ribosomal DNA genes of a limited number of ascomycetes. How the majority of these introns originate remains an open question because few proven cases of recent and pervasive intron origin have been documented. We report here the widespread occurrence of spliceosomal introns (69 introns at 27 different sites) in the small- and large-subunit nuclear-encoded rDNA of lichen-forming and free-living members of the Ascomycota. Our analyses suggest that these spliceosomal introns are of relatively recent origin, i.e., within the Euascomycetes, and have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into rRNAs. The spliceosome itself, and not an external agent (e.g., transposable elements, group II introns), may have given rise to these introns. A nonrandom sequence pattern was found at sites flanking the rRNA spliceosomal introns. This pattern (AG-intron-G) closely resembles the proto-splice site (MAG-intron-R) postulated for intron insertions in pre-mRNA genes. The clustered positions of spliceosomal introns on secondary structures suggest that particular rRNA regions are preferred sites for insertion through reverse-splicing.

Integrating ambiguously aligned regions of DNA sequences in phylogenetic analyses without violating positional homology

Phylogenetic analyses of non-protein-coding nucleotide sequences such as ribosomal RNA genes, internal transcribed spacers, and introns are often impeded by regions of the alignments that are ambiguously aligned. These regions are characterized by the presence of gaps and their uncertain positions, no matter which optimization criteria are used. This problem is particularly acute in large-scale phylogenetic studies and when aligning highly diverged sequences. Accommodating these regions, where positional homology is likely to be violated, in phylogenetic analyses has been dealt with very differently by molecular systematists and evolutionists, ranging from the total exclusion of these regions to the inclusion of every position regardless of ambiguity in the alignment. We present a new method that allows the inclusion of ambiguously aligned regions without violating homology. In this three-step procedure, first homologous regions of the alignment containing ambiguously aligned sequences are delimited. Second, each ambiguously aligned region is unequivocally coded as a new character, replacing its respective ambiguous region. Third, each of the coded characters is subjected to a specific step matrix to account for the differential number of changes (summing substitutions and indels) needed to transform one sequence to another. The optimal number of steps included in the step matrix is the one derived from the pairwise alignment with the greatest similarity and the least number of steps. In addition to potentially enhancing phylogenetic resolution and support, by integrating previously nonaccessible characters without violating positional homology, this new approach can improve branch length estimations when using parsimony.