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.

L1 retrotransposons are an abundant class of transposable elements which pose a threat to genome stability and may play a role in age-related pathologies such as cancer. Recent evidence indicates that L1s become more active in somatic tissues during the course of aging; the mechanisms underlying this phenomenon remain unknown, however. Here we report that the longevity regulating protein, SIRT6, is a powerful repressor of L1 activity. Specifically, SIRT6 binds to the 5′UTR of L1 loci, where it mono-ADP ribosylates the nuclear corepressor protein, KAP1, and facilitates KAP1 interaction with the heterochromatin factor, HP1α, thereby contributing to the packaging of L1 elements into transcriptionally repressive heterochromatin. During the course of aging, and also in response to DNA damage, however, we find that SIRT6 is depleted from L1 loci, allowing for the activation of these previously silenced retroelements.

Bayesian inference using Markov chain Monte Carlo (MCMC) has become one of the primary methods used to infer phylogenies from sequence data. Assessing convergence is a crucial component of these analyses, as it establishes the reliability of the posterior distribution estimates of the tree topology and model parameters sampled from the MCMC. Numerous tests and visualizations have been developed for this purpose, but many of the most popular methods are implemented in ways that make them inconvenient to use for large data sets. RWTY is an R package that implements established and new methods for diagnosing phylogenetic MCMC convergence in a single convenient interface.

Hurricanes are catastrophically destructive. Beyond their toll on human life and livelihoods, hurricanes have tremendous and often long-lasting effects on ecological systems. Despite many examples of mass mortality events following hurricanes, hurricane-induced natural selection has not previously been demonstrated. Immediately after we finished a survey of Anolis scriptus-a common, small-bodied lizard found throughout the Turks and Caicos archipelago-our study populations were battered by Hurricanes Irma and Maria. Shortly thereafter, we revisited the populations to determine whether morphological traits related to clinging capacity had shifted in the intervening six weeks and found that populations of surviving lizards differed in body size, relative limb length and toepad size from those present before the storm. Our serendipitous study, which to our knowledge is the first to use an immediately before and after comparison to investigate selection caused by hurricanes, demonstrates that hurricanes can induce phenotypic change in a population and strongly implicates natural selection as the cause. In the decades ahead, as extreme climate events are predicted to become more intense and prevalent, our understanding of evolutionary dynamics needs to incorporate the effects of these potentially severe selective episodes.

To explain the frequency and distribution of heteromorphic sex chromosomes in the lizard genus Anolis we compared the relative roles of sex chromosome conservation vs. turnover of sex determining mechanisms. We used model based comparative methods to reconstruct karyotype evolution and the presence of heteromorphic sex chromosomes onto a newly generated Anolis phylogeny. We found that heteromorphic sex chromosomes evolved multiple times in the genus. Fluorescent in situ hybridization (FISH) of repetitive DNA showed variable rates of Y chromosome degeneration among Anolis species and identified previously undetected, homomorphic sex chromosomes in two species. We confirmed homology of sex chromosomes in the genus by performing FISH of an X-linked BAC and qPCR of X-linked genes in multiple Anolis species sampled across the phylogeny. Taken together, these results are consistent with long-term conservation of sex chromosomes in the group. Our results pave the way to address additional questions related to Anolis sex chromosome evolution and describe a conceptual framework that can be used to evaluate the origins and evolution of heteromorphic sex chromosomes in other clades.

Secondary contact between divergent populations or incipient species may result in the exchange and introgression of genomic material. We develop a simple DNA sequence measure, called G min, which is designed to identify genomic regions experiencing introgression in a secondary contact model. G min is defined as the ratio of the minimum between-population number of nucleotide differences in a genomic window to the average number of between-population differences. Although it is conceptually simple, one advantage of G min is that it is computationally inexpensive relative to model-based methods for detecting gene flow and it scales easily to the level of whole-genome analysis. We compare the sensitivity and specificity of G min to those of the widely used index of population differentiation, F ST, and suggest a simple statistical test for identifying genomic outliers. Extensive computer simulations demonstrate that G min has both greater sensitivity and specificity for detecting recent introgression than does F ST. Furthermore, we find that the sensitivity of G min is robust with respect to both the population mutation and recombination rates. Finally, a scan of G min across the X chromosome of Drosophila melanogaster identifies candidate regions of introgression between sub-Saharan African and cosmopolitan populations that were previously missed by other methods. These results show that G min is a biologically straightforward, yet powerful, alternative to F ST, as well as to more computationally intensive model-based methods for detecting gene flow.

Drosophila mauritiana is an Indian Ocean island endemic species that diverged from its two sister species, Drosophila simulans and Drosophila sechellia, approximately 240,000 years ago. Multiple forms of incomplete reproductive isolation have evolved among these species, including sexual, gametic, ecological, and intrinsic postzygotic barriers, with crosses among all three species conforming to Haldane’s rule: F1 hybrid males are sterile and F1 hybrid females are fertile. Extensive genetic resources and the fertility of hybrid females have made D. mauritiana, in particular, an important model for speciation genetics. Analyses between D. mauritiana and both of its siblings have shown that the X chromosome makes a disproportionate contribution to hybrid male sterility. But why the X plays a special role in the evolution of hybrid sterility in these, and other, species remains an unsolved problem. To complement functional genetic analyses, we have investigated the population genomics of D. mauritiana, giving special attention to differences between the X and the autosomes. We present a de novo genome assembly of D. mauritiana annotated with RNAseq data and a whole-genome analysis of polymorphism and divergence from ten individuals. Our analyses show that, relative to the autosomes, the X chromosome has reduced nucleotide diversity but elevated nucleotide divergence; an excess of recurrent adaptive evolution at its protein-coding genes; an excess of recent, strong selective sweeps; and a large excess of satellite DNA. Interestingly, one of two centimorgan-scale selective sweeps on the D. mauritiana X chromosome spans a region containing two sex-ratio meiotic drive elements and a high concentration of satellite DNA. Furthermore, genes with roles in reproduction and chromosome biology are enriched among genes that have histories of recurrent adaptive protein evolution. Together, these genome-wide analyses suggest that genetic conflict and frequent positive natural selection on the X chromosome have shaped the molecular evolutionary history of D. mauritiana, refining our understanding of the possible causes of the large X-effect in speciation.

Hybridization is among the evolutionary mechanisms most frequently hypothesized to drive the success of invasive species, in part because hybrids are common in invasive populations. One explanation for this pattern is that biological invasions coincide with a change in selection pressures that limit hybridization in the native range. To investigate this possibility, we studied the introduction of the brown anole (Anolis sagrei) in the southeastern United States. We find that native populations are highly genetically structured. In contrast, all invasive populations show evidence of hybridization among native-range lineages. Temporal sampling in the invasive range spanning 15 y showed that invasive genetic structure has stabilized, indicating that large-scale contemporary gene flow is limited among invasive populations and that hybrid ancestry is maintained. Additionally, our results are consistent with hybrid persistence in invasive populations resulting from changes in natural selection that occurred during invasion. Specifically, we identify a large-effect X chromosome locus associated with variation in limb length, a well-known adaptive trait in anoles, and show that this locus is often under selection in the native range, but rarely so in the invasive range. Moreover, we find that the effect size of alleles at this locus on limb length is much reduced in hybrids among divergent lineages, consistent with epistatic interactions. Thus, in the native range, epistasis manifested in hybrids can strengthen extrinsic postmating isolation. Together, our findings show how a change in natural selection can contribute to an increase in hybridization in invasive populations.