The condition-dependent transcriptional landscape of Burkholderia pseudomallei.

Burkholderia pseudomallei (Bp), the causative agent of the often-deadly infectious disease melioidosis, contains one of the largest prokaryotic genomes sequenced to date, at 7.2 Mb with two large circular chromosomes (1 and 2). To comprehensively delineate the Bp transcriptome, we integrated whole-genome tiling array expression data of Bp exposed to >80 diverse physical, chemical, and biological conditions. Our results provide direct experimental support for the strand-specific expression of 5,467 Sanger protein-coding genes, 1,041 operons, and 766 non-coding RNAs. A large proportion of these transcripts displayed condition-dependent expression, consistent with them playing functional roles. The two Bp chromosomes exhibited dramatically different transcriptional landscapes--Chr 1 genes were highly and constitutively expressed, while Chr 2 genes exhibited mosaic expression where distinct subsets were expressed in a strongly condition-dependent manner. We identified dozens of cis-regulatory motifs associated with specific condition-dependent expression programs, and used the condition compendium to elucidate key biological processes associated with two complex pathogen phenotypes--quorum sensing and in vivo infection. Our results demonstrate the utility of a Bp condition-compendium as a community resource for biological discovery. Moreover, the observation that significant portions of the Bp virulence machinery can be activated by specific in vitro cues provides insights into Bp's capacity as an "accidental pathogen", where genetic pathways used by the bacterium to survive in environmental niches may have also facilitated its ability to colonize human hosts.

Widespread epistasis among beneficial genetic variants revealed by high-throughput genome editing.

The phenotypic effect of any genetic variant can be altered by variation at other genomic loci. Known as epistasis, these genetic interactions shape the genotype-phenotype map of every species, yet their origins remain poorly understood. To investigate this, we employed high-throughput genome editing to measure the fitness effects of 1,826 naturally polymorphic variants in four strains of Saccharomyces cerevisiae. About 31% of variants affect fitness, of which 24% have strain-specific fitness effects indicative of epistasis. We found that beneficial variants are more likely to exhibit genetic interactions and that these interactions can be mediated by specific traits such as flocculation ability. This work suggests that adaptive evolution will often involve trade-offs where a variant is only beneficial in some genetic backgrounds, potentially explaining why many beneficial variants remain polymorphic. In sum, we provide a framework to understand the factors influencing epistasis with single-nucleotide resolution, revealing widespread epistasis among beneficial variants.

Gene-by-environment interactions are pervasive among natural genetic variants.

Gene-by-environment (GxE) interactions, in which a genetic variant's phenotypic effect is condition specific, are fundamental for understanding fitness landscapes and evolution but have been difficult to identify at the single-nucleotide level. Although many condition-specific quantitative trait loci (QTLs) have been mapped, these typically contain numerous inconsequential variants in linkage, precluding understanding of the causal GxE variants. Here, we introduce BARcoded Cas9 retron precise parallel editing via homology (CRISPEY-BAR), a high-throughput precision genome editing strategy, and use it to map GxE interactions of naturally occurring genetic polymorphisms impacting yeast growth. We identified hundreds of GxE variants within condition-specific QTLs, revealing unexpected genetic complexity. Moreover, we found that 93.7% of non-neutral natural variants within ergosterol biosynthesis pathway genes showed GxE interactions, including many impacting antifungal drug resistance through diverse molecular mechanisms. In sum, our results suggest an extremely complex, context-dependent fitness landscape characterized by pervasive GxE interactions while also demonstrating massively parallel genome editing as an effective means for investigating this complexity.

Divergent patterns of selection on metabolite levels and gene expression.

BACKGROUND: Natural selection can act on multiple genes in the same pathway, leading to polygenic adaptation. For example, adaptive changes were found to down-regulate six genes involved in ergosterol biosynthesis-an essential pathway targeted by many antifungal drugs-in some strains of the yeast Saccharomyces cerevisiae. However, the impact of this polygenic adaptation on metabolite levels was unknown. Here, we performed targeted mass spectrometry to measure the levels of eight metabolites in this pathway in 74 yeast strains from a genetic cross. RESULTS: Through quantitative trait locus (QTL) mapping we identified 19 loci affecting ergosterol pathway metabolite levels, many of which overlap loci that also impact gene expression within the pathway. We then used the recently developed v-test, which identified selection acting upon three metabolite levels within the pathway, none of which were predictable from the gene expression adaptation. CONCLUSIONS: These data showed that effects of selection on metabolite levels were complex and not predictable from gene expression data. This suggests that a deeper understanding of metabolism is necessary before we can understand the impacts of even relatively straightforward gene expression adaptations on metabolic pathways.