A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation.

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae.

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

Sex determination through X-Y heterogamety in Salix nigra.

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

Genetic and physiological mechanisms of freezing tolerance in locally adapted populations of a winter annual.

PREMISE: Despite myriad examples of local adaptation, the phenotypes and genetic variants underlying such adaptive differentiation are seldom known. Recent work on freezing tolerance and local adaptation in ecotypes of Arabidopsis thaliana from Italy and Sweden provides an essential foundation for uncovering the genotype-phenotype-fitness map for an adaptive response to a key environmental stress. METHODS: We examined the consequences of a naturally occurring loss-of-function (LOF) mutation in an Italian allele of the gene that encodes the transcription factor CBF2, which underlies a major freezing-tolerance locus. We used four lines with a Swedish genetic background, each containing a LOF CBF2 allele. Two lines had introgression segments containing the Italian CBF2 allele, and two contained deletions created using CRISPR-Cas9. We used a growth chamber experiment to quantify freezing tolerance and gene expression before and after cold acclimation. RESULTS: Freezing tolerance was lower in the Italian (11%) compared to the Swedish (72%) ecotype, and all four experimental CBF2 LOF lines had reduced freezing tolerance compared to the Swedish ecotype. Differential expression analyses identified 10 genes for which all CBF2 LOF lines, and the IT ecotype had similar patterns of reduced cold responsive expression compared to the SW ecotype. CONCLUSIONS: We identified 10 genes that are at least partially regulated by CBF2 that may contribute to the differences in cold-acclimated freezing tolerance between the Italian and Swedish ecotypes. These results provide novel insight into the molecular and physiological mechanisms connecting a naturally occurring sequence polymorphism to an adaptive response to freezing conditions.

Sex-biased gene expression in flowers, but not leaves, reveals secondary sexual dimorphism in Populus balsamifera.

Because sexual dimorphism in plants is often less morphologically conspicuous than in animals, studies of sex-biased gene expression may provide a quantitative metric to better address their commonality, molecular pathways, consistency across tissues and taxa, and evolution. The presence of sex-biased gene expression in tissues other than the androecium or gynoecium, termed secondary sexual characters, suggests that these traits arose after the initial evolution of dioecy. Patterns of sequence evolution may provide evidence of positive selection that drove sexual specialization. We compared gene expression in male and female flowers and leaves of Populus balsamifera to assess the extent of sex-biased expression, and tested whether sex-biased genes exhibit elevated rates of protein evolution. Sex-biased expression was pervasive in floral tissue, but nearly absent in leaf tissue. Female-biased genes in flowers were associated with photosynthesis, whereas male-biased genes were associated with mitochondrial function. Sex-biased genes did not exhibit elevated rates of protein evolution, contrary to results from other studies in animals and plants. Our results suggest that the ecological and physiological constraints associated with the energetics of flowering, rather than sexual conflict, have probably shaped the differences in male and female gene expression in P. balsamifera.

Constraints Imposed by a Natural Landscape Override Offspring Fitness Effects to Shape Oviposition Decisions in Wild Forked Fungus Beetles.

Oviposition site decisions often maximize offspring fitness, but costs constraining choice can cause females to oviposit in poor developmental environments. It is unclear whether these constraints cumulatively outweigh offspring fitness to determine oviposition decisions in wild populations. Understanding how constraints shape oviposition in natural landscapes is a critical step toward revealing how maternal behavior influences fundamental phenomena like the evolution of specialization and the use of sink environments. Here, we used a genetic capture-recapture technique to reconstruct the oviposition decisions of individual females in a natural metapopulation of a beetle (Bolitotherus cornutus) that oviposits on three fungus species. We measured larval fitness-related traits (mass, development time, survival) on each fungus and compared the oviposition preferences of females in laboratory versus field tests. Larval fitness differed substantially among fungi, and females preferred a high-quality (high larval fitness) fungus in laboratory trials. However, females frequently laid eggs on the lowest-quality fungus in the wild. They preferred high-quality fungi when moving between oviposition sites, but this preference disappeared as the distance between sites increased and was inconsistent between study plots. Our results suggest that constraints on oviposition preferences in natural landscapes are sufficiently large to drive oviposition in poor developmental environments even when offspring fitness consequences are severe.

Genetic dissection of the Drosophila melanogaster female head transcriptome reveals widespread allelic heterogeneity.

Modern genetic mapping is plagued by the "missing heritability" problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.

Aneuploidy confers quantitative proteome changes and phenotypic variation in budding yeast.

Aneuploidy, referring here to genome contents characterized by abnormal numbers of chromosomes, has been associated with developmental defects, cancer and adaptive evolution in experimental organisms. However, it remains unresolved how aneuploidy impacts gene expression and whether aneuploidy could directly bring about phenotypic variation and improved fitness over that of euploid counterparts. Here we show, using quantitative mass spectrometry-based proteomics and phenotypic profiling, that levels of protein expression in aneuploid yeast strains largely scale with chromosome copy numbers, following the same trend as that observed for the transcriptome, and that aneuploidy confers diverse phenotypes. We designed a novel scheme to generate, through random meiotic segregation, 38 stable and fully isogenic aneuploid yeast strains with distinct karyotypes and genome contents between 1N and 3N without involving any genetic selection. Through quantitative growth assays under various conditions or in the presence of a panel of chemotherapeutic or antifungal drugs, we found that some aneuploid strains grew significantly better than euploid control strains under conditions suboptimal for the latter. These results provide strong evidence that aneuploidy directly affects gene expression at both the transcriptome and proteome levels and can generate significant phenotypic variation that could bring about fitness gains under diverse conditions. Our findings suggest that the fitness ranking between euploid and aneuploid cells is dependent on context and karyotype, providing the basis for the notion that aneuploidy can directly underlie phenotypic evolution and cellular adaptation.

Genes related to mitochondrial functions are differentially expressed in phosphine-resistant and -susceptible Tribolium castaneum.

BACKGROUND: Phosphine is a valuable fumigant to control pest populations in stored grains and grain products. However, recent studies indicate a substantial increase in phosphine resistance in stored product pests worldwide. RESULTS: To understand the molecular bases of phosphine resistance in insects, we used RNA-Seq to compare gene expression in phosphine-resistant and susceptible laboratory populations of the red flour beetle, Tribolium castaneum. Each population was evaluated as either phosphine-exposed or no phosphine (untreated controls) in triplicate biological replicates (12 samples total). Pairwise analysis indicated there were eight genes differentially expressed between susceptible and resistant insects not exposed to phosphine (i.e., basal expression) or those exposed to phopshine (>8-fold expression and 90 % C.I.). However, 214 genes were differentially expressed among all four treatment groups at a statistically significant level (ANOVA, p < 0.05). Increased expression of 44 cytochrome P450 genes was found in resistant vs. susceptible insects, and phosphine exposure resulted in additional increases of 21 of these genes, five of which were significant among all treatment groups (p < 0.05). Expression of two genes encoding anti-diruetic peptide was 2- to 8-fold reduced in phosphine-resistant insects, and when exposed to phosphine, expression was further reduced 36- to 500-fold compared to susceptible. Phosphine-resistant insects also displayed differential expression of cuticle, carbohydrate, protease, transporter, and many mitochondrial genes, among others. Gene ontology terms associated with mitochondrial functions (oxidation biological processes, monooxygenase and catalytic molecular functions, and iron, heme, and tetrapyyrole binding) were enriched in the significantly differentially expressed dataset. Sequence polymorphism was found in transcripts encoding a known phosphine resistance gene, dihydrolipoamide dehydrogenase, in both susceptible and resistant insects. Phosphine-resistant adults also were resistant to knockdown by the pyrethroid deltamethrin, likely due to the increased cytochrome P450 expression. CONCLUSIONS: Overall, genes associated with the mitochondria were differentially expressed in resistant insects, and these differences may contribute to a reduction in overall metabolism and energy production and/or compensation in resistant insects. These data provide the first gene expression data on the response of phosphine-resistant and -susceptible insects to phosphine exposure, and demonstrate that RNA-Seq is a valuable tool to examine differences in insects that respond differentially to environmental stimuli.

Modification of gel architecture and TBE/TAE buffer composition to minimize heating during agarose gel electrophoresis.

Agarose gel electrophoresis of DNA and RNA is routinely performed using buffers containing either Tris, acetate, and EDTA (TAE) or Tris, borate, and EDTA (TBE). Gels are run at a low, constant voltage (∼10 V/cm) to minimize current and asymmetric heating effects, which can induce band artifacts and poor resolution. In this study, alterations of gel structure and conductive media composition were analyzed to identify factors causing higher electrical currents during horizontal slab gel electrophoresis. Current was reduced when thinner gels and smaller chamber buffer volumes were used, but was not influenced by agarose concentration or the presence of ethidium bromide. Current was strongly dependent on the amount and type of EDTA used and on the concentrations of the major acid-base components of each buffer. Interestingly, resolution and the mobilities of circular versus linear plasmid DNAs were also affected by the chemical form and amount of EDTA. With appropriate modifications to gel structure and buffer constituents, electrophoresis could be performed at high voltages (20-25 V/cm), reducing run times by up to 3-fold. The most striking improvements were observed with small DNAs and RNAs (10-100 bp): high voltages and short run times produced sharper bands and higher resolution.

Tolerability and pharmacokinetics of ACT-280778, a novel nondihydropyridine dual L/T-type calcium channel blocker: early clinical studies in healthy male subjects using adaptive designs.

ACT-280778 is a novel nondihydropyridine dual L/T-type calcium channel blocker. Two clinical studies (AC-067-101 and AC-067-102) were conducted to characterize its safety, tolerability, and pharmacokinetics in healthy male subjects after oral administration of single and multiple doses. Both trials were single-center, randomized, double-blind, placebo-controlled, adaptive design, ascending-dose studies, in which ACT-280778 was administrated as single doses of 2, 5, 15, or 40 mg, or as once-daily doses of 5 or 15 mg for 7 days. Single and multiple doses up to and including 15 mg were well tolerated, and no serious or severe adverse event was reported in either study. A single dose of 40 mg was associated with abnormal electrocardiogram findings resulting in the discontinuation of further treatment at this dose or higher doses. ACT-280778 was rapidly absorbed, and larger than dose-proportional increases of the maximum plasma concentration and area under the plasma concentration-time curve were observed. Food intake delayed the time to maximum plasma concentration and doubled exposure. Urinary excretion of unchanged ACT-280778 was negligible, and accumulation at steady state was modest. Overall, pharmacokinetic and tolerability profiles of ACT-280778 observed in these 2 studies warranted further evaluation of ACT-280778 in a proof-of-concept study in patients with hypertension.

Improving outpatient charge capture.

Hospitals can identify opportunities to enhance revenue collection by closely analyzing outpatient charge-capture data. A hospital can bolster its charge-capture analysis by performing a charge-capture process walk-through and scrutinizing subsystem links, third-party payer contracts, and electronic health record structures. The hospital then can integrate charge-integrity functions into clinical departments as needed by developing charge-reconciliation tools and reports and monitoring their utilization, and incorporating charge-reconciliation responsibilities into clinical department managers' job descriptions and goals.

Acute exposure to mercury drives changes in gene expression in Drosophila melanogaster.

OBJECTIVE: We quantified the effect of acute exposure to a high dosage of inorganic mercury on gene expression in Drosophila melanogaster using RNA-sequencing of whole adult females. RESULTS: We found 119 genes with higher gene expression following treatment (including all 5 Drosophila metallothionine genes and a number of heat shock protein genes), and 31 with lower expression (several of which are involved in egg formation). Our results highlight biological processes and genetic pathways impacted by exposure to this toxic metal, and provide motivation for future studies to understand the genetic basis of response to mercury.

Spatiotemporal profiling defines persistence and resistance dynamics during targeted treatment of melanoma.

Resistance of BRAF-mutant melanomas to targeted therapy arises from the ability of cells to enter a persister state, evade treatment with relative dormancy, and repopulate the tumor when reactivated. Using spatial transcriptomics in patient derived xenograft models, we capture clonal lineage evolution during treatment, finding the persister state to show increased oxidative phosphorylation, decreased proliferation, and increased invasive capacity, with central-to-peripheral gradients. Phylogenetic tracing identifies intrinsic- and acquired-resistance mechanisms (e.g. dual specific phosphatases, Reticulon-4, CDK2) and suggests specific temporal windows of potential therapeutic efficacy. Using deep learning to analyze histopathological slides, we find morphological features of specific cell states, demonstrating that juxtaposition of transcriptomics and histology data enables identification of phenotypically-distinct populations using imaging data alone. In summary, we define state change and lineage selection during melanoma treatment with spatiotemporal resolution, elucidating how choice and timing of therapeutic agents will impact the ability to eradicate resistant clones. Statement of Significance: Tumor evolution is accelerated by application of anti-cancer therapy, resulting in clonal expansions leading to dormancy and subsequently resistance, but the dynamics of this process are incompletely understood. Tracking clonal progression during treatment, we identify conserved, global transcriptional changes and local clone-clone and spatial patterns underlying the emergence of resistance.

A Pan-Cancer Patient-Derived Xenograft Histology Image Repository with Genomic and Pathologic Annotations Enables Deep Learning Analysis.

Patient-derived xenografts (PDX) model human intra- and intertumoral heterogeneity in the context of the intact tissue of immunocompromised mice. Histologic imaging via hematoxylin and eosin (H&E) staining is routinely performed on PDX samples, which could be harnessed for computational analysis. Prior studies of large clinical H&E image repositories have shown that deep learning analysis can identify intercellular and morphologic signals correlated with disease phenotype and therapeutic response. In this study, we developed an extensive, pan-cancer repository of >1,000 PDX and paired parental tumor H&E images. These images, curated from the PDX Development and Trial Centers Research Network Consortium, had a range of associated genomic and transcriptomic data, clinical metadata, pathologic assessments of cell composition, and, in several cases, detailed pathologic annotations of neoplastic, stromal, and necrotic regions. The amenability of these images to deep learning was highlighted through three applications: (i) development of a classifier for neoplastic, stromal, and necrotic regions; (ii) development of a predictor of xenograft-transplant lymphoproliferative disorder; and (iii) application of a published predictor of microsatellite instability. Together, this PDX Development and Trial Centers Research Network image repository provides a valuable resource for controlled digital pathology analysis, both for the evaluation of technical issues and for the development of computational image-based methods that make clinical predictions based on PDX treatment studies. Significance: A pan-cancer repository of >1,000 patient-derived xenograft hematoxylin and eosin-stained images will facilitate cancer biology investigations through histopathologic analysis and contributes important model system data that expand existing human histology repositories.

A multistep genomic screen identifies new genes required for repair of DNA double-strand breaks in Saccharomyces cerevisiae.

BACKGROUND: Efficient mechanisms for rejoining of DNA double-strand breaks (DSBs) are vital because misrepair of such lesions leads to mutation, aneuploidy and loss of cell viability. DSB repair is mediated by proteins acting in two major pathways, called homologous recombination and nonhomologous end-joining. Repair efficiency is also modulated by other processes such as sister chromatid cohesion, nucleosome remodeling and DNA damage checkpoints. The total number of genes influencing DSB repair efficiency is unknown. RESULTS: To identify new yeast genes affecting DSB repair, genes linked to gamma radiation resistance in previous genome-wide surveys were tested for their impact on repair of site-specific DSBs generated by in vivo expression of EcoRI endonuclease. Eight members of the RAD52 group of DNA repair genes (RAD50, RAD51, RAD52, RAD54, RAD55, RAD57, MRE11 and XRS2) and 73 additional genes were found to be required for efficient repair of EcoRI-induced DSBs in screens utilizing both MATa and MATα deletion strain libraries. Most mutants were also sensitive to the clastogenic chemicals MMS and bleomycin. Several of the non-RAD52 group genes have previously been linked to DNA repair and over half of the genes affect nuclear processes. Many proteins encoded by the protective genes have previously been shown to associate physically with each other and with known DNA repair proteins in high-throughput proteomics studies. A majority of the proteins (64%) share sequence similarity with human proteins, suggesting that they serve similar functions. CONCLUSIONS: We have used a genetic screening approach to detect new genes required for efficient repair of DSBs in Saccharomyces cerevisiae. The findings have spotlighted new genes that are critical for maintenance of genome integrity and are therefore of greatest concern for their potential impact when the corresponding gene orthologs and homologs are inactivated or polymorphic in human cells.

Inhibition of Wnt signaling by Wise (Sostdc1) and negative feedback from Shh controls tooth number and patterning.

Mice carrying mutations in Wise (Sostdc1) display defects in many aspects of tooth development, including tooth number, size and cusp pattern. To understand the basis of these defects, we have investigated the pathways modulated by Wise in tooth development. We present evidence that, in tooth development, Wise suppresses survival of the diastema or incisor vestigial buds by serving as an inhibitor of Lrp5- and Lrp6-dependent Wnt signaling. Reducing the dosage of the Wnt co-receptor genes Lrp5 and Lrp6 rescues the Wise-null tooth phenotypes. Inactivation of Wise leads to elevated Wnt signaling and, as a consequence, vestigial tooth buds in the normally toothless diastema region display increased proliferation and continuous development to form supernumerary teeth. Conversely, gain-of-function studies show that ectopic Wise reduces Wnt signaling and tooth number. Our analyses demonstrate that the Fgf and Shh pathways are major downstream targets of Wise-regulated Wnt signaling. Furthermore, our experiments revealed that Shh acts as a negative-feedback regulator of Wnt signaling and thus determines the fate of the vestigial buds and later tooth patterning. These data provide insight into the mechanisms that control Wnt signaling in tooth development and into how crosstalk among signaling pathways controls tooth number and morphogenesis.

Strategies for successful revenue cycle outsourcing.

Revenue cycle outsourcing can offer hospitals and health systems many advantages, including cost savings and revenue gains, but it also carries risks. Some organizations may choose to outsource revenue cycle to third-party service providers; others may opt to develop internal centers of excellence. Hospitals and health systems should consider IT system compatibility, payment arrangements, and incentive and value alignment when selecting an outsourcing partner.

Evolution of a ZW sex chromosome system in willows.

Transitions in the heterogamety of sex chromosomes (e.g., XY to ZW or vice versa) fundamentally alter the genetic basis of sex determination, however the details of these changes have been studied in only a few cases. In an XY to ZW transition, the X is likely to give rise to the W because they both carry feminizing genes and the X is expected to harbour less genetic load than the Y. Here, using a new reference genome for Salix exigua, we trace the X, Y, Z, and W sex determination regions during the homologous transition from an XY system to a ZW system in willow (Salix). We show that both the W and the Z arose from the Y chromosome. We find that the new Z chromosome shares multiple homologous putative masculinizing factors with the ancestral Y, whereas the new W lost these masculinizing factors and gained feminizing factors. The origination of both the W and Z from the Y was permitted by an unexpectedly low genetic load on the Y and this indicates that the origins of sex chromosomes during homologous transitions may be more flexible than previously considered.