Parallel evolution of linezolid-resistant Staphylococcus aureus in patients with cystic fibrosis.

Linezolid is an antibiotic used to treat serious Staphylococcus aureus infections. Resistance to linezolid is considered rare but could emerge with repeated dosing. We recently reported widespread prescription of linezolid for a cohort of patients with cystic fibrosis (CF). The goals of this study were to determine the incidence of linezolid-resistant methicillin-resistant Staphylococcus aureus (MRSA) in CF and identify molecular mechanisms for linezolid resistance. We identified patients who cultured S. aureus resistant to linezolid with minimum inhibitory concentration (MIC) >4 at the University of Iowa CF Center between 2008 and 2018. We obtained isolates from these patients and retested susceptibility to linezolid using broth microdilution. We used whole genome sequencing to perform phylogenetic analysis of linezolid-resistant isolates and examine sequences for mutations or accessory genes that confer linezolid resistance. Between 2008 and 2018, 111 patients received linezolid, and 4 of these patients cultured linezolid-resistant S. aureus. We sequenced 11 resistant and 21 susceptible isolates from these 4 subjects. Phylogenetic analysis indicated that linezolid resistance developed in ST5 or ST105 backgrounds. Three individuals had linezolid-resistant S. aureus with a G2576T mutation in 23S rRNA. One of these subjects additionally had a mutS- mutL- hypermutating S. aureus that produced five resistant isolates with multiple ribosomal subunit mutations. In one subject, the genetic basis for linezolid resistance was unclear. We conclude that linezolid resistant S. aureus can occur through multiple genetic mechanisms in patients with repeated exposure to this antibiotic. IMPORTANCE Patients with cystic fibrosis have persistent lung infections with Staphylococcus aureus that require extensive antibiotic treatments. Linezolid, an antibiotic given by oral or intravenous route, is prescribed repeatedly for patients whose lung disease has progressed. After treatment with linezolid, S. aureus strains can evolve antibiotic resistance through multiple genetic mechanisms. In addition to a common mutation in the 23S ribosomal RNA known to confer linezolid resistance, S. aureus strains can evolve novel resistance based on a combination of mutations affecting the bacterial ribosome. This combination of mutations was observed in a strain that exhibited hypermutation owing to the loss of the DNA repair genes mutS and mutL. In this cohort of patients with cystic fibrosis, linezolid resistance was transient, possibly due to the growth disadvantage of resistant strains. However, ongoing chronic exposure to linezolid may create optimal conditions for the future emergence of resistance to this critical antibiotic.

Parallel Evolution of Linezolid Resistant Staphylococcus aureus in Patients with Cystic Fibrosis.

Background: Linezolid is an antibiotic used to treat serious Staphylococcus aureus infections. Resistance to linezolid is considered rare but could emerge with repeated dosing. We recently reported widespread prescription of linezolid for a cohort of patients with cystic fibrosis (CF). Objectives: The goals of this study were to determine the incidence of linezolid resistance in CF and identify molecular mechanisms for linezolid resistance. Methods: We identified patients with S. aureus resistant to linezolid (MIC > 4) at the University of Iowa CF Center between 2008 and 2018. We obtained isolates from these patients and retested susceptibility to linezolid using broth microdilution. We used whole genome sequencing to perform phylogenetic analysis of linezolid resistant isolates and examine sequences for mutations or accessory genes that confer linezolid resistance. Main Results: Between 2008 and 2018, 111 patients received linezolid and 4 of these patients cultured linezolid resistant S. aureus . We sequenced 11 resistant and 21 susceptible isolates from these 4 subjects. Phylogenetic analysis indicated that linezolid resistance developed in ST5 or ST105 backgrounds. Three individuals had linezolid resistant S. aureus with a G2576T mutation in 23S rRNA. One of these subjects additionally had a mutS - mutL - hypermutating S. aureus that produced 5 resistant isolates with multiple ribosomal subunit mutations. In one subject, the genetic basis for linezolid resistance was unclear. Conclusions: Linezolid resistance evolved in 4 of 111 patients in this study. Linezolid resistance occurred by multiple genetic mechanisms. All resistant strains developed in ST5 or ST105 MRSA backgrounds. Key Point: Linezolid resistance arises through multiple genetic mechanisms and could be facilitated by mutator phenotypes. Linezolid resistance was transient, possibly due to growth disadvantage.

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 (FEV1 ) 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 FEV1 . ST8 MRSA was more common in younger patients and thus has the potential to rise in prevalence as these patients age.

Listeria monocytogenes Occurrence and Adherence to Recommendations: Small and Large Retail Delicatessens in Iowa.

The goal of this research was to determine the prevalence of Listeria monocytogenes in Iowa retail delicatessens and assess environmental aspects that mitigate L. monocytogenes. Fifty-seven small and large retail delicatessens in Iowa were selected randomly. More small operations (n = 43) were included as compared with larger stores, given the higher frequency of violations. An environmental assessment instrument was used to determine environmental factors and practices. At least five microbial samples were collected per site. We collected 286 (74.3%) of 385 microbial samples from small deli operations and 99 (25.7%) of 385 samples from large deli operations. Samples were taken from various zone 1 and 2 areas, such as the slicer, deli case, and meat scale; three (0.08%) samples were positive for L. monocytogenes. Regarding environmental aspects, not preparing, holding or storing ready-to-eat products near raw products (n = 30, 53%) was practiced by the fewest delis. The majority of establishments were observed covering, wrapping, or protecting ready-to-eat products when not in use to prevent contamination (n = 56, 98.2%). Comparisons were made to the U.S. Department of Agriculture Food Safety and Inspection Service Guidance on environmental practices, and 60% of the operations surveyed were in adherence with at least seven of the eight recommendations.

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.

Phylogenetic affiliations of members of the heterogeneous lichen-forming fungi of the genus Lecidea sensu Zahlbruckner (Lecanoromycetes, Ascomycota).

The genus Lecidea Ach. sensu lato (sensu Zahlbruckner) includes almost 1200 species, out of which only 100 species represent Lecidea sensu stricto (sensu Hertel). The systematic position of the remaining species is mostly unsettled but anticipated to represent several unrelated lineages within Lecanoromycetes. This study attempts to elucidate the phylogenetic placement of members of this heterogeneous group of lichen-forming fungi and to improve the classification and phylogeny of Lecanoromycetes. Twenty-five taxa of Lecidea sensu lato and 22 putatively allied species were studied in a broad selection of 268 taxa, representing 48 families of Lecanoromycetes. Six loci, including four ribosomal and two protein-coding genes for 315- and 209-OTU datasets were subjected to maximum likelihood and Bayesian analyses. The resulting well supported phylogenetic relationships within Lecanoromycetes are in agreement with published phylogenies, but the addition of new taxa revealed putative rearrangements of several families (e.g. Catillariaceae, Lecanoraceae, Lecideaceae, Megalariaceae, Pilocarpaceae and Ramalinaceae). As expected, species of Lecidea sensu lato and putatively related taxa are scattered within Lecanoromycetidae and beyond, with several species nested in Lecanoraceae and Pilocarpaceae and others placed outside currently recognized families in Lecanorales and orders in Lecanoromycetidae. The phylogenetic affiliations of Schaereria and Strangospora are outside Lecanoromycetidae, probably with Ostropomycetidae. All species referred to as Lecidea sensu stricto based on morphology (including the type species, Lecidea fuscoatra [L.] Ach.) form, with Porpidia species, a monophyletic group with high posterior probability outside Lecanorales, Peltigerales and Teloschistales, in Lecanoromycetidae, supporting the recognition of order Lecideales Vain. in this subclass. The genus name Lecidea must be redefined to apply only to Lecidea sensu stricto and to include at least some members of the genus Porpidia. Based on morphological and chemical similarities, as well as the phylogenetic relationship of Lecidea pullata sister to Frutidella caesioatra, the new combination Frutidella pullata is proposed here.

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.

Interrelationships of chromalveolates within a broadly sampled tree of photosynthetic protists.

The Chromalveolata "supergroup" is a massive assemblage of single-celled and multicellular protists such as ciliates and kelps that remains to be substantiated in molecular trees. Recent multigene analyses place chromalveolates into two major clades, the SAR (Stramenopiles, Alveolata, and Rhizaria) and the Cryptophyta+Haptophyta. Here we determined 69 new sequences from different chromalveolates to study the interrelationships of its constituent phyla. We included in our trees, the novel groups Telonemia and Katablepharidophyta that have previously been described as chromalvoleate allies. The best phylogenetic resolution resulted from a 6-protein (actin, alpha-tubulin, beta-tubulin, cytosolic HSP70, BIP HSP70, HSP90) and a 5-protein (lacking HSP90) alignment that validated the SAR and cryptophyte+haptophyte clades with the inclusion of telonemids in the former and katablepharids in the latter. We assessed the Plastidophila hypothesis that is based on EF2 data and suggest this grouping may be explained by horizontal gene transfers involving the EF2 gene rather than indicating host relationships.

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.

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.

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.

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.

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.

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).

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.

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.