Are Lethal Alleles Too Abundant in Humans?

Across species, many individuals carry one or more recessive lethal alleles, posing an evolutionary conundrum for their persistence. Using a population genomic approach, Amorim et al. studied the abundance of lethal disease-causing mutations in humans and found that, while appearing more common than expected, most may nonetheless persist at frequencies predicted by mutation-selection balance.

Variation in Recombination Rate: Adaptive or Not?

Rates of meiotic recombination are widely variable both within and among species. However, the functional significance of this variation remains largely unknown. Is the observed within-species variation in recombination rate adaptive? Recent work has revealed new insight into the scale and scope of population-level variation in recombination rate. These data indicate that the magnitude of within-population variation in recombination is similar among taxa. The apparent similarity of the variance in recombination rate among individuals between distantly related species suggests that the relative costs and benefits of recombination that establish the upper and lower bounds may be similar across species. Here we review the current data on intraspecific variation in recombination rate and discuss the molecular and evolutionary costs and benefits of recombination frequency. We place this variation in the context of adaptation and highlight the need for more empirical studies focused on the adaptive value of variation in recombination rate.

Pervasive gene conversion in chromosomal inversion heterozygotes.

Chromosomal inversions shape recombination landscapes, and species differing by inversions may exhibit reduced gene flow in these regions of the genome. Though single crossovers within inversions are not usually recovered from inversion heterozygotes, the recombination barrier imposed by inversions is nuanced by noncrossover gene conversion. Here, we provide a genomewide empirical analysis of gene conversion rates both within species and in species hybrids. We estimate that gene conversion occurs at a rate of 1 × 10-5 to 2.5 × 10-5 converted sites per bp per generation in experimental crosses within Drosophila pseudoobscura and between D. pseudoobscura and its naturally hybridizing sister species D. persimilis. This analysis is the first direct empirical assessment of gene conversion rates within inversions of a species hybrid. Our data show that gene conversion rates in interspecies hybrids are at least as high as within-species estimates of gene conversion rates, and gene conversion occurs regularly within and around inverted regions of species hybrids, even near inversion breakpoints. We also found that several gene conversion events appeared to be mitotic rather than meiotic in origin. Finally, we observed that gene conversion rates are higher in regions of lower local sequence divergence, yet our observed gene conversion rates in more divergent inverted regions were at least as high as in less divergent collinear regions. Given our observed high rates of gene conversion despite the sequence differentiation between species, especially in inverted regions, gene conversion has the potential to reduce the efficacy of inversions as barriers to recombination over evolutionary time.

Gene conversion and linkage: effects on genome evolution and speciation.

Crossing over is well known to have profound effects on patterns of genetic diversity and genome evolution. Far less direct attention has been paid to another distinct outcome of meiotic recombination: noncrossover gene conversion (NCGC). Crossing over and NCGC both shuffle combinations of alleles, and this degradation of linkage disequilibrium (LD) has major evolutionary consequences, ranging from immediate effects on nucleotide diversity to long-term consequences that shape genome evolution, species formation and species persistence. Unlike simple crossing over, NCGC has the potential to alter allele frequencies. Gene conversion can also occur in genomic regions where crossing over does not, and it purportedly exhibits more uniform rates across genomes. Considerable progress has been made towards understanding the mechanisms of gene conversion, and this progress enables us to begin exploring how gene conversion affects processes such as molecular evolution and interspecies gene flow. These topics are timely with the recent shift in focus from a primarily neutral null model of molecular evolution and speciation to one incorporating base levels of selection, making it all the more crucial to understand the basis and evolutionary implications of linkage. Here, we discuss the impact of gene conversion on genome structure and evolution and the current methods for detecting these events. We provide a comprehensive review of how gene conversion breaks down LD and affects both short- and long-term evolutionary processes, and we contrast its impact to that expected from crossing over alone.

In memoriam: Richard G. Harrison - his life and legacy.

Richard G. Harrison passed away unexpectedly on April 12th, 2016. In this memoriam we pay tribute to the life and legacy of an extraordinary scientist, mentor, friend, husband, and father.

Recombination modulates how selection affects linked sites in Drosophila.

One of the most influential observations in molecular evolution has been a strong association between local recombination rate and nucleotide polymorphisms across the genome. This is interpreted as evidence for ubiquitous natural selection. The alternative explanation, that recombination is mutagenic, has been rejected by the absence of a similar association between local recombination rate and nucleotide divergence between species. However, many recent studies show that recombination rates are often very different even in closely related species, questioning whether an association between recombination rate and divergence between species has been tested satisfactorily. To circumvent this problem, we directly surveyed recombination across approximately 43% of the D. pseudoobscura physical genome in two separate recombination maps and 31% of the D. miranda physical genome, and we identified both global and local differences in recombination rate between these two closely related species. Using only regions with conserved recombination rates between and within species and accounting for multiple covariates, our data support the conclusion that recombination is positively related to diversity because recombination modulates Hill-Robertson effects in the genome and not because recombination is predominately mutagenic. Finally, we find evidence for dips in diversity around nonsynonymous substitutions. We infer that at least some of this reduction in diversity resulted from selective sweeps and examine these dips in the context of recombination rate.

Temporal Stability of Molecular Diversity Measures in Natural Populations of Drosophila pseudoobscura and Drosophila persimilis.

Many molecular ecological and evolutionary studies sample wild populations at a single point in time, but that data represents genetic variation from a potentially unrepresentative snapshot in time. Variation across time in genetic parameters may occur quickly in species that produce multiple generations of offspring per year. Here, we compare genetic diversity in wild caught populations of Drosophila persimilis and Drosophila pseudoobscura collected 16 years apart at the same time of year and same site at 4 X-linked and 2 mitochondrial loci to assess genetic stability. We found no major changes in nucleotide diversity in either species, but we observed a drastic shift in Tajima's D between D. pseudoobscura timepoints at 1 locus associated with increased abundance of a set of related haplotypes. Our data also suggests that D. persimilis may have recently accelerated its demographic expansion. While the changes we observed were modest, this study reinforces the importance of considering potential temporal variation in genetic parameters within single populations over short evolutionary timescales.

Interrogating the Roles of Mutation-Selection Balance, Heterozygote Advantage, and Linked Selection in Maintaining Recessive Lethal Variation in Natural Populations.

For nearly a century, evolutionary biologists have observed chromosomes that cause lethality when made homozygous persisting at surprisingly high frequencies (>25%) in natural populations of many species. The evolutionary forces responsible for the maintenance of such detrimental mutations have been heavily debated-are some lethal mutations under balancing selection? We suggest that mutation-selection balance alone cannot explain lethal variation in nature and the possibility that other forces play a role. We review the potential that linked selection in particular may drive maintenance of lethal alleles through associative overdominance or linkage to beneficial mutations or by reducing effective population size. Over the past five decades, investigation into this mystery has tapered. During this time, key scientific advances have provided the ability to collect more accurate data and analyze them in new ways, making the underlying genetic bases and evolutionary forces of lethal alleles timely for study once more.

Effects of premature termination codon polymorphisms in the Drosophila pseudoobscura subclade.

Premature termination codon (PTC) mutations can have dramatic effects--both adaptive and deleterious--on gene expression and function. Here, we examine the number and selective effects of PTC mutations within the Drosophila pseudoobscura subclade using 18 resequenced genomes aligned to the reference genome. We located and characterized 1,679 PTC mutations in 605 genes across each of these genomes relative to the D. pseudoobscura reference genome, and use RT-PCR to confirm transcription of a subset of these genes containing PTC mutations. We confirm previous findings that genes containing PTC mutations are less selectively constrained and less broadly expressed than non-PTC-containing genes, suggesting that the most of these mutations are at least mildly deleterious. Further, we find highly significant codon usage bias in regions downstream of the PTC in 38 of these PTC-containing genes, suggesting that some of these PTC mutations--if not alternatively spliced out of the transcript--have neutral effects. Ultimately, these analyzes support the view that the PTC mutations are mostly detrimental, but are nonetheless common enough in genomes that a subset could be effectively neutral.

Geographic selection in the small heat shock gene complex differentiating populations of Drosophila pseudoobscura.

Environmental temperature plays a crucial role in determining a species distribution and abundance by affecting individual physiological processes, metabolic activities, and developmental rates. Many studies have identified clinal variation in phenotypes associated with response to environmental stresses, but variation in traits associated with climatic adaptation directly attributed to sequence variation within candidate gene regions has been difficult to identify. Insect heat shock genes are possible agents of thermal tolerance because of their involvement in protein folding, traffic, protection, and renaturation at the cellular level in response to temperature stress. Previously, members of the Drosophila small heat shock protein (sHSP) complex (Hsp23, Hsp26, Hsp27, Hsp67Ba) have been implicated as candidate climatic adaptation genes; therefore, this research examines sequence variation at these genes in 2 distant populations of Drosophila pseudoobscura. Flies from Tempe, AZ (n = 30) and Cheney, WA (n = 17) were used in the study. We identify high differentiation in the heat-shock complex (F(ST) : 0.219**, 0.262*, 0.279***, 0.166 not significant) as compared with neighboring genes and Tajima's D values indicative of balancing selection (Mann-Whitney U = 38, n(1) = 10 n(2) = 4, P < 0.05 two-tailed), both of which are suggestive of such climatic adaptation.

The genomics of speciation in Drosophila: diversity, divergence, and introgression estimated using low-coverage genome sequencing.

In nature, closely related species may hybridize while still retaining their distinctive identities. Chromosomal regions that experience reduced recombination in hybrids, such as within inversions, have been hypothesized to contribute to the maintenance of species integrity. Here, we examine genomic sequences from closely related fruit fly taxa of the Drosophila pseudoobscura subgroup to reconstruct their evolutionary histories and past patterns of genic exchange. Partial genomic assemblies were generated from two subspecies of Drosophila pseudoobscura (D. ps.) and an outgroup species, D. miranda. These new assemblies were compared to available assemblies of D. ps. pseudoobscura and D. persimilis, two species with overlapping ranges in western North America. Within inverted regions, nucleotide divergence among each pair of the three species is comparable, whereas divergence between D. ps. pseudoobscura and D. persimilis in non-inverted regions is much lower and closer to levels of intraspecific variation. Using molecular markers flanking each of the major chromosomal inversions, we identify strong crossover suppression in F(1) hybrids extending over 2 megabase pairs (Mbp) beyond the inversion breakpoints. These regions of crossover suppression also exhibit the high nucleotide divergence associated with inverted regions. Finally, by comparison to a geographically isolated subspecies, D. ps. bogotana, our results suggest that autosomal gene exchange between the North American species, D. ps. pseudoobscura and D. persimilis, occurred since the split of the subspecies, likely within the last 200,000 years. We conclude that chromosomal rearrangements have been vital to the ongoing persistence of these species despite recent hybridization. Our study serves as a proof-of-principle on how whole genome sequencing can be applied to formulate and test hypotheses about species formation in lesser-known non-model systems.

Gene flow biases population genetic inference of recombination rate.

Accurate estimates of the rate of recombination are key to understanding a host of evolutionary processes as well as the evolution of the recombination rate itself. Model-based population genetic methods that infer recombination rates from patterns of linkage disequilibrium in the genome have become a popular method to estimate rates of recombination. However, these linkage disequilibrium-based methods make a variety of simplifying assumptions about the populations of interest that are often not met in natural populations. One such assumption is the absence of gene flow from other populations. Here, we use forward-time population genetic simulations of isolation-with-migration scenarios to explore how gene flow affects the accuracy of linkage disequilibrium-based estimators of recombination rate. We find that moderate levels of gene flow can result in either the overestimation or underestimation of recombination rates by up to 20-50% depending on the timing of divergence. We also find that these biases can affect the detection of interpopulation differences in recombination rate, causing both false positives and false negatives depending on the scenario. We discuss future possibilities for mitigating these biases and recommend that investigators exercise caution and confirm that their study populations meet assumptions before deploying these methods.

Translocation of Y-linked genes to the dot chromosome in Drosophila pseudoobscura.

One of the most striking cases of sex chromosome reorganization in Drosophila occurred in the lineage ancestral to Drosophila pseudoobscura, where there was a translocation of Y-linked genes to an autosome. These genes went from being present only in males, never recombining, and having an effective population size of 0.5N to a state of autosomal linkage, where they are passed through both sexes, may recombine, and their effective population size has quadrupled. These genes appear to be functional, and they underwent a drastic reduction in intron size after the translocation. A Y-autosome translocation may pose problems in meiosis if the rDNA locus responsible for X-Y pairing had also moved to an autosome. In this study, we demonstrate that the Y-autosome translocation moved Y-linked genes onto the dot chromosome, a small, mainly heterochromatic autosome with some sex chromosome-like properties. The rDNA repeats occur exclusively on the X chromosome in D. pseudoobscura, but we found that the new Y chromosome of this species harbors four clusters bearing only the intergenic spacer region (IGS) of the rDNA repeats. This arrangement appears analogous to the situation in Drosophila simulans, where X-rDNA to Y-IGS pairing could be responsible for X-Y chromosome pairing. We postulate that the nascent D. pseudoobscura Y chromosome acquired and amplified copies of the IGS, suggesting a potential mechanism for X-Y pairing in D. pseudoobscura.

Fine-scale mapping of recombination rate in Drosophila refines its correlation to diversity and divergence.

Regional rates of recombination often correlate with levels of nucleotide diversity, and either selective or neutral hypotheses can explain this relationship. Regional recombination rates also correlate with nucleotide differences between human and chimpanzee, consistent with models where recombination is mutagenic; however, a lack of correlation is observed in the Drosophila melanogaster group, consistent with models invoking natural selection. Here, we revisit the relationship among recombination, diversity, and interspecies difference by generating empirical estimates of these parameters in Drosophila pseudoobscura. To measure recombination rate, we genotyped 1,294 backcross hybrids at 50 markers across the largest assembled linkage group in this species. Genome-wide diversity was estimated by sequencing a second isolate of D. pseudoobscura at shallow coverage. Alignment to the sequenced genome of the closely related species, Drosophila persimilis, provided nucleotide site orthology. Our findings demonstrate that scale is critical in determining correlates to recombination rate: fine-scale cross-over rate estimates are far stronger predictors of both diversity and interspecies difference than broad-scale estimates. The correlation of fine-scale recombination rate to diversity and interspecies difference appears to be genome-wide, evidenced by examination of an X-linked region in greater detail. Because we observe a strong correlation of cross-over rate with interspecies difference, even after correcting for segregating ancestral variation, we suggest that both mutagenic and selective forces generate these correlations, the latter in regions of low crossing over. We propose that it is not cross-overs per se that are mutagenic, but rather repair of DNA double-strand break precursors via crossing over and gene conversion.

Inversions shape the divergence of Drosophila pseudoobscura and Drosophila persimilis on multiple timescales.

By shaping meiotic recombination, chromosomal inversions can influence genetic exchange between hybridizing species. Despite the recognized importance of inversions in evolutionary processes such as divergence and speciation, teasing apart the effects of inversions over time remains challenging. For example, are their effects on sequence divergence primarily generated through creating blocks of linkage disequilibrium prespeciation or through preventing gene flux after speciation? We provide a comprehensive look into the influence of inversions on gene flow throughout the evolutionary history of a classic system: Drosophila pseudoobscura and Drosophila persimilis. We use extensive whole-genome sequence data to report patterns of introgression and divergence with respect to chromosomal arrangements. Overall, we find evidence that inversions have contributed to divergence patterns between D. pseudoobscura and D. persimilis over three distinct timescales: (1) segregation of ancestral polymorphism early in the speciation process, (2) gene flow after the split of D. pseudoobscura and D. persimilis, but prior to the split of D. pseudoobscura subspecies, and (3) recent gene flow between sympatric D. pseudoobscura and D. persimilis, after the split of D. pseudoobscura subspecies. We discuss these results in terms of our understanding of evolution in this classic system and provide cautions for interpreting divergence measures in other systems.

Intraspecific Genetic Variation for Behavioral Isolation Loci in Drosophila.

Behavioral isolation is considered to be the primary mode of species isolation, and the lack of identification of individual genes for behavioral isolation has hindered our ability to address fundamental questions about the process of speciation. One of the major questions that remains about behavioral isolation is whether the genetic basis of isolation between species also varies within a species. Indeed, the extent to which genes for isolation may vary across a population is rarely explored. Here, we bypass the problem of individual gene identification by addressing this question using a quantitative genetic comparison. Using strains from eight different populations of Drosophila simulans, we genetically mapped the genomic regions contributing to behavioral isolation from their closely related sibling species, Drosophila mauritiana. We found extensive variation in the size of contribution of different genomic regions to behavioral isolation among the different strains, in the location of regions contributing to isolation, and in the ability to redetect loci when retesting the same strain.

PseudoBase: a genomic visualization and exploration resource for the Drosophila pseudoobscura subgroup.

Drosophila pseudoobscura is a classic model system for the study of evolutionary genetics and genomics. Given this long-standing interest, many genome sequences have accumulated for D. pseudoobscura and closely related species D. persimilis, D. miranda, and D. lowei. To facilitate the exploration of genetic variation within species and comparative genomics across species, we present PseudoBase, a database that couples extensive publicly available genomic data with simple visualization and query tools via an intuitive graphical interface, amenable for use in both research and educational settings. All genetic variation (SNPs and indels) within the database is derived from the same workflow, so variants are easily comparable across data sets. Features include an embedded JBrowse interface, ability to pull out alignments of individual genes/regions, and batch access for gene lists. Here, we introduce PseudoBase, and we demonstrate how this resource facilitates use of extensive genomic data from flies of the Drosophila pseudoobscura subgroup.

Homage to Felsenstein 1981, or why are there so few/many species?

If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build-up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between "two-allele" models, which are subject to this effect, and "one-allele" models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.

Genetic and evolutionary correlates of fine-scale recombination rate variation in Drosophila persimilis.

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed "recombinational neighborhoods," where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5-7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.