Complete genome sequence of the Streptomyces bacteriophage Amabiko

Amabiko is a lytic subcluster BE2 bacteriophage that infects Streptomyces scabiei-a bacterium causing common scab in potatoes. Its 131,414 bp genome has a GC content of 49.5% and contains 245 putative protein-coding genes, 45 tRNAs, and one tmRNA. Amabiko is closely related to Streptomyces bacteriophage MindFlayer (gene content similarity: 86.5%).

The Effects of Mutation and Recombination Rate Heterogeneity on the Inference of Demography and the Distribution of Fitness Effects

Disentangling the effects of demography and selection has remained a focal point of population genetic analysis. Knowledge about mutation and recombination is essential in this endeavor; however, despite clear evidence that both mutation and recombination rates vary across genomes, it is common practice to model both rates as fixed. In this study, we quantify how this unaccounted for rate heterogeneity may impact inference using common approaches for inferring selection (DFE-alpha, Grapes, and polyDFE) and/or demography (fastsimcoal2 and δaδi). We demonstrate that, if not properly modeled, this heterogeneity can increase uncertainty in the estimation of demographic and selective parameters and in some scenarios may result in mis-leading inference. These results highlight the importance of quantifying the fundamental evolutionary parameters of mutation and recombination before utilizing population genomic data to quantify the effects of genetic drift (i.e. as modulated by demographic history) and selection; or, at the least, that the effects of uncertainty in these parameters can and should be directly modeled in downstream inference.

Novel Insights into the Landscape of Crossover and Noncrossover Events in Rhesus Macaques (Macaca mulatta)

Meiotic recombination landscapes differ greatly between distantly and closely related taxa, populations, individuals, sexes, and even within genomes; however, the factors driving this variation are yet to be well elucidated. Here, we directly estimate contemporary crossover rates and, for the first time, noncrossover rates in rhesus macaques (Macaca mulatta) from four three-generation pedigrees comprising 32 individuals. We further compare these results with historical, demography-aware, linkage disequilibrium-based recombination rate estimates. From paternal meioses in the pedigrees, 165 crossover events with a median resolution of 22.3 kb were observed, corresponding to a male autosomal map length of 2,357 cM-approximately 15% longer than an existing linkage map based on human microsatellite loci. In addition, 85 noncrossover events with a mean tract length of 155 bp were identified-similar to the tract lengths observed in the only other two primates in which noncrossovers have been studied to date, humans and baboons. Consistent with observations in other placental mammals with PRDM9-directed recombination, crossover (and to a lesser extent noncrossover) events in rhesus macaques clustered in intergenic regions and toward the chromosomal ends in males-a pattern in broad agreement with the historical, sex-averaged recombination rate estimates-and evidence of GC-biased gene conversion was observed at noncrossover sites.

Computational host range prediction-The good, the bad, and the ugly

The rapid emergence and spread of antimicrobial resistance across the globe have prompted the usage of bacteriophages (i.e. viruses that infect bacteria) in a variety of applications ranging from agriculture to biotechnology and medicine. In order to effectively guide the application of bacteriophages in these multifaceted areas, information about their host ranges-that is the bacterial strains or species that a bacteriophage can successfully infect and kill-is essential. Utilizing sixteen broad-spectrum (polyvalent) bacteriophages with experimentally validated host ranges, we here benchmark the performance of eleven recently developed computational host range prediction tools that provide a promising and highly scalable supplement to traditional, but laborious, experimental procedures. We show that machine- and deep-learning approaches offer the highest levels of accuracy and precision-however, their predominant predictions at the species- or genus-level render them ill-suited for applications outside of an ecosystems metagenomics framework. In contrast, only moderate sensitivity (<80 per cent) could be reached at the strain-level, albeit at low levels of precision (<40 per cent). Taken together, these limitations demonstrate that there remains room for improvement in the active scientific field of in silico host prediction to combat the challenge of guiding experimental designs to identify the most promising bacteriophage candidates for any given application.

The effects of mutation and recombination rate heterogeneity on the inference of demography and the distribution of fitness effects

Disentangling the effects of demography and selection has remained a focal point of population genetic analysis. Knowledge about mutation and recombination is essential in this endeavour; however, despite clear evidence that both mutation and recombination rates vary across genomes, it is common practice to model both rates as fixed. In this study, we quantify how this unaccounted for rate heterogeneity may impact inference using common approaches for inferring selection (DFE-alpha, Grapes, and polyDFE) and/or demography (fastsimcoal2 and δaδi). We demonstrate that, if not properly modelled, this heterogeneity can increase uncertainty in the estimation of demographic and selective parameters and in some scenarios may result in mis-leading inference. These results highlight the importance of quantifying the fundamental evolutionary parameters of mutation and recombination prior to utilizing population genomic data to quantify the effects of genetic drift (i.e., as modulated by demographic history) and selection; or, at the least, that the effects of uncertainty in these parameters can and should be directly modelled in downstream inference.

Protective effect of pre-existing natural immunity in a nonhuman primate reinfection model of congenital cytomegalovirus infection

Congenital cytomegalovirus (cCMV) is the leading infectious cause of neurologic defects in newborns with particularly severe sequelae in the setting of primary CMV infection in the first trimester of pregnancy. The majority of cCMV cases worldwide occur after non-primary infection in CMV-seropositive women; yet the extent to which pre-existing natural CMV-specific immunity protects against CMV reinfection or reactivation during pregnancy remains ill-defined. We previously reported on a novel nonhuman primate model of cCMV in rhesus macaques where 100% placental transmission and 83% fetal loss were seen in CD4+ T lymphocyte-depleted rhesus CMV (RhCMV)-seronegative dams after primary RhCMV infection. To investigate the protective effect of preconception maternal immunity, we performed reinfection studies in CD4+ T lymphocyte-depleted RhCMV-seropositive dams inoculated in late first / early second trimester gestation with RhCMV strains 180.92 (n = 2), or RhCMV UCD52 and FL-RhCMVΔRh13.1/SIVgag, a wild-type-like RhCMV clone with SIVgag inserted as an immunological marker, administered separately (n = 3). An early transient increase in circulating monocytes followed by boosting of the pre-existing RhCMV-specific CD8+ T lymphocyte and antibody response was observed in the reinfected dams but not in control CD4+ T lymphocyte-depleted dams. Emergence of SIV Gag-specific CD8+ T lymphocyte responses in macaques inoculated with the FL-RhCMVΔRh13.1/SIVgag virus confirmed reinfection. Placental transmission was detected in only one of five reinfected dams and there were no adverse fetal sequelae. Viral whole genome, short-read, deep sequencing analysis confirmed transmission of both reinfection RhCMV strains across the placenta with ~30% corresponding to FL-RhCMVΔRh13.1/SIVgag and ~70% to RhCMV UCD52, consistent with the mixed human CMV infections reported in infants with cCMV. Our data showing reduced placental transmission and absence of fetal loss after non-primary as opposed to primary infection in CD4+ T lymphocyte-depleted dams indicates that preconception maternal CMV-specific CD8+ T lymphocyte and/or humoral immunity can protect against cCMV infection.

Developing an evolutionary baseline model for humans: jointly inferring purifying selection with population history

Building evolutionarily appropriate baseline models for natural populations is not only important for answering fundamental questions in population genetics - including quantifying the relative contributions of adaptive vs. non-adaptive processes - but it is also essential for identifying candidate loci experiencing relatively rare and episodic forms of selection (e.g., positive or balancing selection). Here, a baseline model was developed for a human population of West African ancestry, the Yoruba, comprising processes constantly operating on the genome (i.e., purifying and background selection, population size changes, recombination rate heterogeneity, and gene conversion). Specifically, to perform joint inference of selective effects with demography, an approximate Bayesian approach was employed that utilizes the decay of background selection effects around functional elements, taking into account genomic architecture. This approach inferred a recent 6-fold population growth together with a distribution of fitness effects that is skewed towards effectively neutral mutations. Importantly, these results further suggest that, while strong and/or frequent recurrent positive selection is inconsistent with observed data, weak to moderate positive selection is consistent but unidentifiable if rare.

Developing an appropriate evolutionary baseline model for the study of human cytomegalovirus

Human cytomegalovirus (HCMV) represents a major threat to human health, contributing to both birth defects in neonates as well as organ transplant failure and opportunistic infections in immunocompromised individuals. HCMV exhibits considerable inter- and intra-host diversity, which likely influences the pathogenicity of the virus. Therefore, understanding the relative contributions of various evolutionary forces in shaping patterns of variation is of critical importance both mechanistically as well as clinically. Herein we present the individual components of an evolutionary baseline model for HCMV, with a particular focus on congenital infections for the sake of illustration - including mutation and recombination rates, the distribution of fitness effects, infection dynamics, as well as compartmentalization - and describe the current state of knowledge of each. By building this baseline model, researchers will be able to better describe the range of possible evolutionary scenarios contributing to observed variation, as well as improve power and reduce false-positive rates when scanning for adaptive mutations in the HCMV genome.

Developing an evolutionary baseline model for humans: jointly inferring purifying selection with population history

Building evolutionarily appropriate baseline models for natural populations is not only important for answering fundamental questions in population genetics - including quantifying the relative contributions of adaptive vs. non-adaptive processes - but it is also essential for identifying candidate loci experiencing relatively rare and episodic forms of selection ( e.g., positive or balancing selection). Here, a baseline model was developed for a human population of West African ancestry, the Yoruba, comprising processes constantly operating on the genome ( i.e. , purifying and background selection, population size changes, recombination rate heterogeneity, and gene conversion). Specifically, to perform joint inference of selective effects with demography, an approximate Bayesian approach was employed that utilizes the decay of background selection effects around functional elements, taking into account genomic architecture. This approach inferred a recent 6-fold population growth together with a distribution of fitness effects that is skewed towards effectively neutral mutations. Importantly, these results further suggest that, while strong and/or frequent recurrent positive selection is inconsistent with observed data, weak to moderate positive selection is consistent but unidentifiable if rare.

Protective effect of pre-existing natural immunity in a nonhuman primate reinfection model of congenital cytomegalovirus infection

Congenital cytomegalovirus (cCMV) is the leading infectious cause of neurologic defects in newborns with particularly severe sequelae in the setting of primary CMV infection in the first trimester of pregnancy. The majority of cCMV cases worldwide occur after non-primary infection in CMV-seropositive women; yet the extent to which pre-existing natural CMV-specific immunity protects against CMV reinfection or reactivation during pregnancy remains ill-defined. We previously reported on a novel nonhuman primate model of cCMV in rhesus macaques where 100% placental transmission and 83% fetal loss were seen in CD4 + T lymphocyte-depleted rhesus CMV (RhCMV)-seronegative dams after primary RhCMV infection. To investigate the protective effect of preconception maternal immunity, we performed reinfection studies in CD4+ T lymphocyte-depleted RhCMV-seropositive dams inoculated in late first / early second trimester gestation with RhCMV strains 180.92 ( n =2), or RhCMV UCD52 and FL-RhCMVΔRh13.1/SIV gag , a wild-type-like RhCMV clone with SIV gag inserted as an immunological marker ( n =3). An early transient increase in circulating monocytes followed by boosting of the pre-existing RhCMV-specific CD8+ T lymphocyte and antibody response was observed in the reinfected dams but not in control CD4+ T lymphocyte-depleted dams. Emergence of SIV Gag-specific CD8+ T lymphocyte responses in macaques inoculated with the FL-RhCMVΔRh13.1/SIV gag virus confirmed reinfection. Placental transmission was detected in only one of five reinfected dams and there were no adverse fetal sequelae. Viral whole genome, short-read, deep sequencing analysis confirmed transmission of both reinfection RhCMV strains across the placenta with ∼30% corresponding to FL-RhCMVΔRh13.1/SIV gag and ∼70% to RhCMV UCD52, consistent with the mixed human CMV infections reported in infants with cCMV. Our data showing reduced placental transmission and absence of fetal loss after non-primary as opposed to primary infection in CD4+ T lymphocyte-depleted dams indicates that preconception maternal CMV-specific CD8+ T lymphocyte and/or humoral immunity can protect against cCMV infection.

Author Summary

Globally, pregnancies in CMV-seropositive women account for the majority of cases of congenital CMV infection but the immune responses needed for protection against placental transmission in mothers with non-primary infection remains unknown. Recently, we developed a nonhuman primate model of primary rhesus CMV (RhCMV) infection in which placental transmission and fetal loss occurred in RhCMV-seronegative CD4+ T lymphocyte-depleted macaques. By conducting similar studies in RhCMV-seropositive dams, we demonstrated the protective effect of pre-existing natural CMV-specific CD8+ T lymphocytes and humoral immunity against congenital CMV after reinfection. A 5-fold reduction in congenital transmission and complete protection against fetal loss was observed in dams with pre-existing immunity compared to primary CMV in this model. Our study is the first formal demonstration in a relevant model of human congenital CMV that natural pre-existing CMV-specific maternal immunity can limit congenital CMV transmission and its sequelae. The nonhuman primate model of non-primary congenital CMV will be especially relevant to studying immune requirements of a maternal vaccine for women in high CMV seroprevalence areas at risk of repeated CMV reinfections during pregnancy.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

The Impact of Sample Size and Population History on Observed Mutational Spectra: A Case Study in Human and Chimpanzee Populations

Recent studies have highlighted variation in the mutational spectra among human populations as well as closely related hominoids-yet little remains known about the genetic and nongenetic factors driving these rate changes across the genome. Pinpointing the root causes of these differences is an important endeavor that requires careful comparative analyses of population-specific mutational landscapes at both broad and fine genomic scales. However, several factors can confound such analyses. Although previous studies have shown that technical artifacts, such as sequencing errors and batch effects, can contribute to observed mutational shifts, other potentially confounding parameters have received less attention thus far. Using population genetic simulations of human and chimpanzee populations as an illustrative example, we here show that the sample size required for robust inference of mutational spectra depends on the population-specific demographic history. As a consequence, the power to detect rate changes is high in certain hominoid populations while, for others, currently available sample sizes preclude analyses at fine genomic scales.

Phylogenomic analyses and host range prediction of cluster P mycobacteriophages

Bacteriophages, infecting bacterial hosts in every environment on our planet, are a driver of adaptive evolution in bacterial communities. At the same time, the host range of many bacteriophages-and thus one of the selective pressures acting on complex microbial systems in nature-remains poorly characterized. Here, we computationally inferred the putative host ranges of 40 cluster P mycobacteriophages, including members from 6 subclusters (P1-P6). A series of comparative genomic analyses revealed that mycobacteriophages of subcluster P1 are restricted to the Mycobacterium genus, whereas mycobacteriophages of subclusters P2-P6 are likely also able to infect other genera, several of which are commonly associated with human disease. Further genomic analysis highlighted that the majority of cluster P mycobacteriophages harbor a conserved integration-dependent immunity system, hypothesized to be the ancestral state of a genetic switch that controls the shift between lytic and lysogenic life cycles-a temperate characteristic that impedes their usage in antibacterial applications.

Complete Genome Sequence of the Cluster P Mycobacteriophage Phegasus

We characterized the complete genome of the cluster P mycobacteriophage Phegasus. Its 47.5-kb genome contains 81 protein-coding genes, 36 of which could be assigned a putative function. Phegasus is most closely related to two subcluster P1 bacteriophages, Mangethe and Majeke, with an average nucleotide identity of 99.63% each.

Comparative Genomics of Closely-Related Gordonia Cluster DR Bacteriophages

Bacteriophages infecting bacteria of the genus Gordonia have increasingly gained interest in the scientific community for their diverse applications in agriculture, biotechnology, and medicine, ranging from biocontrol agents in wastewater management to the treatment of opportunistic pathogens in pulmonary disease patients. However, due to the time and costs associated with experimental isolation and cultivation, host ranges for many bacteriophages remain poorly characterized, hindering a more efficient usage of bacteriophages in these areas. Here, we perform a series of computational genomic inferences to predict the putative host ranges of all Gordonia cluster DR bacteriophages known to date. Our analyses suggest that BiggityBass (as well as several of its close relatives) is likely able to infect host bacteria from a wide range of genera—from Gordonia to Nocardia to Rhodococcus, making it a suitable candidate for future phage therapy and wastewater treatment strategies.

Recommendations for improving statistical inference in population genomics

The field of population genomics has grown rapidly in response to the recent advent of affordable, large-scale sequencing technologies. As opposed to the situation during the majority of the 20th century, in which the development of theoretical and statistical population genetic insights outpaced the generation of data to which they could be applied, genomic data are now being produced at a far greater rate than they can be meaningfully analyzed and interpreted. With this wealth of data has come a tendency to focus on fitting specific (and often rather idiosyncratic) models to data, at the expense of a careful exploration of the range of possible underlying evolutionary processes. For example, the approach of directly investigating models of adaptive evolution in each newly sequenced population or species often neglects the fact that a thorough characterization of ubiquitous nonadaptive processes is a prerequisite for accurate inference. We here describe the perils of these tendencies, present our consensus views on current best practices in population genomic data analysis, and highlight areas of statistical inference and theory that are in need of further attention. Thereby, we argue for the importance of defining a biologically relevant baseline model tuned to the details of each new analysis, of skepticism and scrutiny in interpreting model fitting results, and of carefully defining addressable hypotheses and underlying uncertainties.

The third international hackathon for applying insights into large-scale genomic composition to use cases in a wide range of organisms

In October 2021, 59 scientists from 14 countries and 13 U.S. states collaborated virtually in the Third Annual Baylor College of Medicine & DNANexus Structural Variation hackathon. The goal of the hackathon was to advance research on structural variants (SVs) by prototyping and iterating on open-source software. This led to nine hackathon projects focused on diverse genomics research interests, including various SV discovery and genotyping methods, SV sequence reconstruction, and clinically relevant structural variation, including SARS-CoV-2 variants. Repositories for the projects that participated in the hackathon are available at https://github.com/collaborativebioinformatics.

The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis

Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.

The recombination landscapes of spiny lizards (genus Sceloporus)

Despite playing a critical role in evolutionary processes and outcomes, relatively little is known about rates of recombination in the vast majority of species, including squamate reptiles—the second largest order of extant vertebrates, many species of which serve as important model organisms in evolutionary and ecological studies. This paucity of data has resulted in limited resolution on questions related to the causes and consequences of rate variation between species and populations, the determinants of within-genome rate variation, as well as the general tempo of recombination rate evolution on this branch of the tree of life. In order to address these questions, it is thus necessary to begin broadening our phylogenetic sampling. We here provide the first fine-scale recombination maps for two species of spiny lizards, Sceloporus jarrovii and Sceloporus megalepidurus, which diverged at least 12 Mya. As might be expected from similarities in karyotype, population-scaled recombination landscapes are largely conserved on the broad-scale. At the same time, considerable variation exists at the fine-scale, highlighting the importance of incorporating species-specific recombination maps in future population genomic studies.

The mutationathon highlights the importance of reaching standardization in estimates of pedigree-based germline mutation rates

In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various nonhuman species have been published. However, methodological differences among studies in detecting germline mutations and estimating mutation rates make direct comparisons difficult. Here, we describe the many different steps involved in estimating pedigree-based mutation rates, including sampling, sequencing, mapping, variant calling, filtering, and appropriately accounting for false-positive and false-negative rates. For each step, we review the different methods and parameter choices that have been used in the recent literature. Additionally, we present the results from a ‘Mutationathon,’ a competition organized among five research labs to compare germline mutation rate estimates for a single pedigree of rhesus macaques. We report almost a twofold variation in the final estimated rate among groups using different post-alignment processing, calling, and filtering criteria, and provide details into the sources of variation across studies. Though the difference among estimates is not statistically significant, this discrepancy emphasizes the need for standardized methods in mutation rate estimations and the difficulty in comparing rates from different studies. Finally, this work aims to provide guidelines for computational and statistical benchmarks for future studies interested in identifying germline mutations from pedigrees.

Computational Prediction of Bacteriophage Host Ranges

Increased antibiotic resistance has prompted the development of bacteriophage agents for a multitude of applications in agriculture, biotechnology, and medicine. A key factor in the choice of agents for these applications is the host range of a bacteriophage, i.e., the bacterial genera, species, and strains a bacteriophage is able to infect. Although experimental explorations of host ranges remain the gold standard, such investigations are inherently limited to a small number of viruses and bacteria amendable to cultivation. Here, we review recently developed bioinformatic tools that offer a promising and high-throughput alternative by computationally predicting the putative host ranges of bacteriophages, including those challenging to grow in laboratory environments.

Studying mutation rate evolution in primates-the effects of computational pipelines and parameter choices

This commentary investigates the important role of computational pipeline and parameter choices in performing mutation rate estimation, using the recent article published in this journal by Bergeron et al. entitled "The germline mutational process in rhesus macaque and its implications for phylogenetic dating" as an illustrative example.

Phylogenetic relationships and codon usage bias amongst cluster K mycobacteriophages

Bacteriophages infecting pathogenic hosts play an important role in medical research, not only as potential treatments for antibiotic-resistant infections but also offering novel insights into pathogen genetics and evolution. A prominent example is cluster K mycobacteriophages infecting Mycobacterium tuberculosis, a causative agent of tuberculosis in humans. However, as handling M. tuberculosis as well as other pathogens in a laboratory remains challenging, alternative nonpathogenic relatives, such as Mycobacterium smegmatis, are frequently used as surrogates to discover therapeutically relevant bacteriophages in a safer environment. Consequently, the individual host ranges of the majority of cluster K mycobacteriophages identified to date remain poorly understood. Here, we characterized the complete genome of Stinson, a temperate subcluster K1 mycobacteriophage with a siphoviral morphology. A series of comparative genomic analyses revealed strong similarities with other cluster K mycobacteriophages, including the conservation of an immunity repressor gene and a toxin/antitoxin gene pair. Patterns of codon usage bias across the cluster offered important insights into putative host ranges in nature, highlighting that although all cluster K mycobacteriophages are able to infect M. tuberculosis, they are less likely to have shared an evolutionary infection history with Mycobacterium leprae (underlying leprosy) compared to the rest of the genus’ host species. Moreover, subcluster K1 mycobacteriophages are able to integrate into the genomes of Mycobacterium abscessus and Mycobacterium marinum—two bacteria causing pulmonary and cutaneous infections which are often difficult to treat due to their drug resistance.

Common Polymorphisms in the Glycoproteins of Human Cytomegalovirus and Associated Strain-Specific Immunity

Human cytomegalovirus (HCMV), one of the most prevalent viruses across the globe, is a common cause of morbidity and mortality for immunocompromised individuals. Recent clinical observations have demonstrated that mixed strain infections are common and may lead to more severe disease progression. This clinical observation illustrates the complexity of the HCMV genome and emphasizes the importance of taking a population-level view of genotypic evolution. Here we review frequently sampled polymorphisms in the glycoproteins of HCMV, comparing the variable regions, and summarizing their corresponding geographic distributions observed to date. The related strain-specific immunity, including neutralization activity and antigen-specific cellular immunity, is also discussed. Given that these glycoproteins are common targets for vaccine design and anti-viral therapies, this observed genetic variation represents an important resource for future efforts to combat HCMV infections.

A Fine-Scale Genetic Map for Vervet Monkeys

Despite its important biological role, the evolution of recombination rates remains relatively poorly characterized. This owes, in part, to the lack of high-quality genomic resources to address this question across diverse species. Humans and our closest evolutionary relatives, anthropoid apes, have remained a major focus of large-scale sequencing efforts, and thus recombination rate variation has been comparatively well studied in this group-with earlier work revealing a conservation at the broad- but not the fine-scale. However, in order to better understand the nature of this variation, and the time scales on which substantial modifications occur, it is necessary to take a broader phylogenetic perspective. I here present the first fine-scale genetic map for vervet monkeys based on whole-genome population genetic data from ten individuals and perform a series of comparative analyses with the great apes. The results reveal a number of striking features. First, owing to strong positive correlations with diversity and weak negative correlations with divergence, analyses suggest a dominant role for purifying and background selection in shaping patterns of variation in this species. Second, results support a generally reduced broad-scale recombination rate compared with the great apes, as well as a narrower fraction of the genome in which the majority of recombination events are observed to occur. Taken together, this data set highlights the great necessity of future research to identify genomic features and quantify evolutionary processes that are driving these rate changes across primates.

Recent population genomic insights into the genetic basis of arsenic tolerance in humans: the difficulties of identifying positively selected loci in strongly bottlenecked populations

Recent advances in genomics have enabled researchers to shed light on the evolutionary processes driving human adaptation, by revealing the genetic architectures underlying traits ranging from lactase persistence, to skin pigmentation, to hypoxic response, to arsenic tolerance. Complicating the identification of targets of positive selection in modern human populations is their complex demographic history, characterized by population bottlenecks and expansions, population structure, migration, and admixture. In particular, founder effects and recent strong population size reductions, such as those experienced by the indigenous peoples of the Americas, have severe impacts on genetic variation that can lead to the accumulation of large allele frequency differences between populations due to genetic drift rather than natural selection. While distinguishing the effects of demographic history from selection remains challenging, neglecting neutral processes can lead to the incorrect identification of candidate loci. We here review the recent population genomic insights into the genetic basis of arsenic tolerance in Andean populations, and utilize this example to highlight both the difficulties pertaining to the identification of local adaptations in strongly bottlenecked populations, as well as the importance of controlling for demographic history in selection scans.

The Demographic History of African Drosophila melanogaster

As one of the most commonly utilized organisms in the study of local adaptation, an accurate characterization of the demographic history of Drosophila melanogaster remains as an important research question. This owes both to the inherent interest in characterizing the population history of this model organism, as well as to the well-established importance of an accurate null demographic model for increasing power and decreasing false positive rates in genomic scans for positive selection. Although considerable attention has been afforded to this issue in non-African populations, less is known about the demographic history of African populations, including from the ancestral range of the species. While qualitative predictions and hypotheses have previously been forwarded, we here present a quantitative model fitting of the population history characterizing both the ancestral Zambian population range as well as the subsequently colonized west African populations, which themselves served as the source of multiple non-African colonization events. We here report the split time of the West African population at 72 kya, a date corresponding to human migration into this region as well as a period of climatic changes in the African continent. Furthermore, we have estimated population sizes at this split time. These parameter estimates thus represent an important null model for future investigations in to African and non-African D. melanogaster populations alike.

The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow

The interplay of gene flow, genetic drift, and local selective pressure is a dynamic process that has been well studied from a theoretical perspective over the last century. Wright and Haldane laid the foundation for expectations under an island-continent model, demonstrating that an island-specific beneficial allele may be maintained locally if the selection coefficient is larger than the rate of migration of the ancestral allele from the continent. Subsequent extensions of this model have provided considerably more insight. Yet, connecting theoretical results with empirical data has proven challenging, owing to a lack of information on the relationship between genotype, phenotype, and fitness. Here, we examine the demographic and selective history of deer mice in and around the Nebraska Sand Hills, a system in which variation at the Agouti locus affects cryptic coloration that in turn affects the survival of mice in their local habitat. We first genotyped 250 individuals from 11 sites along a transect spanning the Sand Hills at 660,000 single nucleotide polymorphisms across the genome. Using these genomic data, we found that deer mice first colonized the Sand Hills following the last glacial period. Subsequent high rates of gene flow have served to homogenize the majority of the genome between populations on and off the Sand Hills, with the exception of the Agouti pigmentation locus. Furthermore, mutations at this locus are strongly associated with the pigment traits that are strongly correlated with local soil coloration and thus responsible for cryptic coloration.

Linking a mutation to survival in wild mice

Adaptive evolution in new or changing environments can be difficult to predict because the functional connections between genotype, phenotype, and fitness are complex. Here, we make these explicit connections by combining field and laboratory experiments in wild mice. We first directly estimate natural selection on pigmentation traits and an underlying pigment locus, Agouti, by using experimental enclosures of mice on different soil colors. Next, we show how a mutation in Agouti associated with survival causes lighter coat color through changes in its protein binding properties. Together, our findings demonstrate how a sequence variant alters phenotype and then reveal the ensuing ecological consequences that drive changes in population allele frequency, thereby illuminating the process of evolution by natural selection.

On the Demographic and Selective Forces Shaping Patterns of Human Cytomegalovirus Variation within Hosts

Human cytomegalovirus (HCMV) is a member of the β -herpesvirus subfamily within Herpesviridae that is nearly ubiquitous in human populations, and infection generally results only in mild symptoms. However, symptoms can be severe in immunonaive individuals, and transplacental congenital infection of HCMV can result in serious neurological sequelae. Recent work has revealed much about the demographic and selective forces shaping the evolution of congenitally transmitted HCMV both on the level of hosts and within host compartments, providing insight into the dynamics of congenital infection, reinfection, and evolution of HCMV with important implications for the development of effective treatments and vaccines.

Loci associated with skin pigmentation identified in African populations

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

Direct estimate of the spontaneous germ line mutation rate in African green monkeys

Here, I provide the first direct estimate of the spontaneous mutation rate in an Old World monkey, using a seven individual, three-generation pedigree of African green monkeys. Eight de novo mutations were identified within ∼1.5 Gbp of accessible genome, corresponding to an estimated point mutation rate of 0.94 × 10^-8 per site per generation, suggesting an effective population size of ∼12000 for the species. This estimation represents a significant improvement in our knowledge of the population genetics of the African green monkey, one of the most important nonhuman primate models in biomedical research. Furthermore, by comparing mutation rates in Old World monkeys with the only other direct estimates in primates to date-humans and chimpanzees-it is possible to uniquely address how mutation rates have evolved over longer time scales. While the estimated spontaneous mutation rate for African green monkeys is slightly lower than the rate of 1.2 × 10^-8 per base pair per generation reported in chimpanzees, it is similar to the lower range of rates of 0.96 × 10^-8 -1.28 × 10^-8 per base pair per generation recently estimated from whole genome pedigrees in humans. This result suggests a long-term constraint on mutation rate that is quite different from similar evidence pertaining to recombination rate evolution in primates.

The Demographic and Adaptive History of the African Green Monkey

Relatively little is known about the evolutionary history of the African green monkey (genus Chlorocebus) due to the lack of sampled polymorphism data from wild populations. Yet, this characterization of genetic diversity is not only critical for a better understanding of their own history, but also for human biomedical research given that they are one of the most widely used primate models. Here, I analyze the demographic and selective history of the African green monkey, utilizing one of the most comprehensive catalogs of wild genetic diversity to date, consisting of 1,795,643 autosomal single nucleotide polymorphisms in 25 individuals, representing all five major populations: C. a. aethiops, C. a. cynosurus, C. a. pygerythrus, C. a. sabaeus, and C. a tantalus. Assuming a mutation rate of 5.9 × 10-9 per base pair per generation and a generation time of 8.5 years, divergence time estimates range from 523 to 621 kya for the basal split of C. a. aethiops from the other four populations. Importantly, the resulting tree characterizing the relationship and split-times between these populations differs significantly from that presented in the original genome paper, owing to their neglect of within-population variation when calculating between population-divergence. In addition, I find that the demographic history of all five populations is well explained by a model of population fragmentation and isolation, rather than novel colonization events. Finally, utilizing these demographic models as a null, I investigate the selective history of the populations, identifying candidate regions potentially related to adaptation in response to pathogen exposure.

Characterizing human cytomegalovirus reinfection in congenitally infected infants: an evolutionary perspective

Given the strong selective pressures often faced by populations when colonizing a novel habitat, the level of variation present on which selection may act is an important indicator of adaptive potential. While often discussed in an ecological context, this notion is also highly relevant in our clinical understanding of viral infection, in which the novel habitat is a new host. Thus, quantifying the factors determining levels of variation is of considerable importance for the design of improved treatment strategies. Here, we focus on such a quantification of human cytomegalovirus (HCMV) - a virus which can be transmitted across the placenta, resulting in foetal infection that can potentially cause severe disease in multiple organs. Recent studies using genomewide sequencing data have demonstrated that viral populations in some congenitally infected infants diverge rapidly over time and between tissue compartments within individuals, while in other infants, the populations remain highly stable. Here, we investigate the underlying causes of these extreme differences in observed intrahost levels of variation by estimating the underlying demographic histories of infection. Importantly, reinfection (i.e. population admixture) appears to be an important, and previously unappreciated, player. We highlight illustrative examples likely to represent a single-population transmission from a mother during pregnancy and multiple-population transmissions during pregnancy and after birth.

The ultimate and proximate mechanisms driving the evolution of long tails in forest deer mice

Understanding both the role of selection in driving phenotypic change and its underlying genetic basis remain major challenges in evolutionary biology. Here, we use modern tools to revisit a classic system of local adaptation in the North American deer mouse, Peromyscus maniculatus, which occupies two main habitat types: prairie and forest. Using historical collections, we find that forest‐dwelling mice have longer tails than those from nonforested habitat, even when we account for individual and population relatedness. Using genome‐wide SNP data, we show that mice from forested habitats in the eastern and western parts of their range form separate clades, suggesting that increased tail length evolved independently. We find that forest mice in the east and west have both more and longer caudal vertebrae, but not trunk vertebrae, than nearby prairie forms. By intercrossing prairie and forest mice, we show that the number and length of caudal vertebrae are not correlated in this recombinant population, indicating that variation in these traits is controlled by separate genetic loci. Together, these results demonstrate convergent evolution of the long‐tailed forest phenotype through two distinct genetic mechanisms, affecting number and length of vertebrae, and suggest that these morphological changes—either independently or together—are adaptive.

The Impact of Linked Selection in Chimpanzees: A Comparative Study

Levels of nucleotide diversity vary greatly across the genomes of most species owing to multiple factors. These include variation in the underlying mutation rates, as well as the effects of both direct and linked selection. Fundamental to interpreting the relative importance of these forces is the common observation of a strong positive correlation between nucleotide diversity and recombination rate. While indeed observed in humans, the interpretation of this pattern has been difficult in the absence of high-quality polymorphism data and recombination maps in closely related species. Here, we characterize genetic features driving nucleotide diversity in Western chimpanzees using a recently generated whole genome polymorphism data set. Our results suggest that recombination rate is the primary predictor of nucleotide variation with a strongly positive correlation. In addition, telomeric distance, regional GC-content, and regional CpG-island content are strongly negatively correlated with variation. These results are compared with humans, with both similarities and differences interpreted in the light of the estimated effective population sizes of the two species as well as their strongly differing recent demographic histories.

Inferring the age of a fixed beneficial allele

Estimating the age and strength of beneficial alleles is central to understanding how adaptation proceeds in response to changing environmental conditions. Several haplotype-based estimators exist for inferring the age of segregating beneficial mutations. Here, we develop an approximate Bayesian-based approach that rather estimates these parameters for fixed beneficial mutations in single populations. We integrate a range of existing diversity, site frequency spectrum, haplotype- and linkage disequilibrium-based summary statistics. We show that for strong selective sweeps on de novo mutations the method can estimate allele age and selection strength even in nonequilibrium demographic scenarios. We extend our approach to models of selection on standing variation, and co-infer the frequency at which selection began to act upon the mutation. Finally, we apply our method to estimate the age and selection strength of a previously identified mutation underpinning cryptic colour adaptation in a wild deer mouse population, and compare our findings with previously published estimates as well as with geological data pertaining to the presumed shift in selective pressure.

The population genomics of rapid adaptation: disentangling signatures of selection and demography in white sands lizards

Understanding the process of adaptation during rapid environmental change remains one of the central focal points of evolutionary biology. The recently formed White Sands system of southern New Mexico offers an outstanding example of rapid adaptation, with a variety of species having rapidly evolved blanched forms on the dunes that contrast with their close relatives in the surrounding dark soil habitat. In this study, we focus on two of the White Sands lizard species, Sceloporus cowlesi and Aspidoscelis inornata, for which previous research has linked mutations in the melanocortin-1 receptor gene (Mc1r) to blanched coloration. We sampled populations both on and off the dunes and used a custom sequence capture assay based on probed fosmid libraries to obtain >50 kb of sequence around Mc1r and hundreds of other random genomic locations. We then used model-based statistical inference methods to describe the demographic and adaptive history characterizing the colonization of White Sands. We identified a number of similarities between the two focal species, including strong evidence of selection in the blanched populations in the Mc1r region. We also found important differences between the species, suggesting different colonization times, different genetic architecture underlying the blanched phenotype and different ages of the beneficial alleles. Finally, the beneficial allele is dominant in S. cowlesi and recessive in A. inornata, allowing for a rare empirical test of theoretically expected patterns of selective sweeps under these differing models.