The genome of a giant clam zooxanthella (Cladocopium infistulum) offers few clues to adaptation as an extracellular symbiont with high thermotolerance.

BACKGROUND: Cladocopium infistulum (Symbiodiniaceae) is a dinoflagellate specialized to live in symbiosis with western Pacific giant clams (Tridacnidae). Unlike coral-associated symbionts, which reside within the host cells, C. infistulum inhabits the extracellular spaces of the clam's digestive diverticula. It is phylogenetically basal to a large species complex of stress-tolerant Cladocopium, many of which are associated with important reef-building corals in the genus Porites. This close phylogenetic relationship may explain why C. infistulum exhibits high thermotolerance relative to other tridacnid symbionts. Moreover, past analyses of microsatellite loci indicated that Cladocopium underwent whole-genome duplication prior to the adaptive radiations that led to its present diversity. RESULTS: A draft genome assembly of C. infistulum was produced using long- and short-read sequences to explore the genomic basis for adaptations underlying thermotolerance and extracellular symbiosis among dinoflagellates and to look for evidence of genome duplication. Comparison to three other Cladocopium genomes revealed no obvious over-representation of gene groups or families whose functions would be important for maintaining C. infistulum's unique physiological and ecological properties. Preliminary analyses support the existence of partial or whole-genome duplication among Cladocopium, but additional high-quality genomes are required to substantiate these findings. CONCLUSION: Although this investigation of Cladocopium infistulum revealed no patterns diagnostic of heat tolerance or extracellular symbiosis in terms of overrepresentation of gene functions or genes under selection, it provided a valuable genomic resource for comparative analyses. It also indicates that ecological divergence among Cladocopium species, and potentially among other dinoflagellates, is partially governed by mechanisms other than gene content. Thus, additional high-quality, multiomic data are needed to explore the molecular basis of key phenotypes among symbiotic microalgae.

Environmental specialization and cryptic genetic divergence in two massive coral species from the Florida Keys Reef Tract.

Broadcast-spawning coral species have wide geographical ranges spanning strong environmental gradients, but it is unclear how much spatially varying selection these gradients actually impose. Strong divergent selection might present a considerable barrier for demographic exchange between disparate reef habitats. We investigated whether the cross-shelf gradient is associated with spatially varying selection in two common coral species, Montastraea cavernosa and Siderastrea siderea, in the Florida Keys. To this end, we generated a de novo genome assembly for M. cavernosa and used 2bRAD to genotype 20 juveniles and 20 adults of both species from each of the three reef zones to identify signatures of selection occurring within a single generation. Unexpectedly, each species was found to be composed of four genetically distinct lineages, with gene flow between them still ongoing but highly reduced in 13.0%-54.7% of the genome. Each species includes two sympatric lineages that are only found in the deep (20 m) habitat, while the other lineages are found almost exclusively on the shallower reefs (3-10 m). The two "shallow" lineages of M. cavernosa are also specialized for either nearshore or offshore: comparison between adult and juvenile cohorts indicates that cross-shelf migrants are more than twice as likely to die before reaching adulthood than local recruits. S. siderea and M. cavernosa are among the most ecologically successful species on the Florida Keys Reef Tract, and this work offers important insight into the genomic background of divergent selection and environmental specialization that may in part explain their resilience and broad environmental range.

Ancestral polymorphisms explain the role of chromosomal inversions in speciation.

Understanding the role of chromosomal inversions in speciation is a fundamental problem in evolutionary genetics. Here, we perform a comprehensive reconstruction of the evolutionary histories of the chromosomal inversions in Drosophila persimilis and D. pseudoobscura. We provide a solution to the puzzling origins of the selfish Sex-Ratio arrangement in D. persimilis and uncover surprising patterns of phylogenetic discordance on this chromosome. These patterns show that, contrary to widely held views, all fixed chromosomal inversions between D. persimilis and D. pseudoobscura were already present in their ancestral population long before the species split. Our results suggest that patterns of higher genomic divergence and an association of reproductive isolation genes with chromosomal inversions may be a direct consequence of incomplete lineage sorting of ancestral polymorphisms. These findings force a reconsideration of the role of chromosomal inversions in speciation, not as protectors of existing hybrid incompatibilities, but as fertile grounds for their formation.

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis.

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura.

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement-specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.

Genome-wide analysis of signatures of selection in populations of African honey bees (Apis mellifera) using new web-based tools.

BACKGROUND: With the development of inexpensive, high-throughput sequencing technologies, it has become feasible to examine questions related to population genetics and molecular evolution of non-model species in their ecological contexts on a genome-wide scale. Here, we employed a newly developed suite of integrated, web-based programs to examine population dynamics and signatures of selection across the genome using several well-established tests, including F ST, pN/pS, and McDonald-Kreitman. We applied these techniques to study populations of honey bees (Apis mellifera) in East Africa. In Kenya, there are several described A. mellifera subspecies, which are thought to be localized to distinct ecological regions. RESULTS: We performed whole genome sequencing of 11 worker honey bees from apiaries distributed throughout Kenya and identified 3.6 million putative single-nucleotide polymorphisms. The dense coverage allowed us to apply several computational procedures to study population structure and the evolutionary relationships among the populations, and to detect signs of adaptive evolution across the genome. While there is considerable gene flow among the sampled populations, there are clear distinctions between populations from the northern desert region and those from the temperate, savannah region. We identified several genes showing population genetic patterns consistent with positive selection within African bee populations, and between these populations and European A. mellifera or Asian Apis florea. CONCLUSIONS: These results lay the groundwork for future studies of adaptive ecological evolution in honey bees, and demonstrate the use of new, freely available web-based tools and workflows ( http://usegalaxy.org/r/kenyanbee ) that can be applied to any model system with genomic information.

Genomics of natural populations: gene conversion events reveal selected genes within the inversions of Drosophila pseudoobscura.

When adaptive phenotypic variation or quantitative trait loci map within an inverted segment of a chromosome, researchers often despair because the suppression of crossing over will prevent the discovery of selective target genes that established the rearrangement. If an inversion polymorphism is old enough, then the accumulation of gene conversion tracts offers the promise that quantitative trait loci or selected loci within inversions can be mapped. The inversion polymorphism of Drosophila pseudoobscura is a model system to show that gene conversion analysis is a useful tool for mapping selected loci within inversions. D. pseudoobscura has over 30 different chromosomal arrangements on the third chromosome (Muller C) in natural populations and their frequencies vary with changes in environmental habitats. Statistical tests of five D. pseudoobscura gene arrangements identified outlier genes within inverted regions that had potentially heritable variation, either fixed amino acid differences or differential expression patterns. We use genome sequences of the inverted third chromosome (Muller C) to infer 98,443 gene conversion tracts for a total coverage of 142 Mb or 7.2× coverage of the 19.7 Mb chromosome. We estimated gene conversion tract coverage in the 2,668 genes on Muller C and tested whether gene conversion coverage was similar among arrangements for outlier vs non-outlier loci. Outlier genes had lower gene conversion tract coverage among arrangements than the non-outlier genes suggesting that selection removes exchanged DNA in the outlier genes. These data support the hypothesis that the third chromosome in D. pseudoobscura captured locally adapted combinations of alleles prior to inversion mutation events.

Relating pathogenic loss-of-function mutations in humans to their evolutionary fitness costs.

Causal loss-of-function (LOF) variants for Mendelian and severe complex diseases are enriched in 'mutation intolerant' genes. We show how such observations can be interpreted in light of a model of mutation-selection balance and use the model to relate the pathogenic consequences of LOF mutations at present to their evolutionary fitness effects. To this end, we first infer posterior distributions for the fitness costs of LOF mutations in 17,318 autosomal and 679 X-linked genes from exome sequences in 56,855 individuals. Estimated fitness costs for the loss of a gene copy are typically above 1%; they tend to be largest for X-linked genes, whether or not they have a Y homolog, followed by autosomal genes and genes in the pseudoautosomal region. We compare inferred fitness effects for all possible de novo LOF mutations to those of de novo mutations identified in individuals diagnosed with one of six severe, complex diseases or developmental disorders. Probands carry an excess of mutations with estimated fitness effects above 10%; as we show by simulation, when sampled in the population, such highly deleterious mutations are typically only a couple of generations old. Moreover, the proportion of highly deleterious mutations carried by probands reflects the typical age of onset of the disease. The study design also has a discernible influence: a greater proportion of highly deleterious mutations is detected in pedigree than case-control studies, and for autism, in simplex than multiplex families and in female versus male probands. Thus, anchoring observations in human genetics to a population genetic model allows us to learn about the fitness effects of mutations identified by different mapping strategies and for different traits.

Large-scale multitrait genome-wide association analyses identify hundreds of glaucoma risk loci.

Glaucoma, a leading cause of irreversible blindness, is a highly heritable human disease. Previous genome-wide association studies have identified over 100 loci for the most common form, primary open-angle glaucoma. Two key glaucoma-associated traits also show high heritability: intraocular pressure and optic nerve head excavation damage quantified as the vertical cup-to-disc ratio. Here, since much of glaucoma heritability remains unexplained, we conducted a large-scale multitrait genome-wide association study in participants of European ancestry combining primary open-angle glaucoma and its two associated traits (total sample size over 600,000) to substantially improve genetic discovery power (263 loci). We further increased our power by then employing a multiancestry approach, which increased the number of independent risk loci to 312, with the vast majority replicating in a large independent cohort from 23andMe, Inc. (total sample size over 2.8 million; 296 loci replicated at P < 0.05, 240 after Bonferroni correction). Leveraging multiomics datasets, we identified many potential druggable genes, including neuro-protection targets likely to act via the optic nerve, a key advance for glaucoma because all existing drugs only target intraocular pressure. We further used Mendelian randomization and genetic correlation-based approaches to identify novel links to other complex traits, including immune-related diseases such as multiple sclerosis and systemic lupus erythematosus.

Spatially varying selection between habitats drives physiological shifts and local adaptation in a broadcast spawning coral on a remote atoll in Western Australia.

At the Rowley Shoals in Western Australia, the prominent reef flat becomes exposed on low tide and the stagnant water in the shallow atoll lagoons heats up, creating a natural laboratory for characterizing the mechanisms of coral resilience to climate change. To explore these mechanisms in the reef coral Acropora tenuis, we collected samples from lagoon and reef slope habitats and combined whole-genome sequencing, ITS2 metabarcoding, experimental heat stress, and transcriptomics. Despite high gene flow across the atoll, we identified clear shifts in allele frequencies between habitats at relatively small linked genomic islands. Common garden heat stress assays showed corals from the lagoon to be more resistant to bleaching, and RNA sequencing revealed marked differences in baseline levels of gene expression between habitats. Our results provide new insight into the complex mechanisms of coral resilience to climate change and highlight the potential for spatially varying selection across complex coral reef seascapes to drive pronounced ecological divergence in climate-related traits.

Recovery among post-arrest patients with mild-to-moderate cerebral edema.

BACKGROUND: Cerebral edema after cardiac arrest may be a modifiable cause of secondary brain injury. We aimed to identify processes of care associated with recovery in a cohort of patients with mild to moderate edema. METHODS: We conducted a retrospective cohort study of adults resuscitated from out-of-hospital arrest (OHCA) at a single center from 2010 to 2018. We included those with cerebral edema ranging from mild to moderate (gray to white matter attenuation ratio (GWR) 1.2 to 1.3 on initial brain computerized tomography (CT). We used Pittsburgh Cardiac Arrest Category (PCAC) to adjust for illness severity and considered the following values in the first 24 h of admission as additional predictors: GWR, lab values affecting serum osmolality (sodium, glucose, blood urea nitrogen (BUN)), total osmolality, change in osmolality from 0 to 24 h, cardiac etiology of arrest, targeted temperature to 33 °C (vs 36 °C), time-weighted mean arterial pressure (MAP), partial pressures of arterial oxygen and carbon dioxide and select medications. Our primary outcome was discharge with cerebral performance category 1-3. We used unadjusted and adjusted logistic regression for analysis. RESULTS: We included 214 patients for whom CT was performed median 3.8 [IQR 2.4-5.2] hours after collapse. Median age was 57 [IQR 48-67] years, 82 (38%) were female, and 68 (32%) arrested from ventricular tachycardia or fibrillation. In adjusted models, modifiable processes of care were not associated with outcome. CONCLUSIONS: Illness severity, but not modifiable processes of care, were associated with recovery among post-arrest patients with mild-to-moderate cerebral edema.

Impact of essential workers in the context of social distancing for epidemic control.

New emerging infectious diseases are identified every year, a subset of which become global pandemics like COVID-19. In the case of COVID-19, many governments have responded to the ongoing pandemic by imposing social policies that restrict contacts outside of the home, resulting in a large fraction of the workforce either working from home or not working. To ensure essential services, however, a substantial number of workers are not subject to these limitations, and maintain many of their pre-intervention contacts. To explore how contacts among such "essential" workers, and between essential workers and the rest of the population, impact disease risk and the effectiveness of pandemic control, we evaluated several mathematical models of essential worker contacts within a standard epidemiology framework. The models were designed to correspond to key characteristics of cashiers, factory employees, and healthcare workers. We find in all three models that essential workers are at substantially elevated risk of infection compared to the rest of the population, as has been documented, and that increasing the numbers of essential workers necessitates the imposition of more stringent controls on contacts among the rest of the population to manage the pandemic. Importantly, however, different archetypes of essential workers differ in both their individual probability of infection and impact on the broader pandemic dynamics, highlighting the need to understand and target intervention for the specific risks faced by different groups of essential workers. These findings, especially in light of the massive human costs of the current COVID-19 pandemic, indicate that contingency plans for future epidemics should account for the impacts of essential workers on disease spread.

Extensive Recombination Suppression and Epistatic Selection Causes Chromosome-Wide Differentiation of a Selfish Sex Chromosome in Drosophila pseudoobscura.

Sex-Ratio (SR) chromosomes are selfish X-chromosomes that distort Mendelian segregation and are commonly associated with inversions. These chromosomal rearrangements suppress recombination with Standard (ST) X-chromosomes and are hypothesized to maintain multiple alleles important for distortion in a single large haplotype. Here, we conduct a multifaceted study of the multiply inverted Drosophila pseudoobscura SR chromosome to understand the evolutionary history, genetic architecture, and present-day dynamics that shape this enigmatic selfish chromosome. The D. pseudoobscura SR chromosome has three nonoverlapping inversions of the right arm of the metacentric X-chromosome: basal, medial, and terminal. We find that 23 of 29 Mb of the D. pseudoobscuraX-chromosome right arm is highly differentiated between the Standard and SR arrangements, including a 6.6 Mb collinear region between the medial and terminal inversions. Although crossing-over is heavily suppressed on this chromosome arm, we discover it is not completely eliminated, with measured rates indicating recombination suppression alone cannot explain patterns of differentiation or the near-perfect association of the three SR chromosome inversions in nature. We then demonstrate the ancient basal and medial inversions of the SR chromosome contain genes sufficient to cause weak distortion. In contrast, the younger terminal inversion cannot distort by itself, but contains at least one modifier gene necessary for full manifestation of strong sex chromosome distortion. By parameterizing population genetic models for chromosome-wide linkage disequilibrium with our experimental results, we infer that strong selection acts to maintain the near-perfect association of SR chromosome inversions in present-day populations. Based on comparative genomic analyses, direct recombination experiments, segregation distortion assays, and population genetic modeling, we conclude the combined action of suppressed recombination and strong, ongoing, epistatic selection shape the D. pseudoobscura SR arrangement into a highly differentiated chromosome.

Population genetics of the coral Acropora millepora: Toward genomic prediction of bleaching.

Although reef-building corals are declining worldwide, responses to bleaching vary within and across species and are partly heritable. Toward predicting bleaching response from genomic data, we generated a chromosome-scale genome assembly for the coral Acropora millepora We obtained whole-genome sequences for 237 phenotyped samples collected at 12 reefs along the Great Barrier Reef, among which we inferred little population structure. Scanning the genome for evidence of local adaptation, we detected signatures of long-term balancing selection in the heat-shock co-chaperone sacsin We conducted a genome-wide association study of visual bleaching score for 213 samples, incorporating the polygenic score derived from it into a predictive model for bleaching in the wild. These results set the stage for genomics-based approaches in conservation strategies.

Measuring intolerance to mutation in human genetics.

In numerous applications, from working with animal models to mapping the genetic basis of human disease susceptibility, knowing whether a single disrupting mutation in a gene is likely to be deleterious is useful. With this goal in mind, a number of measures have been developed to identify genes in which protein-truncating variants (PTVs), or other types of mutations, are absent or kept at very low frequency in large population samples-genes that appear 'intolerant' to mutation. One measure in particular, the probability of being loss-of-function intolerant (pLI), has been widely adopted. This measure was designed to classify genes into three categories, null, recessive and haploinsufficient, on the basis of the contrast between observed and expected numbers of PTVs. Such population-genetic approaches can be useful in many applications. As we clarify, however, they reflect the strength of selection acting on heterozygotes and not dominance or haploinsufficiency.

Investigating the viral ecology of global bee communities with high-throughput metagenomics.

Bee viral ecology is a fascinating emerging area of research: viruses exert a range of effects on their hosts, exacerbate impacts of other environmental stressors, and, importantly, are readily shared across multiple bee species in a community. However, our understanding of bee viral communities is limited, as it is primarily derived from studies of North American and European Apis mellifera populations. Here, we examined viruses in populations of A. mellifera and 11 other bee species from 9 countries, across 4 continents and Oceania. We developed a novel pipeline to rapidly and inexpensively screen for bee viruses. This pipeline includes purification of encapsulated RNA/DNA viruses, sequence-independent amplification, high throughput sequencing, integrated assembly of contigs, and filtering to identify contigs specifically corresponding to viral sequences. We identified sequences for (+)ssRNA, (-)ssRNA, dsRNA, and ssDNA viruses. Overall, we found 127 contigs corresponding to novel viruses (i.e. previously not observed in bees), with 27 represented by >0.1% of the reads in a given sample, and 7 contained an RdRp or replicase sequence which could be used for robust phylogenetic analysis. This study provides a sequence-independent pipeline for viral metagenomics analysis, and greatly expands our understanding of the diversity of viruses found in bee communities.

A food-predictive cue attributed with incentive salience engages subcortical afferents and efferents of the paraventricular nucleus of the thalamus.

The paraventricular nucleus of the thalamus (PVT) has been implicated in behavioral responses to reward-associated cues. However, the precise role of the PVT in these behaviors has been difficult to ascertain since Pavlovian-conditioned cues can act as both predictive and incentive stimuli. The "sign-tracker/goal-tracker" rat model has allowed us to further elucidate the role of the PVT in cue-motivated behaviors, identifying this structure as a critical component of the neural circuitry underlying individual variation in the propensity to attribute incentive salience to reward cues. The current study assessed differences in the engagement of specific PVT afferents and efferents in response to presentation of a food-cue that had been attributed with only predictive value or with both predictive and incentive value. The retrograde tracer fluorogold (FG) was injected into the PVT or the nucleus accumbens (NAc) of rats, and cue-induced c-Fos in FG-labeled cells was quantified. Presentation of a predictive stimulus that had been attributed with incentive value elicited c-Fos in PVT afferents from the lateral hypothalamus, medial amygdala (MeA), and the prelimbic cortex (PrL), as well as posterior PVT efferents to the NAc. PVT afferents from the PrL also showed elevated c-Fos levels following presentation of a predictive stimulus alone. Thus, presentation of an incentive stimulus results in engagement of subcortical brain regions; supporting a role for the hypothalamic-thalamic-striatal axis, as well as the MeA, in mediating responses to incentive stimuli; whereas activity in the PrL to PVT pathway appears to play a role in processing the predictive qualities of reward-paired stimuli.

Genomics of Natural Populations: How Differentially Expressed Genes Shape the Evolution of Chromosomal Inversions in Drosophila pseudoobscura.

Chromosomal rearrangements can shape the structure of genetic variation in the genome directly through alteration of genes at breakpoints or indirectly by holding combinations of genetic variants together due to reduced recombination. The third chromosome of Drosophila pseudoobscura is a model system to test hypotheses about how rearrangements are established in populations because its third chromosome is polymorphic for >30 gene arrangements that were generated by a series of overlapping inversion mutations. Circumstantial evidence has suggested that these gene arrangements are selected. Despite the expected homogenizing effects of extensive gene flow, the frequencies of arrangements form gradients or clines in nature, which have been stable since the system was first described >80 years ago. Furthermore, multiple arrangements exist at appreciable frequencies across several ecological niches providing the opportunity for heterokaryotypes to form. In this study, we tested whether genes are differentially expressed among chromosome arrangements in first instar larvae, adult females and males. In addition, we asked whether transcriptional patterns in heterokaryotypes are dominant, semidominant, overdominant, or underdominant. We find evidence for a significant abundance of differentially expressed genes across the inverted regions of the third chromosome, including an enrichment of genes involved in sensory perception for males. We find the majority of loci show additivity in heterokaryotypes. Our results suggest that multiple genes have expression differences among arrangements that were either captured by the original inversion mutation or accumulated after it reached polymorphic frequencies, providing a potential source of genetic variation for selection to act upon. These data suggest that the inversions are favored because of their indirect effect of recombination suppression that has held different combinations of differentially expressed genes together in the various gene arrangement backgrounds.

Evidence for stabilizing selection on codon usage in chromosomal rearrangements of Drosophila pseudoobscura.

There has been a renewed interest in investigating the role of stabilizing selection acting on genome-wide traits such as codon usage bias. Codon bias, when synonymous codons are used at unequal frequencies, occurs in a wide variety of taxa. Standard evolutionary models explain the maintenance of codon bias through a balance of genetic drift, mutation and weak purifying selection. The efficacy of selection is expected to be reduced in regions of suppressed recombination. Contrary to observations in Drosophila melanogaster, some recent studies have failed to detect a relationship between the recombination rate, intensity of selection acting at synonymous sites, and the magnitude of codon bias as predicted under these standard models. Here, we examined codon bias in 2798 protein coding loci on the third chromosome of D. pseudoobscura using whole-genome sequences of 47 individuals, representing five common third chromosome gene arrangements. Fine-scale recombination maps were constructed using more than 1 million segregating sites. As expected, recombination was demonstrated to be significantly suppressed between chromosome arrangements, allowing for a direct examination of the relationship between recombination, selection, and codon bias. As with other Drosophila species, we observe a strong mutational bias away from the most frequently used codons. We find the rate of synonymous and nonsynonymous polymorphism is variable between different amino acids. However, we do not observe a reduction in codon bias or the strength of selection in regions of suppressed recombination as expected. Instead, we find that the interaction between weak stabilizing selection and mutational bias likely plays a role in shaping the composition of synonymous codons across the third chromosome in D. pseudoobscura.