A single amino acid polymorphism in natural Metchnikowin alleles of Drosophila results in systemic immunity and life history tradeoffs.

Antimicrobial peptides (AMPs) are at the interface of interactions between hosts and microbes and are therefore expected to be rapidly evolving in a coevolutionary arms race with pathogens. In contrast, previous work demonstrated that insect AMPs tend to evolve more slowly than the genome average. Metchikowin (Mtk) is a Drosophila AMP that has a single amino acid residue that segregates as either proline (P) or arginine (R) in populations of four different species, some of which diverged more than 10 million years ago. These results suggest that there is a distinct functional importance to each allele. The most likely hypotheses are driven by two main questions: does each allele have a different efficacy against different specific pathogens (specificity hypothesis)? Or, is one allele a more potent antimicrobial, but with a host fitness cost (autoimmune hypothesis)? To assess their functional differences, we created D. melanogaster lines with the P allele, R allele, or Mtk null mutation using CRISPR/Cas9 genome editing and performed a series of life history and infection assays to assess them. In males, testing of systemic immune responses to a repertoire of bacteria and fungi demonstrated that the R allele performs as well or better than the P and null alleles with most infections. Females show some results that contrast with males, with Mtk alleles either not contributing to survival or with the P allele outperforming the R allele. In addition, measurements of life history traits demonstrate that the R allele is more costly in the absence of infection for both sexes. These results are consistent with both the specificity hypothesis (either allele can perform better against certain pathogens depending on context), and the autoimmune hypothesis (the R allele is generally the more potent antimicrobial in males, and carries a fitness cost). These results provide strong in vivo evidence that differential fitness with or without infection and sex-based functional differences in alleles may be adaptive mechanisms of maintaining immune gene polymorphisms in contrast with expectations of rapid evolution. Therefore, a complex interplay of forces including pathogen species and host sex may lead to balancing selection for immune genotypes. Strikingly, this selection may act on even a single amino acid polymorphism in an AMP.

SARS-CoV-2 Mac1 is required for IFN antagonism and efficient virus replication in cell culture and in mice.

Several coronavirus (CoV) encoded proteins are being evaluated as targets for antiviral therapies for COVID-19. Included in these drug targets is the conserved macrodomain, or Mac1, an ADP-ribosylhydrolase and ADP-ribose binding protein encoded as a small domain at the N terminus of nonstructural protein 3. Utilizing point mutant recombinant viruses, Mac1 was shown to be critical for both murine hepatitis virus (MHV) and severe acute respiratory syndrome (SARS)-CoV virulence. However, as a potential drug target, it is imperative to understand how a complete Mac1 deletion impacts the replication and pathogenesis of different CoVs. To this end, we created recombinant bacterial artificial chromosomes (BACs) containing complete Mac1 deletions (ΔMac1) in MHV, MERS-CoV, and SARS-CoV-2. While we were unable to recover infectious virus from MHV or MERS-CoV ΔMac1 BACs, SARS-CoV-2 ΔMac1 was readily recovered from BAC transfection, indicating a stark difference in the requirement for Mac1 between different CoVs. Furthermore, SARS-CoV-2 ΔMac1 replicated at or near wild-type levels in multiple cell lines susceptible to infection. However, in a mouse model of severe infection, ΔMac1 was quickly cleared causing minimal pathology without any morbidity. ΔMac1 SARS-CoV-2 induced increased levels of interferon (IFN) and IFN-stimulated gene expression in cell culture and mice, indicating that Mac1 blocks IFN responses which may contribute to its attenuation. ΔMac1 infection also led to a stark reduction in inflammatory monocytes and neutrophils. These results demonstrate that Mac1 only minimally impacts SARS-CoV-2 replication, unlike MHV and MERS-CoV, but is required for SARS-CoV-2 pathogenesis and is a unique antiviral drug target.

The genetic basis of variation in immune defense against Lysinibacillus fusiformis infection in Drosophila melanogaster.

The genetic causes of phenotypic variation often differ depending on the population examined, particularly if the populations were founded by relatively small numbers of genotypes. Similarly, the genetic causes of phenotypic variation among similar traits (resistance to different xenobiotic compounds or pathogens) may also be completely different or only partially overlapping. Differences in genetic causes for variation in the same trait among populations suggests context dependence for how selection acts on those traits. Similarities in the genetic causes of variation for different traits, on the other hand, suggests pleiotropy which would also influence how natural selection shapes variation in a trait. We characterized immune defense against a natural Drosophila pathogen, the Gram-positive bacterium Lysinibacillus fusiformis, in three different populations and found almost no overlap in the genetic architecture of variation in survival post infection. However, when comparing our results to a similar experiment with the fungal pathogen, B. bassiana, we found a convincing shared QTL peak for both pathogens. This peak contains the Bomanin cluster of Drosophila immune effectors. Loss of function mutants and RNAi knockdown experiments confirms a role of some of these genes in immune defense against both pathogens. This suggests that natural selection may act on the entire cluster of Bomanin genes (and the linked region under the QTL) or specific peptides for specific pathogens.

Positive Selection and Horizontal Gene Transfer in the Genome of a Male-Killing Wolbachia.

Wolbachia are a genus of widespread bacterial endosymbionts in which some strains can hijack or manipulate arthropod host reproduction. Male killing is one such manipulation in which these maternally transmitted bacteria benefit surviving daughters in part by removing competition with the sons for scarce resources. Despite previous findings of interesting genome features of microbial sex ratio distorters, the population genomics of male-killers remain largely uncharacterized. Here, we uncover several unique features of the genome and population genomics of four Arizonan populations of a male-killing Wolbachia strain, wInn, that infects mushroom-feeding Drosophila innubila. We first compared the wInn genome with other closely related Wolbachia genomes of Drosophila hosts in terms of genome content and confirm that the wInn genome is largely similar in overall gene content to the wMel strain infecting D. melanogaster. However, it also contains many unique genes and repetitive genetic elements that indicate lateral gene transfers between wInn and non-Drosophila eukaryotes. We also find that, in line with literature precedent, genes in the Wolbachia prophage and Octomom regions are under positive selection. Of all the genes under positive selection, many also show evidence of recent horizontal transfer among Wolbachia symbiont genomes. These dynamics of selection and horizontal gene transfer across the genomes of several Wolbachia strains and diverse host species may be important underlying factors in Wolbachia's success as a male-killer of divergent host species.

Caenorhabditis elegans ETR-1/CELF has broad effects on the muscle cell transcriptome, including genes that regulate translation and neuroblast migration.

Migration of neuroblasts and neurons from their birthplace is central to the formation of neural circuits and networks. ETR-1 is the Caenorhabditis elegans homolog of the CELF1 (CUGBP, ELAV-like family 1) RNA-processing factor involved in neuromuscular disorders. etr-1 regulates body wall muscle differentiation. Our previous work showed that etr-1 in muscle has a non-autonomous role in neuronal migration, suggesting that ETR-1 is involved in the production of a signal emanating from body wall muscle that controls neuroblast migration and that interacts with Wnt signaling. etr-1 is extensively alternatively-spliced, and we identified the viable etr-1(lq61) mutant, caused by a stop codon in alternatively-spliced exon 8 and only affecting etr-1 isoforms containing exon 8. We took advantage of viable etr-1(lq61) to identify potential RNA targets of ETR-1 in body wall muscle using a combination of fluorescence activated cell sorting (FACS) of body wall muscles from wild-type and etr-1(lq61) and subsequent RNA-seq. This analysis revealed genes whose splicing and transcript levels were controlled by ETR-1 exon 8 isoforms, and represented a broad spectrum of genes involved in muscle differentiation, myofilament lattice structure, and physiology. Genes with transcripts underrepresented in etr-1(lq61) included those involved in ribosome function and translation, similar to potential CELF1 targets identified in chick cardiomyocytes. This suggests that at least some targets of ETR-1 might be conserved in vertebrates, and that ETR-1 might generally stimulate translation in muscles. As proof-of-principle, a functional analysis of a subset of ETR-1 targets revealed genes involved in AQR and PQR neuronal migration. One such gene, lev-11/tropomyosin, requires ETR-1 for alternative splicing, and another, unc-52/perlecan, requires ETR-1 for the production of long isoforms containing 3' exons. In sum, these studies identified gene targets of ETR-1/CELF1 in muscles, which included genes involved in muscle development and physiology, and genes with novel roles in neuronal migration.

Hoisted with his own petard: How sex-ratio meiotic drive in Drosophila affinis creates resistance alleles that limit its spread.

Meiotic drivers are selfish genetic elements that tinker with gametogenesis to bias their own transmission into the next generation of offspring. Such tinkering can have significant consequences on gametogenesis and end up hampering the spread of the driver. In Drosophila affinis, sex-ratio meiotic drive is caused by an X-linked complex that, when in males with a susceptible Y chromosome, results in broods that are typically more than 95% female. Interestingly, D. affinis males lacking a Y chromosome (XO) are fertile and males with the meiotic drive X and no Y produce only sons-effectively reversing the sex-ratio effect. Here, we show that meiotic drive dramatically increases the rate of nondisjunction of the Y chromosome (at least 750X), meaning that the driver is creating resistant alleles through the process of driving. We then model how the O might influence the spread, dynamics and equilibrium of the sex-ratio X chromosome. We find that the O can prevent the spread or reduce the equilibrium frequency of the sex-ratio X chromosome, and it can even lead to oscillations in frequency. Finally, with reasonable parameters, the O is unlikely to lead to the loss of the Y chromosome, but we discuss how it might lead to sex-chromosome turnover indirectly.

The phage gene wmk is a candidate for male killing by a bacterial endosymbiont.

Wolbachia are the most widespread maternally-transmitted bacteria in the animal kingdom. Their global spread in arthropods and varied impacts on animal physiology, evolution, and vector control are in part due to parasitic drive systems that enhance the fitness of infected females, the transmitting sex of Wolbachia. Male killing is one common drive mechanism wherein the sons of infected females are selectively killed. Despite decades of research, the gene(s) underlying Wolbachia-induced male killing remain unknown. Here using comparative genomic, transgenic, and cytological approaches in fruit flies, we identify a candidate gene in the eukaryotic association module of Wolbachia prophage WO, termed WO-mediated killing (wmk), which transgenically causes male-specific lethality during early embryogenesis and cytological defects typical of the pathology of male killing. The discovery of wmk establishes new hypotheses for the potential role of phage genes in sex-specific lethality, including the control of arthropod pests and vectors.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species.

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

The Genome of Drosophila innubila Reveals Lineage-Specific Patterns of Selection in Immune Genes.

Pathogenic microbes can exert extraordinary evolutionary pressure on their hosts. They can spread rapidly and sicken or even kill their host to promote their own proliferation. Because of this strong selective pressure, immune genes are some of the fastest evolving genes across metazoans, as highlighted in mammals and insects. Drosophila melanogaster serves as a powerful model for studying host/pathogen evolution. While Drosophila melanogaster are frequently exposed to various pathogens, little is known about D. melanogaster's ecology, or if they are representative of other Drosophila species in terms of pathogen pressure. Here, we characterize the genome of Drosophila innubila, a mushroom-feeding species highly diverged from D. melanogaster and investigate the evolution of the immune system. We find substantial differences in the rates of evolution of immune pathways between D. innubila and D. melanogaster. Contrasting what was previously found for D. melanogaster, we find little evidence of rapid evolution of the antiviral RNAi genes and high rates of evolution in the Toll pathway. This suggests that, while immune genes tend to be rapidly evolving in most species, the specific genes that are fastest evolving may depend either on the pathogens faced by the host and/or divergence in the basic architecture of the host's immune system.

The Population Genetics of Evolutionary Rescue in Diploids: X Chromosomal versus Autosomal Rescue.

Most population genetic theory assumes that populations adapt to an environmental change without a change in population size. However, environmental changes might be so severe that populations decline in size and, without adaptation, become extinct. This "evolutionary rescue" scenario differs from traditional models of adaptation in that rescue involves a race between adaptation and extinction. While most previous work has usually focused on models of evolutionary rescue in haploids, here we consider diploids. In many species, diploidy introduces a novel feature into adaptation: adaptive evolution might occur either on sex chromosomes or on autosomes. Previous studies of nonrescue adaptation revealed that the relative rates of adaptation on the X chromosome versus autosomes depend on the dominance of beneficial mutations, reflecting differences in effective population size and the efficacy of selection. Here, we extend these results to evolutionary rescue and find that, given equal-sized chromosomes, there is greater parameter space in which the X is more likely to contribute to adaptation than the autosomes relative to standard nonrescue models. We also discuss how subtle effects of dominance can increase the chance of evolutionary rescue in diploids when absolute heterozygote fitness is close to 1. These effects do not arise in standard nonrescue models.

Resistance to natural and synthetic gene drive systems.

Scientists are rapidly developing synthetic gene drive elements intended for release into natural populations. These are intended to control or eradicate disease vectors and pests, or to spread useful traits through wild populations for disease control or conservation purposes. However, a crucial problem for gene drives is the evolution of resistance against them, preventing their spread. Understanding the mechanisms by which populations might evolve resistance is essential for engineering effective gene drive systems. This review summarizes our current knowledge of drive resistance in both natural and synthetic gene drives. We explore how insights from naturally occurring and synthetic drive systems can be integrated to improve the design of gene drives, better predict the outcome of releases and understand genomic conflict in general.

An assessment of the immune costs associated with meiotic drive elements in Drosophila.

Most organisms are constantly adapting to pathogens and parasites that exploit their host for their own benefit. Less studied, but perhaps more ubiquitous, are intragenomic parasites or selfish genetic elements. These include transposable elements, selfish B chromosomes and meiotic drivers that promote their own replication without regard to fitness effects on hosts. Therefore, intragenomic parasites are also a constant evolutionary pressure on hosts. Gamete-killing meiotic drive elements are often associated with large chromosomal inversions that reduce recombination between the drive and wild-type chromosomes. This reduced recombination is thought to reduce the efficacy of selection on the drive chromosome and allow for the accumulation of deleterious mutations. We tested whether gamete-killing meiotic drive chromosomes were associated with reduced immune defence against two bacterial pathogens in three species of Drosophila. We found little evidence of reduced immune defence in lines with meiotic drive. One line carrying the Drosophila melanogaster autosomal Segregation Distorter did show reduced defence, but we were unable to attribute that reduced defence to either genotype or immune gene expression differences. Our results suggest that though gamete-killing meiotic drive chromosomes probably accumulate deleterious mutations, those mutations do not result in reduced capacity for immune defence.

Efflux Pumps in Chromobacterium Species Increase Antibiotic Resistance and Promote Survival in a Coculture Competition Model.

Members of the Chromobacterium genus include opportunistic but often-fatal pathogens and soil saprophytes with highly versatile metabolic capabilities. In previous studies of Chromobacterium subtsugae (formerly C. violaceum) strain CV017, we identified a resistance nodulation division (RND)-family efflux pump (CdeAB-OprM) that confers resistance to several antibiotics, including the bactobolin antibiotic produced by the soil saprophyte Burkholderia thailandensis Here, we show the cdeAB-oprM genes increase C. subtsugae survival in a laboratory competition model with B. thailandensis We also demonstrate that adding sublethal bactobolin concentrations to the coculture increases C. subtsugae survival, but this effect is not through CdeAB-OprM. Instead, the increased survival requires a second, previously unreported pump we call CseAB-OprN. We show that in cells exposed to sublethal bactobolin concentrations, the cseAB-oprN genes are transcriptionally induced, and this corresponds to an increase in bactobolin resistance. Induction of this pump is highly specific and sensitive to bactobolin, while CdeAB-OprM appears to have a broader range of antibiotic recognition. We examine the distribution of cseAB-oprN and cdeAB-oprM gene clusters in members of the Chromobacterium genus and find the cseAB-oprN genes are limited to the nonpathogenic C. subtsugae strains, whereas the cdeAB-oprM genes are more widely distributed among members of the Chromobacterium genus. Our results provide new information on the antibiotic resistance mechanisms of Chromobacterium species and highlight the importance of efflux pumps for saprophytic bacteria existing in multispecies communities.IMPORTANCE Antibiotic efflux pumps are best known for increasing antibiotic resistance of pathogens; however, the role of these pumps in saprophytes is much less well defined. This study describes two predicted efflux pump gene clusters in the Chromobacterium genus, which is comprised of both nonpathogenic saprophytes and species that cause highly fatal human infections. One of the predicted efflux pump clusters is present in every member of the Chromobacterium genus and increases resistance to a broad range of antibiotics. The other gene cluster has more narrow antibiotic specificity and is found only in Chromobacterium subtsugae, a subset of entirely nonpathogenic species. We demonstrate the role of both pumps in increasing antibiotic resistance and demonstrate the importance of efflux-dependent resistance induction for C. subtsugae survival in a dual-species competition model. These results have implications for managing antibiotic-resistant Chromobacterium infections and for understanding the evolution of efflux pumps outside the host.

Baculovirus Molecular Evolution via Gene Turnover and Recurrent Positive Selection of Key Genes.

Hosts and viruses are locked in an evolutionary arms race. Hosts are constantly evolving to suppress virulence and replication, while viruses, which are reliant on host machinery for survival and reproduction, develop counterstrategies to escape this immune defense. Viruses must also adapt to novel conditions while establishing themselves in a host species. Both processes provide strong selection for viral adaptation. Understanding adaptive evolution in insect viruses can help us to better understand adaptive evolution in general and is important due to the use of these viruses as biocontrol agents and for protecting ecologically or economically important species from outbreaks. Here we examine the molecular evolution of baculoviruses and nudiviruses, a group of insect-infecting viruses with key roles in biocontrol. We looked for signatures of selection between genomes of baculoviruses infecting a range of species and within a population of baculoviruses. Both analyses found only a few strong signatures of positive selection, primarily in replication- and transcription-associated genes and several structural protein genes. In both analyses, we detected a conserved complex of genes, including the helicase gene, showing consistently high levels of adaptive evolution, suggesting that they may be key in antagonistic coevolution to escape host suppression. These genes are integral to the baculovirus life cycle and may be good focal genes for developing baculoviruses as effective biocontrol agents or for targeting baculoviruses infecting ecologically relevant species. Recombination and complex genomes make evolution in these double-stranded DNA viruses more efficient than that in smaller RNA viruses with error-prone replication, as seen via signatures of selection in specific genes within a population of baculoviruses.IMPORTANCE Most viral evolutionary studies focus on RNA viruses. While these viruses cause many human and animal diseases, such studies leave us with a lesser understanding of how DNA viruses adapt to hosts and how the host responds to these pathogens. In this paper, we focus on the evolution of baculoviruses, a group of insect-infecting DNA viruses, many of which have been used in biocontrol. We find that most of the genome is under purifying selection, with only a few key genes evolving adaptively. Our results provide a glimpse into how DNA viruses differ from RNA viruses in their evolutionary dynamics and identify genes that are key to DNA virus adaptation, improving our understanding of how this group of pathogens evolves.

A fast-evolving X-linked duplicate of importin-α2 is overexpressed in sex-ratio drive in Drosophila neotestacea.

Selfish genetic elements that manipulate gametogenesis to achieve a transmission advantage are known as meiotic drivers. Sex-ratio X chromosomes (SR) are meiotic drivers that prevent the maturation of Y-bearing sperm in male carriers to result in the production of mainly female progeny. The spread of an SR chromosome can affect host genetic diversity and genome evolution, and can even cause host extinction if it reaches sufficiently high prevalence. Meiotic drivers have evolved independently many times, though only in a few cases is the underlying genetic mechanism known. In this study we use a combination of transcriptomics and population genetics to identify widespread expression differences between the standard (ST) and sex-ratio (SR) X chromosomes of the fly Drosophila neotestacea. We found the X chromosome is enriched for differentially expressed transcripts and that many of these X-linked differentially expressed transcripts had elevated Ka /Ks values between ST and SR, indicative of potential functional differences. We identified a set of candidate transcripts, including a testis-specific, X-linked duplicate of the nuclear transport gene importin-α2 that is overexpressed in SR. We find suggestions of positive selection in the lineage leading to the duplicate and that its molecular evolutionary patterns are consistent with relaxed purifying selection in ST. As these patterns are consistent with involvement in the mechanism of drive in this species, this duplicate is a strong candidate worthy of further functional investigation. Nuclear transport may be a common target for genetic conflict, as the mechanism of the autosomal Segregation Distorter drive system in D. melanogaster involves the same pathway.

The complex contributions of genetics and nutrition to immunity in Drosophila melanogaster.

Both malnutrition and undernutrition can lead to compromised immune defense in a diversity of animals, and "nutritional immunology" has been suggested as a means of understanding immunity and determining strategies for fighting infection. The genetic basis for the effects of diet on immunity, however, has been largely unknown. In the present study, we have conducted genome-wide association mapping in Drosophila melanogaster to identify the genetic basis for individual variation in resistance, and for variation in immunological sensitivity to diet (genotype-by-environment interaction, or GxE). D. melanogaster were reared for several generations on either high-glucose or low-glucose diets and then infected with Providencia rettgeri, a natural bacterial pathogen of D. melanogaster. Systemic pathogen load was measured at the peak of infection intensity, and several indicators of nutritional status were taken from uninfected flies reared on each diet. We find that dietary glucose level significantly alters the quality of immune defense, with elevated dietary glucose resulting in higher pathogen loads. The quality of immune defense is genetically variable within the sampled population, and we find genetic variation for immunological sensitivity to dietary glucose (genotype-by-diet interaction). Immune defense was genetically correlated with indicators of metabolic status in flies reared on the high-glucose diet, and we identified multiple genes that explain variation in immune defense, including several that have not been previously implicated in immune response but which are confirmed to alter pathogen load after RNAi knockdown. Our findings emphasize the importance of dietary composition to immune defense and reveal genes outside the conventional "immune system" that can be important in determining susceptibility to infection. Functional variation in these genes is segregating in a natural population, providing the substrate for evolutionary response to pathogen pressure in the context of nutritional environment.

The population genetics of evolutionary rescue.

Evolutionary rescue occurs when a population that is threatened with extinction by an environmental change adapts to the change sufficiently rapidly to survive. Here we extend the mathematical theory of evolutionary rescue. In particular, we model evolutionary rescue to a sudden environmental change when adaptation involves evolution at a single locus. We consider adaptation using either new mutations or alleles from the standing genetic variation that begin rare. We obtain several results: i) the total probability of evolutionary rescue from either new mutation or standing variation; ii) the conditions under which rescue is more likely to involve a new mutation versus an allele from the standing genetic variation; iii) a mathematical description of the U-shaped curve of total population size through time, conditional on rescue; and iv) the time until the average population size begins to rebound as well as the minimal expected population size experienced by a rescued population. Our analysis requires taking into account a subtle population-genetic effect (familiar from the theory of genetic hitchhiking) that involves "oversampling" of those lucky alleles that ultimately sweep to high frequency. Our results are relevant to conservation biology, experimental microbial evolution, and medicine (e.g., the dynamics of antibiotic resistance).

The Genetic Basis of Natural Variation in Drosophila (Diptera: Drosophilidae) Virgin Egg Retention.

Drosophila melanogaster is able to thrive in harsh northern climates through adaptations in life-history traits and physiological mechanisms that allow for survival through the winter. We examined the genetic basis of natural variation in one such trait, female virgin egg retention, which was previously shown to vary clinally and seasonally. To further our understanding of the genetic basis and evolution of virgin egg retention, we performed a genome-wide association study (GWAS) using the previously sequenced Drosophila Genetic Reference Panel (DGRP) mapping population. We found 29 single nucleotide polymorphisms (SNPs) associated with virgin egg retention and assayed 6 available mutant lines, each harboring a mutation in a candidate gene, for effects on egg retention time. We found that four out of the six mutant lines had defects in egg retention time as compared with the respective controls: mun, T48, Mes-4, and Klp67A Surprisingly, none of these genes has a recognized role in ovulation control, but three of the four genes have known effects on fertility or have high expression in the ovaries. We also found that the SNP set associated with egg retention time was enriched for clinal SNPs. The majority of clinal SNPs had alleles associated with longer egg retention present at higher frequencies in higher latitudes. Our results support previous studies that show higher frequency of long retention times at higher latitude, providing evidence for the adaptive value of virgin egg-retention.

Modeling the Evolution of S. pombe Populations with Multiple Killer Meiotic Drivers.

Meiotic drivers are selfish genetic loci that can be transmitted to more than half of the viable gametes produced by a heterozygote. This biased transmission gives meiotic drivers an evolutionary advantage that can allow them to spread over generations until all members of a population carry the driver. This evolutionary power can also be exploited to modify natural populations using synthetic drivers known as 'gene drives.' Recently, it has become clear that natural drivers can spread within genomes to birth multicopy gene families. To understand intragenomic spread of drivers, we model the evolution of two or more distinct meiotic drivers in a population. We employ the wtf killer meiotic drivers from Schizosaccharomyces pombe, which are multicopy in all sequenced isolates, as models. We find that a duplicate wtf driver identical to the parent gene can spread in a population unless, or until, the original driver is fixed. When the duplicate driver diverges to be distinct from the parent gene, we find that both drivers spread to fixation under most conditions, but both drivers can be lost under some conditions. Finally, we show that stronger drivers make weaker drivers go extinct in most, but not all, polymorphic populations with absolutely linked drivers. These results reveal the strong potential for natural meiotic drive loci to duplicate and diverge within genomes. Our findings also highlight duplication potential as a factor to consider in the design of synthetic gene drives.

Meiotic drive, postzygotic isolation, and the Snowball Effect.

As populations diverge, they accumulate incompatibilities which reduce gene flow and facilitate the formation of new species. Simple models suggest that the genes that cause Dobzhansky-Muller incompatibilities should accumulate at least as fast as the square of the number of substitutions between taxa, the so-called snowball effect. We show, however, that in the special- but possibly common- case in which hybrid sterility is due primarily to cryptic meiotic (gametic) drive, the number of genes that cause postzygotic isolation may increase nearly linearly with the number of substitutions between species.