Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence.

A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the α-subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitution's effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases (See S1 Text for Spanish translation).

Evolution of the tan locus contributed to pigment loss in Drosophila santomea: a response to Matute et al.

We have shown previously that the loss of abdominal pigmentation in D. santomea relative to its sister species D. yakuba resulted, in part, from cis-regulatory mutations at the tan locus. Matute et al. claim, based solely upon extrapolation from genetic crosses of D. santomea and D. melanogaster, a much more divergent species, that at least four X chromosome regions but not tan are responsible for pigmentation differences. Here, we provide additional evidence from introgressions of D. yakuba genes into D. santomea that support a causative role for tan in the loss of pigmentation and present analyses that contradict Matute et al.'s claims. We discuss how the choice of parental species and other factors affect the ability to identify loci responsible for species divergence, and we affirm that all of our previously reported results and conclusions stand.

Epistatic Effects Between Amino Acid Insertions and Substitutions Mediate Toxin resistance of Vertebrate Na+,K+-ATPases.

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino acid substitutions in the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work revealed that hystricognath rodents (e.g., chinchilla) and pterocliform birds (sandgrouse) have convergently evolved amino acid insertions in the H1-H2 loop, but their functional significance was not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance in homologs of each of the two species that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla NKA unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse NKA, the amino acid insertion and substitution Q111R both contribute to an augmented CTS resistance without compromising ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific mutational effects on CTS insensitivity of the enzyme. Our results highlight the diversity of genetic substrates that underlie CTS insensitivity in vertebrate NKA and reveal how amino acid insertions can alter the phenotypic effects of point mutations at key sites in the same protein domain.

The evolution of gene regulation underlies a morphological difference between two Drosophila sister species.

Understanding the mechanisms underlying the morphological divergence of species is one of the central goals of evolutionary biology. Here, we analyze the genetic and molecular bases of the divergence of body pigmentation patterns between Drosophila yakuba and its sister species Drosophila santomea. We found that loss of pigmentation in D. santomea involved the selective loss of expression of the tan and yellow pigmentation genes. We demonstrate that tan gene expression was eliminated through the mutational inactivation of one specific tan cis-regulatory element (CRE) whereas the Tan protein sequence remained unchanged. Surprisingly, we identify three independent loss-of-function alleles of the tan CRE in the young D. santomea lineage. We submit that there is sufficient empirical evidence to support the general prediction that functional evolutionary changes at pleiotropic loci will most often involve mutations in their discrete, modular cis-regulatory elements.

Allele-specific knockouts reveal a role for apontic-like in the evolutionary loss of larval melanin pigmentation in the domesticated silkworm, Bombyx mori.

The domesticated silkworm, Bombyx mori, and its wild progenitor, B. mandarina, are extensively studied as a model case of the evolutionary process of domestication. A conspicuous difference between these species is the dramatic reduction in melanin pigmentation in both larval and adult B. mori. Here we evaluate the efficiency of CRISPR/Cas9-targeted knockouts of pigment-related genes as a tool to understand their potential contributions to domestication-associated melanin pigmentation loss in B. mori. To demonstrate the efficacy of targeted knockouts in B. mandarina, we generated a homozygous CRISPR/Cas9-targeted knockout of yellow-y. In yellow-y knockout mutants, black body colour became lighter throughout the larval, pupal and adult stages, confirming a role for this gene in melanin pigment formation. Further, we performed allele-specific CRISPR/Cas9-targeted knockouts of the pigment-related transcription factor, apontic-like (apt-like) in B. mori × B. mandarina F1 hybrid individuals which exhibit B. mandarina-like larval pigmentation. Knockout of the B. mandarina allele of apt-like in F1 embryos results in white patches on the dorsal integument of larvae, whereas corresponding knockouts of the B. mori allele consistently exhibit normal F1 larval pigmentation. These results demonstrate a contribution of apt-like to the evolution of reduced melanin pigmentation in B. mori. Together, our results demonstrate the feasibility of CRISPR/Cas9-targeted knockouts as a tool for understanding the genetic basis of traits associated with B. mori domestication.

Genomic signatures of spatially divergent selection at clownfish range margins.

Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.

Assortative mating and persistent reproductive isolation in hybrids.

The emergence of new species is driven by the establishment of mechanisms that limit gene flow between populations. A major challenge is reconciling the theoretical and empirical importance of assortative mating in speciation with the ease with which it can fail. Swordtail fish have an evolutionary history of hybridization and fragile prezygotic isolating mechanisms. Hybridization between two swordtail species likely arose via pollution-mediated breakdown of assortative mating in the 1990s. Here we track unusual genetic patterns in one hybrid population over the past decade using whole-genome sequencing. Hybrids in this population formed separate genetic clusters by 2003, and maintained near-perfect isolation over 25 generations through strong ancestry-assortative mating. However, we also find that assortative mating was plastic, varying in strength over time and disappearing under manipulated conditions. In addition, a nearby population did not show evidence of assortative mating. Thus, our findings suggest that assortative mating may constitute an intermittent and unpredictable barrier to gene flow, but that variation in its strength can have a major effect on how hybrid populations evolve. Understanding how reproductive isolation varies across populations and through time is critical to understanding speciation and hybridization, as well as their dependence on disturbance.

A Genomic Map of the Effects of Linked Selection in Drosophila.

Natural selection at one site shapes patterns of genetic variation at linked sites. Quantifying the effects of "linked selection" on levels of genetic diversity is key to making reliable inference about demography, building a null model in scans for targets of adaptation, and learning about the dynamics of natural selection. Here, we introduce the first method that jointly infers parameters of distinct modes of linked selection, notably background selection and selective sweeps, from genome-wide diversity data, functional annotations and genetic maps. The central idea is to calculate the probability that a neutral site is polymorphic given local annotations, substitution patterns, and recombination rates. Information is then combined across sites and samples using composite likelihood in order to estimate genome-wide parameters of distinct modes of selection. In addition to parameter estimation, this approach yields a map of the expected neutral diversity levels along the genome. To illustrate the utility of our approach, we apply it to genome-wide resequencing data from 125 lines in Drosophila melanogaster and reliably predict diversity levels at the 1Mb scale. Our results corroborate estimates of a high fraction of beneficial substitutions in proteins and untranslated regions (UTR). They allow us to distinguish between the contribution of sweeps and other modes of selection around amino acid substitutions and to uncover evidence for pervasive sweeps in untranslated regions (UTRs). Our inference further suggests a substantial effect of other modes of linked selection and of adaptation in particular. More generally, we demonstrate that linked selection has had a larger effect in reducing diversity levels and increasing their variance in D. melanogaster than previously appreciated.

Sexual Dimorphism and Retinal Mosaic Diversification following the Evolution of a Violet Receptor in Butterflies.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Reproductive isolation of hybrid populations driven by genetic incompatibilities.

Despite its role in homogenizing populations, hybridization has also been proposed as a means to generate new species. The conceptual basis for this idea is that hybridization can result in novel phenotypes through recombination between the parental genomes, allowing a hybrid population to occupy ecological niches unavailable to parental species. Here we present an alternative model of the evolution of reproductive isolation in hybrid populations that occurs as a simple consequence of selection against genetic incompatibilities. Unlike previous models of hybrid speciation, our model does not incorporate inbreeding, or assume that hybrids have an ecological or reproductive fitness advantage relative to parental populations. We show that reproductive isolation between hybrids and parental species can evolve frequently and rapidly under this model, even in the presence of substantial ongoing immigration from parental species and strong selection against hybrids. An interesting prediction of our model is that replicate hybrid populations formed from the same pair of parental species can evolve reproductive isolation from each other. This non-adaptive process can therefore generate patterns of species diversity and relatedness that resemble an adaptive radiation. Intriguingly, several known hybrid species exhibit patterns of reproductive isolation consistent with the predictions of our model.

The within-host dynamics of infection in trans-generationally primed flour beetles.

Many taxa exhibit plastic immune responses initiated after primary microbial exposure that provide increased protection against disease-induced mortality and the fitness costs of infection. In several arthropod species, this protection can even be passed from parents to offspring through a phenomenon called trans-generational immune priming. Here, we first demonstrate that trans-generational priming is a repeatable phenomenon in flour beetles (Tribolium castaneum) primed and infected with Bacillus thuringiensis (Bt). We then quantify the within-host dynamics of microbes and host physiological responses in infected offspring from primed and unprimed mothers by monitoring bacterial density and using mRNA-seq to profile host gene expression, respectively, over the acute infection period. We find that priming increases inducible resistance against Bt around a critical temporal juncture where host septicaemic trajectories, and consequently survival, may be determined in unprimed individuals. Our results identify a highly differentially expressed biomarker of priming, containing an EIF4-e domain, in uninfected individuals, as well as several other candidate genes. Moreover, the induction and decay dynamics of gene expression over time suggest a metabolic shift in primed individuals. The identified bacterial and gene expression dynamics are likely to influence patterns of bacterial fitness and disease transmission in natural populations.

A second-generation assembly of the Drosophila simulans genome provides new insights into patterns of lineage-specific divergence.

We create a new assembly of the Drosophila simulans genome using 142 million paired short-read sequences and previously published data for strain w(501). Our assembly represents a higher-quality genomic sequence with greater coverage, fewer misassemblies, and, by several indexes, fewer sequence errors. Evolutionary analysis of this genome reference sequence reveals interesting patterns of lineage-specific divergence that are different from those previously reported. Specifically, we find that Drosophila melanogaster evolves faster than D. simulans at all annotated classes of sites, including putatively neutrally evolving sites found in minimal introns. While this may be partly explained by a higher mutation rate in D. melanogaster, we also find significant heterogeneity in rates of evolution across classes of sites, consistent with historical differences in the effective population size for the two species. Also contrary to previous findings, we find that the X chromosome is evolving significantly faster than autosomes for nonsynonymous and most noncoding DNA sites and significantly slower for synonymous sites. The absence of a X/A difference for putatively neutral sites and the robustness of the pattern to Gene Ontology and sex-biased expression suggest that partly recessive beneficial mutations may comprise a substantial fraction of noncoding DNA divergence observed between species. Our results have more general implications for the interpretation of evolutionary analyses of genomes of different quality.

Population differences in olfaction accompany host shift in Drosophila mojavensis.

Evolutionary shifts in plant-herbivore interactions provide a model for understanding the link among the evolution of behaviour, ecological specialization and incipient speciation. Drosophila mojavensis uses different host cacti across its range, and volatile chemicals emitted by the host are the primary cue for host plant identification. In this study, we show that changes in host plant use between distinct D. mojavensis populations are accompanied by changes in the olfactory system. Specifically, we observe differences in olfactory receptor neuron specificity and sensitivity, as well as changes in sensillar subtype abundance, between populations. Additionally, RNA-seq analyses reveal differential gene expression between populations for members of the odorant receptor gene family. Hence, alterations in host preference are associated with changes in development, regulation and function at the olfactory periphery.

Ancient hybridization and genomic stabilization in a swordtail fish.

A rapidly increasing body of work is revealing that the genomes of distinct species often exhibit hybrid ancestry, presumably due to postspeciation hybridization between closely related species. Despite the growing number of documented cases, we still know relatively little about how genomes evolve and stabilize following hybridization, and to what extent hybridization is functionally relevant. Here, we examine the case of Xiphophorus nezahualcoyotl, a teleost fish whose genome exhibits significant hybrid ancestry. We show that hybridization was relatively ancient and is unlikely to be ongoing. Strikingly, the genome of X. nezahualcoyotl has largely stabilized following hybridization, distinguishing it from examples such as human-Neanderthal hybridization. Hybridization-derived regions are remarkably distinct from other regions of the genome, tending to be enriched in genomic regions with reduced constraint. These results suggest that selection has played a role in removing hybrid ancestry from certain functionally important regions. Combined with findings in other systems, our results raise many questions about the process of genomic stabilization and the role of selection in shaping patterns of hybrid ancestry in the genome.

Landscape of standing variation for tandem duplications in Drosophila yakuba and Drosophila simulans.

We have used whole genome paired-end Illumina sequence data to identify tandem duplications in 20 isofemale lines of Drosophila yakuba and 20 isofemale lines of D. simulans and performed genome wide validation with PacBio long molecule sequencing. We identify 1,415 tandem duplications that are segregating in D. yakuba as well as 975 duplications in D. simulans, indicating greater variation in D. yakuba. Additionally, we observe high rates of secondary deletions at duplicated sites, with 8% of duplicated sites in D. simulans and 17% of sites in D. yakuba modified with deletions. These secondary deletions are consistent with the action of the large loop mismatch repair system acting to remove polymorphic tandem duplication, resulting in rapid dynamics of gain and loss in duplicated alleles and a richer substrate of genetic novelty than has been previously reported. Most duplications are present in only single strains, suggesting that deleterious impacts are common. Drosophila simulans shows larger numbers of whole gene duplications in comparison to larger proportions of gene fragments in D. yakuba. Drosophila simulans displays an excess of high-frequency variants on the X chromosome, consistent with adaptive evolution through duplications on the D. simulans X or demographic forces driving duplicates to high frequency. We identify 78 chimeric genes in D. yakuba and 38 chimeric genes in D. simulans, as well as 143 cases of recruited noncoding sequence in D. yakuba and 96 in D. simulans, in agreement with rates of chimeric gene origination in D. melanogaster. Together, these results suggest that tandem duplications often result in complex variation beyond whole gene duplications that offers a rich substrate of standing variation that is likely to contribute both to detrimental phenotypes and disease, as well as to adaptive evolutionary change.

Genome sequencing reveals complex speciation in the Drosophila simulans clade.

The three species of the Drosophila simulans clade--the cosmopolitan species, D. simulans, and the two island endemic species, D. mauritiana and D. sechellia--are important models in speciation genetics, but some details of their phylogenetic and speciation history remain unresolved. The order and timing of speciation are disputed, and the existence, magnitude, and timing of gene flow among the three species remain unclear. Here we report on the analysis of a whole-genome four-species sequence alignment that includes all three D. simulans clade species as well as the D. melanogaster reference sequence. The alignment comprises novel, paired short-read sequence data from a single highly inbred line each from D. simulans, D. mauritiana, and D. sechellia. We are unable to reject a species phylogeny with a basal polytomy; the estimated age of the polytomy is 242,000 yr before the present. However, we also find that up to 4.6% of autosomal and 2.2% of X-linked regions have evolutionary histories consistent with recent gene flow between the mainland species (D. simulans) and the two island endemic species (D. mauritiana and D. sechellia). Our findings thus show that gene flow has occurred throughout the genomes of the D. simulans clade species despite considerable geographic, ecological, and intrinsic reproductive isolation. Last, our analysis of lineage-specific changes confirms that the D. sechellia genome has experienced a significant excess of slightly deleterious changes and a dearth of presumed favorable changes. The relatively reduced efficacy of natural selection in D. sechellia is consistent with its derived, persistently reduced historical effective population size.

Revisiting an old riddle: what determines genetic diversity levels within species?

Understanding why some species have more genetic diversity than others is central to the study of ecology and evolution, and carries potentially important implications for conservation biology. Yet not only does this question remain unresolved, it has largely fallen into disregard. With the rapid decrease in sequencing costs, we argue that it is time to revive it.

Multiplexed shotgun genotyping for rapid and efficient genetic mapping.

We present a new approach to genotyping based on multiplexed shotgun sequencing that can identify recombination breakpoints in a large number of individuals simultaneously at a resolution sufficient for most mapping purposes, such as quantitative trait locus (QTL) mapping and mapping of induced mutations. We first describe a simple library construction protocol that uses just 10 ng of genomic DNA per individual and makes the approach accessible to any laboratory with standard molecular biology equipment. Sequencing this library results in a large number of sequence reads widely distributed across the genomes of multiplexed bar-coded individuals. We develop a Hidden Markov Model to estimate ancestry at all genomic locations in all individuals using these data. We demonstrate the utility of the approach by mapping a dominant marker allele in D. simulans to within 105 kb of its true position using 96 F1-backcross individuals genotyped in a single lane on an Illumina Genome Analyzer. We further demonstrate the utility of our method by genetically mapping more than 400 previously unassembled D. simulans contigs to linkage groups and by evaluating the quality of targeted introgression lines. At this level of multiplexing and divergence between strains, our method allows estimation of recombination breakpoints to a median of 38-kb intervals. Our analysis suggests that higher levels of multiplexing and/or use of strains with lower levels of divergence are practicable.

Correlated evolution of nearby residues in Drosophilid proteins.

Here we investigate the correlations between coding sequence substitutions as a function of their separation along the protein sequence. We consider both substitutions between the reference genomes of several Drosophilids as well as polymorphisms in a population sample of Zimbabwean Drosophila melanogaster. We find that amino acid substitutions are "clustered" along the protein sequence, that is, the frequency of additional substitutions is strongly enhanced within ≈10 residues of a first such substitution. No such clustering is observed for synonymous substitutions, supporting a "correlation length" associated with selection on proteins as the causative mechanism. Clustering is stronger between substitutions that arose in the same lineage than it is between substitutions that arose in different lineages. We consider several possible origins of clustering, concluding that epistasis (interactions between amino acids within a protein that affect function) and positional heterogeneity in the strength of purifying selection are primarily responsible. The role of epistasis is directly supported by the tendency of nearby substitutions that arose on the same lineage to preserve the total charge of the residues within the correlation length and by the preferential cosegregation of neighboring derived alleles in our population sample. We interpret the observed length scale of clustering as a statistical reflection of the functional locality (or modularity) of proteins: amino acids that are near each other on the protein backbone are more likely to contribute to, and collaborate toward, a common subfunction.