Project ECHO and primary care buprenorphine treatment for opioid use disorder: Implementation and clinical outcomes.

Background: Our rural health system sought to (1) increase the number of primary care clinicians waivered to prescribe buprenorphine for treatment of opioid use disorder (OUD) and (2) consequently increase the number of our patients receiving this treatment. Methods: We used the Project for Extension for Community Health Outcomes (ECHO) tele-education model as an implementation strategy. We examined the number of clinicians newly waivered, the number of patients treated with buprenorphine, the relationship between clinician engagement with ECHO training and rates of buprenorphine prescribing, and treatment retention at 180 days. Results: The number of clinicians with a waiver and number of patients treated increased during and after ECHO training. There was a moderate correlation between the number of ECHO sessions attended by a clinician and number of their buprenorphine prescriptions (r = 0.50, p = 0.01). The 180-day retention rate was 80.7%. Conclusions: Project ECHO was highly effective for increasing access to this evidence-based treatment. The high retention rate in this rural context indicates that most patients are increasing their likelihood of favorable outcomes.

Author Correction: Hybridization is a recurrent evolutionary stimulus in wild yeast speciation.

The original version of the Supplementary Information associated with this Article contained errors in Supplementary Figures 2, 12, 20 and 22. The HTML has been updated to include a corrected version of the Supplementary Information; the original incorrect versions of these Figures can be found as Supplementary Information associated with this Correction.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species.

BACKGROUND: Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time. RESULTS: Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes. CONCLUSIONS: Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Hybridization is a recurrent evolutionary stimulus in wild yeast speciation.

Hybridization can result in reproductively isolated and phenotypically distinct lineages that evolve as independent hybrid species. How frequently hybridization leads to speciation remains largely unknown. Here we examine the potential recurrence of hybrid speciation in the wild yeast Saccharomyces paradoxus in North America, which comprises two endemic lineages SpB and SpC, and an incipient hybrid species, SpC*. Using whole-genome sequences from more than 300 strains, we uncover the hybrid origin of another group, SpD, that emerged from hybridization between SpC* and one of its parental species, the widespread SpB. We show that SpD has the potential to evolve as a novel hybrid species, because it displays phenotypic novelties that include an intermediate transcriptome profile, and partial reproductive isolation with its most abundant sympatric parental species, SpB. Our findings show that repetitive cycles of divergence and hybridization quickly generate diversity and reproductive isolation, providing the raw material for speciation by hybridization.

Persistence of Resident and Transplanted Genotypes of the Undomesticated Yeast Saccharomyces paradoxus in Forest Soil.

One might expect yeasts in soil to be highly dispersed via water or insects, forming ephemeral, genetically heterogeneous populations subject to competition and environmental stochasticity. Here, we report persistence of genotypes of the yeast Saccharomyces paradoxus in space and time. Within 1 km2 in a mixed hardwood forest on scales from centimeters to tens of meters, we detected persistence over 3 years of native genotypes, identified by single nucleotide polymorphisms (SNPs) genome-wide, of the wild yeast Saccharomyces paradoxus growing around Quercus rubra and Quercus alba Yeasts were recovered by enrichment in ethanol-containing medium, which measures only presence or absence, not abundance. Additional transplantation experiments employed strains marked with spontaneous defects in the URA3 gene, which also confer resistance to 5-fluoroorotic acid (5FOA). Plating soil suspensions from transplant sites on 5FOA-containing medium permitted one-step quantification of yeast CFU, with no interference from other unmarked yeasts or microorganisms. After an initial steep decrease in abundance, the yeast densities fluctuated over time, increasing in association with rainfall and decreasing in association with drought. After 18 months, the transplanted yeasts remained in place on the nine sites. In vitro transplantation experiments into nonsterile soil in petri dishes showed similar patterns of persistence and response to moisture and drought. To determine whether Saccharomyces cerevisiae, not previously recovered from soils regionally, can persist in our cold climate sites, we transplanted marked S. cerevisiae alone and in mixture with S. paradoxus in the fall of 2017. Five months later, S. cerevisiae persisted to the same extent as S. paradoxusIMPORTANCESaccharomyces yeasts are intensively studied in biological research and in their domesticated roles in brewing and baking, and yet, remarkably little is known about their mode of life in forest soils. We report here that resident genotypes of the yeast S. paradoxus are persistent on a time scale of years in their microhabitats in forest soils. We also show that resident genotypes can be replaced by transplanted yeast genotypes. The high inoculum levels in experimental transplantations rapidly decreased over time, but the transplanted genotypes persisted at low abundance. We conclude that, in forest soils, Saccharomyces yeasts exist at very low abundance and that dispersal events are rare.

Armillaria.

Lingering in forests around the world, some of the largest and oldest terrestrial organisms on earth hide in plain sight. In this Quick Guide, Sipos et al. shed light on the biology of the Armillaria fungi.

Author Correction: Genome expansion and lineage-specific genetic innovations in the forest pathogenic fungi Armillaria.

In the version of this Article originally published, it was incorrectly stated that "16,687 protein-coding genes were inferred for the most recent common ancestor (MRCA) of Armillaria"; the value was incorrect and it should have read "15,787". This has now been corrected.

Clonal evolution and genome stability in a 2500-year-old fungal individual.

Individuals of the basidiomycete fungus Armillaria are well known for their ability to spread from woody substrate to substrate on the forest floor through the growth of rhizomorphs. Here, we made 248 collections of A. gallica in one locality in Michigan's Upper Peninsula. To identify individuals, we genotyped collections with molecular markers and somatic compatibility testing. We found several different individuals in proximity to one another, but one genetic individual stood out as exceptionally large, covering hundreds of tree root systems over approximately 75 hectares of the forest floor. Based on observed growth rates of the fungus, we estimate the minimum age of the large individual as 2500 years. With whole-genome sequencing and variant discovery, we also found that mutation had occurred within the somatic cells of the individual, reflecting its historical pattern of growth from a single point. The overall rate of mutation over the 90 mb genome, however, was extremely low. This same individual was first discovered in the late 1980s, but its full spatial extent and internal mutation dynamic was unknown at that time. The large individual of A. gallica has been remarkably resistant to genomic change as it has persisted in place.

Genome expansion and lineage-specific genetic innovations in the forest pathogenic fungi Armillaria.

Armillaria species are both devastating forest pathogens and some of the largest terrestrial organisms on Earth. They forage for hosts and achieve immense colony sizes via rhizomorphs, root-like multicellular structures of clonal dispersal. Here, we sequenced and analysed the genomes of four Armillaria species and performed RNA sequencing and quantitative proteomic analysis on the invasive and reproductive developmental stages of A. ostoyae. Comparison with 22 related fungi revealed a significant genome expansion in Armillaria, affecting several pathogenicity-related genes, lignocellulose-degrading enzymes and lineage-specific genes expressed during rhizomorph development. Rhizomorphs express an evolutionarily young transcriptome that shares features with the transcriptomes of both fruiting bodies and vegetative mycelia. Several genes show concomitant upregulation in rhizomorphs and fruiting bodies and share cis-regulatory signatures in their promoters, providing genetic and regulatory insights into complex multicellularity in fungi. Our results suggest that the evolution of the unique dispersal and pathogenicity mechanisms of Armillaria might have drawn upon ancestral genetic toolkits for wood-decay, morphogenesis and complex multicellularity.

Fungus Causing White-Nose Syndrome in Bats Accumulates Genetic Variability in North America with No Sign of Recombination.

Emerging fungal diseases of wildlife are on the rise worldwide, and the white-nose syndrome (WNS) epidemic in North American bats is a catastrophic example. The causal agent of WNS is a single clone of the fungus Pseudogymnoascus destructans. Early evolutionary change in this clonal population has major implications for disease ecology and conservation. Accumulation of variation in the fungus through mutation, and shuffling of variation through recombination, could affect the virulence and transmissibility of the fungus and the durability of what appears to be resistance arising in some bat populations. Our genome-wide analysis shows that the clonal population of P. destructans has expanded in size from a single genotype, has begun to accumulate variation through mutation, and presents no evidence as yet of genetic exchange among individuals. IMPORTANCE Since its discovery in 2006, the emerging infectious disease known as white-nose syndrome has killed millions of bats in North America, making it one of the most devastating wildlife epidemics in recorded history. We demonstrate that there has been as yet only spontaneous mutation across the North American population of P. destructans, and we find no indication of recombination. Thus, selective forces, which might otherwise impact pathogenic virulence, have so far had essentially no genetic variation on which to act. Our study confirmed the time of origin for the first and, thus far, only introduction of P. destructans to North America. This system provides an unprecedented opportunity to follow the evolution of a host-pathogen interaction unfolding in real time.

Yeast Populations Evolve to Resist CdSe Quantum Dot Toxicity.

Engineered nanomaterials are used globally in biomedical, electronic, and optical devices, and are often discarded into the environment. Cell culture experiments have shown that many inorganic nanoparticles are toxic to eukaryotic cells. Here, we show that populations of eukaryotic cells can evolve to survive chronic exposure to toxic CdSe semiconductor quantum dots (QDs). We grew yeast Saccharomyces cerevisiae for 24 days in liquid medium containing QDs prepared daily at half the minimum inhibitory concentration (MIC50) of the progenitor yeast cells. After 24 days, the cells grew normally under constant exposure to QDs. We concluded that these cells evolved to resist QD toxicity. Surprisingly, when we removed QDs from the growth medium, some of the evolved cells grew poorly, i.e., they grew better in the presence of QDs. Finally, genetic analysis confirmed that the ubiquitin ligase gene bul1 was mutated in the evolved cells, which suggests that this gene may be implicated in increased CdSe QD tolerance. This study shows that chronic exposure to QDs can exert selective pressure causing irreversible genetic changes leading to adaptation.

Population genomics reveals structure at the individual, host-tree scale and persistence of genotypic variants of the undomesticated yeast Saccharomyces paradoxus in a natural woodland.

Genetic diversity in experimental, domesticated and wild populations of the related yeasts, Saccharomyces cerevisiae and Saccharomyces paradoxus, has been well described at the global scale. We investigated the population genomics of a local population on a small spatial scale to address two main questions. First, is there genomic variation in a S. paradoxus population at a spatial scale spanning centimetres (microsites) to tens of metres? Second, does the distribution of genomic variants persist over time? Our sample consisted of 42 S. paradoxus strains from 2014 and 43 strains from 2015 collected from the same 72 microsites around four host trees (Quercus rubra and Quercus alba) within 1 km2 in a mixed hardwood forest in southern Ontario. Six additional S. paradoxus strains recovered from adjacent maple and beech trees in 2015 are also included in the sample. Whole-genome sequencing and genomic SNP analysis revealed five differentiated groups (clades) within the sampled area. The signal of persistence of genotypes in their microsites from 2014 to 2015 was highly significant. Isolates from the same tree tended to be more related than strains from different trees, with limited evidence of dispersal between trees. In growth assays, one genotype had a significantly longer lag phase than the other strains. Our results indicate that different clades coexist at fine spatial scale and that population structure persists over at least a one-year interval in these wild yeasts, suggesting the efficacy of yearly sampling to follow longer term genetic dynamics in future studies.

Note: "Exact" quantum Monte Carlo calculations of the barrier for the H + H2 reaction at the sub-microhartree level.

We report the third in a series of 'exact' quantum Monte Carlo calculations for the potential energy of the saddle point of the barrier for the reaction H + H2 → H2 + H. The barrier heights determined are 9.61 ± 0.01 in 1992/94, 9.608 ± 0.001 in 2003, and 9.6089 ± 0.0001 in 2016 (this work), all in kcal/mole and successively a factor of ten more accurate. The new value is below the lowest value from explicitly correlated Gaussian calculations and within the estimated limits of extrapolated multireference configuration calculations.

A Genetic Incompatibility Accelerates Adaptation in Yeast.

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker's yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

The underlying structure of adaptation under strong selection in 12 experimental yeast populations.

The aims of this study were to determine (i) whether adaptation under strong selection occurred through mutations in a narrow target of one or a few nucleotide sites or a broad target of numerous sites and (ii) whether the programs of adaptation previously observed from three experimental populations were unique or shared among populations that underwent parallel evolution. We used archived population samples from a previous study, representing 500 generations of experimental evolution in 12 populations under strong selection, 6 populations in a high-salt environment and 6 populations in a low-glucose environment. Each set of six populations included four with sexual reproduction and two with exclusively asexual reproduction. Populations were sampled as resequenced genomes of 115 individuals and as bulk samples from which frequencies of mutant alleles were estimated. In a high-salt environment, a broad target of 11 mutations within the proton exporter, PMA1, was observed among the six populations, in addition to expansions of the ENA gene cluster. This pattern was shared among populations that underwent parallel evolution. In a low-glucose environment, two programs of adaptation were observed. The originally observed pattern of mutation in MDS3/MKT1 in population M8 was a narrow target of a single nucleotide, unique to this population. Among the other five populations, the three mutations were shared in a broad target, sensing/signaling genes RAS1 and RAS2. RAS1/RAS2 mutations were not observed in the high-salt populations; PMA1 mutations were observed only in a high-salt environment.

Genomewide mutation dynamic within a long-lived individual of Armillaria gallica.

Mutation is the ultimate source of all genetic variation in populations and yet they remain unobservable and buried in the past. Long-lived individuals of Armillaria gallica, a common opportunistic fungal pathogen of tree roots in temperate forests of the northern hemisphere, provide a spatial context for examining the mutational dynamic. Each individual of A. gallica arises in a single mating between two haploid gametes and the resulting diploid then grows vegetatively to occupy a discrete spatial territory often including many adjacent tree root systems. In effect, this leaves a spatial record of growth over time within which mutations can be localized. To identify mutations, the entire genomes of three spatially separated samples of one individual of A. gallica approximately 200 × 60 m were sequenced and compared. In this comparison, mutations and regions of loss of heterozygosity (LOH) were identified then assayed in another 20 isolates from the same individual by conventional PCR and Sanger sequencing. The genotype network of all mutations and LOH were without internal conflict. Further, the spatial pattern of genotypes was nonrandom and appeared to reflect the vegetative expansion leading to the present-day individual. The results reflect the spectrum of spontaneous mutation in nature and provide insight into cellular generation times.

Cellular effects and epistasis among three determinants of adaptation in experimental populations of Saccharomyces cerevisiae.

Epistatic interactions in which the phenotypic effect of an allele is conditional on its genetic background have been shown to play a central part in various evolutionary processes. In a previous study (J. B. Anderson et al., Curr. Biol. 20:1383-1388, 2010; J. R. Dettman, C. Sirjusingh, L. M. Kohn, and J. B. Anderson, Nature 447:585-588, 2007), beginning with a common ancestor, we identified three determinants of fitness as mutant alleles (each designated with the letter "e") that arose in replicate Saccharomyces cerevisiae populations propagated in two different environments, a low-glucose and a high-salt environment. In a low-glucose environment, MDS3e and MKT1e interacted positively to confer a fitness advantage. Also, PMA1e from a high-salt environment interacted negatively with MKT1e in a low-glucose environment, an example of a Dobzhansky-Muller incompatibility that confers reproductive isolation. Here we showed that the negative interaction between PMA1e and MKT1e is mediated by alterations in intracellular pH, while the positive interaction between MDS3e and MKT1e is mediated by changes in gene expression affecting glucose transporter genes. We specifically addressed the evolutionary significance of the positive interaction by showing that the presence of the MDS3 mutation is a necessary condition for the spread and fixation of the new mutations at the identical site in MKT1. The expected mutations in MKT1 rose to high frequencies in two of three experimental populations carrying MDS3e but not in any of three populations carrying the ancestral allele. These data show how positive and negative epistasis can contribute to adaptation and reproductive isolation.

Determinants of divergent adaptation and Dobzhansky-Muller interaction in experimental yeast populations.

Divergent adaptation can be associated with reproductive isolation in speciation [1]. We recently demonstrated the link between divergent adaptation and the onset of reproductive isolation in experimental populations of the yeast Saccharomyces cerevisiae evolved from a single progenitor in either a high-salt or a low-glucose environment [2]. Here, whole-genome resequencing and comparative genome hybridization of representatives of three populations revealed 17 mutations, six of which explained the adaptive increases in mitotic fitness. In two populations evolved in high salt, two different mutations occurred in the proton efflux pump gene PMA1 and the global transcriptional repressor gene CYC8; the ENA genes encoding sodium efflux pumps were overexpressed once through expansion of this gene cluster and once because of mutation in the regulator CYC8. In the population from low glucose, one mutation occurred in MDS3, which modulates growth at high pH, and one in MKT1, a global regulator of mRNAs encoding mitochondrial proteins, the latter recapitulating a naturally occurring variant. A Dobzhansky-Muller (DM) incompatibility between the evolved alleles of PMA1 and MKT1 strongly depressed fitness in the low-glucose environment. This DM interaction is the first reported between experimentally evolved alleles of known genes and shows how reproductive isolation can arise rapidly when divergent selection is strong.

Monte Carlo simulations of single- and multistep enzyme-catalyzed reaction sequences: effects of diffusion, cell size, enzyme fluctuations, colocalization, and segregation.

Following the discovery of slow fluctuations in the catalytic activity of an enzyme in single-molecule experiments, it has been shown that the classical Michaelis-Menten (MM) equation relating the average enzymatic velocity and the substrate concentration may hold even for slowly fluctuating enzymes. In many cases, the average velocity is that given by the MM equation with time-averaged values of the fluctuating rate constants and the effect of enzyme fluctuations is simply averaged out. The situation is quite different for a sequence of reactions. For colocalization of a pair of enzymes in a sequence to be effective in promoting reaction, the second must be active when the first is active or soon after. If the enzymes are slowly varying and only rarely active, the product of the first reaction may diffuse away before the second enzyme is active, and colocalization may have little value. Even for single-step reactions the interplay of reaction and diffusion with enzyme fluctuations leads to added complexities, but for multistep reactions the interplay of reaction and diffusion, cell size, compartmentalization, enzyme fluctuations, colocalization, and segregation is far more complex than for single-step reactions. In this paper, we report the use of stochastic simulations at the level of whole cells to explore, understand, and predict the behavior of single- and multistep enzyme-catalyzed reaction systems exhibiting some of these complexities. Results for single-step reactions confirm several earlier observations by others. The MM relationship, with altered constants, is found to hold for single-step reactions slowed by diffusion. For single-step reactions, the distribution of enzymes in a regular grid is slightly more effective than a random distribution. Fluctuations of enzyme activity, with average activity fixed, have no observed effects for simple single-step reactions slowed by diffusion. Two-step sequential reactions are seen to be slowed by segregation of the enzymes for each step, and results of the calculations suggest limits for cell size. Colocalization of enzymes for a two-step sequence is seen to promote reaction, and rates fall rapidly with increasing distance between enzymes. Low frequency fluctuations of the activities of colocalized enzymes, with average activities fixed, can greatly reduce reaction rates for sequential reactions.

Genome-wide investigation of reproductive isolation in experimental lineages and natural species of Neurospora: identifying candidate regions by microarray-based genotyping and mapping.

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe-dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky-Muller model of genic incompatibility. We developed microarray-based, restriction-site associated DNA (RAD) markers that identified approximately 1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk-segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only approximately 2% of the nuclear genome. This study demonstrates how advances in high-throughput, genome-wide genotyping can be applied to mapping reproductive isolation genes and speciation research.