Humans have generated ecological and environmental disturbances, such as vineyards, across the globe. Disturbed environments create widespread and repeated selective pressures that can drive colonization and local adaptation in microbial species. We investigated the distribution of fermentative yeast species in vineyards compared to nearby arboreal habitats and measured their resistance to two commonly used vineyard antimicrobials, copper and sulfite. We analyzed 4 101 strains, representing 70 species, collected from grapevine- and oak-associated substrates at 17 vineyard and 20 non-vineyard sites in the USA and Slovenia. Species frequency varied with geography and substrate, but the majority of species commonly present in vineyards were also found in non-vineyard arboreal environments, representing a potential source for vineyard colonization and exploitation of sugar from grapes. Species varied in both copper and sulfite resistance, but only Saccharomyces cerevisiae showed elevated resistance in vineyard compared to non-vineyard samples. Our results indicate that S. cerevisiae has uniquely taken advantage of vineyard environments through adaptations that appear either unnecessary or inaccessible to other yeast species present in vineyards.

Humans have generated ecological and environmental disturbances, such as vineyards, across the globe. Disturbed environments create widespread and repeated selective pressures that can drive colonization and local adaptation in microbial species. We investigated the distribution of fermentative yeast species in vineyards compared to nearby arboreal habitats and measured their resistance to two commonly used vineyard antimicrobials, copper and sulfite. We analyzed 4,101 strains, representing 70 species, collected from grapevine- and oak-associated substrates at 17 vineyard and 20 non-vineyard sites in the USA and Slovenia. Species frequency varied with geography and substrate, but the majority of species commonly present in vineyards were also found in non-vineyard arboreal environments, representing a potential source for vineyard colonization and exploitation of sugar from grapes. Species varied in both copper and sulfite resistance, but only Saccharomyces cerevisiae showed elevated resistance in vineyard compared to non-vineyard samples. Our results indicate that S. cerevisiae has uniquely taken advantage of vineyard environments through adaptations that appear either unnecessary or inaccessible to other yeast species present in vineyards.

ABSTRACT Evolutionary compromises are thought to be common under fluctuating selection because the mutations that best enable adaptation to one environmental context can often be detrimental in others. Yet, prior experimental work has shown that generalists can sometimes perform as well as specialists in their own environments. Here we use a highly replicated evolutionary experiment (N = 448 asexual lineages of the brewer’s yeast) to show that even though fluctuation between two environmental conditions often induces evolutionary compromises (at least early on), it can also help reveal difficult to reach adaptive outcomes that ultimately improve performance in both environments. Specifically, we begin by showing that yeast adaptation to chemical stress can involve fitness tradeoffs with stress-free environments and that, accordingly, lineages that are repeatedly exposed to occasional stress tend to respond by trading performance for breadth of adaptation. We then show that on rare occasions, fluctuating selection leads to the evolution of no-cost generalists that can even outcompete constant selection specialists in their own environments. We propose that the discovery of these broader and more effective adaptive outcomes under fluctuating selection could be partially facilitated by changes in the adaptive landscape that result from having to deal with fitness tradeoffs across different environmental conditions. Overall, our findings indicate that reconciling the short- and long-term evolutionary consequences of fluctuating selection could significantly improve our understanding of the evolution of specialization and generalism.

Head and neck squamous cell carcinoma (HNSCC) includes a subset of cancers driven by human papillomavirus (HPV). Here we use single-cell RNA-seq to profile both HPV-positive and HPV-negative oropharyngeal tumors, uncovering a high level of cellular diversity within and between tumors. First, we detect diverse chromosomal aberrations within individual tumors, suggesting genomic instability and enabling the identification of malignant cells even at pathologically negative margins. Second, we uncover diversity with respect to HNSCC subtypes and other cellular states such as the cell cycle, senescence and epithelial-mesenchymal transitions. Third, we find heterogeneity in viral gene expression within HPV-positive tumors. HPV expression is lost or repressed in a subset of cells, which are associated with a decrease in HPV-associated cell cycle phenotypes, decreased response to treatment, increased invasion and poor prognosis. These findings suggest that HPV expression diversity must be considered during diagnosis and treatment of HPV-positive tumors, with important prognostic ramifications.

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

Background: Improved prognostication of oropharyngeal squamous cell carcinoma (OPSCC) may facilitate individualized patient management. The goal of this study was to develop and validate a prognostic signature based on microRNA sequencing (miRNA-seq) analysis. Methods: We collected tumor specimens for miRNA-seq analysis from OPSCC patients treated at Washington University in St Louis (n = 324) and Vanderbilt University (n = 130). OPSCC patients (n = 79) from The Cancer Genome Atlas Program were also included for independent validation. Univariate and multivariable Cox regression analyses were performed to identify miRNAs associated with disease outcomes. All statistical tests were 2-sided. Results: By miRNA-seq profiling analysis, we identified a 26-miRNA signature. Based on computed risk scores of the signature, we classified the patients into low- and high-risk groups. In the training cohort, the high-risk group had much shorter overall survival compared with the low-risk group (hazard ratio [HR] = 3.80, 95% confidence interval [CI] = 2.37 to 6.10, P < .001). Subgroup analysis further revealed that the signature was prognostic for HPV-positive OPSCCs (HR = 3.07, 95% CI = 1.65 to 5.71, P < .001). Multivariable analysis indicated that the signature was independent of common clinicopathologic factors for OPSCCs. Importantly, the miRNA signature was a statistically significant predictor of overall survival in independent validation cohorts (The Cancer Genome Atlas Program cohort: HR = 6.05, 95% CI = 2.10 to 17.37, P < .001; Vanderbilt cohort: HR = 7.98, 95% CI = 3.99 to 15.97, P < .001; Vanderbilt HPV-positive cohort: HR = 8.71, 95% CI = 2.70 to 28.14, P < .001). Conclusions: The miRNA signature is a robust and independent prognostic tool for risk stratification of OPSCCs including HPV-positive OPSCCs.

Increasing evidence has elucidated the clinicopathological significance of tumor microenvironment (TME) cells. However, TME differences associated with human papillomavirus (HPV) infection in oropharyngeal squamous cell carcinoma (OPSCC) have not been well characterized. In our study, we comprehensively determined the TME infiltration patterns in 315 OPSCC patients, and systematically correlated the TME phenotypes with genomic characteristics and clinical features of OPSCCs. In this way, we observed the enrichment of high endothelial cells and adaptive immune cells in HPV-positive (HPV+) OPSCCs, in contrast to the enrichment of fibroblasts and capillary endothelial cells in HPV- negative (HPV-) OPSCCs. By focusing on immune checkpoint genes, we constructed a coexpression network using genes that were differentially expressed between HPV+ and HPV- OPSCCs. Functional analysis of the network indicated that HPV+ OPSCCs had elevated immune activities by promoting adaptive immune response and suppressing activities related to extracellular matrix organization. Subsequently, clinical analysis showed that identified TME-relevant genes were closely associated with the prognosis and therapy response in OPSCC. Importantly, results from the TME analysis were further validated using an independent OPSCC cohort.

Oropharyngeal squamous cell carcinoma (SCC) is increasing in incidence and, in Western countries, strongly associated with transcriptionally-active high-risk human papillomavirus (HPV). Within HPV-positive tumors, there is wide morphologic diversity with numerous histologic subtypes of SCC. There are also variable degrees of keratinization, anaplasia, stromal fibrosis, and maturing squamous differentiation. Unlike in the uterine cervix, where associations between HPV types and lineages/sublineages within types have been investigated with some clear correlations identified, little to no data exists for oropharyngeal SCC. In this study, for a large cohort of oropharyngeal SCC patients, we performed RTPCR for high-risk HPV. For the HPV positive patients, we sequenced the DNA of the entire HPV16 genome and determined lineages and sublineages, correlating HPV status, genotype, and HPV16 lineages/sublineages with SCC subtype and various histologic features. Of the 259 patients, 224 (86.5%) were high-risk HPV positive, of which 210/224 (93.8%) were HPV type 16 and 6/224 (2.7%) HPV type 33. Of the four HPV16 lineages, A was the most frequent (192/214 or 89.8%) and of the HPV16 A sublineages, A1 was the most frequent (112/210 or 53.3%). Patients with HPV negative tumors were more often keratinizing vs other types (23/35 or 65.7%) and thus more likely to have more maturing squamous differentiation and stromal desmoplasia. There was no significant correlation between HPV type (16 versus other), between HPV16 lineage (A versus others), or HPV16 A sublineages (A1 or A2 versus others) and morphologic type of SCC nor the various morphologic features of anaplasia/multinucleation, degree of keratinization, nor amount of stromal desmoplasia. In summary, in our cohort, there was no correlation between the type of HPV, the HPV 16 lineage or sublineage, and any of the histologic features or morphologic SCC subtypes.

Early studies estimate that 5% to 10% of oropharyngeal squamous cell carcinomas overexpress p16 but are unassociated with transcriptionally-active high-risk human papillomavirus (HPV). Patients with discordant HPV testing may experience clinical outcomes that differ from traditional expectations. To document the rate of p16 and HPV mRNA positivity, characterize patients with discordant testing, and identify features that may warrant selective use of HPV-specific testing after p16 IHC, a multi-institutional, retrospective review of oropharyngeal squamous cell carcinoma patients with p16 IHC and HPV mRNA testing by reverse transcriptase polymerase chain reaction was performed. Of the 467 patients, most had T1 or T2 tumors (71%), 82% were p16 positive, and 84% were HPV mRNA positive. Overall, most tumors were nonkeratinizing (378, 81%), which was strongly associated with p16 and HPV positivity (93% and 95%, respectively). Overall, 81% of patients were double positive, 14% double negative, and 4.9% discordant (3.4% p16 negative/HPV mRNA positive and 1.5% p16 positive/HPV mRNA negative). The survival rates of these discordant patient groups fell squarely between the 2 concordant groups, although in multivariate analysis for both disease-free survival and overall survival, discordant patients were not found to have statistically significantly different outcomes. Reclassifying patients by applying HPV mRNA testing when p16 results and morphology do not match, or when p16 results are equivocal, improved prognostication slightly over p16 or HPV mRNA testing alone. Patients with discordant testing demonstrate a borderline significant trend toward survival differences from those with concordant tests. When evaluated independently, patients who were p16 negative but HPV mRNA positive had a prognosis somewhat closer to double-positive patients, while those who were p16 positive, but HPV mRNA negative had a prognosis closer to that of double-negative patients. We suggest an algorithm whereby confirmatory HPV mRNA testing is performed in patients where p16 status is not consistent with tumor morphology. This captures a majority of discordant patients and improves, albeit modestly, the prognostication.

Extracellular RNAs (exRNAs) have attracted great attention due to their essential role in cell-to-cell communication as well as their potential as non-invasive disease biomarkers. However, at present, there is no consensus on the best method to profile exRNA expression, which leads to significant variability across studies. To address this issue, we established an experimental pipeline for comprehensive profiling of small exRNAs isolated from cell culture. By evaluating six RNA extraction protocols, we developed an improved method for robust recovery of vesicle-bound exRNAs. With this method, we performed small RNA sequencing of exosomes (EXOs), microvesicles (MVs) and source cells from 14 cancer cell lines. Compared to cells, EXOs and MVs were similarly enriched in tRNAs and rRNAs, but depleted in snoRNAs. By miRNA profiling analysis, we identified a subset of miRNAs, most noticeably miR-122-5p, that were significantly over-represented in EXOs and MVs across all 14 cell lines. In addition, we also identified a subset of EXO miRNAs associated with cancer type or human papillomavirus (HPV) status, suggesting their potential roles in HPV-induced cancers. In summary, our work has laid a solid foundation for further standardization on exRNA analysis across various cellular systems.

Background: Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue.Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies.We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge. Results:We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays.Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations.In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics. Conclusions:Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking.These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution.The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Background Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue. Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies. We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge. Results We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays. Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations. In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics. Conclusions Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking. These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution. The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Background: Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue.Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies.We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge. Results:We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays.Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations.In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics. Conclusions:Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking.These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution.The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Strains of Saccharomyces cerevisiae used to make beer, bread, and wine are genetically and phenotypically distinct from wild populations associated with trees. The origins of these domesticated populations are not always clear; human-associated migration and admixture with wild populations have had a strong impact on S. cerevisiae population structure. We examined the population genetic history of beer strains and found that ale strains and the S. cerevisiae portion of allotetraploid lager strains were derived from admixture between populations closely related to European grape wine strains and Asian rice wine strains. Similar to both lager and baking strains, ale strains are polyploid, providing them with a passive means of remaining isolated from other populations and providing us with a living relic of their ancestral hybridization. To reconstruct their polyploid origin, we phased the genomes of two ale strains and found ale haplotypes to both be recombinants between European and Asian alleles and to also contain novel alleles derived from extinct or as yet uncharacterized populations. We conclude that modern beer strains are the product of a historical melting pot of fermentation technology.

Migraine and other primary headache disorders affect a large population and cause debilitating pain. Establishing animal models that display behavioral correlates of long-lasting and ongoing headache, the most common and disabling symptom of migraine, is vital for the elucidation of disease mechanisms and identification of drug targets. We have developed a mouse model of headache, using dural application of capsaicin along with a mixture of inflammatory mediators (IScap) to simulate the induction of a headache episode. This elicited intermittent head-directed wiping and scratching as well as the phosphorylation of c-Jun N-terminal kinase in trigeminal ganglion neurons. Interestingly, dural application of IScap preferentially induced FOS protein expression in the excitatory but not inhibitory cervical/medullary dorsal horn neurons. The duration of IScap-induced behavior and the number of FOS-positive neurons correlated positively in individual mice; both were reduced to the control level by the pretreatment of antimigraine drug sumatriptan. Dural application of CGRP(8-37), the calcitonin gene-related peptide (CGRP) receptor antagonist, also effectively blocked IScap-induced behavior, which suggests that the release of endogenous CGRP in the dura is necessary for IScap-induced nociception. These data suggest that dural IScap-induced nocifensive behavior in mice may be mechanistically related to the ongoing headache in humans. In addition, dural application of IScap increased resting time in female mice. Taken together, we present the first detailed study using dural application of IScap in mice. This headache model can be applied to genetically modified mice to facilitate research on the mechanisms and therapeutic targets for migraine headache.

Saccharomyces cerevisiae and its sibling species S. paradoxus are known to inhabit temperate arboreal habitats across the globe. Despite their sympatric distribution in the wild, S. cerevisiae is predominantly associated with human fermentations. The apparent ecological differentiation of these species is particularly striking in Europe where S. paradoxus is abundant in forests and S. cerevisiae is abundant in vineyards. However, ecological differences may be confounded with geographic differences in species abundance. To compare the distribution and abundance of these two species we isolated Saccharomyces strains from over 1,200 samples taken from vineyard and forest habitats in Slovenia. We isolated numerous strains of S. cerevisiae and S. paradoxus as well as small number of S. kudriavzevii strains from both vineyard and forest environments. We find S. cerevisiae less abundant than S. paradoxus on oak trees within and outside the vineyard, but more abundant on grapevines and associated substrates. Analysis of the uncultured microbiome shows that both S. cerevisiae and S. paradoxus are rare species in soil and bark samples, but can be much more common in grape must. In contrast to S. paradoxus, European strains of S. cerevisiae have acquired multiple traits thought to be important for life in the vineyard and dominance of wine fermentations. We conclude that S. cerevisiae and S. paradoxus currently share both vineyard and non-vineyard habitats in Slovenia and we discuss factors relevant to their global distribution and relative abundance.

In budding yeasts, fermentation in the presence of oxygen evolved around the time of a whole genome duplication (WGD) and is thought to confer dominance in high-sugar environments because ethanol is toxic to many species. While there are many fermentative yeast species, only Saccharomyces cerevisiae consistently dominates wine fermentations. In this study, we use co-culture experiments and intrinsic growth rate assays to examine the relative fitness of non-WGD and WGD yeast species across environments to assess when S. cerevisiae's ability to dominate high-sugar environments arose. We show that S. cerevisiae dominates nearly all other non-WGD and WGD species except for its sibling species S. paradoxus in both grape juice and a high-sugar rich medium. Of the species we tested, S. cerevisiae and S. paradoxus have evolved the highest ethanol tolerance and intrinsic growth rate in grape juice. However, the ability of S. cerevisiae and S. paradoxus to dominate certain species depends on the temperature and the type of high-sugar environment. Our results indicate that dominance of high-sugar environments evolved much more recently than the WGD, most likely just prior to or during the differentiation of Saccharomyces species, and that evolution of multiple traits contributes to S. cerevisiae's ability to dominate wine fermentations.

Recent genetic studies suggest that dysfunction of ion channels and transporters may contribute to migraine pathophysiology. A migraine-associated frameshift mutation in the TWIK-related spinal cord K+ (TRESK) channel results in nonfunctional channels. Moreover, mutant TRESK subunits exert a dominant-negative effect on whole cell TRESK currents and result in hyperexcitability of small-diameter trigeminal ganglion (TG) neurons, suggesting that mutant TRESK may increase the gain of the neuronal circuit underlying migraine headache. However, the nonmigraine-associated TRESK C110R variant exhibits the same effect on TRESK currents as the mutant subunits in Xenopus oocytes, suggesting that dysfunction of TRESK is not sufficient to cause migraine. Here, we confirmed that the C110R variant formed nonfunctional channels and exerted a dominant-negative effect on TRESK currents in HEK293T cells, similar to the migraine-associated mutant TRESK. To compare the functional consequences of TRESK mutations/variants in a more physiological setting, we expressed the mutant TRESK and the C110R variant in cultured mouse TG neurons and investigated their effects on background K+ currents and neuronal excitability. Both mutant TRESK and the C110R variant reduced the endogenous TRESK currents in TG neurons, but the effect of the C110R variant was significantly smaller. Importantly, only TG neurons expressing mutant TRESK subunits, but not those expressing the C110R variant, exhibited a significant increase in excitability. Thus only the migraine-associated TRESK mutation, but not the C110R variant, reduces the endogenous TRESK currents to a degree that affects TG excitability. Our results support a potential causal relationship between the frameshift TRESK mutation and migraine susceptibility.