Genetic architecture, spatial heterogeneity, and the coevolutionary arms race between newts and snakes.

Coevolution between two species can lead to exaggerated phenotypes that vary in a correlated manner across space. However, the conditions under which we expect such spatially varying coevolutionary patterns in polygenic traits are not well-understood. We investigate the coevolutionary dynamics between two species undergoing reciprocal adaptation across space and time, using simulations inspired by the Taricha newt - Thamnophis garter snake system. One striking observation from this system is that newts in some areas carry much more tetrodotoxin than in other areas, and garter snakes that live near more toxic newts tend to be more resistant to this toxin, a correlation seen across several broad geographic areas. Furthermore, snakes seem to be "winning" the coevolutionary arms race, i.e., having a high level of resistance compared to local newt toxicity, despite substantial variation in both toxicity and resistance across the range. We explore how possible genetic architectures of the toxin and resistance traits would affect the coevolutionary dynamics by manipulating both mutation rate and effect size of mutations across many simulations. We find that coevolutionary dynamics alone were not sufficient in our simulations to produce the striking mosaic of levels of toxicity and resistance observed in nature, but simulations with ecological heterogeneity (in trait costliness or interaction rate) did produce such patterns. We also find that in simulations, newts tend to "win" across most combinations of genetic architectures, although the species with higher mutational genetic variance tends to have an advantage.

Expanding the stdpopsim species catalog, and lessons learned for realistic genome simulations.

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.

Ten simple rules for an inclusive summer coding program for non-computer-science undergraduates.

Since 2015, we have run a free 9-week summer program that provides non-computer science (CS) undergraduates at San Francisco State University (SFSU) with experience in coding and doing research. Undergraduate research experiences remain very limited at SFSU and elsewhere, so the summer program provides opportunities for many more students beyond the mentoring capacity of our university laboratories. In addition, we were concerned that many students from historically underrepresented (HU) groups may be unable to take advantage of traditional summer research programs because these programs require students to relocate or be available full time, which is not feasible for students who have family, work, or housing commitments. Our program, which is local and part-time, serves about 5 times as many students as a typical National Science Foundation (NSF) Research Experiences for Undergraduates (REU) program, on a smaller budget. Based on our experiences, we present 10 simple rules for busy faculty who want to create similar programs to engage non-CS HU undergraduates in computational research. Note that while some of the strategies we implement are based on evidence-based publications in the social sciences or education research literature, the original suggestions we make here are based on our trial-and-error experiences, rather than formal hypothesis testing.

CpG-creating mutations are costly in many human viruses.

Mutations can occur throughout the virus genome and may be beneficial, neutral or deleterious. We are interested in mutations that yield a C next to a G, producing CpG sites. CpG sites are rare in eukaryotic and viral genomes. For the eukaryotes, it is thought that CpG sites are rare because they are prone to mutation when methylated. In viruses, we know less about why CpG sites are rare. A previous study in HIV suggested that CpG-creating transition mutations are more costly than similar non-CpG-creating mutations. To determine if this is the case in other viruses, we analyzed the allele frequencies of CpG-creating and non-CpG-creating mutations across various strains, subtypes, and genes of viruses using existing data obtained from Genbank, HIV Databases, and Virus Pathogen Resource. Our results suggest that CpG sites are indeed costly for most viruses. By understanding the cost of CpG sites, we can obtain further insights into the evolution and adaptation of viruses.