• The Interactive Fly

    Mechanisms of Evolution

    Extensive local adaptation within the chemosensory system following Drosophila melanogaster's global expansion

    How organisms adapt to new environments is of fundamental biological interest, but poorly understood at the genetic level. Chemosensory systems provide attractive models to address this problem, because they lie between external environmental signals and internal physiological responses. To investigate how selection has shaped the well-characterized chemosensory system of Drosophila melanogaster, this study analysed genome-wide data from five diverse populations. By couching population genomic analyses of chemosensory protein families, including odorant receptors, gustatory receptors, and odorant-binding proteins, within parallel analyses of other large families, it was demonstrated that chemosensory proteins are not outliers for adaptive divergence between species. However, chemosensory families often display the strongest genome-wide signals of recent selection within D. melanogaster. Recent adaptation has operated almost exclusively on standing variation, and patterns of adaptive mutations predict diverse effects on protein function. Finally, evidence is provided that chemosensory proteins have experienced relaxed constraint, and it is argued that this has been important for their rapid adaptation over short timescales (Arguello, 2016).

    Gene expression clines reveal local adaptation and associated trade-offs at a continental scale

    Local adaptation, where fitness in one environment comes at a cost in another, should lead to spatial variation in trade-offs between life history traits and may be critical for population persistence. Recent studies have sought genomic signals of local adaptation, but often have been limited to laboratory populations representing two environmentally different locations of a species' distribution. This study measured gene expression, as a proxy for fitness, in males of Drosophila subobscura, occupying a 20 degrees latitudinal and 11 ° C thermal range. Uniquely, six populations were sampled, and both common garden and semi-natural responses to identify signals of local adaptation were identified. Contrasting patterns of investment were found: transcripts with expression positively correlated to latitude were enriched for metabolic processes, expressed across all tissues whereas negatively correlated transcripts were enriched for reproductive processes, expressed primarily in testes. When using only the end populations, to compare the results to previous studies, it was found that locally adaptive patterns were obscured. While phenotypic trade-offs between metabolic and reproductive functions across widespread species are well-known, the results identify underlying genetic and tissue responses at a continental scale that may be responsible for this. This may contribute to understanding population persistence under environmental change (Porcelli, 2016).

    Experimental evolution of gene expression and plasticity in alternative selective regimes

    Little is known of how gene expression and its plasticity evolves as populations adapt to different environmental regimes. Expression is expected to evolve adaptively in all populations but only those populations experiencing environmental heterogeneity are expected to show adaptive evolution of plasticity. This study measured the transcriptome in a cadmium-enriched diet and a salt-enriched diet for experimental populations of Drosophila melanogaster that evolved for ~130 generations in one of four selective regimes: two constant regimes maintained in either cadmium or salt diets and two heterogeneous regimes that varied either temporally or spatially between the two diets. For populations evolving in constant regimes, a strong signature of counter-gradient evolution was found; the evolved expression differences between populations adapted to alternative diets is opposite to the plastic response of the ancestral population that is naive to both diets. Based on expression patterns in the ancestral populations, a set of genes was identified for which selection in heterogeneous regimes was predicted to result in increases in plasticity, and the expected pattern was found. In contrast, a set of genes where reduced plasticity was predicted did not follow expectation. Nonetheless, both gene sets showed a pattern consistent with adaptive expression evolution in heterogeneous regimes, highlighting the difference between observing 'optimal' plasticity and improvements in environment-specific expression. Looking across all genes, there is evidence in all regimes of differences in biased allele expression across environments ('allelic plasticity') and this is more common among genes with plasticity in total expression (Huang, 2016).

    Reasons for success: rapid evolution for desiccation resistance and life-history changes in the polyphagous fly Anastrepha ludens

    Species that exhibit broad ranges of distribution may successfully navigate environmental changes by modifying some of their life history traits. Environmental humidity imposes a critical stress that organisms may overcome by increasing their resistance to desiccation. This study used experimental evolution to investigate adaptation to desiccation in the tephritid Anastrepha ludens, a species with high fecundity, late maturation and long lifespan. This study measured morphological, physiological, developmental as well as demographic changes involved in the adaptation to desiccation. Notwithstanding a low heritability (h2 = 0.237), desiccation resistance evolved extremely rapidly and few negative trade-offs were detected. Selected flies exhibited correlated increases in longevity, body size, the amount of body lipids and bulk water content, and in the duration of the pupal stage. Females further delayed sexual maturation, decreased daily fecundity but retained high lifetime reproductive potential. No differences in male mating competitiveness were found. Selected and control lines differed in longevity but not in total female fecundity, demonstrating that A. ludens flies have the capability for fast adaptation to desiccation without loosing their reproductive capability. Thus, it seems that a rapid evolutionary response to desiccation in this polyphagous insect works as a buffer for environmental variation and reduces the strength of selection on reproductive traits (Tejeda, 2016).

    Cold adaptation increases rates of nutrient flow and metabolic plasticity during cold exposure in Drosophila melanogaster

    Metabolic flexibility is an important component of adaptation to stressful environments, including thermal stress and latitudinal adaptation. The direct relationship between selection on thermal stress hardiness and metabolic flux has not previously been tested. This study investigated flexibility of nutrient catabolism during cold stress in Drosophila artificially selected for fast or slow recovery from chill coma (i.e. cold-hardy or -susceptible), specifically testing the hypothesis that stress adaptation increases metabolic turnover. Using 13C-labelled glucose, this study first showed that cold-hardy flies more rapidly incorporate ingested carbon into amino acids and newly synthesized glucose, permitting rapid synthesis of proline, a compound shown elsewhere to improve survival of cold stress. Second, using glucose and leucine tracers cold-hardy flies were shown to have higher oxidation rates than cold-susceptible flies before cold exposure, similar oxidation rates during cold exposure, and returned to higher oxidation rates during recovery. Additionally, cold-hardy flies transferred compounds among body pools more rapidly during cold exposure and recovery. Increased metabolic turnover may allow cold-adapted flies to better prepare for, resist and repair/tolerate cold damage. This work illustrates for the first time differences in nutrient fluxes associated with cold adaptation, suggesting that metabolic costs associated with cold hardiness could invoke resource-based trade-offs that shape life histories (Williams, 2016).

    Colder environments did not select for a faster metabolism during experimental evolution of Drosophila melanogaster

    The effect of temperature on the evolution of metabolism has been the subject of debate for a century; however, no consistent patterns have emerged from comparisons of metabolic rate within and among species living at different temperatures. This study used experimental evolution to determine how metabolism evolves in populations of Drosophila melanogaster exposed to one of three selective treatments: a constant 16 ° C, a constant 25 ° C, or temporal fluctuations between 16 and 25 ° C. August Krogh's controversial hypothesis was tested that colder environments select for a faster metabolism. Given that colder environments also experience greater seasonality, the hypothesis was also tested that temporal variation in temperature may be the factor that selects for a faster metabolism. The metabolic rate of flies from each selective treatment was measured at 16, 20.5, and 25 ° C. Although metabolism was faster at higher temperatures, flies from the selective treatments had similar metabolic rates at each measurement temperature. Based on variation among genotypes within populations, heritable variation in metabolism was likely sufficient for adaptation to occur. It is concluded that colder or seasonal environments do not necessarily select for a faster metabolism. Rather, other factors besides temperature likely contribute to patterns of metabolic rate over thermal clines in nature (Alton, 2016).

    Adaptive patterns of phenotypic plasticity in laboratory and field environments in Drosophila melanogaster

    Identifying mechanisms of adaptation to variable environments is essential in developing a comprehensive understanding of evolutionary dynamics in natural populations. Phenotypic plasticity allows for phenotypic change in response to changes in the environment, and as such may play a major role in adaptation to environmental heterogeneity. This study examined the plasticity of stress response in D. melanogaster originating from two distinct geographic regions and ecological habitats. Adults were given a short-term, 5-day exposure to combinations of temperature and photoperiod to elicit a plastic response for three fundamental aspects of stress tolerance that vary adaptively with geography. This was replicated in both the laboratory and in outdoor enclosures in the field. In the laboratory, geographic origin was found to be the primary determinant of the stress response. Temperature and the interaction between temperature and photoperiod were also found to significantly affect stress resistance. In the outdoor enclosures, plasticity was distinct among traits and between geographic regions. These results demonstrate that short-term exposure of adults to ecologically relevant environmental cues results in predictable effects on multiple aspects of fitness. These patterns of plasticity vary among traits and are highly distinct between the two examined geographic regions, consistent with patterns of local adaptation to climate and associated environmental parameters (Mathur, 2016).

    Genomic trajectories to desiccation resistance: Convergence and divergence among replicate selected Drosophila lines

    Adaptation to environmental stress is critical for long-term species persistence. With climate change and other anthropogenic stressors compounding natural selective pressures, understanding the nature of adaptation is as important as ever in evolutionary biology. This study investigated this issue in a set of replicated Drosophila lines selected for increased desiccation resistance, a classical physiological trait that has been closely linked to Drosophila species distributions. Pooled whole-genome sequencing was used to compare the genetic basis of their selection responses. While selected SNPs in replicates of the same treatment (desiccation-selection or lab adaptation) tended to change frequency in the same direction, suggesting some commonality in the selection response, candidate SNP and gene lists often differed among replicates. Three of the five desiccation-selection replicates showed significant overlap at the gene and network level. All five replicates showed enrichment for ovary-expressed genes, suggesting maternal effects on the selected trait. Divergence between pairs of replicate lines for desiccation-candidate SNPs was greater than between pairs of control lines. This difference also far exceeded the divergence between pairs of replicate lines for neutral SNPs. Overall, while there was overlap in the direction of allele frequency changes and the network and functional categories affected by desiccation selection, replicates showed unique responses at all levels likely reflecting hitchhiking effects, and highlighting the challenges in identifying candidate genes from these types of experiments when traits are likely to be polygenic (Griffin, 2016).

    Independent natural genetic variation of punishment- versus relief-memory

    A painful event establishes two opponent memories: cues that are associated with pain onset are remembered negatively, whereas cues that coincide with the relief at pain offset acquire positive valence. Such punishment- versus relief-memories are conserved across species, including humans, and the balance between them is critical for adaptive behaviour with respect to pain and trauma. In the fruit fly, Drosophila melanogaster as a study case, this study found that both punishment- and relief-memories display natural variation across wild-derived inbred strains, but they do not covary, suggesting a considerable level of dissociation in their genetic effectors. This provokes the question whether there may be heritable inter-individual differences in the balance between these opponent memories in man, with potential psycho-clinical implications (Appel, 2016).

    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. This 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 understanding of the genetic basis and evolution of virgin egg retention, a genome-wide association study (GWAS) was performed using the previously sequenced Drosophila Genetic Reference Panel (DGRP) mapping population. Twenty-nine single nucleotide polymorphisms (SNPs) associated with virgin egg retention were found, and six available mutant lines, each harboring a mutation in a candidate gene, were examined for effects on egg retention time. 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. 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. These 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 (Akhund-Zade, 2016).

    Testing for local adaptation and evolutionary potential along altitudinal gradients in rainforest Drosophila: beyond laboratory estimates

    Predicting how species will respond to the rapid climatic changes predicted this century is an urgent task. Species distribution models (SDMs) use the current relationship between environmental variation and species' abundances to predict the effect of future environmental change on their distributions. However, two common assumptions of SDMs are likely to be violated in many cases: (i) that the relationship of environment with abundance or fitness is constant throughout a species' range and will remain so in future and (ii) that abiotic factors (e.g. temperature, humidity) determine species' distributions. These assumptions were tested by relating field abundance of the rainforest fruit fly Drosophila birchii to ecological change across gradients that include its low and high altitudinal limits. Then, how such ecological variation affects the fitness of 35 D. birchii families transplanted in 591 cages to sites along two altitudinal gradients, was tested to determine whether genetic variation in fitness responses could facilitate future adaptation to environmental change. Overall, field abundance was highest at cooler, high-altitude sites, and declined towards warmer, low-altitude sites. By contrast, cage fitness (productivity) increased towards warmer, lower-altitude sites, suggesting that biotic interactions (absent from cages) drive ecological limits at warmer margins. In addition, the relationship between environmental variation and abundance varied significantly among gradients, indicating divergence in ecological niche across the species' range. However, there was no evidence for local adaptation within gradients, despite greater productivity of high-altitude than low-altitude populations when families were reared under laboratory conditions. Families also responded similarly to transplantation along gradients, providing no evidence for fitness trade-offs that would favour local adaptation. These findings highlight the importance of (i) measuring genetic variation in key traits under ecologically relevant conditions, and (ii) considering the effect of biotic interactions when predicting species' responses to environmental change (O'Brien, 2017)

    Genomics of parallel experimental evolution in Drosophila

    What are the genomic foundations of adaptation in sexual populations? This question was addressed using fitness-character and whole-genome sequence data from 30 Drosophila laboratory populations. These 30 populations are part of a nearly forty-year laboratory radiation featuring three selection regimes, each shared by ten populations for up to 837 generations, with moderately large effective population sizes. Each of three sets of ten populations that shared a selection regime consist of five populations that have long been maintained under that selection regime, paired with five populations that had only recently been subjected to that selection regime. A high degree of evolutionary parallelism in fitness phenotypes was found when most-recent selection regimes are shared, as in previous studies from this laboratory. Genomic parallelism was also found with respect to the frequencies of single-nucleotide polymorphisms, transposable elements, insertions, and structural variants, which was expected. Entirely unexpected was a high degree of parallelism for linkage disequilibrium. The evolutionary genetic changes among these sexual populations are rapid and genomically extensive. This pattern may be due to segregating functional genetic variation that is abundantly maintained genome-wide by selection, variation that responds immediately to changes of selection regime (Graves, 2017).

    Geographical analysis of diapause inducibility in European Drosophila melanogaster populations

    Seasonal overwintering in insects represents an adaptation to stressful environments and in European Drosophila melanogaster females, low temperatures and short photoperiods can induce an ovarian diapause. Diapause may represent a recent (<15Ky) adaptation to the colonisation of temperate Europe by D. melanogaster from tropical sub-Saharan Africa, because African D. melanogaster and the sibling species D. simulans, have been reported to fail to undergo diapause. Over the past few centuries, D. melanogaster have also invaded North America and Australia, and eastern populations on both continents show a predictable latitudinal cline in diapause induction. In Europe however, a new diapause-enhancing timeless allele, ls-tim, is observed at high levels in southern Italy ( approximately 80%), where it appears to have arisen and has spread throughout the continent with a frequency of approximately 20% in Scandinavia. Given the phenotype of ls-tim and its geographical distribution, it was predicted that it would work against any latitudinal cline in diapause induction within Europe. Indeed this study revealed that any latitudinal cline for diapause in Europe is very weak, as predicted by ls-tim frequencies. In contrast, ls-tim frequencies were determined in North America and it was observed that they would be expected to strengthen the latitudinal pattern of diapause. The results reveal how a newly arisen mutation, can, via the stochastic nature of where it initially arose, blur an otherwise adaptive geographical pattern (Pegoraro, 2017)

    Evolution of circadian rhythms in Drosophila melanogaster populations reared in constant light and dark regimes for over 330 generations

    Organisms are believed to have evolved circadian clocks as adaptations to deal with cyclic environmental changes, and therefore it has been hypothesized that evolution in constant environments would lead to regression of such clocks. This study examined whether circadian clocks and the associated properties evolve differently under constant light and constant darkness. Activity-rest, adult emergence and oviposition rhythms were measured of D. melanogaster populations that have been maintained for over 19 years (~330 generations) under three different light regimes - constant light (LL), light-dark cycles of 12:12 h (LD) and constant darkness (DD). While circadian rhythms in all the three behaviors persist in both LL and DD stocks with no differences in circadian period, they differed in certain aspects of the entrained rhythms when compared to controls reared in rhythmic environment (LD). Interestingly, it was also observed that DD stocks have evolved significantly higher robustness or power of free-running activity-rest and adult emergence rhythms compared to LL stocks. Thus, this study, in addition to corroborating previous results of circadian clock evolution in constant light, also highlights that, contrary to the expected regression of circadian clocks, rearing in constant darkness leads to the evolution of more robust circadian clocks which may be attributed to an intrinsic adaptive advantage of circadian clocks and/or pleiotropic functions of clock genes in other traits (Shindey, 2017).

    Transcriptional polymorphism of piRNA regulatory genes underlies the mariner activity in D. simulans testes

    During colonization of new areas, natural populations have to deal with changing environments, and transposable elements (TEs) can be useful "tools" in the adaptation process since they are major contributor to the structural and functional evolution of genomes. In this general context, the activity (copy number, transcriptional and excision rate) of the mariner mos1 element was estimated in 19 natural populations of D. simulans. It is shown (1) that mos1 expression is always higher and more variable in testes than in ovaries; (2) that mos1 activity is higher in colonizing populations compared to the sub-Saharan African ones (ancestral populations); (3) that mos1 variations in transcript levels and copy number are negatively correlated to transcriptional variations of piRNA genes, aubergine and argonaute3. Furthermore, mos1 levels of expression in testes highly contrast with the low expression patterns of ago3. These results strongly suggest that the expression polymorphism of piRNA genes could be responsible for the mos1 variations, first between male and female germlines and second, according to the status of natural populations (colonizing or not). These results provide new perspectives about TEs and piRNA genes co-evolution in Drosophila germlines (Saint-Leandre, 2017).

    A test for gene flow among sympatric and allopatric Hawaiian picture-winged Drosophila

    The Hawaiian Drosophila are one of the most species-rich endemic groups in Hawaii and a spectacular example of adaptive radiation. Drosophila silvestris and D. heteroneura are two closely related picture-winged Drosophila species that occur sympatrically on Hawaii Island and are known to hybridize in nature, yet exhibit highly divergent behavioral and morphological traits driven largely through sexual selection. Their closest-related allopatric species, D. planitibia from Maui, exhibits hybrid male sterility and reduced behavioral reproductive isolation when crossed experimentally with D. silvestris or D. heteroneura. A modified four-taxon test for gene flow was applied to recently obtained genomes of the three Hawaiian Drosophila species. The analysis indicates recent gene flow in sympatry, but also, although less extensive, between allopatric species. This study underscores the prevalence of gene flow, even in taxonomic groups considered classic examples of allopatric speciation on islands. The potential confounding effects of gene flow in phylogenetic and population genetics inference are discussed, as well as the implications for conservation (Kang, 2017).

    Experimental test and refutation of a classic case of molecular adaptation in Drosophila melanogaster

    Identifying the genetic basis for adaptive differences between species requires explicit tests of historical hypotheses concerning the effects of past changes in gene sequence on molecular function, organismal phenotype and fitness. This challenge was addressed by combining ancestral protein reconstruction with biochemical experiments and physiological analysis of transgenic animals that carry ancestral genes. A widely held hypothesis of molecular adaptation was tested in this study-that changes in the alcohol dehydrogenase protein (ADH) along the lineage leading to Drosophila melanogaster increased the catalytic activity of the enzyme and thereby contributed to the ethanol tolerance and adaptation of the species to its ethanol-rich ecological niche. These experiments strongly refute the predictions of the adaptive ADH hypothesis and caution against accepting intuitively appealing accounts of historical molecular adaptation that are based on correlative evidence. The experimental strategy employed can be used to decisively test other adaptive hypotheses and the claims they entail about past biological causality (Siddiq, 2017).

    Adaptive evolution of gene expression in Drosophila

    Gene expression levels are important quantitative traits that link genotypes to molecular functions and fitness. In Drosophila, population-genetic studies have revealed substantial adaptive evolution at the genomic level, but the evolutionary modes of gene expression remain controversial. This study presents evidence that adaptation dominates the evolution of gene expression levels in flies. 64% of the observed expression divergence across seven Drosophila species are adaptive changes driven by directional selection. The results are derived from time-resolved data of gene expression divergence across a family of related species, using a probabilistic inference method for gene-specific selection. Adaptive gene expression is stronger in specific functional classes, including regulation, sensory perception, sexual behavior, and morphology. Moreover, a large group of genes was identifed with sex-specific adaptation of expression, which predominantly occurs in males. This analysis opens an avenue to map system-wide selection on molecular quantitative traits independently of their genetic basis (Nourmohammad, 2017).

    Adaptive evolution of gene expression in Drosophila

    Gene expression levels are important quantitative traits that link genotypes to molecular functions and fitness. In Drosophila, population-genetic studies have revealed substantial adaptive evolution at the genomic level, but the evolutionary modes of gene expression remain controversial. This study presents evidence that adaptation dominates the evolution of gene expression levels in flies. 64% of the observed expression divergence across seven Drosophila species are adaptive changes driven by directional selection. The results are derived from time-resolved data of gene expression divergence across a family of related species, using a probabilistic inference method for gene-specific selection. Adaptive gene expression is stronger in specific functional classes, including regulation, sensory perception, sexual behavior, and morphology. Moreover, a large group of genes was identified with sex-specific adaptation of expression, which predominantly occurs in males. This analysis opens an avenue to map system-wide selection on molecular quantitative traits independently of their genetic basis (Nourmohammad, 2017).

    The inference of adaptation exploits the complex dependence of the expression divergence on the evolutionary distance between species. It reflects two fundamental evolutionary features of quantitative traits. First, such traits generate a divergence pattern with two distinct molecular clocks: at a short evolutionary distance, the divergence is always near the expected value under neutrality; at a longer distance, it depends jointly on stabilizing and directional selection. This feature reconciles seemingly contradictory results of previous studies: analysis of closely related species produces a signal of neutral evolution, whereas evolutionary constraint becomes apparent for more distant species. Second, the phenotypic evolution of gene expression decouples from details of its genetic basis. This explains why overall strong selection on gene expression levels was found even though selection on individual QTLs is often weak. The probabilistic extension of the curren inference scheme, which is based on gene-specific expression divergence, identifies functional gene classes associated with adaptive evolution of regulation (Nourmohammad, 2017).

    The selection model underlying this analysis is a single-peak fitness seascape, which contains components of stabilizing and directional selection on a quantitative trait. These components are well-established notions of quantitative genetics on micro-evolutionary timescales. Each of them can provide a snapshot of the predominant selection pressure in a population. However, the description of selection remains incomplete as a description of selection over macro-evolutionary periods. If selection on a trait is directional at a given evolutionary time, will that selection relax after the trait value has significantly adapted in the direction of selection? If selection is stabilizing, can it be assumed that the optimal trait value will remain invariant in the context of a different species? To address these questions, a conceptual and quantitative synthesis of stabilizing and directional selection is need. The single-peak seascape model arguably provides the simplest such synthesis. It also provides a simple picture of continual adaptation over macro-evolutionary periods: a species follows a moving fitness peak, and this process generates positive fitness flux but no net increase in fitness (Nourmohammad, 2017).

    This method of selection inference can be applied to a spectrum of molecular quantitative traits with a complex genetic basis, provided that comparative data from multiple, sufficiently diverged species are available. Such traits include genome-wide protein levels, protein-DNA binding interactions, and enzymatic activities. For most of these traits, only partial knowledge is available of the underlying genetic loci and their effects on trait and fitness. This method complements QTL studies and opens a way to infer quantitative phenotype-fitness maps at the systems level (Nourmohammad, 2017).

    Regulation of gene expression and RNA editing in Drosophila adapting to divergent microclimates

    Determining the mechanisms by which a species adapts to its environment is a key endeavor in the study of evolution. In particular, relatively little is known about how transcriptional processes are fine-tuned to adjust to different environmental conditions. Here we study Drosophila melanogaster from 'Evolution Canyon' in Israel, which consists of two opposing slopes with divergent microclimates. Several hundred differentially expressed genes and dozens of differentially edited sites were identified between flies from each slope; these changes were correlate with genetic differences, and CRISPR mutagenesis was used to validate that an intronic SNP in prominin regulates its editing levels. It was also demonstrated that while temperature affects editing levels at more sites than genetic differences, genetically regulated sites tend to be less affected by temperature. This work shows the extent to which gene expression and RNA editing differ between flies from different microclimates, and provides insights into the regulation responsible for these differences (Yablonovitch, 2017).

    Inter- and intra-specific genomic divergence in Drosophila montana shows evidence for cold adaptation

    The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. This study describes the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species. Branch tests were used to identify genes showing accelerated divergence in contrasts between cold- and warm adapted species, and about 250 genes were identified that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. Evidence was sought of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, and no strong evidence was found for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Three populations of D. montana were resequenced from its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales (Parker, 2018).

    Effects of evolutionary history on genome wide and phenotypic convergence in Drosophila populations

    Studies combining experimental evolution and next-generation sequencing have found that adaptation in sexually reproducing populations is primarily fueled by standing genetic variation. Consequently, the response to selection is rapid and highly repeatable across replicate populations. Some studies suggest that the response to selection is highly repeatable at both the phenotypic and genomic levels, and that evolutionary history has little impact. Other studies suggest that even when the response to selection is repeatable phenotypically, evolutionary history can have significant impacts at the genomic level. This study tests two hypotheses that may explain this discrepancy. Hypothesis 1: Past intense selection reduces evolutionary repeatability at the genomic and phenotypic levels when conditions change. Hypothesis 2: Previous intense selection does not reduce evolutionary repeatability, but other evolutionary mechanisms may. These hypotheses were tested using D. melanogaster populations that were subjected to 260 generations of intense selection for desiccation resistance and have since been under relaxed selection for the past 230 generations. It was found that, with the exception of longevity and to a lesser extent fecundity, 230 generations of relaxed selection has erased the extreme phenotypic differentiation previously found. No signs were found of genetic fixation, and only limited evidence of genetic differentiation between previously desiccation resistance selected populations and their controls. These findings suggest that evolution in this system is highly repeatable even when populations have been previously subjected to bouts of extreme selection. It is therefore concluded that evolutionary repeatability can overcome past bouts of extreme selection in Drosophila experimental evolution, provided experiments are sufficiently long and populations are not inbred (Phillips, 2018).

    Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

    Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. This study aimed to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae. A phylogenetic and population genetic analysis framework was used to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. Strong evidence was found of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains. The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts (Dizz, 2018).

    Ultrastructural variation and adaptive evolution of the ovipositor in the endemic Hawaiian Drosophilidae

    Ecological diversification of the endemic Hawaiian Drosophilidae has been accompanied by striking divergence in egg morphology, and ovarian structure and function. To determine how these flies successfully oviposit in a variety of breeding substrates, Scanning Electron Microscopy was used to examine the ultrastructure of the ovipositor of a sample of 65 Drosophila species and five Scaptomyza species of this hyperdiverse monophyletic group. The Drosophila species analyzed included representatives of the fungus-breeding haleakalae group, the leaf-breeding antopocerus and modified tarsus groups, the modified mouthparts species group, the nudidrosophila, and the picture wing clade; the latter sample of 41 species from four species groups included stem- and bark-breeders, as well as tree sap flux-breeders. Ovipositor length was found to vary more than 12-fold among Hawaiian drosophilids, with the longest ovipositors observed in the bark-breeding species and the shortest among the Scaptomyza and fungus-breeders. More noteworthy is the striking variation in overall shape and proportions of the ovipositor, in the shape of the apical region, and in the pattern of sensory structures or ovisensilla. Ultrastructural observations of the pair of long subapical sensilla on the ventral side identify these as taste bristles. Ovipositor form correlates strongly with the oviposition substrate used by the species, being of a distinctive shape and size in each case. It is inferred that the observed morphological divergence in the ovipositor is adaptive and the product of natural selection for successful reproduction in alternate microhabitats (Craddock, 2018).

    Estimating the timing of multiple admixture pulses during local ancestry inference

    Admixture, the mixing of genetically distinct populations, is increasingly recognized as a fundamental biological process. One major goal of admixture analyses is to estimate the timing of admixture events. Whereas most methods today can only detect the most recent admixture event, this study presents coalescent theory and associated software that can be used to estimate the timing of multiple admixture events in an admixed population. This approach was extensively validated and the conditions under which it can successfully distinguish one from two-pulse admixture models was validated. This approach was applied to real and simulated data of Drosophila melanogaster. Evidence was found of a single very recent pulse of cosmopolitan ancestry contributing to African populations as well as evidence for more ancient admixture among genetically differentiated populations in sub-Saharan Africa. These results suggest this method can quantify complex admixture histories involving genetic material introduced by multiple discrete admixture pulses. The new method facilitates the exploration of admixture and its contribution to adaptation, ecological divergence, and speciation (Medina, 2018).

    The making of a pest: Insights from the evolution of chemosensory receptor families in a pestiferous and invasive fly

    Drosophila suzukii differs from other melanogaster group members in their proclivity for laying eggs in fresh fruit rather than in fermenting fruits. Earlier work has revealed how the olfactory landscape of D. suzukii is dominated by volatiles derived from its unique niche. This study annotated the Olfactory receptors and Gustatory Receptors in D. suzukii and two close relatives, D. biarmipes and D. takahashii, to identify candidate chemoreceptors associated with D. suzukii's unusual niche utilization. A total of 71 Or genes were annotated in D. suzukii, with nine of those being pseudogenes (12.7 %). Alternative splicing of two genes brings the total to 62 genes encoding 66 Ors. Duplications of Or23a and Or67a expanded D. suzukii's Or repertoire, while pseudogenization of Or74a, Or85a, and Or98b reduced the number of functional Ors to roughly the same as other annotated species in the melanogaster group. Seventy-one intact Gr genes and three pseudogenes were annotated in D. suzukii. Alternative splicing in three genes brings the total number of Grs to 81. Signatures of positive selection were identified in two Ors and three Grs at nodes leading to D. suzukii, while three copies in the largest expanded Or lineage, Or67a, also showed signs of positive selection at the external nodes. This analysis of D. suzukii's chemoreceptor repertoires in the context of nine melanogaster group drosophilids, including two of its closest relatives (D. biarmipes and D. takahashii), revealed several candidate receptors associated with the adaptation of D. suzukii to its unique ecological niche (Hickner, 2016).

    A genome-wide scan for genes under balancing selection in Drosophila melanogaster

    In the history of population genetics balancing selection has been considered as an important evolutionary force, yet until today little is known about its abundance and its effect on patterns of genetic diversity. Several well-known examples of balancing selection have been reported from humans, mice, plants, and parasites. However, only very few systematic studies have been carried out to detect genes under balancing selection. This study carried out a genome scan in Drosophila melanogaster to find signatures of balancing selection in a derived (European) and an ancestral (African) population. A total of 34 genomes were scanned, searching for regions of high genetic diversity and an excess of SNPs with intermediate frequency. In total, 183 candidate genes were found: 141 in the European population and 45 in the African one, with only three genes shared between both populations. Most differences between both populations were observed on the X chromosome, though this might be partly due to false positives. Functionally, an overrepresentation of genes involved in neuronal development and circadian rhythm were found. Furthermore, some of the top genes identified are involved in innate immunity. These results revealed evidence of genes under balancing selection in European and African populations. More candidate genes have been found in the European population. They are involved in several different functions (Croze, 2017).

    Nucleotide diversity inflation as a genome-wide response to experimental lifespan extension in Drosophila melanogaster

    Evolutionary theory predicts that antagonistically selected alleles, such as those with divergent pleiotropic effects in early and late life, may often reach intermediate population frequencies due to balancing selection, an elusive process when sought out empirically. Alternatively, genetic diversity may increase as a result of positive frequency-dependent selection and genetic purging in bottlenecked populations. While experimental evolution systems with directional phenotypic selection typically result in at least local heterozygosity loss, this study reports that selection for increased lifespan in Drosophila melanogaster leads to an extensive genome-wide increase of nucleotide diversity in the selected lines compared to replicate control lines, pronounced in regions with no or low recombination, such as chromosome 4 and centromere neighborhoods. These changes, particularly in coding sequences, are most consistent with the operation of balancing selection and the antagonistic pleiotropy theory of aging and life history traits that tend to be intercorrelated. Genes involved in antioxidant defenses, along with multiple lncRNAs, were among those most affected by balancing selection. Despite the overwhelming genetic diversification and the paucity of selective sweep regions, two genes with functions important for central nervous system and memory, Ptp10D and Ank2, evolved under positive selection in the longevity lines. Overall, the 'evolve-and-resequence' experimental approach proves successful in providing unique insights into the complex evolutionary dynamics of genomic regions responsible for longevity (Michalak, 2017).

    Genetic convergence in the evolution of male-specific color patterns in Drosophila

    Convergent evolution provides a type of natural replication that can be exploited to understand the roles of contingency and constraint in the evolution of phenotypes and the gene networks that control their development. For sex-specific traits, convergence offers the additional opportunity for testing whether the same gene networks follow different evolutionary trends in males versus females. Thus study used an unbiased, systematic mapping approach to compare the genetic basis of evolutionary changes in male-limited pigmentation in several pairs of Drosophila species that represent independent evolutionary transitions. A strong evidence for repeated recruitment of the same genes to specify similar pigmentation in different species was found. At one of these genes, ebony, convergent evolution of sexually dimorphic and monomorphic expression through cis-regulatory changes was observed. However, this functional convergence has a different molecular basis in different species, reflecting both parallel fixation of ancestral alleles and independent origin of distinct mutations with similar functional consequences. These results show that a strong evolutionary constraint at the gene level is compatible with a dominant role of chance at the molecular level (Signor, 2016).

    Recent studies have led to an emerging consensus that convergent phenotypes often reflect evolutionary changes in the same genes (the 'evolutionary hotspot' model). However, most case studies supporting the hotspot model relied heavily on candidate gene approaches, which results in a positive ascertainment bias and some difficulty in interpreting the results. Although genetic mapping is more laborious than candidate gene analysis, it offers several key advantages: it is unbiased in that the loci implicated in other species are no more likely to be discovered than any other genes, it produces direct evidence of a gene's causative role in trait evolution (and, equally importantly, allows other genes to be ruled out), and it provides a quantitative estimate of the relative importance of each gene to the overall genetic architecture of a phenotype. This study applied this approach to the ananassae species subgroup, in which multiple species have independently evolved similar male-specific color patterns. The convergent phenotypes were found ti be controlled by overlapping, but not identical, sets of genes in different evolutionary contrasts. The only gene that was implicated in all three high-resolution mapping crosses and could not be ruled out in the fourth, low-resolution cross (D. m. malerkotliana/D. bipectinata) was ebony. Moreover, ebony is the strongest QTL in all contrasts, contributing 35% to 80% of overall divergence. ebony has previously been implicated in both intraspecific and interspecific differences in pigmentation in several other Drosophila species. Clearly, ebony fits the definition of an evolutionary hotspot. ebony encodes an enzyme, β-alanyl-dopamine synthase, that synthesizes light pigment precursors so that higher ebony expression causes lighter pigmentation. Like other pigmentation genes, it has additional roles in other tissues, such as control of circadian locomotor activity in brain glial cells. Every time ebony has been implicated in the evolution of pigmentation, cis-regulatory rather than coding mutations were involved. Presumably, this pattern reflects the ability of cis-regulatory mutations to overcome pleiotropic constraints by uncoupling gene functions in different cell types (Signor, 2016).

    Repeated involvement of the same gene in multiple phenotypic transitions could potentially result from different evolutionary processes: recurrent de novo mutation, lineage sorting of ancestral variation, or interspecific introgression. Adaptation from standing variation is likely to be faster than awaiting a new mutation because potentially beneficial alleles are available immediately and are likely to be present at higher frequencies than alleles arising de novo. Cases of parallel fixation appear to be fairly common, most notably when a reservoir of newly adaptive alleles is available to colonizing populations. Similarly, introgressive hybridization has been implicated, among other traits, in the evolution of wing color patterns in Heliconius butterflies and high-altitude adaptation in humans. Genome-wide analyses suggest that interspecific introgression may play a more important role in evolution than previously thought. This study found that in D. malerkotliana, D. bipectinata, and D. parabipectinata, at least some of the putative causative SNPs at the ebony locus are shared across species, most likely reflecting ancestral lineage sorting. However, these variants are not found either in D. pseudoananassae or in the ercepeae species complex, suggesting a role for independent ebony mutations in these taxa. The region of the ebony locus associated with the evolution of color patterns in the bipectinata complex evolves so rapidly that it cannot be aligned with other taxa. Outside of the ananassae subgroup, evolutionary changes in ebony expression that lead to divergent pigmentation map to a different, more upstream region of the gene; population-genetic analysis rules out this region in the bipectinata complex. It can be concluded that the widely convergent involvement of ebony in the evolution of color patterns is not due solely to fixation of pre-existing variation, but reflects independent origin of distinct, though functionally similar, cis-regulatory mutations (Signor, 2016).

    Cuticle color depends not just on the level of ebony, but on the balance between ebony, tan, yellow, and potentially other enzymes. For example, tan encodes a β-alanyl-dopamine hydrolase, which reverses the chemical reaction catalyzed by ebony. Although tan is implicated in the evolution of pigmentation in several species, it has been ruled out in many other studies including this one (Table S5). Interestingly, the evolutionary changes that were found to involve tan are never male limited; for example, tan controls a female-limited color polymorphism in D. erecta, and the secondary loss of pigmentation in D. santomea affects both sexes. It is possible that ebony could be favored as the male hotspot due to its dosage sensitivity and chromosomal location. ebony, but not tan or yellow, has a semi-dominant loss-of-function phenotype, suggesting that cis-regulatory mutations in ebony could be more readily visible to directional selection. At the same time, yellow and tan are X-linked, whereas ebony is autosomal, suggesting that it could harbor higher levels of standing genetic variation when not under directional selection. Interestingly, in D. melanogaster, D. americana, and the bipectinata species complex, the role of ebony in phenotypic evolution appears to derive from a combination of pre-existing and de novo mutations. Thus, chromosomal sex, the topology of the regulatory network, the kinetics of the pigment synthesis pathway, and population-genetic factors may all contribute to the evolutionary hotspot status of ebony (Signor, 2016).

    In an accompanying paper (Yassin, 2016), a similarly unbiased and systematic approach was taken to map the genetic basis of natural variation in female-specific abdominal pigmentation in multiple species of the Drosophila montium species subgroup, which is closely related to the ananassae lineage. This variation was found to map to the pdm3 transcription factor in several distantly related species. Moreover, convergent involvement of pdm3 appears to reflect independent mutations in this gene in different species. In contrast, ebony does not contribute to color pattern variation in any of the four montium-subgroup species examined. Why do evolutionary hotspots differ in closely related lineages? Is this a matter of historical contingency or different gene network topology in different clades? Or are different genes within a shared network favored for cis-regulatory evolution in different sexes, resulting in sex-specific evolutionary hotspots? Intriguingly, in the only montium subgroup species in which ebony was implicated in color pattern variation, the pigmentation phenotype is male specific, supporting the latter hypothesis (Signor, 2016).

    One can envision two principal mechanisms for the gain and loss of sex-specific traits. First (the instructive model), the genes responsible for phenotypic changes may be the same genes that are regulated in a dimorphic manner to generate sex-specific phenotypes, as was observed for ebony. Alternatively (the permissive model), the causative genes could be monomorphic, while sexual dimorphism is encoded in parallel to or downstream of these genes. For example, ebony could have been expressed equally between sexes in all species, while a different, 'gatekeeper' gene is sexually dimorphic in all species. In the latter scenario, high levels of ebony expression in both sexes would mask the dimorphic phenotypes promoted by the sex-specific gatekeeper gene, while low ebony levels would uncover the underlying dimorphism. Thus, ebony would be the causative gene responsible for differences between lineages, while the gatekeeper gene is responsible for sexual dimorphism (Signor, 2016).

    These two models suggest very different mechanisms for the evolution of sexual dimorphism. Under the instructive model, gain and loss of sex-specific traits is caused by frequent changes in sex-specific gene regulation. Under the permissive model, the targets of the sex determination pathway can remain static over long evolutionary distances, while the underlying sexual dimorphism is revealed or obscured by sexually monomorphic genetic changes that occur elsewhere in the developmental pathway (Signor, 2016).

    The current results argue for the instructive model: in all evolutionary contrasts, sex-specific pigmentation is associated with sex-specific ebony expression, and sexually monomorphic pigmentation is associated with monomorphic expression. Thus, sex-specific transcriptional regulation of ebony has been gained and/or lost several times within the ananassae species subgroup. A similar pattern has been observed in the expression of desat-F, a hydrocarbon desaturase enzyme involved in the synthesis of cuticular pheromones. It appears that sex-specific gene regulation can be gained and lost quite easily over short evolutionary timescales and that the evolution of sexually dimorphic traits is more likely to follow the instructive model (Signor, 2016).

    Although the study of color pattern evolution in Drosophila has largely been dominated by candidate gene analyses, it has now been enriched by several unbiased, high-resolution genetic mapping studies, including this work. These studies have gone beyond spotlighting individual genes to provide a more holistic picture of the genetic architecture of evolutionary changes and have confirmed the predominance of cis-regulatory mutations in phenotypic evolution. Although no single gene is involved in all cases, the number of players appears to be limited, and most genes have been implicated repeatedly in multiple taxa. Collectively, parallel genetic analyses in multiple species suggest a 'toolkit model' of convergent evolution. For any trait, there are a limited number of genes that can potentially evolve to produce phenotypic changes. Within that toolkit, the relative likelihood of each gene's involvement may depend on its position in the regulatory network that controls the development of that trait, on the historical contingencies specific to each evolving lineage, and, potentially, on sex. Together, these trends result in a pattern where convergent phenotypes have distinct yet overlapping genetic basis in different species (Signor, 2016).

    Patterning the insect eye: From stochastic to deterministic mechanisms

    While most processes in biology are highly deterministic, stochastic mechanisms are sometimes used to increase cellular diversity. In human and Drosophila eyes, photoreceptors sensitive to different wavelengths of light are distributed in stochastic patterns, and one such patterning system has been analyzed in detail in the Drosophila retina. Interestingly, some species in the dipteran family Dolichopodidae (the "long legged" flies, or "Doli") instead exhibit highly orderly deterministic eye patterns. In these species, alternating columns of ommatidia (unit eyes) produce corneal lenses of different colors. Occasional perturbations in some individuals disrupt the regular columns in a way that suggests that patterning occurs via a posterior-to-anterior signaling relay during development, and that specification follows a local, cellular-automaton-like rule. It is hypothesized that the regulatory mechanisms that pattern the eye are largely conserved among flies and that the difference between unordered Drosophila and ordered dolichopodid eyes can be explained in terms of relative strengths of signaling interactions rather than a rewiring of the regulatory network itself. A simple stochastic model is presented that is capable of explaining both the stochastic Drosophila eye and the striped pattern of Dolichopodidae eyes and thereby characterize the least number of underlying developmental rules necessary to produce both stochastic and deterministic patterns. Only small changes to model parameters are needed to also reproduce intermediate, semi-random patterns observed in another Doli species, and quantification of ommatidial distributions in these eyes suggests that their patterning follows similar rules (Ebadi, 2018).

    Evolution of herbivory in Drosophilidae linked to loss of behaviors, antennal responses, odorant receptors, and ancestral diet

    Herbivory is a key innovation in insects, yet has only evolved in one-third of living orders. The evolution of herbivory likely involves major behavioral changes mediated by remodeling of canonical chemosensory modules. Herbivorous flies in the genus Scaptomyza (Drosophilidae) are compelling species in which to study the genomic architecture linked to the transition to herbivory because they recently evolved from microbe-feeding ancestors and are closely related to Drosophila melanogaster. This study found that Scaptomyza flava, a leaf-mining specialist on plants in the family (Brassicaceae), was not attracted to yeast volatiles in a four-field olfactometer assay, whereas D. melanogaster was strongly attracted to these volatiles. Yeast-associated volatiles, especially short-chain aliphatic esters, elicited strong antennal responses in D. melanogaster, but weak antennal responses in electroantennographic recordings from S. flava. The genome of S. flava was sequenced, and this species' odorant receptor repertoire was characterized. Orthologs of odorant receptors, which detect yeast volatiles in D. melanogaster and mediate critical host-choice behavior, were deleted or pseudogenized in the genome of S. flava. These genes were lost step-wise during the evolution of Scaptomyza. Additionally, Scaptomyza has experienced gene duplication and likely positive selection in paralogs of Or67b in D. melanogaster. Olfactory sensory neurons expressing Or67b are sensitive to green-leaf volatiles. Major trophic shifts in insects are associated with chemoreceptor gene loss as recently evolved ecologies shape sensory repertoires (Goldman-Huertas, 2005).

    Understanding the origins and consequences of trophic shifts, especially the transition to herbivory, has been a central problem in evolutionary biology. The paleontological record suggests that evolutionary transitions to herbivory have been rare in insects, and the first transitions to herbivory in vertebrates occurred long after the colonization of land. However, species radiations result from herbivorous transitions in insects and vertebrates, suggesting that herbivory is a key innovation. Identifying functional genomic changes associated with the evolutionary transition to herbivory could yield insight into the mechanisms that have driven their success. However, the origins of the most diverse clades of herbivorous insects are ancient and date to the Jurassic or earlier, limiting meaningful genomic comparisons. In contrast, herbivory has evolved more times in Diptera than in any other order. The Drosophilidae is an excellent system to study the evolution of herbivory from a functional genomic perspective because it includes several transitions to herbivory, and the genomic model Drosophila melanogaster (Goldman-Huertas, 2005).

    The transition to herbivory involves adaptations in physiology, morphology, and behavior. The evolution of sensory repertoires could reinforce or even precipitate these adaptations through adaptive loss or relaxation of functional constraint subsequent to a trophic shift. Adaptive loss of chemoreceptors has been rarely shown but occurs in nematodes, although their olfactory systems are distinct from insects. Families of mammalian olfactory receptor proteins have been remodeled during transitions to flight, aquatic lifestyles, and frugivory. Similarly, the evolution of diet specialization in Drosophila species correlates with chemoreceptor gene losses, and hematophagous flies have lost gustatory receptors that detect sweet compounds. More profound changes such as the evolution of new protein families are associated with major evolutionary transitions such as the evolution of flight in insects. Although gene loss is unlikely to be a driving force of innovation, loss-of-function mutations may be exeptations that allow novel behaviors to evolve by disrupting ancestral attractions. If detection of different chemical cues becomes selected in a novel niche, then loss through relaxed constraint may indicate which chemical cues have changed during a trophic shift (Goldman-Huertas, 2005).

    The chemosensory repertoires of many drosophilid species have been functionally annotated. The genus Drosophila includes 23 species with published genome sequences, and D. melanogaster presents the most fully characterized insect olfactory system, allowing potential linkage of receptor remodeling to a mechanistic understanding of behavioral change (Goldman-Huertas, 2005).

    Most drosophilids feed on yeast and other microbes growing on decaying plant tissues. Adult female D. melanogaster and distantly related species innately prefer yeast chemical cues to those produced by the fruit on which they oviposit. D. melanogaster detects volatiles with chemoreceptors of several different protein families, but especially receptors from the odorant receptor (OR) gene family, some of which, such as Or42b, are highly conserved across species. Or42b is necessary for attraction and orientation to vinegar and aliphatic esters. Similar compounds activate Or42b across many Drosophila species, suggesting that volatile cues for yeast, and the associated receptors, are conserved across the Drosophilidae (Goldman-Huertas, 2005).

    The ancestral feeding niche for the genus Scaptomyza (Drosophilidae) is microbe-feeding, but Scaptomyza use decaying leaves and stems rather than the fermenting fruit used by D. melanogaster and other members of the subgenus Sophophora. The close association of Scaptomyza with decaying plant tissues may have precipitated the evolution of herbivory <20 MyBP. Adult females of the species S. flava feed and oviposit on living leaves of many cruciferous plants (Brassicales) including Arabidopsis thaliana. Females puncture leaves with serrated ovipositors to create feeding and oviposition sites, and larvae mine and pupate within the living leaves (Goldman-Huertas, 2005).

    This study used Scaptomyza as a model to test the hypothesis that functional loss of chemosensory genes has played a role in a major ecological transition to herbivory in insects. It was hypothesized that the conserved detection of yeast volatiles would be lost in the herbivorous Scaptomyza lineage. This loss was tested by comparing D. melanogaster and S. flava at behavioral, physiological, and genetic levels. First, it was hypothesized that gravid ovipositing S. flava females would not be attracted to yeast volatiles. Second, it was hypothesized that the olfactory sensory organs of S. flava would have a decreased ability to detect individual yeast volatiles and volatile mixtures. Third, chemoreceptor genes from the OR gene family implicated in detection of yeast volatiles would be lost in the S. flava genome. Finally, it was predicted that chemoreceptor genes potentially mediating detection of plant volatiles would show evidence of positive selection and possibly, neofunctionalization (Goldman-Huertas, 2005).

    Olfaction is used by insects to find resources, mates, and oviposition substrates. This study tested the hypothesis that S. flava is not attracted to yeast volatiles, whereas D. melanogaster is attracted to yeast volatiles. A four-field olfactometer assay was used in which filtered air blown through four corners of a diamond-shaped arena establishes four independent airfields. Two of the four fields were exposed to yeast volatiles from Saccharomyces cerevisiae cultures. The presence of gravid adult females of both species in either yeast or control fields was recorded every 6 s for 10 min. D. melanogaster flies spent 82.4 ± 18.2% SD of the assay time in yeast-volatile fields and more time in yeast-volatile fields than S. flava. S. flava did not spend more time in yeast-volatile fields and divided residence time evenly between yeast and control fields, consistent with a loss of attraction to yeast volatiles in S. flava flies (Goldman-Huertas, 2005).

    Because S. flava flies failed to increase their residence time in olfactometer quadrants exposed to yeast volatiles, it was hypothesized that S. flava antennal olfactory sensory neurons (OSNs) were deficient in their ability to detect yeast volatiles. This hypothesis was tested by conducting electroantennogram (EAG) measurements in adult D. melanogaster and S. flava flies of both sexes 4-20 d after eclosion, exposed to the same yeast volatiles used in the olfactometer assays and to crushed rosette leaves of the host plant of S. flava flies in laboratory colonies (Arabidopsis thaliana accession Col-0). EAG responses are driven by the aggregate depolarization of OSNs in the antennae and scale with the concentration and identity of stimulants. No difference were found between sexes and data for male and female flies were combined. Consistently lower EAG signals were recorded in S. flava flies compared with D. melanogaster, preventing interspecific comparisons of signal amplitude, possibly due to differences in electrical properties of antennae (Goldman-Huertas, 2005).

    The antennae of S. flava were more strongly stimulated by Arabidopsis volatiles than by yeast, whereas the antennae of D. melanogaster were more responsive to volatiles from yeast than those from Arabidopsis. Responses were recorded to a small panel of three volatiles associated with A. thaliana [(Z)-3-hexenol, myrcene, phenethyl isothiocyanate] and two with S. cerevisiae (2-phenylethanol, ethyl acetate). Antennae of both species detected all volatiles compared with a negative control. The antennae of S. flava were most responsive to (Z)-3-hexenol, a volatile produced by damaged leaves of many plant species, and were also highly attuned to phenethyl isothiocyanate, a hydrolyzed product of glucosinolates, which are the major defensive compound in host plants of S. flava. Responses to myrcene and 2-phenylethanol were not in the expected direction, although 2-phenylethanol, as a widespread floral volatile, may remain an important chemical cue for Scaptomyza adults (Goldman-Huertas, 2005).

    Antennae of S. flava were less responsive to yeast and the yeast-associated volatile ethyl acetate than to plant-related volatiles, but these relative comparisons were insufficient to prove that the detection threshold for yeast volatiles had decreased in Scaptomyza. Therefore the sensitivity of S. flava and D. melanogaster flies to this and other short-chain aliphatic esters was tested by exposing females to half-log dilution series of ethyl acetate, ethyl propionate, and isobutyl acetate. Sensitivity was defined as the first concentration increase that generated an increased antennal response. S. flava was insensitive to ethyl acetate at the concentrations tested. S. flava was also less responsive to ethyl propionate and isobutyl acetate compared with D. melanogaster. Scaptomyza is considerably less sensitive to short aliphatic esters, which may account for differences in signal strength in response to plant and yeast volatile mixtures and the lack of attraction to yeast volatiles by S. flava. This unresponsiveness is consistent with the fact that deficits in the production of aliphatic esters in a yeast strain decreased attractiveness to D. melanogaster flies (Goldman-Huertas, 2005).

    The lack of attraction and minimal EAG response to yeast volatiles in S. flava suggested that chemosensory genes have been lost or changed in herbivorous Scaptomyza species. ORs are expressed in the dendrites of OSNs in the antennae and maxillary palps and are the primary receptors by which most neopteran insects detect odors in their environments. The OR family has been functionally annotated in D. melanogaster, and members of subfamily H OR genes in particular are highly conserved and enriched in receptors for aliphatic esters, a group of compounds S. flava detected poorly (Goldman-Huertas, 2005).

    To characterize changes in the OR gene repertoire in S. flava associated with the olfactory phenotypes, the genome of S. flava and annotated OR genes were annotated by using reciprocal tBLASTn searches of previously annotated Drosophila OR protein sequences against this de novo S. flava genome assembly. 65 full-length ORFs for OR genes were found in S. flava. Consistent with previous OR gene-naming conventions, ORs were named after the D. melanogaster ortholog or the most closely related gene, with the exception of OrN1 and OrN2 orthologs, which are not present in D. melanogaster (Goldman-Huertas, 2005).

    Protein translations of S. flava genes were included in a phylogeny of D. melanogaster, Drosophila virilis, Drosophila mojavensis, and Drosophila grimshawi OR protein sequences to assess homology. The latter three species are the closest relatives of Scaptomyza with fully sequenced genomes (Goldman-Huertas, 2005).

    S. flava retains duplicates of Or42a, Or67a, Or74a, Or83c, Or98a, and OrN2 found in other sequenced Drosophila species. Scaptomyza also has duplications not shared with close relatives, although nine of these genes are pseudogenized. The majority of paralogs (56%) were found on the same scaffold in tandem arrays. The functional significance of these gene duplications is not yet clear, but it is suggestive that Or67b, with three copies in S. flava, is in single copy in nearly all sequenced Drosophila. In D. melanogaster, neurons expressing Or67b respond to green leaf volatiles such as (Z)-3-hexenol, to which S. flava also has a robust antennal response (Goldman-Huertas, 2005).

    Only four widely conserved ORs were uniquely lost (Or22a and Or85d) or pseudogenized (Or9a, Or42b) in the Scaptomyza lineage. Syntenic regions flanking OR losses were recovered in the genome assembly. Orthologs of Or9a, Or22a, and Or42b are intact in 23 Drosophila species with genome sequences, and Or85d is missing only in the Drosophila albomicans and Drosophila rhopaloa genome assemblies. As predicted, orthologs of ORs that persist in microbe-feeding Drosophila species and are lost in S. flava, function in yeast-volatile detection. Or42b is highly conserved in sequence among Drosophila species, and the receptor is highly attuned to aliphatic esters at low concentrations. Knockouts of Or42b in adult D. melanogaster result in failure to orient in flight toward aliphatic ester odor plumes , and rescuing these neurons restores attraction to yeast volatiles. Similarly, no sequences similar to Or22a were present in the S. flava assembly, although conserved intergenic regions were found in S. flava that flank Or22a in other Drosophila species. Or22a also detects aliphatic esters and in the specialist species Drosophila erecta and Drosophila sechellia, Or22a detects volatiles produced by host fruit. Both Or22a and Or42b are activated by floral volatiles of Arum palestinum, which mimics yeast fermentation volatiles and attracts a diversity of drosophilids. Finally, Or85d orthologs were not detected in the S. flava genome by BLAST or by inspection of genome regions flanking Or85d in other species. Or85d is expressed in the maxillary palps and in D. melanogaster is responsive to the yeast metabolites 2-heptanol, ethyl acetate, and isoamyl acetate. Or85d is highly sensitive to phenethyl acetate, a common volatile of many yeast species. In D. melanogaster, Or9a is activated by a broad range of ketone-, alcohol-, and carboxylic acid-containing ligands. Some of these ligands, such as acetoin, are common yeast volatiles and strong attractants. The consequences of Or9a pseudogenization will require further study (Goldman-Huertas, 2005).

    A time-calibrated phylogeny of the family Drosophilidae suggests that herbivory evolved in Scaptomyza ca.13.5 million years ago (95% highest posterior density 10.02–17.48 million years ago), overlapping with age ranges inferred from previous analyses. Ancestral state reconstructions were performed in the APE package by using an equal rates model. This analysis indicated that microbe feeding is ancestral in Drosophila and Scaptomyza (99.7% probability) and that herbivory evolved once within the genus Scaptomyza (Goldman-Huertas, 2005).

    It was hypothesized that OR gene losses would coincide with the evolution of herbivory. Degenerate PCR primers were developed from genomes of multiple Scaptomyza and Drosophila species that targeted exonic sequences of Or22a and Or9a, and conserved, flanking, intergenic sequences of Or42b and Or85d (Goldman-Huertas, 2005).

    Gene losses in S. flava were confirmed by PCR screen in three natural populations, with the exception of SflaOr9a-1, which appeared to be present in a functional copy in a population from Arizona. A preliminary genome assembly of Scaptomyza pallida was consistent with PCR screening results for OR loss patterns in this species. The presence/absence of S. flava gene losses was reconstructed along ancestral nodes and found that three of the four OR gene losses in S. flava (Or22a, Or85d, Or42b) coincided with or preceded the evolution of herbivory in Scaptomyza. Losses were shared by herbivorous congeners. Or22a, while lost in S. flava, is intact in the microbe-feeding species Scaptomyza apicata and S. pallida and is also lost in two other herbivorous species, Scaptomyza nigrita and Scaptomyza graminum (Goldman-Huertas, 2005).

    Specialist, microbe-feeding Drosophila species, such as D. sechellia and D. erecta have an accelerated rate of chemoreceptor gene loss, but this pattern could also be due to nearly neutral processes. S. flava feeds almost exclusively on plants within the Brassicales, and it was hypothesized that this species has experienced an accelerated rate of chemosensory gene loss compared with other microbe-feeding Drosophila species. This hypothesis was tested by coding homologous groups of ORs as present or absent in S. flava, D. virilis, D. mojavensis and D. grimshawi (the closest Drosophila relatives of Scaptomyza), and two models of gene loss were inferred in the Brownie software package. No evidence was found for the alternative model of increased rate of loss in Scaptomyza, but it cannot be ruled out that there were insufficient loss events to parameterize the more complex model or that other chemoreceptor gene families have undergone accelerated loss in S. flava. Also, S. flava is oligophagous, feeding on many plant species in the Brassicales, and it is less specialized than D. sechellia and D. erecta (Goldman-Huertas, 2005).

    Because the shift to herbivory in Scaptomyza likely involved many changes in olfactory cues, it was hypothesized that some S. flava OR genes should bear signatures of episodic positive selection, as flies adapted to a novel environment. To test this hypothesis, null and alternative (branch-site) models were inferred in PAML 4.7a where subsets of codons in extant S. flava ORs could evolve under (i) purifying or neutral selection or (ii) purifying, neutral, or positive selection, relative to 12 Drosophila species. A phylogeny-aware alignment program, PRANK, was used to identify regions where indels were probable while minimizing sensitivity to alignment errors. Alignments where more than one taxon had an inferred indel in greater than two regions were trimmed by using Gblocks to remove columns with ambiguous homology(Goldman-Huertas, 2005).

    After correcting for false discovery, two ORs were found in which the branch-site model consistent with episodic positive selection was more likely than the null model. Or88a had the strongest statistical support for the branch-site model. In D. melanogaster, Or88a functions in recognition of male and virgin female conspecifics . Two other branches among the S. flava Or67b paralogs also supported the branch-site model: an ancestral branch preceding a Scaptomyza-specific duplication event and a branch leading to Or67b-3. Homologs of this gene in D. melanogaster encode ORs that respond to the green-leaf volatile (Z)-3-hexenol, one of the most salient ligands found in EAG studies of S. flava. Experimental, functional, and population-based tests are needed to verify whether positive selection has fixed amino acid changes in the Scaptomyza lineage (Goldman-Huertas, 2005).

    It is concluded that trophic transitions in the history of animal life, such as herbivory, may be mediated by genetic changes in chemosensory repertoires. The majority of Drosophilidae feed on microbes, and distantly related drosophilid lineages are attracted by the same yeast-mimicking floral scent produced by A. palestinum. A subset of the ORs stimulated by this scent are highly conserved in other drosophilids, which may be part of a homologous and conserved olfactory circuit used to find fermenting host substrates across the family. It was hypothesized that mutations disrupting the function of OR homologs in this conserved olfactory circuit could mediate the evolution of herbivory or other novel food preferences (Goldman-Huertas, 2005).

    S. flava, an herbivorous drosophilid, has lost orthologs of ORs involved in this generalized yeast olfactory circuit. Consistent with these findings, S. flava did not respond to yeast volatiles in a behavioral assay. Antennae of S. flava were weakly activated by active yeast cultures and short-chain aliphatic esters, key compounds found in yeast volatile blends and known ligands of ORs in D. melanogaster lost in S. flava. However, retention of some ORs implicated in yeast-volatile detection, such as Or92a and Or59b, implies that S. flava may retain the ability to detect some untested yeast compounds (Goldman-Huertas, 2005).

    It is hypothesized that OR genes would be intact in nonherbivorous Scaptomyza and gene losses would coincide with the transition to herbivory. Or22a loss did coincide with the evolution of herbivory, but losses of Or42b and Or85d likely predate the evolution of plant feeding. These more ancient losses of conserved yeast-volatile receptors suggest ancestral Scaptomyza may have already evolved novel olfactory pathways that were later co-opted by herbivorous lineages, and in fact, many Scaptomyza species feed on microbes living within decaying leaves or in leaf mines produced by other insects. Sister groups of many major herbivorous insect lineages also feed on detritus and fungi, suggesting that the transition from microbe feeding to herbivory may be common. The genetic changes that underlie host-finding remain to be identified, but recently duplicated ORs, such as the unique triplication of Or67b in Scaptomyza, are likely candidates for further functional study. Subtle, targeted remodeling of chemoreceptor repertoires may be a general mechanism driving changes in behavior, facilitating trophic shifts and ultimately diversification in animals (Goldman-Huertas, 2005).

    Drosophila adaptation to viral infection through defensive symbiont evolution

    Microbial symbionts can modulate host interactions with biotic and abiotic factors. Such interactions may affect the evolutionary trajectories of both host and symbiont. Wolbachia protects Drosophila melanogaster against several viral infections and the strength of the protection varies between variants of this endosymbiont. Since Wolbachia is maternally transmitted, its fitness depends on the fitness of its host. Therefore, Wolbachia populations may be under selection when Drosophila is subjected to viral infection. This study shows that in D. melanogaster populations selected for increased survival upon infection with Drosophila C virus there is a strong selection coefficient for specific Wolbachia variants, leading to their fixation. Flies carrying these selected Wolbachia variants have higher survival and fertility upon viral infection when compared to flies with the other variants. These findings demonstrate how the interaction of a host with pathogens shapes the genetic composition of symbiont populations. Furthermore, host adaptation can result from the evolution of its symbionts, with host and symbiont functioning as a single evolutionary unit (Faria, 2016).

    Generality of toxins in defensive symbiosis: Ribosome-inactivating proteins and defense against parasitic wasps in Drosophila

    While it has become increasingly clear that multicellular organisms often harbor microbial symbionts that protect their hosts against natural enemies, the mechanistic underpinnings underlying most defensive symbioses are largely unknown. Spiroplasma bacteria are widespread associates of terrestrial arthropods, and include strains that protect diverse Drosophila flies against parasitic wasps and nematodes. Recent work implicated a ribosome-inactivating protein (RIP) encoded by Spiroplasma, and related to Shiga-like toxins in enterohemorrhagic Escherichia coli, in defense against a virulent parasitic nematode in the woodland fly, Drosophila neotestacea. This study tested the generality of RIP-mediated protection by examining whether Spiroplasma RIPs also play a role in wasp protection, in D. melanogaster and D. neotestacea. Strong evidence was found for a major role of RIPs, with ribosomal RNA (rRNA) from the larval endoparasitic wasps, Leptopilina heterotoma and Leptopilina boulardi, exhibiting the hallmarks of RIP activity. In Spiroplasma-containing hosts, parasitic wasp ribosomes show abundant site-specific depurination in the alpha-sarcin/ricin loop of the 28S rRNA, with depurination occurring soon after wasp eggs hatch inside fly larvae. Interestingly, the pupal ectoparasitic wasp, Pachycrepoideus vindemmiae, was found to escape protection by Spiroplasma, and its ribosomes do not show high levels of depurination. It was also shown that fly ribosomes show little evidence of targeting by RIPs. Finally, the genome of D. neotestacea's defensive Spiroplasma was found to encode a diverse repertoire of RIP genes, which are differ in abundance. This work suggests that specificity of defensive symbionts against different natural enemies may be driven by the evolution of toxin repertoires, and that toxin diversity may play a role in shaping host-symbiont-enemy interactions (Ballinger, 2017).

    An innovative ovipositor for niche exploitation impacts genital coevolution between sexes in a fruit-damaging Drosophila

    Limited attention has been given to ecological factors influencing the coevolution of male and female genitalia. The innovative ovipositor of Drosophila suzukii, an invading fruit pest, represents an appealing case to document this phenomenon. The serrated saw-like ovipositor is used to pierce the hard skin of ripening fruits that are not used by other fruit flies that prefer soft decaying fruits. This study highlights another function of the ovipositor related to its involvement in genital coupling during copulation. The morphology and coupling of male and female genitalia in this species were compared to its sibling species, Drosophila subpulchrella, and to an outgroup species, Drosophila biarmipes. These comparisons and a surgical manipulation indicated that the shape of male genitalia in D. suzukii has had to be adjusted to ensure tight coupling, despite having to abandon the use of a hook-like structure, paramere, because of the more linearly elongated ovipositor. This phenomenon demonstrates that ecological niche exploitation can directly affect the mechanics of genital coupling and potentially cause incompatibility among divergent forms. This model case provides new insights towards elucidating the importance of the dual functions of ovipositors in other insect species that potentially induce genital coevolution and ecological speciation (Muto, 2018).

    Evolutionary History of the Hymenoptera

    Hymenoptera (sawflies, wasps, ants, and bees) are one of four mega-diverse insect orders, comprising more than 153,000 described and possibly up to one million undescribed extant species. As parasitoids, predators, and pollinators, Hymenoptera play a fundamental role in virtually all terrestrial ecosystems and are of substantial economic importance. To understand the diversification and key evolutionary transitions of Hymenoptera, most notably from phytophagy to parasitoidism and predation (and vice versa) and from solitary to eusocial life, this study inferred the phylogeny and divergence times of all major lineages of Hymenoptera by analyzing 3,256 protein-coding genes in 173 insect species. These analyses suggest that extant Hymenoptera started to diversify around 281 million years ago (mya). The primarily ectophytophagous sawflies are found to be monophyletic. The species-rich lineages of parasitoid wasps constitute a monophyletic group as well. The little-known, species-poor Trigonaloidea are identified as the sister group of the stinging wasps (Aculeata). Finally, the evolutionary root of bees were located within the apoid wasp family "Crabronidae." These results reveal that the extant sawfly diversity is largely the result of a previously unrecognized major radiation of phytophagous Hymenoptera that did not lead to wood-dwelling and parasitoidism. They also confirm that all primarily parasitoid wasps are descendants of a single endophytic parasitoid ancestor that lived around 247 mya. These findings provide the basis for a natural classification of Hymenoptera and allow for future comparative analyses of Hymenoptera, including their genomes, morphology, venoms, and parasitoid and eusocial life styles (Peters, 2017).

    Phylogeography of the Subgenus Drosophila (Diptera: Drosophilidae): Evolutionary history of faunal divergence between the old and the new worlds

    The current subgenus Drosophila (the traditional immigrans-tripunctata radiation) includes major elements of temperate drosophilid faunas in the northern hemisphere. Despite previous molecular phylogenetic analyses, the phylogeny of the subgenus Drosophila has not fully been resolved: the resulting trees have more or less varied in topology. One possible factor for such ambiguous results is taxon-sampling that has been biased towards New World species in previous studies. In this study, taxon sampling was balanced between Old and New World species, and phylogenetic relationships among 45 ingroup species selected from ten core species groups of the subgenus Drosophila were analyzed using nucleotide sequences of three nuclear and two mitochondrial genes. Based on the resulting phylogenetic tree, ancestral distributions and divergence times were estimated for each clade to test Throckmorton's hypothesis that there was a primary, early-Oligocene disjunction of tropical faunas and a subsequent mid-Miocene disjunction of temperate faunas between the Old and the New Worlds that occurred in parallel in separate lineages of the Drosophilidae. The results of this study substantially support Throckmorton's hypothesis of ancestral migrations via the Bering Land Bridge mainly from the Old to the New World, and subsequent vicariant divergence of descendants between the two Worlds occurred in parallel among different lineages of the subgenus Drosophila. However, the results also indicate that these events took place multiple times over a wider time range than Throckmorton proposed, from the late Oligocene to the Pliocene (Izumitani, 2016).

    Adaptive protein divergence of BMP ligands takes place under developmental and evolutionary constraints

    The bone morphogenetic protein (BMP) signaling network, comprising evolutionary conserved BMP2/4/Decapentaplegic (Dpp) and Chordin/Short gastrulation (Sog), is widely utilized for dorsal-ventral (DV) patterning during animal development. A similar network is required for posterior crossvein (PCV) formation in the Drosophila pupal wing. Although both transcriptional and post-transcriptional regulation of co-factors in the network appears to give rise to tissue-specific and species-specific properties, their mechanisms are incompletely understood. In Drosophila, BMP5-8 type ligands, Screw (Scw) and /aGlass bottom boat (Gbb), form heterodimers with Dpp for DV patterning and PCV development, respectively. Sequence analysis indicates that the Scw ligand contains two N-glycosylation motifs; one being highly conserved between BMP2/4 and BMP5-8 type ligands, and the other being Scw ligand-specific. The data reveal that N-glycosylation of the Scw ligand boosts BMP signaling both in cell culture and in the embryo. In contrast, N-glycosylation modifications of Gbb or Scw ligands reduce the consistency of PCV development. These results suggest that tolerance for structural changes of BMP5-8 type ligands is dependent on developmental constraints. Furthermore, gain and loss of N-glycosylation motifs in conserved signaling molecules under evolutionary constraints appear to constitute flexible modules to adapt to developmental processes (Tauscher, 2016).

    Genomics of parallel adaptation at two timescales in Drosophila

    Two interesting unanswered questions are the extent to which both the broad patterns and genetic details of adaptive divergence are repeatable across species, and the timescales over which parallel adaptation may be observed. Drosophila melanogaster is a key model system for population and evolutionary genomics. Findings from genetics and genomics suggest that recent adaptation to latitudinal environmental variation (on the timescale of hundreds or thousands of years) associated with Out-of-Africa colonization plays an important role in maintaining biological variation in the species. Additionally, studies of interspecific differences between D. melanogaster and its sister species D. simulans have revealed that a substantial proportion of proteins and amino acid residues exhibit adaptive divergence on a roughly few million years long timescale. This study used population genomic approaches to attack the problem of parallelism between D. melanogaster and a highly diverged conger, D. hydei, on two timescales. D. hydei, a member of the repleta group of Drosophila, is similar to D. melanogaster, in that it too appears to be a recently cosmopolitan species and recent colonizer of high latitude environments. Parallelism was observed both for genes exhibiting latitudinal allele frequency differentiation within species and for genes exhibiting recurrent adaptive protein divergence between species. Greater parallelism was observed for long-term adaptive protein evolution and this parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences. Thus, despite the roughly 50 million years of time separating D. melanogaster and D. hydei, and despite their considerably divergent biology, they exhibit substantial parallelism, suggesting the existence of a fundamental predictability of adaptive evolution in the genus (Zhao, 2017).

    Inferring the distribution of selective effects from a time inhomogeneous model

    A Poisson random field model has been developed for estimating the distribution of selective effects of newly arisen nonsynonymous mutations that could be observed as polymorphism or divergence in samples of two related species under the assumption that the two species populations are not at mutation-selection-drift equilibrium. The model is applied to 91 Drosophila genes by comparing levels of polymorphism in an African population of D. melanogaster with divergence to a reference strain of D. simulans. Based on the difference of gene expression level between testes and ovaries, the 91 genes were classified as 33 male-biased, 28 female-biased, and 30 sex-unbiased genes. Under a Bayesian framework, Markov chain Monte Carlo simulations are implemented to the model in which the distribution of selective effects is assumed to be Gaussian with a mean that may differ from one gene to the other to sample key parameters. Based on estimates, the majority of newly-arisen nonsynonymous mutations that could contribute to polymorphism or divergence in Drosophila species are mildly deleterious with a mean scaled selection coefficient of -2.81, while almost 86% of the fixed differences between species are driven by positive selection. There are only 16.6% of the nonsynonymous mutations observed in sex-unbiased genes that are under positive selection in comparison to 30% of male-biased and 46% of female-biased genes that are beneficial. It was also estimated that D. melanogaster and D. simulans may have diverged 1.72 million years ago (Amei, 2019).

    The ability of Drosophila hybrids to locate food declines with parental divergence

    Hybrids are generally less fit than their parental species, and the mechanisms underlying their fitness reductions can manifest through different traits. For example, hybrids can have physiological, behavioral, or ecological defects, and these defects can generate reproductive isolation between their parental species. However, the rate that mechanisms of postzygotic isolation other than hybrid sterility and inviability evolve has remained largely uninvestigated, despite isolated studies showing that behavioral defects in hybrids are not only possible but might be widespread.This work studied a fundamental animal behavior - the ability of individuals to find food - and tested the rate at which it breaks down in hybrids. The ability of hybrids from 94 pairs of Drosophila species to find food was measured, and this ability was shown to decrease with increasing genetic divergence between the parental species and that male hybrids are more strongly (and negatively) affected than females. These findings quantify the rate that hybrid dysfunction evolves across the diverse radiation of Drosophila and highlights the need for future investigations of the genetic and neurological mechanisms that affect a hybrid's ability to find a suitable substrate on which to feed and breed (Turissini, 2017).

    The Drosophila speciation factor HMR localizes to genomic insulator sites

    Hybrid incompatibility between Drosophila melanogaster and D. simulans is caused by a lethal interaction of the proteins encoded by the Hmr (Hybrid male rescue) and Lhr (Lethal hybrid rescue) genes. In D. melanogaster the loss of HMR results in mitotic defects, an increase in transcription of transposable elements and a deregulation of heterochromatic genes. To better understand the molecular mechanisms that mediate HMR's function, this study measured genome-wide localization of HMR in D. melanogaster tissue culture cells by chromatin immunoprecipitation. Interestingly, HMR was found to localize to genomic insulator sites that can be classified into two groups. One group belongs to gypsy insulators and another one borders HP1a bound regions at active genes. The transcription of the latter group genes is strongly affected in larvae and ovaries of Hmr mutant flies. These data suggest a novel link between HMR and insulator proteins, a finding that implicates a potential role for genome organization in the formation of species (Gerland, 2017).

    Biodiversity is the result of the emergence and the extinction of species. New species form by pre- and post-zygotic isolation mediated by genetic incompatibility. One of the best characterized examples of hybrid incompatibility is the gene pair Hybrid male rescue (Hmr) and Lethal hybrid rescue (Lhr). Hmr and Lhr cause hybrid incompatibility between the closely related fly species Drosophila melanogaster and D. simulans. Hmr diverged in both Drosophila sibling species under positive selection. HMR and LHR from both species interact physically and localize predominantly to centromeric regions. A reduction of HMR expression results in a misregulation of transposable elements, satellite DNAs and heterochromatic genes. The major difference between HMR and LHR in D. melanogaster and D. simulans is their substantial difference in protein amounts which has been proposed to result in a lethal gain of function in male hybrids. High levels of HMR and LHR in hybrids and overexpression of these proteins in pure species lead to an increased number of binding sites of the complex. Such spreading phenomena based on protein amount have been observed for several chromatin-associated complexes such as the dosage compensation complex, the polycomb complex or components of pericentromeric heterochromatin. In most cases, the precise mechanisms for targeting and spreading are not fully understood. Interestingly, several of the components involved in these processes show signs of adaptive evolution and differ substantially even in very closely related organisms. This observation has spurred a model of a dynamic genome that drives the adaptive evolution of chromatin-associated factors (Gerland, 2017 and references therein).

    Eukaryotic genomes of closely related species differ mostly in the amount and sequence of repetitive DNA. This DNA is often derived from transposable elements, which are highly mutagenic and are therefore under tight transcriptional control by the cellular machinery. During evolution transposons or transposon-derived sequences occasionally adopted structural or novel cis-regulatory functions, thereby contributing to the evolution of new, species-specific, phenotypic traits. Genomic insulators are a particular class of such novel, fast evolving, cis-regulatory elements that show signs of transposon ancestry. A strong expansion of these elements is observed in arthropods, which also experienced a successive gain in the number of insulator binding proteins during evolution. In fact, the Drosophila genome harbours a large variety of insulator proteins such as CTCF, BEAF-32, Su(Hw), Mod(mdg4) and CP190, which all affect nuclear architecture. Different Drosophila species underwent multiple genomic rearrangements and transposon invasions, which presumably resulted in an adaptive response of regulatory DNA binding factors to maintain spatial and temporal gene expression. For example, binding sites for the insulator proteins BEAF-32 and CTCF show a high degree of variability when compared among very closely related species. The gain of new insulator sites is associated with chromosome rearrangements, new born genes and species-specific transcription regulation. Similar to insulator proteins, which tend to cluster in specific nuclear regions, the speciation factor HMR clusters at centromeres or pericentromeric regions in diploid cells but is also detected at distinct euchromatic regions along the chromosome arms in polytene chromosomes. A unifying feature for many of these sites is their close proximity to binding sits of the Heterochromatin Protein 1 (HP1a), a HMR interactor and a well-characterized heterochromatic mark (Gerland, 2017).

    Various studies describe HMR's localization to heterochromatin, but the molecular details on HMR's binding sites and its recruitment to these sites are not well understood. To get new insights into HMR's association to chromatin, this study measured HMR's genome-wide localization by chromatin immunoprecipitation (ChIP) in the D. melanogaster embryonic S2 cell line. This analysis demonstrated an extensive colocalization of HMR with a subset of insulator sites across the genome. HMR's binding to genomic gypsy insulators, which constitute the major group of its binding sites, is dependent on the residing insulator protein complex. In a second group, HMR borders heterochromatin together with the insulator protein BEAF-32. In agreement with previous low-resolution techniques in cell lines and fly tissue, these binding sites are enriched at pericentromeric regions, the cytological region 31 on the 2nd chromosome and the entire 4th chromosome. At most of these sites, HMR associates to the promoters of actively transcribed genes. Interestingly, these genes code for transcripts that have been reported to be downregulated in Hmr mutant larvae and ovaries. Altogether, these data provide evidence for a functional link between HMR and insulator proteins, which potentially results in hybrid incompatibilities due to the adaptive evolution of these genome-organizing complexes (Gerland, 2017).

    HMR localizes to centromeric and pericentromeric regions in D. melanogaster cell lines as well as in mitotically dividing embryonic cells where it has been suggested to act as a repressor of transposable elements. Mutations in Hmr lead to overexpression of satellite DNA and transposable elements in ovaries and larvae. Such a derepression is also observed in hybrid flies, where HMR and LHR levels are higher than the ones in pure species and result in a widespread distribution of the HMR/LHR complex. To better understand the targeting principles that mediate HMR binding within the D. melanogaster genome, it was asked whether it is possible to identify HMR binding sites by applying ChIP-Seq in the D. melanogaster S2 cell line. Combining this approach with RNAi mediated knockdown experiments this study uncovered a strong colocalization of HMR with gypsy insulator binding sites and demonstrated that HMR binding to these sites depends on the presence of the residing insulator protein complex. Notably, HMR associates only with a subset of all Su(Hw) binding sites, but almost all those sites can be classified as gypsy-like sites bound by CP190 and mod(mdg4) in addition to Su(Hw) (Gerland, 2017).

    Besides dispersed binding of HMR at genomic gypsy insulator sites along the chromosome arms, dense clusters of HMR binding sites were observed around the centromere and on the 4th chromosome where it potentially serves to separate HP1a binding domains from highly active genes. This dense clustering of binding sites around the centromere correlates well with the strong colocalization of HMR signals with the centromeric H3 variant CID in immunolocalization experiments. Due to its biochemical interaction and partial colocalization with the heterochromatin protein HP1a in Drosophila embryos, HMR has been suggested to be a bona-fide heterochromatin component. However, in contrast to HP1a, this study detected very distinct HMR binding sites within the genome. When HMR is found close to an HP1a binding domain, it rather borders it than covering the whole domain. The sharp HMR binding signals and the fact that almost all euchromatic HMR binding sites contain putative insulator elements, suggest a role of HMR in separating chromatin domains. A distinct boundary that separates constitutive heterochromatin from the core centromere has also been postulated by Olszak and colleagues who suggest that transition zones between heterochromatin and euchromatin are hotspots for sites of CID misincorporation. Unfortunately, centromeres are notoriously difficult to study by next generation sequencing due to their highly repetitive nature. In addition, the microscopic resolution is not sufficiently high to allow a distinction between a binding to the core centromere chromatin and the chromatin immediately adjacent to it. Therefore, it cannot be ruled that HMR binds large domains at the central region of the Drosophila centromere. However, the fact that the purification of chromatin containing the centromeric H3 variant CID did not identify HMR, suggests that it may very well also form a boundary between pericentromeric heterochromatin and the core centromere. To which extent and by which mechanism HMR fulfils a functional role at these genomic sites remains to be elucidated (Gerland, 2017).

    The genomic sites, where HMR was found bound next to an HP1a domain, are highly enriched for recognition sites of the insulator protein BEAF-32. Interestingly, a depletion of BEAF-32 in S2 cells results in an increased rate of mitotic defects, which is very reminiscent of the phenotype detected when HMR is depleted. Similarly to flies carrying a mutation in the Hmr gene, flies in which BEAF-32 is only contributed maternally have defects in female fertility. BEAF-32's role in maintaining associated promoter regions in an environment that facilitates high transcription levels has been suggested to be functionally relevant for this phenotype. Strikingly, this study found most HMR/BEAF-32 binding sites located between HP1a containing heterochromatin and the transcription start site of a highly active gene. HP1a chromatin might fulfil a repressive function at these genomic regions and HMR might block this repressive impact on the neighbouring gene body. However, no extensive spreading of HP1a or H3K9me3 is seen upon HMR knockdown suggesting that the repressive effect is not directly mediated by HP1a binding or the HMR knock down not efficient enough. As there is evidence that HP1a can also promote gene transcription, HMR may also function as a co-activator for HP1a. Currently, HMR binding next to HP1a containing chromatin is considered as a unifying feature of transcriptionally affected genes but the potential mechanism by which HMR exerts its function is as yet unknown (Gerland, 2017).

    Although HMR depletion has a substantial effect on the transcription of multiple transposons, HMR was found enriched only at the 5' insulator region of the gypsy or gtwin retrotransposons and to similar sites within the genome that are presumably derived from these elements. These sites are occupied by insulator proteins Su(Hw), CP190 and Mod(mdg4) and often display enhancer blocking activity in transgenic assays. Artificial targeting of HMR to DNA placed between an enhancer and a promoter of a reporter gene can block the transcription activity, suggesting that HMR may indeed play a role in setting up endogenous boundary elements. Similar to what is known for Su(Hw), HMR binding to this class of binding sites is dependent on the presence of the structural protein CP190, which has a key function in the stabilization of insulator protein complexes. However, as no strong physical interaction is observed between CP190 and HMR, the loss of HMR binding upon a reduction of CP190 levels may also be the result of increased nucleosome occupancy. Such increase in Histone H3 binding cannot be observed upon HMR removal suggesting that HMR acts downstream of CP190. Interestingly, CP190 loss impairs HMR binding to gypsy-like insulator sites but has weak effect on HMR binding to sites containing BEAF-32 recognition motifs. Notably, in contrast to BEAF-32, CP190 is not required for oogenesis, suggesting that the lack of HMR binding to the class 1 sites may be responsible for the female sterility phenotype observed in Hmr mutant flies (Gerland, 2017).

    How can the current findings be integrated with the lethal phenotype of increased HMR/LHR levels in male hybrids? It is tempting to speculate that multiple additional binding sites that are observed in hybrids and on polytene chromosomes of fly strains over-expressing HMR constitute boundary regions. An increased binding to such boundaries, which have been shown to cluster and form aggregates in vivo, may trigger a massive change in nuclear architecture. In turn, this could indirectly activate multiple transposable elements similar to what is observed when centromere clustering is disturbed. Such a disturbed nuclear architecture may then trigger the activation of a cell cycle checkpoint which has been previously suggested to be a major cause of hybrid lethality (Gerland, 2017).

    Altogether, these data provide a novel link between HMR and cis-regulatory elements bound by insulator proteins. It is speculated that divergent evolution of such genomic elements and their corresponding binding factors in sibling species is triggering hybrid incompatibilities (Gerland, 2017).

    Persistent postmating, prezygotic reproductive isolation between populations

    Studying reproductive barriers between populations of the same species is critical to understand how speciation may proceed. Growing evidence suggests postmating, prezygotic (PMPZ) reproductive barriers play an important role in the evolution of early taxonomic divergence. However, the contribution of PMPZ isolation to speciation is typically studied between species in which barriers that maintain isolation may not be those that contributed to reduced gene flow between populations. Moreover, in internally fertilizing animals, PMPZ isolation is related to male ejaculate-female reproductive tract incompatibilities but few studies have examined how mating history of the sexes can affect the strength of PMPZ isolation and the extent to which PMPZ isolation is repeatable or restricted to particular interacting genotypes. This study addressed these outstanding questions using multiple populations of Drosophila montana. A recurrent pattern of PMPZ isolation is shown, with flies from one population exhibiting reproductive incompatibility in crosses with all three other populations, while those three populations were fully fertile with each other. Reproductive incompatibility is due to lack of fertilization and is asymmetrical, affecting female fitness more than males. There was no effect of male or female mating history on reproductive incompatibility, indicating that PMPZ isolation persists between populations. No evidence was found of variability in fertilization outcomes attributable to different female x male genotype interactions, and in combination with other results, suggests that PMPZ isolation is not driven by idiosyncratic genotype x genotype interactions. These results show PMPZ isolation as a strong, consistent barrier to gene flow early during speciation and suggest several targets of selection known to affect ejaculate-female reproductive tract interactions within species that may cause this PMPZ isolation (Garlovsky, 2018).

    The hpRNA/RNAi pathway is essential to resolve intragenomic conflict in the Drosophila male germline

    Intragenomic conflicts are fueled by rapidly evolving selfish genetic elements, which induce selective pressures to innovate opposing repressive mechanisms. This is patently manifest in sex-ratio (SR) meiotic drive systems, in which distorter and suppressor factors bias and restore equal transmission of X and Y sperm. This study reveals that multiple SR suppressors in Drosophila simulans (Nmy and Tmy) encode related hairpin RNAs (hpRNAs), which generate endo-siRNAs that repress the paralogous distorters Dox and MDox. All components in this drive network are recently evolved and largely testis restricted. To connect SR hpRNA function to the RNAi pathway, D. simulans null mutants of Dcr-2 and AGO2 were generated. Strikingly, these core RNAi knockouts massively derepress Dox and MDox and are in fact completely male sterile and exhibit highly defective spermatogenesis. Altogether, these data reveal how the adaptive capacity of hpRNAs is critically deployed to restrict selfish gonadal genetic systems that can exterminate a species (Lin, 2018).

    RNA interference (RNAi) has long been recognized as a versatile experimental technique, but its endogenous biological utilities have been less tangible. This topic is in principle more accessible in invertebrates, several of which express diverse endogenous siRNAs (endo-siRNAs) via dedicated RNAi machinery that is distinct from the related miRNA pathway. However, while RNAi mutants in nematodes and flies are compromised at defending viruses and affected by certain extreme environmental perturbations, RNAi mutants generally exhibit few overt phenotypes under non-sensitized conditions (Lin, 2018).

    The biological logic of the Drosophila hairpin RNA (hpRNA) pathway has been described, in which inverted repeat transcripts preferentially generate endo-siRNAs in the testis and repress specific highly complementary mRNAs. Although all known hpRNAs are recently evolved, clear evidence is observed for siRNA:target co-evolution, indicating adaptive properties of this regulatory network. While ovaries detectably express hpRNAs and endo-siRNAs, RNAi mutants have relatively little consequence in females. Instead, genetic ablation of RNAi causes spermatogenesis defects and male subfertility. Nevertheless, Drosophila melanogaster RNAi mutant males are fertile, suggesting this species can formally cope without siRNAs, at least within the laboratory setting (Lin, 2018).

    In searching for other manifestations of the hpRNA pathway, the Winters sex-ratio (SR) system of Drosophila simulans was investigated. This meiotic drive system is absent from D. melanogaster and was born within D. simulans subclade species that diverged ~240,000 years ago. Despite its recent de novo appearance, Winters SR factors have profound activities. The Distorter on X (Dox) promotes X chromosome transmission by suppressing Y-bearing sperm, a patently undesirable 'wild-type' gene activity that must be silenced in order to maintain the D. simulans species. An antidote is encoded by autosomal Not much yang (Nmy), to which an inverted repeat with a sequence similarity to Dox was mapped. Signatures consistent with positive selection on Winters factors have been detected within D. simulans populations, indicating the system is actively evolving under an 'arms-race' scenario. While the relationship of Dox and Nmy was evocative of homology-dependent silencing, there is currently no evidence (1) for molecular species constituting the active output of Nmy, (2) that Nmy directly or indirectly regulates the expression of Dox, (3) whether Winters factors are truly distinct from other SR systems, or (4) that the RNAi pathway participates in SR control (Lin, 2018).

    This study provides first molecular evidence that hpRNA-siRNAs are functional mediators of son protection in the Winters SR system. Moreover, this study reveals that a second, previously uncloned SR suppressor in this species, known as the Durham SR system, involves a previously unknown hpRNA-siRNA locus termed Tmy. Although defined as genetically separable SR systems, this study also shows that Nmy and Tmy are paralogous and have partially overlapping capacity to suppress both Dox as well as its progenitor locus MDox. To demonstrate a connection to the RNAi pathway, CRISPR/Cas9 was used to engineer dcr-2 and ago2 null mutants in this non-model fruit fly species. Remarkably, these exhibit profound, testis-specific phenotypes that are much more severe than their well-studied D. melanogaster counterparts, in that they are completely male sterile due to profound defects in spermatogenesis progression and harbor massive synergistic derepression of Dox and MDox transcripts, consistent with loss of collaborative suppression by Nmy and Tmy hpRNAs (Lin, 2018).

    Altogether, these data demonstrate unanticipated complexity of sex-distorting factors and hpRNA-suppressing loci in D. simulans, all of which are rapidly evolving and none of which exist in D. melanogaster. Thus, RNAi is a key pathway that resolves intragenomic conflict that ensures species survival and fulfills roles in adaptive gonadal gene regulation that are more commonly attributed to the piRNA pathway. This highlights the need to explore a wider range of species to more fully appreciate the evolving functions of germline small RNA pathways (Lin, 2018).

    This study provides critical linkages among the RNAi pathway, hpRNA biogenesis and function, and suppression of SR bias. The evolutionary behavior of SR systems conforms to a 'Red Queen' effect, in which a seemingly static outcome (equal transmission of X and Y sperm) actually involves intense, opposing and rapidly evolving genetic programs. These studies provide striking evidence that endogenous RNAi is a central molecular pathway that resolves SR distortion and may potentially have an impact on hybrid sterility. Molecular and genetic evidence are provided of a potential hierarchy, in that Dox is a prime direct target of Nmy based on the observation that it supplies the majority of targeting siRNAs. Nevertheless, Tmy provides a secondary defense: since Tmy exhibits extensive complementarity to Dox that is non-overlapping with Nmy, Tmy siRNAs maintain targeting to Dox even in nmy mutants. Tmy can directly repress Dox in sensor assays, and most importantly, RNAi mutants exhibit strongly elevated Dox transcripts in testis, consistent with co-targeting endogenous action of both hpRNAs on Dox. Moreover, these principles are shown to be true for MDox, strongly implying this locus as a functional distorter. Overall, this study reveals that multiple evolutionary nascent hpRNAs (Nmy/Tmy) are individually required for species preservation via suppression of Winters and Durham SR (Dox/MDox) distorters in D. simulans (Lin, 2018).

    In the future, further dissection of the genetic contributions of the individual hpRNAs and distorters will shed light on their relative contributions to Winters and Durham SR systems and the extent to which they are distinct systems or partially overlapping as indicated by our studies. This will be a challenge in a non-model fly that lacks the genetic tools available in D. melanogaster but will be important to understand how these newly emerging factors are endowed with such powerful activities. For example, compelling hypotheses to test include whether specific deletions of the Tmy hairpin alone can recapitulate SR bias, whether Tmy exhibits derepression of other distorter factors, and whether nmy/tmy double mutant might exhibit only SR or may prove to recapitulate sterility found in RNAi mutants. Thus far, Nmy and Tmy loci have proven recalcitrant to repeated attempts for CRISPR/Cas9 targeting, and it is not clear whether something about their repeat structure affects this endeavor. Moreover, the close linkage of these loci will be a challenge for any efforts to generate recombinants. Still, it will be worthwhile to pursue the generation of new genetic tools (Lin, 2018).

    Remarkably, there is a third SR distortion system in D. simulans ('Paris'). While it is genetically complex, it was recently shown to depend on HP1D2, a recently evolved paralog of the piRNA factor Rhino. Thus, there are apparently molecular linkages of SR systems with small RNA systems, on both driving and suppressing sides. Although the mechanism of Paris SR remains to be determined, HP1D2 protein localizes to the heterochromatic Y chromosome, which provides a connection to the observation that driving Paris alleles prevent segregation of Y chromatids during meiosis II. On the other hand, the defect in Winters SR appears to be post-meiotic, indicating mechanistic diversity in how depletion of male sperm might be achieved by de novo genes (Lin, 2018).

    On the other hand, the sister species D. melanogaster appears to lack SR systems, testament to the extremely rapid rise and fall of SR systems during evolution. The adaptive properties of hpRNAs make them ideal genetic elements to tame such selfish meiotic drive elements, which are theorized to be under constant cycles of emergence, suppression, and disappearance. There is mounting evidence that genetic systems that manifest in SR defects are often associated with sterility. The findings of this study support the notion that the success or failure to resolve intragenomic conflicts using RNAi would be intimately connected to speciation (Lin, 2018).

    While these roles were unexpected in light of the fact that metazoan RNAi biology has been so challenging to appreciate, the situation would have been different had a species only slightly diverged from D. melanogaster initially been selected as a genetic model. This parallels the inference that RNAi might have been recognized earlier had certain budding yeasts other than Saccharomyces cerevisiae been studied earlier. Indeed, functional studies across a broader phylogeny will be necessary to appreciate the evolving requirements of small RNA regulation, beyond standard model organisms. The availability of the D. simulans RNAi mutants of this study opens the door to molecular identification of novel selfish genetic elements that induce SR bias and/or hybrid sterility, in both the Winters and Durham systems (Lin, 2018).

    Identification of misexpressed genetic elements in hybrids between Drosophila-related species

    Crosses between close species can lead to genomic disorders, often considered to be the cause of hybrid incompatibility, one of the initial steps in the speciation process. How these incompatibilities are established and what are their causes remain unclear. To understand the initiation of hybrid incompatibility, reciprocal crosses were performed between two species of Drosophila (D. mojavensis and D. arizonae) that diverged less than 1 Mya. A genome-wide transcriptomic analysis was performed on ovaries from parental lines and on hybrids from reciprocal crosses. Using an innovative procedure of co-assembling transcriptomes, it was shown that parental lines differ in the expression of their genes and transposable elements. Reciprocal hybrids presented specific gene categories and few transposable element families misexpressed relative to the parental lines. Because TEs are mainly silenced by piwi-interacting RNAs (piRNAs), it was hypothesized that in hybrids the deregulation of specific TE families is due to the absence of such small RNAs. Small RNA sequencing confirmed the hypothesis, and it is therefore proposed that TEs can indeed be major players of genome differentiation and be implicated in the first steps of genomic incompatibilities through small RNA regulation (Lopez-Maestre, 2017).

    Introduction of a male-harming mitochondrial haplotype via 'Trojan Females' achieves population suppression in fruit flies

    Pests are a global threat to biodiversity, ecosystem function, and human health. Pest control approaches are thus numerous, but their implementation costly, damaging to non-target species, and ineffective at low population densities. The Trojan Female Technique (TFT) is a prospective self-perpetuating control technique that is species-specific and predicted to be effective at low densities. The goal of the TFT is to harness naturally-occurring mutations in the mitochondrial genome that impair male fertility while having no effect on females. This study provides proof-of-concept for the TFT, by showing that introduction of a male fertility-impairing mtDNA haplotype into replicated populations of Drosophila melanogaster causes numerical population suppression, with the magnitude of effect positively correlated with its frequency at trial inception. Further development of the TFT could lead to establishing a control strategy that overcomes limitations of conventional approaches, with broad applicability to invertebrate and vertebrate species, to control environmental and economic pests (Wolff, 2017).

    Male mate choice via cuticular hydrocarbon pheromones drives reproductive isolation between Drosophila species

    Mate discrimination is a key mechanism restricting gene flow between species. While studied extensively with respect to female mate choice, mechanisms of male mate choice between species are far less studied. Thus, there is little knowledge of the relative frequency, importance, or overall contribution of male mate discrimination to reproductive isolation. This study estimated the relative contributions of male and female choice to reproductive isolation between Drosophila simulans and D. sechellia, and showed that male mate discrimination accounts for the majority of the current isolation between these species. It was further demonstrate that males discriminate based on female cuticular hydrocarbon pheromones, and collect evidence supporting the hypothesis that male mate discrimination may alleviate the costs associated with heterospecific courtship and mating. These findings highlight the potentially significant contribution of male mate choice to the formation of reproductive isolating barriers, and thus the speciation process (Shahandeh, 2017).

    The rate of evolution of postmating-prezygotic reproductive isolation in Drosophila

    Reproductive isolation is an intrinsic aspect of species formation. For that reason, the identification of the precise isolating traits, and the rates at which they evolve, is crucial to understanding how species originate and persist. Previous work has measured the rates of evolution of prezygotic and postzygotic barriers to gene flow, yet no systematic analysis has studied the rates of evolution of postmating-prezygotic (PMPZ) barriers. This study measured the magnitude of two barriers to gene flow that act after mating occurs but before fertilization. The magnitude of a premating barrier (female mating rate in nonchoice experiments) and two postzygotic barriers (hybrid inviability and hybrid sterility) was also measured for all pairwise crosses of all nine known extant species within the melanogaster subgroup. The results indicate that PMPZ isolation evolves faster than hybrid inviability but slower than premating isolation. Next, postzygotic isolation was partitioned into different components; as expected, hybrid sterility evolves faster than hybrid inviability. These results lend support for the hypothesis that, in Drosophila, reproductive isolation mechanisms (RIMs) that act early in reproduction (or in development) tend to evolve faster than those that act later in the reproductive cycle. Finally, whether there was evidence for reinforcing selection at any RIM was tested. No evidence was found for generalized evolution of reproductive isolation via reinforcement which indicates that there is no pervasive evidence of this evolutionary process. These results indicate that PMPZ RIMs might have important evolutionary consequences in initiating speciation and in the persistence of new species (Turissini, 2017).

    Genomic changes following the reversal of a Y chromosome to an autosome in Drosophila pseudoobscura

    Robertsonian translocations resulting in fusions between sex chromosomes and autosomes shape Karyotype evolution by creating new sex chromosomes from autosomes. These translocations can also reverse sex chromosomes back into autosomes, which is especially intriguing given the dramatic differences between autosomes and sex chromosomes. To study the genomic events following a Y chromosome reversal, this study investigated an autosome-Y translocation in Drosophila pseudoobscura. The ancestral Y chromosome fused to a small autosome (the dot chromosome) approximately 10-15 Mya. Single molecule real-time sequencing reads were used to assemble the D. pseudoobscura dot chromosome, including the Y-to-dot translocation. It was found that the intervening sequence between the ancestral Y and the rest of the dot chromosome is only ∼78 Kb and is not repeat-dense, suggesting that the centromere now falls outside, rather than between, the fused chromosomes. The Y-to-dot region is 100 times smaller than the D. melanogaster Y chromosome, owing to changes in repeat landscape. However, a consistent reduction in intron sizes across the Y-to-dot region was not found. Instead, deletions in intergenic regions and possibly a small ancestral Y chromosome size may explain the compact size of the Y-to-dot translocation (Chang, 2017).

    Linked genetic variation and not genome structure causes widespread differential expression associated with chromosomal inversions

    Chromosomal inversions are widely thought to be favored by natural selection because they suppress recombination between alleles that have higher fitness on the same genetic background or in similar environments. Nonetheless, few selected alleles have been characterized at the molecular level. Gene expression profiling provides a powerful way to identify functionally important variation associated with inversions and suggests candidate phenotypes. This study characterized differential expression patterns associated with two chromosomal inversions found in natural Drosophila melanogaster populations. To isolate the impacts of genome structure, synthetic chromosomal inversions were engineered on controlled genetic backgrounds with breakpoints that closely match each natural inversion. Synthetic inversions have negligible effects on gene expression. Nonetheless, natural inversions have broad-reaching regulatory impacts in cis and trans Furthermore, differentially expressed genes associated with both natural inversions are enriched for loci associated with immune response to bacterial pathogens. The results support the idea that inversions in D. melanogaster experience natural selection to maintain associations between functionally related alleles to produce complex phenotypic outcomes (Said, 2018).

    An investigation of Y chromosome incorporations in 400 species of Drosophila and related genera

    Y chromosomes are widely believed to evolve from a normal autosome through a process of massive gene loss (with preservation of some male genes), shaped by sex-antagonistic selection and complemented by occasional gains of male-related genes. The net result of these processes is a male-specialized chromosome. This might be expected to be an irreversible process, but it was found in 2005 that the Drosophila pseudoobscura Y chromosome was incorporated into an autosome. Y chromosome incorporations have important consequences: a formerly male-restricted chromosome reverts to autosomal inheritance, and the species may shift from an XY/XX to X0/XX sex-chromosome system. In order to assess the frequency and causes of this phenomenon Y chromosome incorporations was sought in 400 species from Drosophila and related genera. One additional large scale event of Y chromosome incorporation was found, affecting the whole montium subgroup (40 species in our sample); overall 13% of the sampled species (52/400) have Y incorporations. While previous data indicated that after the Y incorporation the ancestral Y disappeared as a free chromosome, the much larger data set analyzed here indicates that a copy of the Y survived as a free chromosome both in montium and pseudoobscura species, and that the current Y of the pseudoobscura lineage results from a fusion between this free Y and the neoY. The 400 species sample also showed that the previously suggested causal connection between X-autosome fusions and Y incorporations is, at best, weak: the new case of Y incorporation (montium) does not have X-autosome fusion, whereas nine independent cases of X-autosome fusions were not followed by Y incorporations. Y incorporation is an underappreciated mechanism affecting Y chromosome evolution; these results show that at least in Drosophila it plays a relevant role and highlight the need of similar studies in other groups (Dupim, 2018).

    A hidden markov model approach for simultaneously estimating local ancestry and admixture time using next generation sequence data in samples of arbitrary ploidy

    Admixture-the mixing of genomes from divergent populations-is increasingly appreciated as a central process in evolution. This study introduced a novel hidden Markov model for estimating local ancestry that models the read pileup data, rather than genotypes, is generalized to arbitrary ploidy, and can estimate the time since admixture during local ancestry inference. This method can simultaneously estimate the time since admixture and local ancestry with good accuracy, and it performs well on samples of high ploidy-i.e. 100 or more chromosomes. As this method is very general, it will be useful for local ancestry inference in a wider variety of populations than what previously has been possible. The method was applied to pooled sequencing data derived from populations of Drosophila melanogaster on an ancestry cline on the east coast of North America. Regions of local recombination rates were found to be negatively correlated with the proportion of African ancestry, suggesting that selection against foreign ancestry is the least efficient in low recombination regions. Finally it was shown that clinal outlier loci are enriched for genes associated with gene regulatory functions, consistent with a role of regulatory evolution in ecological adaptation of admixed D. melanogaster populations. These results illustrate the potential of local ancestry inference for elucidating fundamental evolutionary processes (Corbett-Detig, 2017).

    Symbiont strain is the main determinant of variation in Wolbachia-mediated protection against viruses across Drosophila species

    Wolbachia is a common heritable bacterial symbiont in insects. Its evolutionary success lies in the diverse phenotypic effects it has on its hosts coupled to its propensity to move between host species over evolutionary timescales. In a survey of natural host-symbiont associations in a range of Drosophila species, this study found that 10 of 16 Wolbachia strains protected their hosts against viral infection. By moving Wolbachia strains between host species, the symbiont genome was found to have a much greater influence on the level of antiviral protection than the host genome. The reason for this was that the level of protection depended on the density of the symbiont in host tissues, and Wolbachia rather than the host controlled density. The finding that virus resistance and symbiont density are largely under the control of symbiont genes in this system has important implications both for the evolution of these traits and for public health programs using Wolbachia to prevent mosquitoes from transmitting disease (Martinez, 2017).

    Inter- and intra-species variation in genome-wide gene expression of Drosophila in response to parasitoid wasp attack

    A study of inter- and intra-species variation in resistance to parasitoid attack used RNA-seq after parasitization in lines experimentally selected for increased resistance. A core set of genes was found that are consistently up-regulated after parasitoid attack. Another set showed no up-regulation or expression in D. sechellia, the species unable to raise an immune response against parasitoids. This set consists largely of genes that are lineage-restricted to the melanogaster subgroup. Artificially selected lines did not show significant differences in gene expression with respect to non-selected lines, but several genes showed differential exon usage. This study has shown substantial similarities, but also notable differences, in the transcriptional responses to parasitoid attack among four closely related Drosophila species. In contrast, within D. melanogaster, the responses were remarkably similar. It was confirmed that in the short-term, selection does not act on a pre-activation of the immune response. Instead it may target alternative mechanisms such as differential exon usage. In the long-term, support was found for the hypothesis that the ability to immunologically resist parasitoid attack is contingent on new genes that are restricted to the melanogaster subgroup (Salazar-Jaramillo, 2017).

    Geographic variation in the spotted-wing drosophila, Drosophila suzukii (Diptera: Drosophilidae), based on mitochondrial DNA sequences

    The spotted-wing drosophila (SWD), Drosophila suzukii, is an economically damaging pest that was originally native to a few Asian countries, including Korea, but is now found in North America and Europe. Portions of the mitochondrial (mt) COI and ND4 genes were sequenced from a total of 195 individuals collected mainly from Korea. GenBank-registered COI sequences were combined with the COI data to assess the worldwide diversity, divergence, and relatedness of SWD haplotypes. A total of 139 haplotypes were obtained from the concatenated COI and ND4 sequences. Most haplotypes were confined to single localities, but 12 of them were found in more than two localities, and one haplotype (SWDCN61) was found from Korea to Canada. A dataset combining GenBank sequences with the current data identified a total of 94 worldwide COI haplotypes with a maximum sequence divergence (MSD) of 5.433% (32 bp). Although most haplotypes were found in only a single country, a few haplotypes were found commonly in China, Korea, and Japan; these occurred at a higher frequency and were often involved in introductions. A rough estimate of genetic diversity in each country showed higher diversity in ancestral distributional ranges, but the invasion over Asian countries seems to have been substantial because haplotype diversity was only 2.35 to 3.97-fold lower in the U.S.A, Canada, and Italy than that in the populations' ancestral ranges (Choi, 2017).

    Rapidly evolving Toll-3/4 genes encode male-specific Toll-like receptors in Drosophila

    Animal Toll-like receptors (TLRs) have evolved through a pattern of duplication and divergence. Whereas mammalian TLRs directly recognize microbial ligands, Drosophila Tolls bind endogenous ligands downstream of both developmental and immune signaling cascades. This study found that most Toll genes in Drosophila evolve slowly with little gene turnover (gains/losses), consistent with their important roles in development and indirect roles in microbial recognition. In contrast, the Toll-3/4 genes were found to have experienced an unusually rapid rate of gene gains and losses, resulting in lineage-specific Toll-3/4s and vastly different gene repertoires among Drosophila species, from zero copies (e.g., D. mojavensis) to nineteen copies (e.g., D. willistoni). In D. willistoni, strong evidence was found for positive selection in Toll-3/4 genes, localized specifically to an extracellular region predicted to overlap with the binding site of Spatzle, the only known ligand of insect Tolls. However, because Spatzle genes are not experiencing similar selective pressures, it was hypothesize that Toll-3/4s may be rapidly evolving because they bind to a different ligand, akin to TLRs outside of insects. Unlike other Toll genes in D. melanogaster, Toll-3 and Toll-4 have apparently escaped from essential developmental roles. It is proposed that the Toll-3/4 genes represent an exceptionally rapidly evolving lineage of Drosophila Toll genes, which play an unusual, as-yet-undiscovered role in the male germline (Levin, 2017).

    A fitness trade-off between seasons causes multigenerational cycles in phenotype and population size

    Although seasonality is widespread and can cause fluctuations in the intensity and direction of natural selection, there is little information about the consequences of seasonal fitness trade-offs for population dynamics. This study exposed populations of Drosophila melanogaster to repeated seasonal changes in resources across 58 generations and used experimental and mathematical approaches to investigate how viability selection on body size in the non-breeding season could affect demography. Opposing seasonal episodes of natural selection on body size were shown to interact with both direct and delayed density dependence to cause populations to undergo predictable multigenerational density cycles. These results provide evidence that seasonality can set the conditions for life-history trade-offs and density dependence, which can, in turn, interact to cause multigenerational population cycles (Betini, 2017).

    Manipulation of feeding regime alters sexual dimorphism for lifespan and reduces sexual conflict in Drosophila melanogaster

    Sexual dimorphism for lifespan (SDL) is widespread, but poorly understood. A leading hypothesis, which was tested in this study, is that strong SDL can reduce sexual conflict by allowing each sex to maximize its sex-specific fitness. Replicated experimental evolution lines of the fruit fly, Drosophila melanogaster, were used that had been maintained for over 360 generations on either unpredictable 'Random' or predictable 'Regular' feeding regimes. This evolutionary manipulation of feeding regime led to robust, enhanced SDL in Random over control, Regular lines. Enhanced SDL was associated with a significant increase in the fitness of focal males, tested with wild-type (WT) females. This was due to sex-specific changes to male life history, manifested as increased early reproductive output and reduced survival. In contrast, focal female fitness, tested with WT males, did not differ across regimes. Hence increased SDL was associated with a reduction in sexual conflict, which increased male fitness and maintained fitness in females. Differences in SDL were not associated with developmental time or developmental survival. Overall, the results showed that the expression of enhanced SDL, resulting from experimental evolution of feeding regimes, was associated with male-specific changes in life history, leading to increased fitness and reduced sexual conflict (Duxbury, 2017).

    Fear creates an Allee effect: experimental evidence from seasonal populations

    Allee effects, a decline in individual fitness at low population size or density, driven by predation can play a strong role in the decline of small populations but are conventionally thought to occur when generalist predators target specific prey (i.e. type II functional response). However, aside from direct consumption, fear of predators could also increase vigilance and reduce time spent foraging as population size decreases, as has been observed in wild mammals living in social groups. To investigate the role of fear on fitness in relation to population density in a species with limited sociality, varying densities of Drosophila melanogaster were exposed to mantid predators either during an experimental breeding season or non-breeding season. The presence of mantids in either season decreased the reproductive performance of individuals but only at low breeding densities, providing evidence for an Allee effect. The experimental results were used to parametrize a mathematical model to examine the population consequences of fear at low densities. Fear tended to destabilize population dynamics and increase the risk of extinction up to sevenfold. The study provides unique experimental evidence that the indirect effects of the presence of predators can cause an Allee effect and has important consequences for understanding of the dynamics of small populations (Elliott, 2017).

    Pleiotropic effects of DDT resistance on male size and behaviour

    Understanding the evolution and spread of insecticide resistance requires knowing the relative fitness of resistant organisms. In the absence of insecticides, resistance is predicted to be costly. The Drosophila melanogaster DDT resistance allele (DDT-R) is associated with a male mating cost. This could be because resistant males are generally smaller, but DDT-R may also alter courtship behaviours. This study tested for body size and courtship effects of DDT-R on mating success in competitive and non-competitive mating trials respectively. Relative aggression was assessed in resistant and susceptible males because aggression can also influence mating success. While the effect of DDT-R on male size partly contributed to reduced mating success, resistant males also had lower rates of courtship and were less aggressive than susceptible males. These differences contribute to the observed DDT-R mating costs. Additionally, these pleiotropic effects of DDT-R are consistent with the history and spread of resistance alleles in nature (Rostant, 2017).

    Perceptive costs of reproduction drive ageing and physiology in male Drosophila

    Costs of reproduction are thought to result from natural selection optimizing organismal fitness within putative physiological constraints. Phenotypic and population genetic studies of reproductive costs are plentiful across taxa, but an understanding of their mechanistic basis would provide important insight into the diversity in life-history traits, including reproductive effort and ageing. This study dissected the causes and consequences of specific costs of reproduction in male Drosophila melanogaster. Key survival and physiological costs of reproduction arise from perception of the opposite sex, and they are reversed by the act of mating. In the absence of pheromone perception, males are free from reproductive costs on longevity, stress resistance and fat storage. The costs of perception and the benefits of mating are both mediated by evolutionarily conserved neuropeptidergic signalling molecules, as well as the transcription factor dFoxo. These results provide a molecular framework in which certain costs of reproduction arise as a result of self-imposed 'decisions' in response to perceptive neural circuits, which then orchestrate the control of life-history traits independently of physical or energetic effects associated with mating itself (Harvanek, 2017).

    One prophage WO gene rescues cytoplasmic incompatibility in Drosophila melanogaster

    Wolbachia are maternally inherited, intracellular bacteria at the forefront of vector control efforts to curb arbovirus transmission. In international field trials, the cytoplasmic incompatibility (CI) drive system of wMel Wolbachia is deployed to replace target vector populations, whereby a Wolbachia-induced modification of the sperm genome kills embryos. However, Wolbachia in the embryo rescue the sperm genome impairment, and therefore CI results in a strong fitness advantage for infected females that transmit the bacteria to offspring. The two genes responsible for the wMel-induced sperm modification of CI, cifA and cifB, were recently identified in the eukaryotic association module of prophage WO, but the genetic basis of rescue is unresolved. This study used transgenic and cytological approaches to demonstrate that maternal cifA expression independently rescues CI and nullifies embryonic death caused by wMel Wolbachia in Drosophila melanogaster. Discovery of cifA as the rescue gene and previously one of two CI induction genes establishes a "Two-by-One" model that underpins the genetic basis of CI. Results highlight the central role of prophage WO in shaping Wolbachia phenotypes that are significant to arthropod evolution and vector control (Shropshire, 2018).

    Repeated duplication of Argonaute2 is associated with strong selection and testis specialization in Drosophila

    Argonaute2 (Ago2) is a rapidly evolving nuclease in the Drosophila melanogaster RNA interference (RNAi) pathway that targets viruses and transposable elements in somatic tissues. This study reconstruct the history of Ago2 duplications across the Drosophila obscura group, and patterns of gene expression were used to infer new functional specialization. Some duplications were shown to be old, shared by the entire species group, and losses may be common, including previously undetected losses in the lineage leading to D. pseudoobscura. While the original (syntenic) gene copy has generally retained the ancestral ubiquitous expression pattern, most of the novel Ago2 paralogues have independently specialized to testis-specific expression. Using population genetic analyses, it was shown that most testis-specific paralogues have significantly lower genetic diversity than the genome-wide average. This suggests recent positive selection in three different species, and model-based analyses provide strong evidence of recent hard selective sweeps in or near four of the six D. pseudoobscura Ago2 paralogues. It is speculated that the repeated evolution of testis-specificity in obscura group Ago2 genes, combined with their dynamic turnover and strong signatures of adaptive evolution, may be associated with highly derived roles in the suppression of transposable elements or meiotic drive. This study highlights the lability of RNAi pathways, even within well-studied groups such as Drosophila, and suggests that strong selection may act quickly after duplication in RNAi pathways, potentially giving rise to new and unknown RNAi functions in non-model species (Lewis, 2016).

    Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

    Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster, however only the Cyp6g1 locus has previously been shown to be relevant to field populations. This study used genome-wide association studies to identify DDT-associated polygenes and used selective sweep analyses to assess their adaptive significance. Two candidate DDT resistance loci were identified and verified. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection (Schmidt, 2017).

    Structural variants and selective sweep foci contribute to insecticide resistance in the Drosophila genetic reference panel

    Patterns of nucleotide polymorphism within populations of Drosophila melanogaster suggest that insecticides have been the selective agents driving the strongest recent bouts of positive selection. However, there is a need to explicitly link selective sweeps to the particular insecticide phenotypes that could plausibly account for the drastic selective responses that are observed in these non-target insects. This study screened the Drosophila Genetic Reference Panel with two common insecticides; malathion (an organophosphate) and permethrin (a pyrethroid). Genome-wide association studies map survival on malathion to the two of the largest sweeps in the D. melanogaster genome; Ace and Cyp6g1. Malathion survivorship also correlates with lines which have high levels of Cyp12d1, Jheh1 and Jheh2 transcript abundance. Permethrin phenotypes map to the largest cluster of P450 genes in the Drosophila genome, however in contrast to a selective sweep driven by insecticide use, the derived allele seems to be associated with susceptibility. These results underscore previous findings that highlight the importance of structural variation to insecticide phenotypes: Cyp6g1 exhibits copy number variation and transposable element insertions, Cyp12d1 is tandemly duplicated, the Jheh loci are associated with a Bari1 transposable element insertion, and a Cyp6a17 deletion is associated with susceptibility (Battlay, 2018).

    Directional selection reduces developmental canalization against genetic and environmental perturbations in Drosophila wings

    Natural selection may enhance or weaken the robustness of phenotypes against genetic or environmental perturbations. However, important aspects of the relationship between adaptive evolution and canalization remain unclear. Recent work showed that the evolution of larger wing size in a high altitude natural population of Drosophila melanogaster was accompanied by decanalized wing development--specifically a loss of robustness to genetic perturbation. But that study did not address environmental robustness, and it compared populations that may have numerous biological differences. This study performed artificial selection on this same trait in D. melanogaster (larger wing length) and directly test whether this directional selection resulted in decanalization. In general, size-selected replicates show greater frequencies of wing defects than control replicates both after mutagenesis (genetic perturbation) and when subjected to high temperature stress (environmental perturbation), although the increase in defect frequency varies importantly among replicates. These results support the hypothesis that directional selection may result in the loss of both genetic and environmental robustness-offering a rare window into the relationship between adaptation and canalization (Groth, 2018).

    Genetic trade-offs between male reproductive traits in Drosophila melanogaster

    In Drosophila melanogaster, males engage in both extensive pre- and post-copulatory competition for the opportunity to mate with females and subsequently sire offspring. The selection pressure for increased male reproductive success has resulted in the evolution of a wide diversity of sexual traits. However, despite strong selection, individuals often exhibit considerable phenotypic variation in the expression of these traits, and it is unclear if any of this variation is owing to underlying genetic trade-offs. Using hemiclonal flies this study examined how male reproductive success covaries with their ability to induce long-term stimulation of oogenesis and oviposition in their mates, and how this relationship may change over time. It was found that males from hemiclone lines with phenotypes that were more successful in a short-term reproductive 'scramble' environment were less effective at stimulating long-term fecundity in females. Furthermore, it was observed that males from hemiclone lines which showed the most improvement over a longer reproductive interaction period also tended to stimulate higher long-term fecundity in females. Together, these results indicate the presence of genetic trade-offs between different male reproductive traits and offer insights into the maintenance of their variation (Filice, 2018).

    X chromosome drive in a widespread Palearctic woodland fly, Drosophila testacea

    Selfish genes that bias their own transmission during meiosis can spread rapidly in populations, even if they contribute negatively to the fitness of their host. Driving X chromosomes provide a clear example of this type of selfish propagation. These chromosomes have important evolutionary and ecological consequences, and can be found in a broad range of taxa including plants, mammals, and insects. This study reports a new case of X chromosome drive (X drive) in a widespread woodland fly, Drosophila testacea. Males carrying the driving X (SR males) sire 80-100% female offspring, and possess a diagnostic X chromosome haplotype that is perfectly associated with the sex ratio distortion phenotype. The majority of sons produced by SR males are sterile and appear to lack a Y chromosome, suggesting that meiotic defects involving the Y chromosome may underlie X drive in this species. Abnormalities in sperm cysts of SR males reflect that some spermatids are failing to develop properly, confirming that drive is acting during gametogenesis. By screening wild-caught flies using progeny sex ratios and a diagnostic marker, it was demonstrated that the driving X is present in wild populations at a frequency of ~10% and that suppressors of drive are segregating in the same population. The testacea species group appears to be a hotspot for X drive, and D. testacea is a promising model to compare driving X chromosomes in closely related species, some of which may even be younger than the chromosomes themselves (Keais, 2017).

    A Pooled Sequencing Approach Identifies a Candidate Meiotic Driver in Drosophila

    Meiotic drive occurs when a selfish element increases its transmission frequency above the Mendelian ratio by hijacking the asymmetric divisions of female meiosis. New methods to reliably detect meiotic drive are therefore needed, particularly for discovering moderate-strength drivers that are likely to be more prevalent in natural populations than strong drivers. This study reports an efficient method that uses sequencing of large pools of backcross (BC1) progeny to test for deviations from Mendelian segregation genome-wide of single-nucleotide polymorphisms (SNPs) that distinguish the parental strains. Meiotic drive can be detected by a characteristic pattern of decay in distortion of SNP frequencies, caused by recombination unlinking the driver from distal loci. Control crosses allow allele-frequency distortion caused by meiotic drive to be distinguished from distortion resulting from developmental effects. This approach was used to test whether chromosomes with extreme telomere-length differences segregate at Mendelian ratios, as telomeric regions are a potential hotspot for meiotic drive due to their roles in meiotic segregation and multiple observations of high rates of telomere sequence evolution. Using four different pairings of long and short telomere strains, this study found no evidence that extreme telomere-length variation causes meiotic drive in Drosophila. However, one candidate meiotic driver was identified in a centromere-linked region that shows an ~8% increase in transmission frequency, corresponding to a ~54:46 segregation ratio. These results show that candidate meiotic drivers of moderate strength can be readily detected and localized in pools of F1 progeny (Wei, 2017).

    TP53 copy number expansion is associated with the evolution of increased body size and an enhanced DNA damage response in elephants
    A major constraint on the evolution of large body sizes in animals is an increased risk of developing cancer. There is no correlation, however, between body size and cancer risk. This lack of correlation is often referred to as 'Peto's Paradox'. This study showed that the elephant genome encodes 20 copies of the tumor suppressor gene TP53 (see Drosophila p53) and that the increase in TP53 copy number occurred coincident with the evolution of large body sizes, the evolution of extreme sensitivity to genotoxic stress, and a hyperactive TP53 signaling pathway in the elephant (Proboscidean) lineage. Furthermore several of the TP53 retrogenes (TP53RTGs) were shown to be transcribed and likely translated. While TP53RTGs do not appear to directly function as transcription factors, they do contribute to the enhanced sensitivity of elephant cells to DNA damage and the induction of apoptosis by regulating activity of the TP53 signaling pathway. These results suggest that an increase in the copy number of TP53 may have played a direct role in the evolution of very large body sizes and the resolution of Peto's paradox in Proboscideans (Sulak, 2016).

    Gene duplications circumvent trade-offs in enzyme function: Insect adaptation to toxic host plants

    Herbivorous insects and their adaptations against plant toxins provide striking opportunities to investigate the genetic basis of traits involved in coevolutionary interactions. Target site insensitivity to cardenolides has evolved convergently across six orders of insects, involving identical substitutions in the Na,K-ATPase gene and repeated convergent gene duplications. The large milkweed bug, Oncopeltus fasciatus, has three copies of the Na,K-ATPase alpha-subunit gene that bear differing numbers of amino acid substitutions in the binding pocket for cardenolides. To analyze the effect of these substitutions on cardenolide resistance and to infer possible trade-offs in gene function, the cardenolide-sensitive Na,K-ATPase of Drosophila melanogaster was expressed in vitro and four distinct combinations were introduced of substitutions observed in the three gene copies of O. fasciatus. With an increasing number of substitutions, the sensitivity of the Na,K-ATPase to a standard cardenolide decreased in a stepwise manner. At the same time, the enzyme's overall activity decreased significantly with increasing cardenolide resistance and only the least substituted mimic of the Na,K-ATPase α1C copy maintained activity similar to the wild-type enzyme. These results suggest that the Na,K-ATPase copies in O. fasciatus have diverged in function, enabling specific adaptations to dietary cardenolides while maintaining the functionality of this critical ion carrier (Dalla, 2016).

    Rapid functional and sequence differentiation of a tandemly-repeated species-specific multigene family in Drosophila

    Gene clusters of recently duplicated genes are hotbeds for evolutionary change. However, understanding of how mutational mechanisms and evolutionary forces shape the structural and functional evolution of these clusters is hindered by the high sequence identity among the copies, which typically results in their inaccurate representation in genome assemblies. The presumed testis-specific, chimeric gene Sdic originated and tandemly expanded in Drosophila melanogaster, contributing to increased male-male xion. Using various types of massively parallel sequencing data, the organization, sequence evolution, and functional attributes of the different Sdic copies were examined. By leveraging long-read sequencing data, both copy number and order differences were uncovered from the currently accepted annotation for the Sdic region. Despite evidence for pervasive gene conversion affecting the Sdic copies, signatures of two episodes of diversifying selection were uncovered, that have contributed to the evolution of a variety of C-termini and miRNA binding site compositions. Expression analyses involving RNA-seq datasets from 59 different biological conditions revealed distinctive expression breadths among the copies, with three copies being transcribed in females, opening the possibility to a sexually antagonistic effect. Phenotypic assays using Sdic knock-out strains indicated that should this antagonistic effect exist, it does not compromise female fertility. These results strongly suggest that the genome consolidation of the Sdic gene cluster is more the result of a quick exploration of different paths of molecular tinkering by different copies than a mere dosage increase, which could be a recurrent evolutionary outcome in the presence of persistent sexual selection (Clifton, 2016).

    Transcriptional interference promotes rapid expression divergence of Drosophila nested genes

    Nested genes are the most common form of protein-coding overlap in eukaryotic genomes. Previous studies have shown that nested genes accumulate rapidly over evolutionary time, typically via the insertion of short young duplicate genes into long introns. However, the evolutionary relationship between nested genes remains unclear. This study compare RNA-seq expression profiles of nested, proximal intra-chromosomal, intermediate intra-chromosomal, distant intra-chromosomal, and inter-chromosomal gene pairs in two Drosophila species. Expression profiles of nested genes were found to be more divergent than those of any other class of genes, supporting the hypothesis that concurrent expression of nested genes is deleterious due to transcriptional interference. Further analysis reveals that expression profiles of derived nested genes are more divergent than those of their ancestral un-nested orthologs, which are more divergent than those of un-nested genes with similar genomic features. Thus, gene expression divergence between nested genes is likely caused by selection against nesting of genes with insufficiently divergent expression profiles, as well as by continued expression divergence after nesting. Moreover, expression divergence and sequence evolutionary rates are elevated in young nested genes and reduced in old nested genes, indicating that a burst of rapid evolution occurs after nesting. Together, these findings suggest that similarity between expression profiles of nested genes is deleterious due to transcriptional interference, and that natural selection addresses this problem both by eradicating highly deleterious nestings and by enabling rapid expression divergence of surviving nested genes, thereby quickly limiting or abolishing transcriptional interference (Assis, 2016).

    The goddard and saturn genes are essential for Drosophila male fertility and may have arisen de novo

    New genes arise through a variety of mechanisms, including the duplication of existing genes and the de novo birth of genes from non-coding DNA sequences. While there are numerous examples of duplicated genes with important functional roles, the functions of de novo genes remain largely unexplored. Many newly evolved genes are expressed in the male reproductive tract, suggesting that these evolutionary innovations may provide advantages to males experiencing sexual selection. Using testis-specific RNA interference, 11 putative de novo genes in Drosophila melanogaster for effects on male fertility, and two, goddard and saturn, were identified that are essential for spermatogenesis and sperm function. goddard knockdown males fail to produce mature sperm, while saturn knockdown males produce fewer sperm that function inefficiently once transferred to females. Consistent with a de novo origin, both genes are identifiable only in Drosophila and are predicted to encode proteins with no sequence similarity to any annotated protein. However, since high levels of divergence prevented the unambiguous identification of the non-coding sequences from which each gene arose, saturn and saturn are considered to be putative de novo genes. Within Drosophila, both genes have been lost in certain lineages, but show conserved, male-specific patterns of expression in the species in which they are found. Goddard is consistently found in single-copy and evolves under purifying selection. In contrast, saturn has diversified through gene duplication and positive selection. These data suggest that de novo genes can evolve essential roles in male reproduction (Gubala, 2017).

    Recurrent gene duplication leads to diverse repertoires of centromeric histones in Drosophila species

    Despite their essential role in the process of chromosome segregation in most eukaryotes, centromeric histones show remarkable evolutionary lability. Not only have they been lost in multiple insect lineages, but they have also undergone gene duplication in multiple plant lineages. Based on detailed study of a handful of model organisms including Drosophila melanogaster, centromeric histone duplication is considered to be rare in animals. Using a detailed phylogenomic study, this study found that Cid, the centromeric histone gene, has undergone at least four independent gene duplications during Drosophila evolution. Duplicate Cid genes were found in D. eugracilis (Cid2), in the montium species subgroup (Cid3, Cid4) and in the entire Drosophila subgenus (Cid5). Cid3, Cid4, Cid5 all localize to centromeres in their respective species. Some Cid duplicates are primarily expressed in the male germline. With rare exceptions, Cid duplicates have been strictly retained after birth, suggesting that they perform non-redundant centromeric functions, independent from the ancestral Cid. Indeed, each duplicate encodes a distinct N-terminal tail, which may provide the basis for distinct protein-protein interactions. Finally, it was shown some Cid duplicates evolve under positive selection whereas others do not. Taken together, these results support the hypothesis that Drosophila Cid duplicates have subfunctionalized. Thus, these gene duplications provide an unprecedented opportunity to dissect the multiple roles of centromeric histones (Kursel, 2017).

    Tandem duplications lead to novel expression patterns through exon shuffling in Drosophila yakuba

    One common hypothesis to explain the impacts of tandem duplications is that whole gene duplications commonly produce additive changes in gene expression due to copy number changes. This study used genome wide RNA-seq data from a population sample of Drosophila yakuba to test this 'gene dosage' hypothesis. Little evidence was observed of expression changes in response to whole transcript duplication capturing 5' and 3' UTRs. Among whole gene duplications, evidence was observed that dosage sharing across copies is likely to be common. The lack of expression changes after whole gene duplication suggests that the majority of genes are subject to tight regulatory control and therefore not sensitive to changes in gene copy number. Rather, changes were observed in expression level due to both shuffling of regulatory elements and the creation of chimeric structures via tandem duplication. Additionally, 30 de novo gene structures were observed arising from tandem duplications, 23 of which form with expression in the testes. Thus, the value of tandem duplications is likely to be more intricate than simple changes in gene dosage. The common regulatory effects from chimeric gene formation after tandem duplication may explain their contribution to genome evolution (Rogers, 2017).

    Diverse cis-regulatory mechanisms contribute to expression evolution of tandem gene duplicates

    Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. This study directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 million years ago. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, they were shown to diverge from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. This analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. It is argued that these findings will encourage studying cis-regulatory evolution at the whole-locus level in order to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression (Baudouin-Gonzalez, 2017).

    Genes relocated between Drosophila chromosome arms evolve under relaxed selective constraints relative to non-relocated genes

    Gene duplication creates a second copy of a gene either in tandem to the ancestral locus or dispersed to another chromosomal location. Gene relocations may be as common as canonical dispersed duplications in which both the ancestral and derived copies are retained. Relocated genes appear to be under more selective constraints than the derived copies of canonical duplications, and they are possibly as conserved as single-copy non-relocated genes. To test this hypothesis, comparative genomics, population genetics, gene expression, and functional analyses were combined to assess the selection pressures acting on relocated, duplicated, and non-relocated single-copy genes in Drosophila genomes. Relocated genes were found to evolve faster than single-copy non-relocated genes, and there is no evidence that this faster evolution is driven by positive selection. In addition, relocated genes are less essential for viability and male fertility than single-copy non-relocated genes, suggesting that relocated genes evolve fast because of relaxed selective constraints. However, relocated genes evolve slower than the derived copies of canonical dispersed duplicated genes. We therefore conclude that relocated genes are under more selective constraints than canonical duplicates, but are not as conserved as single-copy non-relocated genes (Hart, 2018).

    Origin, composition, and structure of the supernumerary B chromosome of Drosophila melanogaster

    The number of chromosomes carried by an individual species is one of its defining characteristics. Some species, however, can also carry supernumerary chromosomes referred to as B chromosomes. B chromosomes were recently identified in a laboratory stock of Drosophila melanogaster enabling them to be subjected to extensive molecular analysis. The B chromosomes by pulsed-field gel electrophoresis and determined their composition through next-generation sequencing. Although these B chromosomes carry no known euchromatic sequence, they are rich in transposable elements and long arrays of short nucleotide repeats, the most abundant being the uncharacterized AAGAT satellite repeat. Fluorescent in-situ hybridization on metaphase chromosome spreads revealed this repeat is located on Chromosome 4, strongly suggesting the origin of the B chromosomes is Chromosome 4. Cytological and quantitative comparisons of signal intensity between Chromosome 4 and the B chromosomes supports the hypothesis that the structure of the B chromosome is an isochromosome. Also, the identification is reported of a new B chromosome variant in a related laboratory stock. This B chromosome has a similar repeat signature as the original but is smaller and much less prevalent. Additional stocks with similar genotypes were examined and B chromosomes were found, but these stocks lacked the AAGAT satellite repeat. This molecular characterization of D. melanogaster B chromosomes is the first step towards understanding how supernumerary chromosomes arise from essential chromosomes and what may be necessary for their stable inheritance (Hanlon, 2018).

    Functional interplay between ribosomal protein paralogues in the eRpL22 family in Drosophila melanogaster

    Duplicated ribosomal protein (RP) genes in the Drosophila melanogaster eRpL22 family encode structurally-divergent and differentially-expressed rRNA-binding RPs. eRpL22 is expressed ubiquitously and eRpL22-like expression is tissue-restricted with highest levels in the adult male germline. Paralogue functional equivalence was explored using the GAL4-UAS system for paralogue knockdown or overexpression and a conditional eRpL22-like knockout in a heat- shock flippase/FRT line. Ubiquitous eRpL22 knockdown with Actin-GAL4 resulted in embryonic lethality, confirming eRpL22 essentiality. eRpL22-like knockdown (60%) was insufficient to cause lethality; yet, conditional eRpL22-like knockout at one hour following egg deposition caused lethality within each developmental stage. Therefore, each paralogue is essential. Variation in timing of heat-shock-induced eRpL22-like knockout highlighted early embryogenesis as the critical period where eRpL22-like expression (not compensated for by eRpL22) is required for normal development of several organ systems, including testis development and subsequent sperm production. To determine if eRpL22-like can substitute for eRpL22, Actin-GAL4 for ubiquitous eRpL22 knockdown and eRpL22-like-FLAG (or FLAG-eRpL22: control) overexpression. Emergence of adults demonstrated that ubiquitous eRpL22-like-FLAG or FLAG-eRpL22 expression eliminates embryonic lethality resulting from eRpL22 depletion. Adults rescued by eRpL22-like-FLAG (but not by FLAG-eRpL22) overexpression had reduced fertility and longevity. It is concluded that eRpL22 paralogue roles are not completely interchangeable and include functionally-diverse roles in development and spermatogenesis (Mageeney, 2018).

    Rapid Global Spread of wRi-like Wolbachia across Multiple Drosophila

    Maternally transmitted Wolbachia, Spiroplasma, and Cardinium bacteria are common in insects, but their interspecific spread is poorly understood. Because Wolbachia cannot survive outside host cells, spread between distantly related host species requires horizontal transfers that are presumably rare. This study documents spread of wRi-like Wolbachia among eight highly diverged Drosophila hosts (10-50 million years) over only about 14,000 years (5,000-27,000). Comparing 110 wRi-like genomes, it was found </=0.02% divergence from the wRi variant that spread rapidly through California populations of D. simulans. The hosts include both globally invasive species (D. simulans, D. suzukii, and D. ananassae) and narrowly distributed Australian endemics (D. anomalata and D. pandora). Phylogenetic analyses that include mtDNA genomes indicate introgressive transfer of wRi-like Wolbachia between closely related species D. ananassae, D. anomalata, and D. pandora but no horizontal transmission within species. These analyses suggest D. ananassae as the Wolbachia source for the recent wRi invasion of D. simulans and D. suzukii as the source of Wolbachia in its sister species D. subpulchrella. Although six of these wRi-like variants cause strong cytoplasmic incompatibility, two cause no detectable reproductive effects, indicating that pervasive mutualistic effects complement the reproductive manipulations for which Wolbachia are best known. "Super spreader" variants like wRi may be particularly useful for controlling insect pests and vector-borne diseases with Wolbachia transinfections (Turelli, 2018).

    Fine scale mapping of genomic introgressions within the Drosophila yakuba clade

    The process of speciation involves populations diverging over time until they are genetically and reproductively isolated. Hybridization between nascent species was long thought to directly oppose speciation. A natural place to look for individuals with admixed ancestry (indicative of introgression) is in regions where species co-occur. In west Africa, D. santomea and D. yakuba hybridize on the island of Sao Tome, while D. yakuba and D. teissieri hybridize on the nearby island of Bioko. This report quantifies the genomic extent of introgression between the three species of the Drosophila yakuba clade (D. yakuba, D. santomea), D. teissieri). The genomes of 86 individuals were sequenced from all three species. A new statistical framework was developed and applied, using a hidden Markov approach, to identify introgression. Introgression was found to have occurred between both species pairs but most introgressed segments are small (on the order of a few kilobases). After ruling out the retention of ancestral polymorphism as an explanation for these similar regions, this study found that the sizes of introgressed haplotypes indicate that genetic exchange is not recent (>1,000 generations ago). It was additionally shown that in both cases, introgression was rarer on X chromosomes than on autosomes which is consistent with sex chromosomes playing a large role in reproductive isolation. Even though the two species pairs have stable contemporary hybrid zones, providing the opportunity for ongoing gene flow, the results indicate that genetic exchange between these species is currently rare (Turissini, 2017).

    A Maladaptive Combination of Traits Contributes to the Maintenance of a Drosophila Hybrid Zone

    Drosophila teissieri and D. yakuba diverged approximately 3 mya and are thought to share a large, ancestral, African range. These species now co-occur in parts of continental Africa and in west Africa on the island of Bioko. While D. yakuba is a human commensal, D. teissieri seems to be associated with Parinari fruits, restricting its range to forests. Genome data indicate introgression, despite no evidence of contemporary hybridization. This study reports the discovery of D. yakuba-D. teissieri hybrids at the interface of secondary forests and disturbed, open habitats on Bioko. Hybrids are the F1 progeny of D. yakuba females and D. teissieri males. At high temperatures like those found on Bioko, D. teissieri females are generally less receptive to mating, and in combination with temperature effects on egg lay and egg-to-adult viability, this decreases the potential for gene flow between female D. teissieri and male D. yakuba relative to the reciprocal cross. Field and laboratory experiments demonstrate that F1 hybrids have a maladaptive combination of D. yakuba behavior and D. teissieri physiology, generating additional barriers to gene flow. Nevertheless, analysis of introgressed and non-introgressed regions of the genome indicate that, while rare, gene flow is relatively recent. These observations identify precise intrinsic and extrinsic factors that, along with hybrid male sterility, limit gene flow and maintain these species. These data contribute to a growing body of literature that suggests the Gulf of Guinea may be a hotspot for hybridization (Cooper, 2018).

    Comparative genomic analyses provide new insights into the evolutionary dynamics of heterochromatin in Drosophila

    The term heterochromatin has been long considered synonymous with gene silencing, but it is now clear that the presence of transcribed genes embedded in pericentromeric heterochromatin is a conserved feature in the evolution of eukaryotic genomes. Several studies have addressed the epigenetic changes that enable the expression of genes in pericentric heterochromatin, yet little is known about the evolutionary processes through which this has occurred. By combining genome annotation analysis and high-resolution cytology, this study has identified and mapped 53 orthologs of D. melanogaster heterochromatic genes in the genomes of two evolutionarily distant species, D. pseudoobscura and D. virilis. The results show that the orthologs of the D. melanogaster heterochromatic genes are clustered at three main genomic regions in D. virilis and D. pseudoobscura. In D. virilis, the clusters lie in the middle of euchromatin, while those in D. pseudoobscura are located in the proximal portion of the chromosome arms. Some orthologs map to the corresponding Muller C element in D. pseudoobscura and D. virilis, while others localize on the Muller B element, suggesting that chromosomal rearrangements that have been instrumental in the fusion of two separate elements involved the progenitors of genes currently located in D. melanogaster heterochromatin. These results demonstrate an evolutionary repositioning of gene clusters from ancestral locations in euchromatin to the pericentromeric heterochromatin of descendent D. melanogaster chromosomes. Remarkably, in both D. virilis and D. pseudoobscura the gene clusters show a conserved association with the HP1a protein, one of the most highly evolutionarily conserved epigenetic marks. In light of these results, a new scenario is suggested whereby ancestral HP1-like proteins (and possibly other epigenetic marks) may have contributed to the evolutionary repositioning of gene clusters into heterochromatin (Caizzi, 2016).

    Genomics of natural populations: How differentially expressed genes shape the evolution of chromosomal inversions in Drosophila pseudoobscura

    The third chromosome of Drosophila pseudoobscura is a model system to test hypotheses about how rearrangements are established in populations because its third chromosome is polymorphic for > 30 gene arrangements that were generated by a series of overlapping inversion mutations. Circumstantial evidence has suggested that these gene arrangements are selected. Despite the expected homogenizing effects of extensive gene flow, the frequencies of arrangements form gradients or clines in nature, which have been stable since the system was first described more than 80 years ago. Furthermore, multiple arrangements exist at appreciable frequencies across several ecological niches providing the opportunity for heterokaryotypes to form. This study tested whether genes are differentially expressed among chromosome arrangements in first instar larvae, adult females and males. In addition, it was asked whether transcriptional patterns in heterokaryotypes are dominant, semidominant, overdominant, or underdominant. Evidence was found for a significant abundance of differentially expressed genes across the inverted regions of the third chromosome, including an enrichment of genes involved in sensory perception for males. The majority of loci show additivity in heterokaryotypes. These results suggest that multiple genes have expression differences among arrangements that were either captured by the original inversion mutation or accumulated after it reached polymorphic frequencies, providing a potential source of genetic variation for selection to act upon. These data suggest that the inversions are favored because of their indirect effect of recombination suppression that has held different combinations of differentially expressed genes together in the various gene arrangement backgrounds (Fuller, 2016).

    The origin of chromosomal inversions as a source of segmental duplications in the Sophophora subgenus of Drosophila

    Chromosomal inversions can contribute to the adaptation of organisms to their environment by capturing particular advantageous allelic combinations of a set of genes included in the inverted fragment and also by advantageous functional changes due to the inversion process itself that might affect not only the expression of flanking genes but also their dose and structure. Of the two mechanisms originating inversions -ectopic recombination, and staggered double-strand breaks and subsequent repair- only the latter confers the inversion the potential to have dosage effects and/or to generate advantageous chimeric genes. In Drosophila subobscura, there is ample evidence for the adaptive character of its chromosomal polymorphism, with an important contribution of some warm-climate arrangements such as E1+2+9+12. This study has characterized the breakpoints of inversion E12 and established that it originated through the staggered-break mechanism like four of the five inversions of D. subobscura previously studied. This mechanism that also predominates in the D. melanogaster lineage might be prevalent in the Sophophora subgenus and contribute to the adaptive character of the polymorphic and fixed inversions of its species. Finally, the D. subobscura inversion breakpoint regions were shown to have generally been disrupted by additional structural changes that occurred at different time scales (Puerma, 2016).

    Life history evolution and cellular mechanisms associated with increased size in high-altitude Drosophila

    Understanding the physiological and genetic basis of growth and body size variation has wide-ranging implications, from cancer and metabolic disease to the genetics of complex traits. This study examined the evolution of body and wing size in high-altitude Drosophila melanogaster from Ethiopia, flies with larger size than any previously known population. Specifically, attempts were made to identify life history characteristics and cellular mechanisms that may have facilitated size evolution. The large-bodied Ethiopian flies laid significantly fewer but larger eggs relative to lowland, smaller-bodied Zambian flies. The highland flies were found to achieve larger size in a similar developmental period, potentially aided by a reproductive strategy favoring greater provisioning of fewer offspring. At the cellular level, cell proliferation was a strong contributor to wing size evolution, but both thorax and wing size increases involved important changes in cell size. Nuclear size measurements were consistent with elevated somatic ploidy as an important mechanism of body size evolution. The significance of these results for the genetic basis of evolutionary changes in body and wing size in Ethiopian D. melanogaster is discussed (Lack, 2016).

    Accelerated pseudogenization on the neo-X chromosome in Drosophila miranda

    Y chromosomes often degenerate via the accumulation of pseudogenes and transposable elements. By contrast, little is known about X-chromosome degeneration. This study compared the pseudogenization process between genes on the neo-sex chromosomes in Drosophila miranda and their autosomal orthologues in closely related species. The pseudogenization rate on the neo-X is much lower than the rate on the neo-Y, but appears to be higher than the rate on the orthologous autosome in D. pseudoobscura. Genes under less functional constraint and/or genes with male-biased expression tend to become pseudogenes on the neo-X, indicating the accumulation of slightly deleterious mutations and the feminization of the neo-X. A weak trend was found that the genes with female-benefit/male-detriment effects identified in D. melanogaster are pseudogenized on the neo-X, implying the masculinization of the neo-X. These observations suggest that both X and Y chromosomes can degenerate due to a complex suite of evolutionary forces (Nozawa, 2016).

    The landscape of A-to-I RNA editome is shaped by both positive and purifying selection

    The hydrolytic deamination of adenosine to inosine (A-to-I editing) in precursor mRNA induces variable gene products at the post-transcription level. How and to what extent A-to-I RNA editing diversifies transcriptome is not fully characterized in the evolution, and very little is known about the selective constraints that drive the evolution of RNA editing events. This study on A-to-I RNA editing, by generating a global profile of A-to-I editing for a phylogeny of seven Drosophila species, presents a model system spanning an evolutionary timeframe of approximately 45 million years. Of totally 9281 editing events identified, 5150 (55.5%) are located in the coding sequences (CDS) of 2734 genes. Phylogenetic analysis places these genes into 1,526 homologous families, about 5% of total gene families in the fly lineages. Based on conservation of the editing sites, the editing events in CDS are categorized into three distinct types, representing events on singleton genes (type I), and events not conserved (type II) or conserved (type III) within multi-gene families. While both type I and II events are subject to purifying selection, notably type III events are positively selected, and highly enriched in the components and functions of the nervous system. The tissue profiles are documented for three editing types, and their critical roles are further implicated by their shifting patterns during holometabolous development and in post-mating response. In conclusion, three A-to-I RNA editing types are found to have distinct evolutionary dynamics. It appears that nervous system functions are mainly tested to determine if an A-to-I editing is beneficial for an organism. The coding plasticity enabled by A-to-I editing creates a new class of binary variations, which is a superior alternative to maintain heterozygosity of expressed genes in a diploid mating system (Yu, 2016).

    Weak polygenic selection drives the rapid adaptation of the chemosensory system: lessons from the upstream regions of the major gene families

    The animal chemosensory system is involved in essential biological processes, most of them mediated by proteins encoded in multigene families. These multigene families have been fundamental for the adaptation to new environments, significantly contributing to phenotypic variation. This adaptive potential contrasts, however, with the lack of studies at their upstream regions, especially taking into account the evidence linking their transcriptional changes to certain phenotypic effects. This study explicitly characterised the contribution of the upstream sequences of the major chemosensory gene families to rapid adaptive processes. For that, the genome sequences of 158 lines were analyzed from a population of Drosophila melanogaster that recently colonised North America, and functional and transcriptional data available for this species were integrated. Both, strong negative and strong positive selection were found to shape transcriptional evolution at the genome-wide level. The chemosensory upstream regions, however, exhibit a distinctive adaptive landscape, including multiple mutations of small beneficial effect and a reduced number of cis-regulatory elements. Together, these results suggest that the promiscuous and partially redundant transcription and function of the chemosensory genes provide evolutionarily opportunities for rapid adaptive episodes through weak polygenic selection (Librado, 2016).

    Evolutionary dynamics of abundant stop codon readthrough

    This study leveraged comparative genomic evidence across 21 Anopheles mosquitoes to systematically annotate translational stop codon readthrough genes in the malaria vector Anopheles gambiae, and to provide the first study of abundant readthrough evolution, by comparison with 20 Drosophila species. Evolutionary signatures were identified of conserved, functional readthrough of 353 stop codons in the malaria vector, Anopheles gambiae, and of 51 additional Drosophila melanogaster stop codons. Most differences between the readthrough repertoires of the two species arose from readthrough gain or loss in existing genes, rather than birth of new genes or gene death; that readthrough-associated RNA structures are sometimes gained or lost while readthrough persists; that readthrough is more likely to be lost at TAA and TAG stop codons; and that readthrough is under continued purifying evolutionary selection in mosquito, based on population genetic evidence. Readthrough-associated gene properties were determined that predate readthrough, and differences were identified in the characteristic properties of readthrough genes between clades. More than 600 functional readthrough stop codons were identified in mosquito and 900 in fruit fly, provide evidence of readthrough control of peroxisomal targeting, and the phylogenetic extent of abundant readthrough as following divergence from centipede was refined (Jungreis, 2016).

    Molecular population genetics of the Polycomb genes in Drosophila subobscura
    Polycomb group (PcG) proteins are important regulatory factors that modulate the chromatin state. They form protein complexes that repress gene expression by the introduction of posttranslational histone modifications. The study of PcG proteins divergence in Drosophila revealed signals of coevolution among them and an acceleration of the nonsynonymous evolutionary rate in the lineage ancestral to the obscura group species, mainly in subunits of the Pcl-PRC2 complex. This study examined the nucleotide polymorphism of PcG genes in a natural population of D. subobscura to detect whether natural selection has also modulated the evolution of these important regulatory genes in a more recent time scale. Results show that most genes are under the action of purifying selection and present a level and pattern of polymorphism consistent with predictions of the neutral model, the exceptions being Su(z)12 and Pho. MK tests indicate an accumulation of adaptive changes in the SU(Z)12 protein during the divergence of D. subobscura and D. guanche. In contrast, the HKA test shows a deficit of polymorphism at Pho. The most likely explanation for this reduced variation is the location of this gene in the dot-like chromosome and would indicate that this chromosome also has null or very low recombination in D. subobscura, as reported in D. melanogaster (Calvo-Martin, 2017).

    Translating natural genetic variation to gene expression in a computational model of the Drosophila gap gene regulatory network
    This study applied a sequence-level model of gap gene expression in the early development of Drosophila to analyze single nucleotide polymorphisms (SNPs) in a panel of natural sequenced D. melanogaster lines. Using a thermodynamic modeling framework, both analytical and computational descriptions are provided of how single-nucleotide variants affect gene expression. The analysis reveals that the sequence variants increase (decrease) gene expression if located within binding sites of repressors (activators). The sign of SNP influence (activation or repression) may change in time and space. The thermodynamic modeling approach predicts non-local and non-linear effects arising from SNPs, and combinations of SNPs, in individual fly genotypes. Simulation of individual fly genotypes using this model reveals that this non-linearity reduces to almost additive inputs from multiple SNPs. Further, signatures are seen of the action of purifying selection in the gap gene regulatory regions. To infer the specific targets of purifying selection, the patterns of polymorphism in the data were analyzed at two phenotypic levels: the strengths of binding and expression. Combinations of SNPs show evidence of being under selective pressure, while individual SNPs do not. The model predicts that SNPs appear to accumulate in the genotypes of the natural population in a way biased towards small increases in activating action on the expression pattern (Gursky, 2017).

    Structure of the Transcriptional Regulatory Network Correlates with Regulatory Divergence in Drosophila

    Transcriptional control of gene expression is regulated by biochemical interactions between cis-regulatory DNA sequences and trans-acting factors that form complex regulatory networks. Genetic changes affecting both cis- and trans-acting sequences in these networks have been shown to alter patterns of gene expression as well as higher-order organismal phenotypes. This study investigated how the structure of these regulatory networks relates to patterns of polymorphism and divergence in gene expression. To do this, a transcriptional regulatory network inferred for Drosophila melanogaster was compared to differences in gene regulation observed between two strains of D. melanogaster as well as between two pairs of closely related species: Drosophila sechellia and Drosophila simulans, and D. simulans and D. melanogaster. The number of transcription factors predicted to directly regulate a gene ("in-degree") was negatively correlated with divergence in both gene expression (mRNA abundance) and cis-regulation. This observation suggests that the number of transcription factors directly regulating a gene's expression affects the conservation of cis-regulation and gene expression over evolutionary time. The hypothesis that transcription factors regulating more target genes (higher "out-degree") are less likely to evolve changes in their cis-regulation and expression (presumably due to increased pleiotropy) was also tested, but little support was found for this predicted relationship. Taken together, these data show how the architecture of regulatory networks can influence regulatory evolution (Yang, 2017).

    Cis- and trans-regulatory effects on gene expression in a natural population of Drosophila melanogaster

    Cis- and trans-regulatory mutations are important contributors to transcriptome evolution. Quantifying their relative contributions to intraspecific variation in gene expression is essential for understanding the population genetic processes that underlie evolutionary changes in gene expression. This study has examined this issue by quantifying genome-wide allele specific expression (ASE) variation using a crossing scheme that produces F1 hybrids between 18 different Drosophila melanogaster strains sampled from the Drosophila Genetic Reference Panel (DGRP) and a reference strain from another population. Head and body samples from F1 adult females were subjected to RNA-seq and the subsequent ASE quantification. Cis- and trans-regulatory effects on expression variation were estimated from these data. A higher proportion of genes showed significant cis-regulatory variation (~28%) than those showed significant trans-regulatory variation (~9%). The sizes of cis-regulatory effects on expression variation were 1.98 and 1.88 times larger than trans-regulatory effects in heads and bodies, respectively. A generalized linear model analysis revealed that both cis- and trans-regulated expression variation was strongly associated with nonsynonymous nucleotide diversity and tissue specificity. Interestingly, trans-regulated variation showed a negative correlation with local recombination rate. Also, analysis on proximal transposon element (TE) insertions suggested that they affect transcription levels of ovary-expressed genes more pronouncedly than genes not expressed in the ovary, possibly due to defense mechanisms against TE mobility in the germline. Collectively, this detailed quantification of ASE variations from a natural population has revealed a number of new relationships between genomic factors and the effects of cis- and trans-regulatory factors on expression variation (Osada, 2017).

    Pigmentation pattern and developmental constraints: flight muscle attachment sites delimit the thoracic trident of Drosophila melanogaster

    In their seminal paper published in 1979, Gould and Lewontin argued that some traits arise as by-products of the development of other structures and not for direct utility in themselves. This study shows that this applies to the trident, a pigmentation pattern observed on the thorax of Drosophila melanogaster. Using reporter constructs, it was shown that the expression domain of several genes encoding pigmentation enzymes follows the trident shape. This domain is complementary to the expression pattern of stripe (sr), which encodes an essential transcription factor specifying flight muscle attachment sites. sr limits the expression of these pigmentation enzyme genes to the trident by repressing them in its own expression domain, i.e. at the flight muscle attachment sites. Evidence is given that repression of not only yellow but also other pigmentation genes, notably tan, is involved in the trident shape. The flight muscle attachment sites and sr expression patterns are remarkably conserved in dipterans reflecting the essential role of sr. The data suggest that the trident is a by-product of flight muscle attachment site patterning that arose when sr was co-opted for the regulation of pigmentation enzyme coding genes (Gibert, 2018).

    Holding it together: rapid evolution and positive selection in the synaptonemal complex of Drosophila

    The synaptonemal complex (SC) is a highly conserved meiotic structure that functions to pair homologs and facilitate meiotic recombination in most eukaryotes. Five Drosophila SC proteins have been identified and localized within the complex: C(3)G, C(2)M, CONA, ORD, and the newly identified Corolla. The SC is required for meiotic recombination in Drosophila and absence of these proteins leads to reduced crossing over and chromosomal nondisjunction. The proteins display little apparent sequence conservation outside the genus. To identify factors that explain this lack of apparent conservation, a molecular evolutionary analysis of these genes was performed across the Drosophila genus. For the five SC components, gene sequence similarity declines rapidly with increasing phylogenetic distance and only ORD and C(2)M are identifiable outside of the Drosophila genus. SC gene sequences have a higher dN/dS (omega) rate ratio than the genome wide average and this can in part be explained by the action of positive selection in almost every SC component. Omega estimates for the five SC components are in accordance with their physical position within the SC. Components interacting with chromatin evolve slowest and components comprising the central elements evolve the most rapidly. Thus, the Drosophila SC is proposed to be evolving rapidly due to two combined effects: 1) a high rate of evolution can be partly explained by low purifying selection on protein components whose function is to simply hold chromosomes together, 2) positive selection in the SC is driven by its sex-specificity combined with its role in facilitating both recombination and centromere clustering in the face of recurrent bouts of drive in female meiosis (Hemmer, 2016).

    Genomic and transcriptomic associations identify a new insecticide resistance phenotype for the selective sweep at the Cyp6g1 locus of Drosophila melanogaster

    Scans of the Drosophila melanogaster genome have identified organophosphate resistance loci among those with the most pronounced signature of positive selection. In this study the molecular basis of resistance to the organophosphate insecticide azinphos-methyl was investigated using the Drosophila Genetic Reference Panel and genome-wide association. Recently released full transcriptome data was used to extend the utility of the Drosophila Genetic Reference Panel resource beyond traditional genomewide association studies to allow systems genetics analyses of phenotypes. Both genomic and transcriptomic associations independently identified Cyp6g1, a gene involved in resistance to DDT and neonicotinoid insecticides, as the top candidate for azinphos-methyl resistance. This was verified by transgenically overexpressing Cyp6g1 using natural regulatory elements from a resistant allele, resulting in a 6.5-fold increase in resistance. Four novel candidate genes associated with azinphos-methyl resistance were found, all of which are involved in either regulation of fat storage or nervous system development. In Cyp6g1, a demonstrable resistance locus was found, a verification that transcriptome data can be used to identify variants associated with insecticide resistance, and an overlap was found between peaks of a genome-wide association study and a genome-wide selective sweep analysis (Battlay, 2015).

    Rapid divergence and convergence of life-history in experimentally evolved Drosophila melanogaster

    selection experiments are alluring in their simplicity, power, and ability to inform about how evolution works. A longstanding challenge facing evolution experiments with metazoans is that significant generational turnover takes a long time. This study presents data from a unique system of experimentally evolved laboratory populations of Drosophila melanogaster that have experienced three distinct life-history selection regimes. The goal of the study was to determine how quickly populations of a certain selection regime diverge phenotypically from their ancestors, and how quickly they converge with independently derived populations that share a selection regime. Results indicate that phenotypic divergence from an ancestral population occurs rapidly, within dozens of generations, regardless of that population's evolutionary history. Similarly, populations sharing a selection treatment converge on common phenotypes in this same time frame, regardless of selection pressures those populations may have experienced in the past. These patterns of convergence and divergence were found to emerge much faster than expected, suggesting that intermediate evolutionary history has transient effects in this system. These results are applicable to other experimental evolution projects, and suggest that many relevant questions can be sufficiently tested on shorter timescales than previously thought (Burke, 2016).

    Survey of global genetic diversity within the Drosophila immune system

    Numerous studies across a wide range of taxa have demonstrated that immune genes are routinely among the most rapidly evolving genes in the genome. This observation, however, does not address what proportion of immune genes undergo strong selection during adaptation to novel environments. This study determined the extent of very recent divergence in genes with immune function across five populations of Drosophila melanogaster;t immune genes do not show an overall trend of recent rapid adaptation. The population-based approach uses a set of carefully matched control genes to account for the effects of demography and local recombination rate, allowing identification of whether specific immune functions are putative targets of strong selection. Evidence was found that viral defense genes are rapidly evolving in Drosophila at multiple time scales. Local adaptation to bacteria and fungi is less extreme and primarily occurs through changes in recognition and effector genes rather than large-scale changes to the regulation of the immune response. Surprisingly, genes in the Toll pathway, which show a high rate of adaptive substitutions between the D. melanogaster and D. simulans lineages, show little population differentiation. Quantifying the flies for resistance to a generalist Gram-positive bacterial pathogen, it was found that this genetic pattern of low population differentiation was recapitulated at the phenotypic level. In sum, these results highlight the complexity of immune evolution and suggest that Drosophila immune genes do not follow a uniform trajectory of strong directional selection as flies encounter new environments (Early, 2016).

    Genome-wide analysis of long-term evolutionary domestication in Drosophila melanogaster

    Experimental evolutionary genomics now allows biologists to test fundamental theories concerning the genetic basis of adaptation. This laboratory conducted one of the longest laboratory evolution experiments with any sexually-reproducing metazoan, Drosophila melanogaster. Next-generation resequencing data from this experiment was conducted to examine genome-wide patterns of genetic variation over an evolutionary time-scale that approaches 1,000 generations. Measures of variation within and differentiation between populations were compared to simulations based on a variety of evolutionary scenarios. This analysis yielded no clear evidence of hard selective sweeps, whereby natural selection acts to increase the frequency of a newly-arising mutation in a population until it becomes fixed. Evidence was found for selection acting on standing genetic variation, as independent replicate populations exhibit similar population-genetic dynamics, without obvious fixation of candidate alleles under selection. A hidden-Markov model test for selection also found widespread evidence for selection. More genetic variation was found genome-wide, and less differentiation was found between replicate populations genome-wide, than arose in any of the simulated evolutionary scenarios (Phillips, 2016).

    Variation in the intensity of selection on codon bias over time causes contrasting patterns of base composition evolution in Drosophila

    Four-fold degenerate coding sites form a major component of the genome, and are often used to make inferences about selection and demography, so that understanding their evolution is important. Despite previous efforts, many questions regarding the causes of base composition changes at these sites in Drosophila remain unanswered. To shed further light on this issue, a new whole-genome polymorphism dataset was obtained from D. simulans. Samples were analyzed from the putatively ancestral range of D. simulans, as well as an existing polymorphism dataset from an African population of D. melanogaster. By using D. yakuba as an outgroup, clear evidence was found for selection on 4-fold sites along both lineages over a substantial period, with the intensity of selection increasing with GC content. Based on an explicit model of base composition evolution, it is suggested that the observed AT-biased substitution pattern in both lineages is probably due to an ancestral reduction in selection intensity, and is unlikely to be the result of an increase in mutational bias towards AT alone. By using two polymorphism-based methods for estimating selection coefficients over different timescales, it was shown that the selection intensity on codon usage has been rather stable in D. simulans in the recent past, but the long-term estimates in D. melanogaster are much higher than the short-term ones, indicating a continuing decline in selection intensity, to such an extent that the short-term estimates suggest that selection is only active in the most GC-rich parts of the genome. Finally, evidence is provided for complex evolutionary patterns in the putatively neutral short introns, which cannot be explained by the standard GC-biased gene conversion model. These results reveal a dynamic picture of base composition evolution (Jackson, 2017).

    A variable genetic architecture of melanic evolution in Drosophila melanogaster

    Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. The genetic basis of this melanism was investigated in three populations using a new bulk segregant analysis mapping method. Nineteen distinct QTL regions were identified from 9 mapping crosses, with several QTL peaks overlapping between two or all populations, and yet different crosses involving the same melanic population commonly yielded distinct QTLs. The strongest QTLs often overlapped well-known pigmentation genes, but wide signals of genetic differentiation (FST) was typically not found between lightly and darkly pigmented populations at these genes. Instead, small numbers of highly differentiated SNPs were found at the probable causative genes. A simulation analysis showed that these patterns of polymorphism were consistent with selection on standing genetic variation. Overall, these results suggest that even for potentially simpler traits like pigmentation, the complexity of adaptive trait evolution poses important challenges for QTL mapping and population genetic analysis (Bastide, 2016).

    Promoter shape varies across populations and affects promoter evolution and expression noise

    Animal promoters initiate transcription either at precise positions (narrow promoters) or dispersed regions (broad promoters), a distinction referred to as promoter shape. Although highly conserved, the functional properties of promoters with different shapes and the genetic basis of their evolution remain unclear. This study used natural genetic variation across a panel of 81 Drosophila lines to measure changes in transcriptional start site (TSS) usage, identifying thousands of genetic variants affecting transcript levels (strength) or the distribution of TSSs within a promoter (shape). The results identify promoter shape as a molecular trait that can evolve independently of promoter strength. Broad promoters typically harbor shape-associated variants, with signatures of adaptive selection. Single-cell measurements demonstrate that variants modulating promoter shape often increase expression noise, whereas heteroallelic interactions with other promoter variants alleviate these effects. These results uncover new functional properties of natural promoters and suggest the minimization of expression noise as an important factor in promoter evolution (Schor, 2016).

    Predictable phenotypic, but not karyotypic, evolution of populations with contrasting initial history

    The relative impact of selection, chance and history will determine the predictability of evolution. There is a lack of empirical research on this subject, particularly in sexual organisms. This study used experimental evolution to test the predictability of evolution. The real-time evolution was analyzed of Drosophila subobscura populations derived from contrasting European latitudes placed in a novel laboratory environment. Each natural population was sampled twice within a three-year interval. Evolutionary responses at both phenotypic (life-history, morphological and physiological traits) and karyotypic levels were studied for around 30 generations of laboratory culture. The results show (1) repeatable historical effects between years in the initial state, at both phenotypic and karyotypic levels; (2) predictable phenotypic evolution with general convergence except for body size; and (3) unpredictable karyotypic evolution. It is concluded that the predictability of evolution is contingent on the trait and level of organization, highlighting the importance of studying multiple biological levels with respect to evolutionary patterns (Simoes, 2017).

    Naturally-segregating variation at Ugt86Dd contributes to nicotine resistance in Drosophila melanogaster

    Identifying the sequence polymorphisms underlying complex trait variation is a key goal of genetics research, since knowing the precise causative molecular events allows insight into the pathways governing trait variation. Genetic analysis of complex traits in model systems regularly starts by constructing QTL maps, but generally fails to identify causative sequence polymorphisms. Previous studies mapped a series of QTL contributing to resistance to nicotine in a Drosophila melanogaster multiparental mapping resource, and this study use a battery of functional tests to resolve QTL to the molecular level. One large-effect QTL resided over a cluster of UDP-glucuronosyltransferases, and quantitative complementation tests using deficiencies eliminating subsets of these detoxification genes revealed allelic variation impacting resistance. RNAseq showed that Ugt86Dd had significantly higher expression in genotypes that are more resistant to nicotine, and anterior midgut-specific RNAi of this gene reduced resistance. A segregating 22-bp frameshift deletion in Ugt86Dd, and accounting for the InDel during mapping largely eliminates the QTL, implying the event explains the bulk of the effect of the mapped locus. CRISPR/Cas9 editing a relatively resistant genotype to generate lesions in Ugt86Dd that recapitulate the naturally-occurring putative loss-of-function allele leads to a large reduction in resistance. Despite this major effect of the deletion, the allele appears to be very rare in wild-caught populations, and likely explains only a small fraction of the natural variation for the trait. Nonetheless, this putatively causative coding InDel can be a launchpad for future mechanistic exploration of xenobiotic detoxification (Highfill, 2017).

    Multiple P450s and variation in neuronal genes underpins the response to the insecticide Imidacloprid in a population of Drosophila melanogaster

    Insecticide resistance is an economically important example of evolution in response to intense selection pressure. In this study, the genetics of resistance to the neonicotinoid insecticide imidacloprid is explored using the Drosophila Genetic Reference Panel, a collection of inbred Drosophila melanogaster genotypes derived from a single population in North Carolina. Imidacloprid resistance varied substantially among genotypes, and more resistant genotypes tended to show increased capacity to metabolize and excrete imidacloprid. Variation in resistance level was then associated with genomic and transcriptomic variation, implicating several candidate genes involved in central nervous system function and the cytochrome P450s Cyp6g1 and Cyp6g2. CRISPR-Cas9 mediated removal of Cyp6g1 suggested that it contributed to imidacloprid resistance only in backgrounds where it was already highly expressed. Cyp6g2, previously implicated in juvenile hormone synthesis via expression in the ring gland, was shown to be expressed in metabolically relevant tissues of resistant genotypes. Cyp6g2 overexpression was shown to both metabolize imidacloprid and confer resistance. These data collectively suggest that imidacloprid resistance is influenced by a variety of previously known and unknown genetic factors (Denecke, 2017).

    Nearly neutral evolution across the Drosophila melanogaster genome

    Under the nearly neutral theory of molecular evolution, the proportion of effectively neutral mutations is expected to depend upon the effective population size (Ne). This study investigate whether this is the case across the genome of Drosophila melanogaster using polymorphism data from North American and African lines. This study shows that the ratio of the number of nonsynonymous and synonymous polymorphisms is negatively correlated to the number of synonymous polymorphisms, even when the nonindependence is accounted for. The relationship is such that the proportion of effectively neutral nonsynonymous mutations increases by approximately 45% as Ne is halved. However, it was also shown that this relationship is steeper than expected from an independent estimate of the distribution of fitness effects from the site frequency spectrum. A number of potential explanations for this was investigated, and it was shown, using simulation, that this is consistent with a model of genetic hitchhiking: Genetic hitchhiking depresses diversity at neutral and weakly selected sites, but has little effect on the diversity of strongly selected sites (Castellano, 2018).

    Pleiotropy modulates the efficacy of selection in Drosophila melanogaster

    Pleiotropy is the well-established idea that a single mutation affects multiple phenotypes. If a mutation has opposite effects on fitness when expressed in different contexts, then genetic conflict arises. Pleiotropic conflict is expected to reduce the efficacy of selection by limiting the fixation of beneficial mutations through adaptation, and the removal of deleterious mutations through purifying selection. While this has been widely discussed, in particular in the context of a putative "gender load", it has yet to be systematically quantified. This work empirically estimates to which extent different pleiotropic regimes impede the efficacy of selection in Drosophila melanogaster. Whole-genome polymorphism data from a single African population and divergence data from D. simulans were used to estimate the fraction of adaptive fixations (alpha), the rate of adaptation (omegaA) and the direction of selection (DoS). After controlling for confounding covariates, it was found that the different pleiotropic regimes have a relatively small, but significant, effect on selection efficacy. Specifically, the results suggest that pleiotropic sexual antagonism may restrict the efficacy of selection, but that this conflict can be resolved by limiting the expression of genes to the sex where they are beneficial. Intermediate levels of pleiotropy across tissues and life stages can also lead to maladaptation in D. melanogaster, due to inefficient purifying selection combined with low frequency of mutations that confer a selective advantage. Thus, this study highlights the need to consider the efficacy of selection in the context of antagonistic pleiotropy, and of genetic conflict in general (Fraisse, 2018).

    Selection on inversion breakpoints favors proximity to pairing sensitive sites in Drosophila melanogaster

    Chromosomal inversions are widespread among taxa, and have been implicated in a number of biological processes including adaptation, sex chromosome evolution, and segregation distortion. Consistent with selection favoring linkage between loci, it is well established that length is a selected trait of inversions. However, the factors that affect the distribution of inversion breakpoints remain poorly understood. 'Sensitive sites' have been mapped on all euchromatic chromosome arms in Drosophila melanogaster, and may be a source of natural selection on inversion breakpoint positions. Briefly, sensitive sites are genomic regions wherein proximal structural rearrangements result in large reductions in local recombination rates in heterozygotes. This study shows that breakpoints of common inversions are significantly more likely to lie within a cytological band containing a sensitive site than are breakpoints of rare inversions. Furthermore, common inversions for which neither breakpoint intersects a sensitive site are significantly longer than rare inversions, but common inversions whose breakpoints intersect a sensitive site show no evidence for increased length. These results are interpreted to mean that selection favors inversions whose breakpoints disrupt synteny near to sensitive sites, possibly because these inversions suppress recombination in large genomic regions. This is the first evidence consistent with positive selection acting on inversion breakpoint positions (Corbett-Detig, 2016).

    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. This study introduced 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 this approach, it was applied to genome-wide resequencing data from 125 lines in Drosophila melanogaster and diversity levels was reliably predicted at the 1Mb scale. The results corroborate estimates of a high fraction of beneficial substitutions in proteins and untranslated regions (UTR). They allow distinguishing between the contribution of sweeps and other modes of selection around amino acid substitutions and uncovered evidence for pervasive sweeps in untranslated regions (UTRs). This inference further suggests a substantial effect of other modes of linked selection and of adaptation in particular. More generally, it was demonstrated that linked selection has had a larger effect in reducing diversity levels and increasing their variance in D. melanogaster than previously appreciated (Elyashiv, 2016).

    Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster

    Highly repetitive satellite DNA (satDNA) repeats are found in most eukaryotic genomes. SatDNAs are rapidly evolving and have roles in genome stability and chromosome segregation. Their repetitive nature poses a challenge for genome assembly and makes progress on the detailed study of satDNA structure difficult. This study used single-molecule sequencing long reads from Pacific Biosciences (PacBio) to determine the detailed structure of all major autosomal complex satDNA loci in Drosophila melanogaster, with a particular focus on the 260-bp and Responder satellites. The optimal de novo assembly methods and parameter combinations were determined that were required to produce a high-quality assembly of these previously unassembled satDNA loci and validate this assembly using molecular and computational approaches. It was determined that the computationally intensive PBcR-BLASR assembly pipeline yielded better assemblies than the faster and more efficient pipelines based on the MHAP hashing algorithm, and it is essential to validate assemblies of repetitive loci. The assemblies reveal that satDNA repeats are organized into large arrays interrupted by transposable elements. The repeats in the center of the array tend to be homogenized in sequence, suggesting that gene conversion and unequal crossovers lead to repeat homogenization through concerted evolution, although the degree of unequal crossing over may differ among complex satellite loci. Evidence was found for higher-order structure within satDNA arrays that suggest recent structural rearrangements. These assemblies provide a platform for the evolutionary and functional genomics of satDNAs in pericentric heterochromatin (Khost, 2017).

    High rate of translocation-based gene birth on the Drosophila Y chromosome

    The Y chromosome is a unique genetic environment defined by a lack of recombination and male-limited inheritance. The Drosophila Y chromosome has been gradually acquiring genes from the rest of the genome, with only seven Y-linked genes being gained over the past 63 million years (0.12 gene gains per million years). Using a next-generation sequencing (NGS)-powered genomic scan, this study shows that gene transfers to the Y chromosome are much more common than previously suspected: at least 25 have arisen across three Drosophila species over the past 5.4 million years (1.67 per million years for each lineage). The gene transfer rate is significantly lower in Drosophila melanogaster than in the Drosophila simulans clade, primarily due to Y-linked retrotranspositions being significantly more common in the latter. Despite all Y-linked gene transfers being evolutionarily recent (<1 million years old), only three showed evidence for purifying selection (omega </= 0.14). Thus, although the resulting Y-linked functional gene acquisition rate (0.25 new genes per million years) is double the longer-term estimate, the fate of most new Y-linked genes is defined by rapid degeneration and pseudogenization. These results show that Y-linked gene traffic, and the molecular mechanisms governing these transfers, can diverge rapidly between species, revealing the Drosophila Y chromosome to be more dynamic than previously appreciated. This analytical method provides a powerful means to identify Y-linked gene transfers and will help illuminate the evolutionary dynamics of the Y chromosome in Drosophila and other species (Tobler, 2017).

    The pdm3 locus is a hotspot for recurrent evolution of female-limited color dimorphism in Drosophila

    Sex-limited polymorphisms are an intriguing form of sexual dimorphism that offer unique opportunities to reconstruct the evolutionary changes that decouple male and female traits encoded by a shared genome. This study investigated the genetic basis of a Mendelian female-limited color dimorphism (FLCD) that segregates in natural populations of more than 20 species of the Drosophila montium subgroup. In these species, females have alternative abdominal color morphs, light and dark, whereas males have only one color morph in each species. A comprehensive molecular phylogeny of the montium subgroup supports multiple origins of FLCD. Despite this, FLCD mapped to the same locus in four distantly related species-the transcription factor POU domain motif 3 (pdm3), which acts as a repressor of abdominal pigmentation in D. melanogaster. In D. serrata, FLCD maps to a structural variant in the first intron of pdm3; however, this variant is not found in the three other species-D. kikkawai, D. leontia, and D. burlai-and sequence analysis strongly suggests the pdm3 alleles responsible for FLCD originated independently at least three times. It is proposed that cis-regulatory changes in pdm3 form sexually dimorphic and monomorphic alleles that segregate within species and are preserved, at least in one species, by structural variation. Surprisingly, pdm3 has not been implicated in the evolution of sex-specific pigmentation outside the montium subgroup, suggesting that the genetic paths to sexual dimorphism may be constrained within a clade but variable across clades (Yassin, 2016).

    Mapping QTL contributing to variation in posterior lobe morphology between strains of Drosophila melanogaster

    Closely-related, and otherwise morphologically similar insect species frequently show striking divergence in the shape and/or size of male genital structures, a phenomenon thought to be driven by sexual selection. Comparative interspecific studies can help elucidate the evolutionary forces acting on genital structures to drive this rapid differentiation. However, genetic dissection of sexual trait divergence between species is frequently hampered by the difficulty generating interspecific recombinants. Intraspecific variation can be leveraged to investigate the genetics of rapidly-evolving sexual traits; this study carried out out a genetic analysis of variation in the posterior lobe within D. melanogaster. The lobe is a male-specific process emerging from the genital arch of D. melanogaster and three closely-related species, is essential for copulation, and shows radical divergence in form across species. There is also abundant variation within species in the shape and size of the lobe, and while this variation is considerably more subtle than that seen among species, it nonetheless provides the raw material for QTL mapping. An advanced intercross population was created from a pair of phenotypically-different inbred strains, and after phenotyping and genotyping-by-sequencing the recombinants, several QTL contributing to various measures of lobe morphology were mapped. The additional generations of crossing over in the mapping population led to QTL intervals that are smaller than is typical for an F2 mapping design. The intervals that were mapped overlap with a pair of lobe QTL that was previously identified in an independent mapping cross, potentially suggesting a level of shared genetic control of trait variation. The QTL additionally implicate a suite of genes that have been shown to contribute to the development of the posterior lobe. These loci are strong candidates to harbor naturally-segregating sites contributing to phenotypic variation within D. melanogaster, and may also be those contributing to divergence in lobe morphology between species (Hackett, 2006).

    Sex-specific selection and sex-biased gene expression in humans and flies

    Sexual dimorphism results from sex-biased gene expression, which evolves when selection acts differently on males and females. While there is an intimate connection between sex-biased gene expression and sex-specific selection, few empirical studies have studied this relationship directly. This study compare the two on a genome-wide scale in humans and flies. A distinctive "Twin Peaks" pattern in humans that relates the strength of sex-specific selection, quantified by genetic divergence was yound between male and female adults at autosomal loci, to the degree of sex-biased expression. Genes with intermediate degrees of sex-biased expression show evidence of ongoing sex-specific selection, while genes with either little or completely sex-biased expression do not. This pattern apparently results from differential viability selection in males and females acting in the current generation. The Twin Peaks pattern is also found in Drosophila using a different measure of sex-specific selection acting on fertility. A simple model was developed that successfully recapitulates the Twin Peaks. These results suggest that many genes with intermediate sex-biased expression experience ongoing sex-specific selection in humans and flies (Cheng, 2016).

    Sexual selection shapes development and maturation rates in Drosophila

    Explanations for the evolution of delayed maturity usually invoke trade-offs mediated by growth, but processes of reproductive maturation continue long after growth has ceased. This study tested whether sexual selection shapes the rate of posteclosion maturation in the fruit fly Drosophila melanogaster. Populations maintained for more than 100 generations under a short generation time and polygamous mating system evolved faster posteclosion maturation and faster egg-to-adult development of males, when compared to populations kept under short generations and randomized monogamy that eliminated sexual selection. An independent assay demonstrated that more mature males have higher fitness under polygamy, but this advantage disappears under monogamy. In contrast, for females greater maturity was equally advantageous under polygamy and monogamy. Furthermore, monogamous populations evolved faster development and maturation of females relative to polygamous populations, with no detectable trade-offs with adult size or egg-to-adult survival. These results suggest that a major aspect of male maturation involves developing traits that increase success in sexual competition, whereas female maturation is not limited by investment in traits involved in mate choice or defense against male antagonism. Moreover, rates of juvenile development and adult maturation can readily evolve in opposite directions in the two sexes, possibly implicating polymorphisms with sexually antagonistic pleiotropy (Hollis, 2016).

    Sex-biased transcriptome divergence along a latitudinal gradient

    Sex-dependent gene expression is likely an important genomic mechanism that allows sex-specific adaptation to environmental changes. Among Drosophila species, sex-biased genes display remarkably consistent evolutionary patterns; male-biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex-biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex-specific selection and the evolution of sex-biased gene expression. This study assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11,631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male-biased genes, there was no overrepresentation of X-linked genes in males. By contrast, X-linked divergence was elevated in females, especially for female-biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. These results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro- and micro-ecological spatial scales (Allen, 2017).

    The large X-effect on secondary sexual characters and the genetics of variation in sex comb tooth number in Drosophila subobscura

    This study examined the genetics of a secondary sexual trait, male sex comb size in Drosophila subobscura, to evaluate the amount of variation attributable to the X-chromosome. This species bears unusually large sex combs for its species group, and therefore, this trait may be a good candidate for having been affected by natural or sexual selection. Significant heritable variation was observed in number of teeth of the distal sex comb across strains. While reciprocal F1 crosses seemed to implicate a disproportionate X-chromosome effect, further examination in the F2 progeny showed that transgressive autosomal effects inflated the estimate of variation associated with the X-chromosome in the F1. Instead, the X-chromosome appears to confer the smallest contribution of all major chromosomes to the observed phenotypic variation. Further, effects on copulation latency or duration associated with the observed phenotypic variation were not observed. Overall, this study presents an examination of the genetics underlying segregating phenotypic variation within species and illustrates two common pitfalls associated with some past studies of the genetic basis of secondary sexual traits (Mittleman, 2017).

    Direct benefits of choosing a high fitness mate can offset the indirect costs associated with intralocus sexual conflict

    Intralocus sexual conflict generates a cost to mate choice: high fitness partners transmit genetic variation that confers lower fitness to offspring of the opposite sex. Earlier work in the fruit fly, Drosophila melanogaster, revealed that these indirect genetic costs were sufficient to reverse potential "good genes" benefits of sexual selection. However, mate choice can also confer direct fitness benefits by inducing larger numbers of progeny. This study considers whether direct benefits through enhanced fertility could offset the costs associated with intralocus sexual conflict in D. melanogaster. Using hemiclonal analysis, it was found that females mated to high fitness males produced 11% more offspring compared to those mated to low fitness males, and high fitness females produced 37% more offspring than low fitness females. These direct benefits more than offset the reduction in offspring fitness caused by intralocus sexual conflict, creating a net fitness benefit for each sex to pairing with a high fitness partner. These findings highlight the need to consider both direct and indirect effects when investigating the fitness impacts of mate choice. Direct fitness benefits may shelter sexually antagonistic alleles from selection, suggesting a novel mechanism for the maintenance of fitness variation (Pischedda, 2017).

    Evolutionary dynamics of male reproductive genes in the Drosophila virilis subgroup

    Postcopulatory sexual selection (PCSS) is a potent evolutionary force that can drive rapid changes of reproductive genes within species, and thus has the potential to generate reproductive incompatibilities between species. Male seminal fluid proteins (SFPs) are major players in postmating interactions, and are important targets of PCSS in males. The virilis subgroup of Drosophila exhibits strong interspecific gametic incompatibilities, and can serve as a model to study the genetic basis of PCSS and gametic isolation. However, reproductive genes in this group have not been characterized. This study utilized short-read RNA sequencing of male reproductive organs to examine the evolutionary dynamics of reproductive genes in members of the virilis subgroup: D. americana, D. lummei, D. novamexicana, and D. virilis. The majority of male reproductive transcripts are testes-biased, accounting for ~15% of all annotated genes. Ejaculatory bulb-biased transcripts largely code for lipid metabolic enzymes, and contain orthologs of the D. melnaogaster ejaculatory bulb protein, Peb-me, which is involved in mating-plug formation. In addition, 71 candidate SFPs were identified, and this gene set was show to have the highest rate of nonsynonymous codon substitution relative to testes- and ejaculatory bulb-biased genes. Furthermore, orthologs were identified of 35 D. melanogaster SFPs that have conserved accessory gland expression in the virilis group. Finally, several of the SFPs that have the highest rate of nonsynonymous codon substitution were shown to reside on chromosomal regions that contribute to paternal gametic incompatibility between species. These results show that SFPs rapidly diversify in the virilis group, and suggest that they likely play a role in PCSS and/or gametic isolation (Ahmed-Braimah, 2017).

    Reproductive isolation through experimental manipulation of sexually antagonistic coevolution in Drosophila melanogaster

    Promiscuity can drive the evolution of sexual conflict before and after mating occurs. Post mating, the male ejaculate can selfishly manipulate female physiology, leading to a chemical arms race between the sexes. Theory suggests that drift and sexually antagonistic coevolution can cause allopatric populations to evolve different chemical interactions between the sexes, thereby leading to postmating reproductive barriers and speciation. There is, however, little empirical evidence supporting this form of speciation. This theory was tested by creating an experimental evolutionary model of Drosophila melanogaster populations undergoing different levels of interlocus sexual conflict. Allopatric populations under elevated sexual conflict were shown to exhibit assortative mating, indicating premating reproductive isolation. Further, these allopatric populations also show reduced copulation duration and sperm defense ability when mating happens between individuals across populations compared to that within the same population, indicating postmating prezygotic isolation. Sexual conflict can cause reproductive isolation in allopatric populations through the coevolution of chemical (postmating prezygotic) as well as behavioural (premating) interactions between the sexes. Thus it study provides the first comprehensive evidence of postmating (as well as premating) reproductive isolation due to sexual conflict (Syed, 2017).

    Male relatedness and familiarity are required to modulate male-induced harm to females in Drosophila

    Males compete over mating and fertilization, and often harm females in the process. Inclusive fitness theory predicts that increasing relatedness within groups of males may relax competition and discourage male harm of females as males gain indirect benefits. Recent studies in Drosophila melanogaster are consistent with these predictions, and have found that within-group male relatedness increases female fitness, though others have found no effects. Importantly, these studies did not fully disentangle male genetic relatedness from larval familiarity, so the extent to which modulation of harm to females is explained by male familiarity remains unclear. This study performed a fully factorial design, isolating the effects of male relatedness and larval familiarity on female harm. While no differences were found in male courtship or aggression, there was a significant interaction between male genetic relatedness and familiarity on female reproduction and survival. Relatedness among males increased female lifespan, reproductive lifespan and overall reproductive success, but only when males were familiar. By showing that both male relatedness and larval familiarity are required to modulate female harm, these findings reconcile previous studies, shedding light on the potential role of indirect fitness effects on sexual conflict and the mechanisms underpinning kin recognition in fly populations (Le Page, 2017).

    A species-specific multigene family mediates differential sperm displacement in Drosophila melanogaster

    Sperm competition is a postcopulatory selection mechanism in species in which females mate with multiple males. Despite its evolutionary relevance in shaping male traits, the genetic mechanisms underlying sperm competition are poorly understood. A recently originated multigene family specific to Drosophila melanogaster, Sdic, is important for the outcome of sperm competition in doubly mated females, although the mechanistic nature of this phenotype remained unresolved. This study compared doubly mated females, second mated to either Sdic knockout or nonknockout males, and directly visualize sperm dynamics in the female reproductive tract. A less effective removal of first-to-mate male's sperm within the female's sperm storage organs is consistent with a reduced sperm competitive ability of the Sdic knockout males. These results highlight the role young genes can play in driving the evolution of sperm competition (Jayaswal, 2018).

    Precopulatory but not postcopulatory male reproductive traits diverge in response to mating system manipulation in Drosophila melanogaster

    Competition between males creates potential for pre- and postcopulatory sexual selection and conflict. Theory predicts that males facing risk of sperm competition should evolve traits to secure their reproductive success. If those traits are costly to females, the evolution of such traits may also increase conflict between the sexes. Conversely, under the absence of sperm competition, one expectation is for selection on male competitive traits to relax thereby also relaxing sexual conflict. Experimental evolution studies are a powerful tool to test this expectation. Studies in multiple insect species have yielded mixed and partially conflicting results. This study evaluated male competitive traits and male effects on female costs of mating in Drosophila melanogaster after replicate lines evolved for more than 50 generations either under enforced monogamy or sustained polygamy, thus manipulating the extent of intrasexual competition between males. In a setting where males competed directly with a rival male for access to a female and fertilization of her ova, polygamous males had superior reproductive success compared to monogamous males. When comparing reproductive success solely in double mating standard sperm competition assays, however, no difference was found in male sperm defense competitiveness between the different selection regimes. Instead, monogamous males were found to be inferior in precopulatory competition, which indicates that in this system, enforced monogamy relaxed selection on traits important in precopulatory rather than postcopulatory competition. These findings are discussed in the context of findings from previous experimental evolution studies in Drosophila and other invertebrate species (Wensing, 2017).

    Evolving doublesex expression correlates with the origin and diversification of male sexual ornaments in the Drosophila immigrans species group

    Male ornaments and other sex-specific traits present some of the most dramatic examples of evolutionary innovations. Comparative studies of similar but independently evolved traits are particularly important for identifying repeated patterns in the evolution of these traits. Male-specific modifications of the front legs have evolved repeatedly in Drosophilidae and other Diptera. The best understood of these novel structures is the sex comb of Drosophila melanogaster and its close relatives. This study examined the evolution of another male foreleg modification, the sex brush, found in the distantly related Drosophila immigrans species group. Similar to the sex comb, it was found that the origin of the sex brush correlates with novel, spatially restricted expression of the doublesex (dsx) transcription factor, the primary effector of the Drosophila sex determination pathway. The diversity of Dsx expression patterns in the immigrans species group closely reflects the differences in the presence, position, and size of the sex brush. Together with previous work on sex comb evolution, these observations suggest that tissue-specific activation of dsx expression may be a common mechanism responsible for the evolution of sexual dimorphism and particularly for the origin of novel male-specific ornaments (Rice, 2018).

    Evolutionary divergence in competitive mating success through female mating bias for good genes

    Despite heritable variation for univariate sexually selected traits, recent analyses exploring multivariate traits find evidence consistent with the lek paradox in showing no genetic variation available to choosy females, and therefore no genetic benefits of choice. This study used the preferences of Drosophila melanogaster females to exert bidirectional selection on competitive male mating success to test for the presence and nature of genetic variation underlying this multivariate trait. Male mating success diverged between selection regimens, and flies from success-selected lines had a smaller burden of deleterious, recessive mutations that affect egg-to-adult viability, were better sperm competitors (sperm offence), and did not demonstrate reduced desiccation resistance or components of female fitness (traits thought to trade off with attractiveness) relative to flies from failure-selected populations. Mating success remained subject to inbreeding depression in success-selected lines, suggesting that variation in mating success remains, thanks to numerous genes of small effect. Together, these results provide unique evidence for the evolutionary divergence in male mating success, demonstrating that genetic variation is not exhausted along the axis of precopulatory sexual selection and that female mating biases align with the avoidance of bad genes (Dugand, 2018).

    Evolution of sex-biased gene expression and dosage compensation in the eye and brain of Heliconius butterflies

    Differences in behavior and life history traits between females and males are the basis of divergent selective pressures between sexes. It has been suggested that a way for the two sexes to deal with different life history requirements is through sex-biased gene expression. A comparative sex-biased gene expression analysis was performed of the combined eye and brain transcriptome from five Heliconius species, H. charithonia, H. sara, H. erato, H. melpomene and H. doris, representing five of the main clades from the Heliconius phylogeny. The degree of sexual dimorphism in gene expression is not conserved across Heliconius. Most of the sex-biased genes identified in each species are not sex-biased in any other, suggesting that sexual selection might have driven sexually dimorphic gene expression. Only three genes shared sex-biased expression across multiple species: ultraviolet opsin UVRh1 and orthologs of Drosophila Kruppel-homolog 1 and CG9492. It was also observed that in some species female-biased genes have higher evolutionary rates, but in others, male-biased genes show the fastest rates when compared with unbiased genes, suggesting that selective forces driving sex-biased gene evolution in Heliconius act in a sex- and species-specific manner. Furthermore, dosage compensation was found in all the Heliconius tested, providing additional evidence for the conservation of dosage compensation across Lepidoptera. Finally, sex-biased genes are significantly enriched on the Z, a pattern that could be a result of sexually antagonistic selection (Catalan, 2018).

    Evolution of a central neural circuit underlies Drosophila mate preferences

    Courtship rituals serve to reinforce reproductive barriers between closely related species. Drosophila melanogaster and Drosophila simulans exhibit reproductive isolation, owing in part to the fact that D. melanogaster females produce 7,11-heptacosadiene, a pheromone that promotes courtship in D. melanogaster males but suppresses courtship in D. simulans males. This study compared pheromone-processing pathways in D. melanogaster and D. simulans males to define how these sister species endow 7,11-heptacosadiene with the opposite behavioural valence to underlie species discrimination. Males of both species were shown to detect 7,11-heptacosadiene using homologous peripheral sensory neurons, but this signal is differentially propagated to P1 neurons, which control courtship behaviour. A change in the balance of excitation and inhibition onto courtship-promoting neurons transforms an excitatory pheromonal cue in D. melanogaster into an inhibitory cue in D. simulans. These results reveal how species-specific pheromone responses can emerge from conservation of peripheral detection mechanisms and diversification of central circuitry, and demonstrate how flexible nodes in neural circuits can contribute to behavioural evolution (Seeholzer, 2018).

    The sensory periphery has been proposed to be the most evolutionarily labile element of the nervous system, as changes in the expression or tuning of sensory receptors can allow for the emergence of species-specific behaviours without necessitating potentially more complex developmental rewiring of central pathways in the brain. By contrast, the current results suggest that species-specific behavioural responses to 7,11-HD arise through functional alterations in how pheromone signals are propagated through a structurally conserved central circuit that consists of parallel excitatory and feedforward inhibitory branches. By reweighting the balance of excitatory vAB3 and inhibitory mAL signalling to P1 neurons, 7,11-HD is transformed from an excitatory signal that promotes courtship in D. melanogaster males into an inhibitory signal that suppresses courtship in D. simulans males. Although the analysis highlights the importance of this pheromone pathway in shaping species-specific mate preferences, the possibility cannot be excluded that additional inputs to P1 neurons or other targets of mAL and vAB3 neurons also contribute to divergent courtship decisions. Nevertheless, the data suggest that the branched architecture of pheromone-processing pathways serves as a substrate for the evolution of mate preferences, pointing to the potential existence of favourable sites within neural circuits to instantiate adaptive behavioural changes, analogous to how specific nodes within developmental regulatory networks contribute to morphological diversity (Seeholzer, 2018).

    The conserved tuning of Ppk23+ neurons suggests that D. simulans males dedicate this sensory pathway to the detection of a female pheromones, rather than sensing the chemical cues carried by conspecific females. Consistent with this idea, although P1 neurons are sufficient to elicit courtship in D. simulans males, they are not excited by the taste of a D. simulans female. Given that D. simulans cuticular pheromones are sexually monomorphic and offer ambiguous signals for mate recognition, males probably rely on additional sensory inputs for their arousal. Indeed, the fervent courtship exhibited by D. simulans ppk23 mutants towards females of different species demonstrates that D. simulans males can be aroused in the absence of any species-specific excitatory cue. Dienes, such as 7,11-HD, actively suppress this arousal, presumably through recruitment of strong mAL-mediated inhibition via Ppk23+ pathways. Together, these observations reinforce the notion that pheromone communication in Drosophila serves to focus a male's desire, such that flies lacking cuticular pheromones can be inherently attractive and appropriate mate choices are honed by specific inhibitory chemical cues (Seeholzer, 2018).

    Peripheral adaptations are likely to play an important role in the evolution of novel chemical sensitivities. In this case, however, preserving the sensory periphery while varying central circuits provides a mechanism to alter the behavioural valence of a single chemical cue. As D. melanogaster and D. simulans diverged, their reproductive isolation was probably strengthened by the ability of both species to detect the same pheromone but assign it a different meaning through these central circuit modifications. Notably, in D. melanogaster, P1 neuron excitability is regulated by the social history of a male. This suggests that both experience-dependent and evolutionary adaptations may act on the same neural substrate to modify sensory integration and mate choices, similar to how phenotypic plasticity may facilitate morphological evolution35. Therefore, functional reweighting of sensory inputs at flexible nodes in the nervous system, shaped by evolutionary selection or individual experience, may allow for alternative behavioural responses to the same sensory signal (Seeholzer, 2018).

    Phylogenomic Insights into the Evolution of Stinging Wasps and the Origins of Ants and Bees

    The stinging wasps (Hymenoptera: Aculeata) are an extremely diverse lineage of hymenopteran insects, encompassing over 70,000 described species and a diversity of life history traits, including ectoparasitism, cleptoparasitism, predation, pollen feeding (bees [Anthophila] and Masarinae), and eusociality (social vespid wasps, ants, and some bees). The most well-studied lineages of Aculeata are the ants, which are ecologically dominant in most terrestrial ecosystems, and the bees, the most important lineage of angiosperm-pollinating insects. Establishing the phylogenetic affinities of ants and bees helps in understanding and reconstruction of patterns of social evolution as well as leading to full appreciation of the biological implications of the switch from carnivory to pollen feeding (pollenivory). Despite recent advancements in aculeate phylogeny, considerable uncertainty remains regarding higher-level relationships within Aculeata, including the phylogenetic affinities of ants and bees. Ultraconserved element (UCE) phylogenomics was used to resolve relationships among stinging-wasp families, gathering sequence data from >800 UCE loci and 187 samples, including 30 out of 31 aculeate families. The 187-taxon dataset was analyzed using multiple analytical approaches, and several alternative taxon sets were analyzed. Alternative hypotheses for the phylogenetic positions of ants and bees were tested. The results present a highly supported phylogeny of the stinging wasps. Most importantly, it was found unequivocal evidence that ants are the sister group to bees+apoid wasps (Apoidea) and that bees are nested within a paraphyletic Crabronidae. It was also demonstrated that taxon choice can fundamentally impact tree topology and clade support in phylogenomic inference (Branstetter, 2017).

    Interactions between the developmental and adult social environments mediate group dynamics and offspring traits in Drosophila melanogaster

    Developmental conditions can strongly influence adult phenotypes and social interactions, which in turn affect key evolutionary processes such as sexual selection and sexual conflict. While the implications of social interactions in phenotypically mixed populations at the individual level are increasingly well known, how these effects influence the fate of groups remains poorly understood, which limits understanding of the broader ecological implications. To address this problem this study manipulated adult phenotypes and social composition in Drosophila melanogaster - by experimentally manipulating the larval density of the group-members - and measured a range of group-level outcomes across the lifespan of groups. Adult groups composed of exclusively low larval-density individuals showed high courtship levels, and low early reproductive rates, group growth rates, offspring mass and offspring eclosion success, relative to high larval-density or mixed larval-density groups. Furthermore, high larval-density groups had lower survival. Offspring mass increased with time, but at a reduced rate in groups when male group members (but not females) were from a mixture of larval-densities; peak reproductive rates were also earlier in these groups. These results suggest that that variation in developmental conditions experienced by adult group members can modify the reproductive output of groups (Morimoto, 2017).

    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. This study shows that changes in host plant use between distinct D. mojavensis populations are accompanied by changes in the olfactory system. Specifically, differences were observed 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 (Crowley-Gall, 2016).

    Shedding light on the grey zone of speciation along a continuum of genomic divergence

    Speciation results from the progressive accumulation of mutations that decrease the probability of mating between parental populations or reduce the fitness of hybrids-the so-called species barriers. Of primary importance is the prevalence of gene flow between diverging entities, which is central in most species concepts and has been widely discussed in recent years. This study explored the continuum of speciation thanks to a comparative analysis of genomic data from 61 pairs of populations/species of animals with variable levels of divergence. The intermediate "grey zone" of speciation, in which taxonomy is often controversial, spans from 0.5% to 2% of net synonymous Divergence, irrespective of species life history traits or ecology. Thanks to appropriate modeling of among-locus variation in genetic drift and introgression rate, this study clarifies the status of the majority of ambiguous cases and uncovers a number of cryptic species. This analysis also reveals the high incidence in animals of semi-isolated species (when some but not all loci are affected by barriers to gene flow) and highlights the intrinsic difficulty, both statistical and conceptual, of delineating species in the grey zone of speciation (Roux, 2016).

    Genetic differentiation and adaptive evolution at reproductive loci in incipient Drosophila species

    Accessory gland proteins (Acps; see Drosophila Sex Peptide and 0vulin) are part of the seminal fluid of male Drosophila flies. Some Acps have exceptionally high evolutionary rates and evolve under positive selection. Proper interactions between Acps and female reproductive molecules are essential for fertilization. These observations lead to suggestions that fast evolving Acps could be involved in speciation by promoting reproductive incompatibilities between emerging species. To test this hypothesis, population genetics data were used for three sibling species: D. mayaguana, D. parisiena, and D. straubae. The latter two species are morphologically very similar and show only incipient reproductive isolation. This system allowed examination of Acp evolution at different time frames in respect to speciation and reproductive isolation. Comparing data of 14 Acp loci with data obtained for other genomic regions, it was found that some Acps show extraordinarily high levels of divergence between D. mayaguana and its two sister species D. parisiena, and D. straubae. This divergence was likely driven by adaptive evolution at several loci. No fixed nucleotide differences were found between D. parisiena and D. straubae, however. Nevertheless, some Acp loci did show significant differentiation between these species associated with signs of positive selection; these loci may be involved in this early phase of the speciation process (Almeida, 2016).

    Early events in speciation: Cryptic species of Drosophila aldrichi

    Understanding the earliest events in speciation remains a major challenge in evolutionary biology. Thus identifying species whose populations are beginning to diverge can provide useful systems to study the process of speciation. Drosophila aldrichi, a cactophilic fruit fly species with a broad distribution in North America, has long been assumed to be a single species owing to its morphological uniformity. While previous reports either of genetic divergence or reproductive isolation among different D. aldrichi strains have hinted at the existence of cryptic species, the evolutionary relationships of this species across its range have not been thoroughly investigated. This study shows that D. aldrichi actually is paraphyletic with respect to its closest relative, Drosophila wheeleri, and that divergent D. aldrichi lineages show complete hybrid male sterility when crossed. The data support the interpretation that there are at least two species of D. aldrichi, making these flies particularly attractive for studies of speciation in an ecological and geographical context (Castro Vargas, 2017).

    Genetic divergence and the number of hybridizing species affect the path to homoploid hybrid speciation

    Hybridization is often maladaptive and in some instances has led to the loss of biodiversity. However, hybridization can also promote speciation, such as during homoploid hybrid speciation, thereby generating biodiversity. Despite examples of homoploid hybrid species, the importance of hybridization as a speciation mechanism is still widely debated, and a general understanding is lacking of the conditions most likely to generate homoploid hybrid species. This study shows that the level of genetic divergence between hybridizing species has a large effect on the probability that their hybrids evolve reproductive isolation. Populations of hybrids formed by parental species with intermediate levels of divergence were more likely to mate assortatively, and discriminate against their parental species, than those generated from weakly or strongly diverged parental species. Reproductive isolation was also found between hybrid populations, suggesting differential sorting of parental traits across populations. Finally, hybrid populations derived from three species were more likely to evolve reproductive isolation than those derived from two species, supporting arguments that hybridization-supplied genetic diversity can lead to the evolution of novel "adaptive systems" and promote speciation. These results illustrate when hybridization and admixture is expected to promote hybrid speciation. Whether homoploid hybrid speciation is a common speciation mechanism in general remains an outstanding empirical question (Comeault, 2018).

    Tissue-specific cis-regulatory divergence implicates eloF in inhibiting interspecies mating in Drosophila

    Reproductive isolation is a key component of speciation. In many insects, a major driver of this isolation is cuticular hydrocarbon pheromones, which help to identify potential intraspecific mates. When the distributions of related species overlap, there may be strong selection on mate choice for intraspecific partners because interspecific hybridization carries significant fitness costs. Drosophila has been a key model for the investigation of reproductive isolation; although both male and female mate choices have been extensively investigated, the genes underlying species recognition remain largely unknown. To explore the molecular mechanisms underlying Drosophila speciation, tissue-specific cis-regulatory divergence was identifed using RNA sequencing (RNA-seq) in D. simulans x D. sechellia hybrids. By focusing on cis-regulatory changes specific to female oenocytes, the tissue that produces cuticular hydrocarbons, this study identified a small number of candidate genes were identified. One of these, the fatty acid elongase eloF, broadly affects the hydrocarbons present on D. sechellia and D. melanogaster females, as well as the propensity of D. simulans males to mate with them. Therefore, cis-regulatory changes in eloF may be a major driver in the sexual isolation of D. simulans from multiple other species. The RNA-seq approach proved to be far more efficient than quantitative trait locus (QTL) mapping in identifying candidate genes; the same framework can be used to pinpoint candidate drivers of cis-regulatory divergence in traits differing between any interfertile species (Combs, 2018).

    Gene flow mediates the role of sex chromosome meiotic drive during complex speciation

    During speciation, sex chromosomes often accumulate interspecific genetic incompatibilities faster than the rest of the genome. The drive theory posits that sex chromosomes are susceptible to recurrent bouts of meiotic drive and suppression, causing the evolutionary build-up of divergent cryptic sex-linked drive systems and, incidentally, genetic incompatibilities. To assess the role of drive during speciation, this study combined high-resolution genetic mapping of X-linked hybrid male sterility with population genomics analyses of divergence and recent gene flow between the fruitfly species, Drosophila mauritiana and D. simulans. The findings reveal a high density of genetic incompatibilities and a corresponding dearth of gene flow on the X chromosome. Surprisingly, a known drive element recently migrated between species and, rather than contributing to interspecific divergence, caused a strong reduction in local sequence divergence, undermining the evolution of hybrid sterility. Gene flow can therefore mediate the effects of selfish genetic elements during speciation (Meiklejohn, 2018).

    Pigmentation and fitness trade-offs through the lens of artificial selection.

    Pigmentation is a classic phenotype that varies widely and adaptively in nature both within and among taxa. Genes underlying pigmentation phenotype are highly pleiotropic, creating the potential for functional trade-offs. However, the basic tenets of this trade-off hypothesis with respect to life-history traits have not been directly addressed. In natural populations of Drosophila melanogaster, the degree of melanin pigmentation covaries with fecundity and several other fitness traits. To examine correlations and potential trade-offs associated with variation in pigmentation, replicate outbred populations were selected for extreme pigmentation phenotypes. Replicate populations responded rapidly to the selection regime and after 100 generations of artificial selection were phenotyped for pigmentation as well as the two basic fitness parameters of fecundity and longevity. The data demonstrate that selection on pigmentation resulted in a significant shift in both fecundity and longevity profiles. Selection for dark pigmentation resulted in greater fecundity and no pronounced change in longevity, whereas selection for light pigmentation decreased longevity but did not affect fecundity. These results indicate the pleiotropic nature of alleles underlying pigmentation phenotype and elucidate possible trade-offs between pigmentation and fitness traits that may shape patterns of phenotypic variation in natural populations (Rajpurohit, 2017).

    Evolutionary history of LTR-retrotransposons among 20 Drosophila species

    The presence of transposable elements (TEs) in genomes is known to explain in part the variations of genome sizes among eukaryotes. Even among closely related species, the variation of TE amount may be striking, as for example between the two sibling species, Drosophila melanogaster and D. simulans. However, not much is known concerning the TE content and dynamics among other Drosophila species. The sequencing of several Drosophila genomes, covering the two subgenus Sophophora and Drosophila, revealed a large variation of the repeat content among these species but no much information is known concerning their precise TE content. The identification of some consensus sequences of TEs from the various sequenced Drosophila species allowed to get an idea concerning their variety in term of diversity of superfamilies but the used classification remains very elusive and ambiguous. This study focused on LTR-retrotransposons because they represent the most widely represented class of TEs in the Drosophila genomes, describing the phylogenetic relationship of each LTR-retrotransposon family described in 20 Drosophila species, computing their proportion in their respective genomes and identifying several new cases of horizontal transfers. CAll these results give a clearer view on the evolutionary history of LTR retrotransposons among Drosophila that seems to be mainly driven by vertical transmissions although the implications of horizontal transfers, losses and intra-specific diversification are clearly also at play (Bargues, 2017).

    Retrotransposons are the major contributors to the expansion of the Drosophila ananassae Muller F element

    The discordance between genome size and the complexity of eukaryotes can partly be attributed to differences in repeat density. The Muller F element (~5.2 Mb) is the smallest chromosome in Drosophila melanogaster, but it is substantially larger (>18.7 Mb) in Drosophila ananassae. To identify the major contributors to the expansion of the F element and to assess their impact, the genome sequence was improved and the genes in a 1.4 Mb region of the D. ananassae F element, and a 1.7 Mb region from the D element were annotated for comparison. Transposons (particularly LTR and LINE retrotransposons) were found to be major contributors to this expansion (78.6%), while Wolbachia sequences integrated into the D. ananassae genome are minor contributors (0.02%). Both D. melanogaster and D. ananassae F element genes exhibit distinct characteristics compared to D element genes (e.g., larger coding spans, larger introns, more coding exons, lower codon bias), but these differences are exaggerated in D. ananassae. Compared to D. melanogaster, the codon bias observed in D. ananassae F element genes can primarily be attributed to mutational biases instead of selection. The 5' ends of F element genes in both species are enriched in H3K4me2 while the coding spans are enriched in H3K9me2. Despite differences in repeat density and gene characteristics, D. ananassae F element genes show a similar range of expression levels compared to genes in euchromatic domains. This study improves understanding of how transposons can affect genome size and how genes can function within highly repetitive domains (Leung, 2017).

    Dissecting the satellite DNA landscape in three cactophilic Drosophila sequenced genomes

    Eukayote genomes are replete with repetitive DNAs. This class includes tandemly repeated satellite DNAs (satDNA) which are among the most abundant, fast evolving (yet poorly studied) genomic components. This study used high throughput sequencing data from three cactophilic Drosophila species, D. buzzatii, D. seriema and D. mojavensis, to access and study their whole satDNA landscape. Five satDNAs were identified, three previously described (pBuM, DBC-150 and CDSTR198) and two novel ones (CDSTR138 and CDSTR130). Only pBuM is shared among all three species. The satDNA repeat length falls within only two classes, between 130-200bp or between 340-390bp. FISH on metaphase and polytene chromosomes revealed the presence of satDNA arrays in at least one of the following genomic compartments: centromeric, telomeric, subtelomeric or dispersed along euchromatin. The chromosomal distribution ranges from a single chromosome to almost all chromosomes of the complement. Interspersion were revealed between pBuM and CDSTR130 in the microchromosomes of D. mojavensis. Phylogenetic analyses showed that the pBuM satDNA underwent concerted evolution at both interspecific and intraspecific levels. Based on RNAseq data, transcription activity was found for pBuM (in D. mojavensis) and CDSTR198 (in D. buzzatii) in all five analyzed developmental stages, most notably in pupae and adult males. These data revealed that cactophilic Drosophila present the lowest amount of satDNAs (1.9% to 2.9%) within the Drosophila genus reported so far (de Lima, 2017).

    Pervasive epigenetic effects of Drosophila euchromatic transposable elements impact their evolution

    Transposable elements (TEs) are widespread genomic parasites, and their evolution has remained a critical question in evolutionary genomics. This paper describes a study of the relatively unexplored epigenetic impacts of TEs and provides the first genome-wide quantification of such effects in D. melanogaster and D. simulans. Surprisingly, the spread of repressive epigenetic marks (histone H3K9me2) to nearby DNA occurs at >50% of euchromatic TEs, and can extend up to 20 kb. This results in differential epigenetic states of genic alleles and, in turn, selection against TEs. Interestingly, the lower TE content in D. simulans compared to D. melanogaster correlates with stronger epigenetic effects of TEs and higher levels of host genetic factors known to promote epigenetic silencing. This study demonstrates that the epigenetic effects of euchromatic TEs, and host genetic factors modulating such effects, play a critical role in the evolution of TEs both within and between species (Lee, 2017).

    Genomic analysis of P elements in natural populations of Drosophila melanogaster

    The Drosophila melanogaster P transposable element provides one of the best cases of horizontal transfer of a mobile DNA sequence in eukaryotes. Invasion of natural populations by the P element has led to a syndrome of phenotypes known as P-M hybrid dysgenesis that emerges when strains differing in their P element composition mate and produce offspring. This study compared estimates of genomic P element content with gonadal dysgenesis phenotypes for isofemale strains obtained from three worldwide populations of D. melanogaster to illuminate the molecular basis of natural variation in cytotype status. P element abundance estimated from genome sequences of isofemale strains is shown to be highly correlated across different bioinformatics approaches, but abundance estimates are sensitive to method and filtering strategies as well as incomplete inbreeding of isofemale strains. P element content was found to vary significantly across populations, with strains from a North American population having fewer P elements but a higher proportion of full-length elements than strains from populations sampled in Europe or Africa. Despite these geographic differences in P element abundance and structure, neither the number of P elements nor the ratio of full-length to internally-truncated copies is strongly correlated with the degree of gonadal dysgenesis exhibited by an isofemale strain. Thus, variation in P element abundance and structure across different populations does not necessarily lead to corresponding geographic differences in gonadal dysgenesis phenotypes. Finally, it was confirmed that population differences in the abundance and structure of P elements that are observed from isofemale lines can also be observed in pool-seq samples from the same populations. This work supports the view that genomic P element content alone is not sufficient to explain variation in gonadal dysgenesis across strains of D. melanogaster, and informs future efforts to decode the genomic basis of geographic and temporal differences in P element induced phenotypes (Bergman, 2017).

    Phenotypic plasticity through transcriptional regulation of the evolutionary hotspot gene tan in Drosophila melanogaster

    Phenotypic plasticity is the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. Phenotypic plasticity can be adaptive. Furthermore, it is thought to facilitate evolution. Although phenotypic plasticity is a widespread phenomenon, its molecular mechanisms are only beginning to be unravelled. Environmental conditions can affect gene expression through modification of chromatin structure, mainly via histone modifications, nucleosome remodelling or DNA methylation, suggesting that phenotypic plasticity might partly be due to chromatin plasticity. As a model of phenotypic plasticity, abdominal pigmentation was studied of Drosophila melanogaster females, which is temperature sensitive. Abdominal pigmentation is indeed darker in females grown at 18 degrees C than at 29 degrees C. This phenomenon is thought to be adaptive as the dark pigmentation produced at lower temperature increases body temperature. This study showed that temperature modulates the expression of tan (t), a pigmentation gene involved in melanin production. t is expressed 7 times more at 18 ° C than at 29 &176; C in female abdominal epidermis. Genetic experiments show that modulation of t expression by temperature is essential for female abdominal pigmentation plasticity. Temperature modulates the activity of an enhancer of t without modifying compaction of its chromatin or level of the active histone mark H3K27ac. By contrast, the active mark H3K4me3 on the t promoter is strongly modulated by temperature. The H3K4 methyl-transferase involved in this process is likely Trithorax, since it regulates t expression and the H3K4me3 level on the t promoter and also participates in female pigmentation and its plasticity. Interestingly, t was previously shown to be involved in inter-individual variation of female abdominal pigmentation in Drosophila melanogaster, and in abdominal pigmentation divergence between Drosophila species. Sensitivity of t expression to environmental conditions might therefore give more substrate for selection, explaining why this gene has frequently been involved in evolution of pigmentation (Gibert, 2016).

    On the long-term stability of clines in some metabolic genes in Drosophila melanogaster

    Very little information exists for long-term changes in genetic variation in natural populations. This study took the unique opportunity to compare a set of data for SNPs in 15 metabolic genes from eastern US collections of Drosophila melanogaster that span a large latitudinal range and represent two collections separated by 12 to 13 years. This was expanded to a 22-year interval for the Adh gene and approximately 30 years for the G6pd and Pgd genes. During these intervals, five genes showed a statistically significant change in average SNP allele frequency corrected for latitude. While much remains unchanged, five genes were seen that wx latitudinal clines have been lost or gained and two where the slope significantly changes. The long-term frequency shift towards a southern favored Adh S allele reported in Australia populations is not observed in the eastern US over a period of 21 years. There is no general pattern of southern-favored or northern-favored alleles increasing in frequency across the genes. This observation points to the fluid nature of some allelic variation over this time period and the action of selective responses or migration that may be more regional than uniformly imposed across the cline (Cogni, 2017).


    Ahmed-Braimah, Y. H., Unckless, R. L. and Clark, A. G. (2017). Evolutionary dynamics of male reproductive genes in the Drosophila virilis subgroup. G3 (Bethesda) [Epub ahead of print]. PubMed ID: 28739599

    Akhund-Zade, J., Bergland, A. O., Crowe, S. O. and Unckless, R. L. (2016). The genetic basis of natural variation in Drosophila (Diptera: Drosophilidae) virgin egg retention. J Insect Sci [Epub ahead of print]. PubMed ID: 28042107

    Allen, S. L., Bonduriansky, R., Sgro, C. M. and Chenoweth, S. F. (2017). Sex-biased transcriptome divergence along a latitudinal gradient. Mol Ecol [Epub ahead of print]. PubMed ID: 28100025

    Almeida, F. C. and DeSalle, R. (2016). Genetic differentiation and adaptive evolution at reproductive loci in incipient Drosophila species. J Evol Biol. PubMed ID: 27883252

    Alton, L. A., Condon, C., White, C. R. and Angilletta, M. J. (2016). Colder environments did not select for a faster metabolism during experimental evolution of Drosophila melanogaster. Evolution [Epub ahead of print]. PubMed ID: 27757954

    Amei, A. and Zhou, S. (2019). Inferring the distribution of selective effects from a time inhomogeneous model. PLoS One 14(1): e0194709. PubMed ID: 30657757

    Appel, M., Scholz, C. J., Kocabey, S., Savage, S., Konig, C. and Yarali, A. (2016). Independent natural genetic variation of punishment- versus relief-memory. Biol Lett 12(12). PubMed ID: 28003518

    Arguello, J. R., Cardoso-Moreira, M., Grenier, J. K., Gottipati, S., Clark, A. G. and Benton, R. (2016). Extensive local adaptation within the chemosensory system following Drosophila melanogaster's global expansion. Nat Commun 7: ncomms11855. PubMed ID: 27292132

    Assis, R. (2016). Transcriptional interference promotes rapid expression divergence of Drosophila nested genes. Genome Biol Evol 8(10):3149-3158. PubMed ID: 27664180

    Ballinger, M. J. and Perlman, S. J. (2017). Generality of toxins in defensive symbiosis: Ribosome-inactivating proteins and defense against parasitic wasps in Drosophila. PLoS Pathog 13(7): e1006431. PubMed ID: 28683136

    Bargues, N. and Lerat, E. (2017). Evolutionary history of LTR-retrotransposons among 20 Drosophila species. Mob DNA 8: 7. PubMed ID: 28465726

    Bastide, H., Lange, J. D., Lack, J. B., Yassin, A. and Pool, J. E. (2016). A variable genetic architecture of melanic evolution in Drosophila melanogaster. Genetics [Epub ahead of print]. PubMed ID: 27638419

    Battlay, P., Schmidt, J. M., Fournier-Level, A. and Robin, C. (2016). Genomic and transcriptomic associations identify a new insecticide resistance phenotype for the selective sweep at the Cyp6g1 locus of Drosophila melanogaster. G3 (Bethesda) [Epub ahead of print]. PubMed ID: 27317781

    Battlay, P., Leblanc, P. B., Green, L., Garud, N. R., Schmidt, J. M., Fournier-Level, A. and Robin, C. (2018). Structural variants and selective sweep foci contribute to insecticide resistance in the Drosophila genetic reference panel. G3 (Bethesda). PubMed ID: 30190421

    Baudouin-Gonzalez, L., Santos, M. A., Tempesta, C., Sucena, E., Roch, F. and Tanaka, K. (2017). Diverse cis-regulatory mechanisms contribute to expression evolution of tandem gene duplicates. Mol Biol Evol [Epub ahead of print]. PubMed ID: 28961967

    Bergman, C. M., Han, S., Nelson, M. G., Bondarenko, V. and Kozeretska, I. (2017). Genomic analysis of P elements in natural populations of Drosophila melanogaster. PeerJ 5: e3824. PubMed ID: 28929030

    Betini, G. S., McAdam, A. G., Griswold, C. K. and Norris, D. R. (2017). A fitness trade-off between seasons causes multigenerational cycles in phenotype and population size. Elife 6 [Epub ahead of print]. PubMed ID: 28164780

    Branstetter, M. G., Danforth, B. N., Pitts, J. P., Faircloth, B. C., Ward, P. S., Buffington, M. L., Gates, M. W., Kula, R. R. and Brady, S. G. (2017). Phylogenomic Insights into the Evolution of Stinging Wasps and the Origins of Ants and Bees. Curr Biol 27(7): 1019-1025. PubMed ID: 28376325

    Burke, M.K., Barter, T.T., Cabral, L.G., Kezos, J.N., Phillips, M.A., Rutledge, G.A., Phung, K.H., Chen, R.H., Nguyen, H.D., Mueller, L.D. and Rose, M.R. (2016). Rapid divergence and convergence of life-history in experimentally evolved Drosophila melanogaster. Evolution [Epub ahead of print]. PubMed ID: 27431916

    Caizzi, R., Moschetti, R., Piacentini, L., Fanti, L., Marsano, R. M. and Dimitri, P. (2016). Comparative genomic analyses provide new insights into the evolutionary dynamics of heterochromatin in Drosophila. PLoS Genet 12: e1006212. PubMed ID: 27513559

    Calvo-Martin, J. M., Papaceit, M. and Segarra, C. (2017). Molecular population genetics of the Polycomb genes in Drosophila subobscura. PLoS One 12(9): e0185005. PubMed ID: 28910411

    Castro Vargas, C., Richmond, M. P., Ramirez Loustalot Laclette, M. and Markow, T. A. (2017). Early events in speciation: Cryptic species of Drosophila aldrichi. Ecol Evol 7(12): 4220-4228. PubMed ID: 28649335

    Catalan, A., Macias-Munoz, A. and Briscoe, A. D. (2018). Evolution of sex-biased gene expression and dosage compensation in the eye and brain of Heliconius butterflies. Mol Biol Evol. PubMed ID: 29931127

    Castellano, D., James, J. and Eyre-Walker, A. (2018). Nearly neutral evolution across the Drosophila melanogaster genome. Mol Biol Evol 35(11): 2685-2694. PubMed ID: 30418639

    Chang, C.H. and Larracuente, A.M. (2017). Genomic changes following the reversal of a Y chromosome to an autosome in Drosophila pseudoobscura. Evolution [Epub ahead of print]. PubMed ID: 28322435

    Cheng, C. and Kirkpatrick, M. (2016). Sex-specific selection and sex-biased gene expression in humans and flies. PLoS Genet 12: e1006170. PubMed ID: 27658217

    Choi, D. S., Park, J. S., Kim, M. J., Kim, J. S., Jeong, S. Y., Jeong, J. S., Park, J. and Kim, I. (2017). Geographic variation in the spotted-wing drosophila, Drosophila suzukii (Diptera: Drosophilidae), based on mitochondrial DNA sequences. Mitochondrial DNA A DNA Mapp Seq Anal: 1-11. PubMed ID: 28129731

    Clifton, B. D., Librado, P., Yeh, S. D., Solares, E., Real, D., Jayasekera, S., Zhang, W., Shi, M., Park, R., Magie, R., Ma, H. C., Xia, X. Q., Marco, A., Rozas, J. and Ranz, J. M. (2016). Rapid functional and sequence differentiation of a tandemly-repeated species-specific multigene family in Drosophila. Mol Biol Evol [Epub ahead of print]. PubMed ID: 27702774

    Cogni, R., Kuczynski, K., Koury, S., Lavington, E., Behrman, E. L., O'Brien, K. R., Schmidt, P. S. and Eanes, W. F. (2017). On the long-term stability of clines in some metabolic genes in Drosophila melanogaster. Sci Rep 7: 42766. PubMed ID: 28220806

    Comeault, A. A. and Matute, D. R. (2018). Genetic divergence and the number of hybridizing species affect the path to homoploid hybrid speciation. Proc Natl Acad Sci U S A 115(39): 9761-9766. PubMed ID: 30209213

    Combs, P. A., Krupp, J. J., Khosla, N. M., Bua, D., Petrov, D. A., Levine, J. D. and Fraser, H. B. (2018). Tissue-specific cis-regulatory divergence implicates eloF in inhibiting interspecies mating in Drosophila. Curr Biol 28(24):3969-3975. PubMed ID: 30503619

    Cooper, B. S., Sedghifar, A., Nash, W. T., Comeault, A. A. and Matute, D. R. (2018). A Maladaptive Combination of Traits Contributes to the Maintenance of a Drosophila Hybrid Zone. Curr Biol. PubMed ID: 30174184

    Corbett-Detig, R. and Nielsen, R. (2017). A hidden markov model approach for simultaneously estimating local ancestry and admixture time using next generation sequence data in samples of arbitrary ploidy. PLoS Genet 13(1): e1006529. PubMed ID: 28045893

    Corbett-Detig, R. (2016). Selection on inversion breakpoints favors proximity to pairing sensitive sites in Drosophila melanogaster. Genetics [Epub ahead of print]. PubMed ID: 27343234

    Craddock, E. M., Kambysellis, M. P., Franchi, L., Francisco, P., Grey, M., Hutchinson, A., Nanhoo, S. and Antar, S. (2018). Ultrastructural variation and adaptive evolution of the ovipositor in the endemic Hawaiian Drosophilidae. J Morphol. PubMed ID: 30397938

    Crowley-Gall, A., Date, P., Han, C., Rhodes, N., Andolfatto, P., Layne, J. E. and Rollmann, S. M. (2016). Population differences in olfaction accompany host shift in Drosophila mojavensis. Proc Biol Sci 283. PubMed ID: 27581882

    Croze, M., Wollstein, A., Bozicevic, V., Zivkovic, D., Stephan, W. and Hutter, S. (2017). BMC Evol Biol 17(1): 15. PubMed ID: 28086750

    Dalla, S. and Dobler, S. (2016). Gene duplications circumvent trade-offs in enzyme function: Insect adaptation to toxic host plants. Evolution [Epub ahead of print]. PubMed ID: 27683239

    de Lima, L. G., Svartman, M. and Kuhn, G. C. S. (2017). Dissecting the satellite DNA landscape in three cactophilic Drosophila sequenced genomes. G3 (Bethesda) [Epub ahead of print]. PubMed ID: 28659292

    Denecke, S., Fusetto, R., Martelli, F., Giang, A., Battlay, P., Fournier-Level, A., RA, O. H. and Batterham, P. (2017). Multiple P450s and variation in neuronal genes underpins the response to the insecticide Imidacloprid in a population of Drosophila melanogaster. Sci Rep 7(1): 11338. PubMed ID: 28900129

    Diaz, F., Allan, C. W. and Matzkin, L. M. (2018). Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila. BMC Evol Biol 18(1): 144. PubMed ID: 30236055

    Dugand, R. J., Kennington, W. J. and Tomkins, J. L. (2018). Evolutionary divergence in competitive mating success through female mating bias for good genes. Sci Adv 4(5): eaaq0369. PubMed ID: 29806021

    Dupim, E. G., Goldstein, G., Vanderlinde, T., Vaz, S. C., Krsticevic, F., Bastos, A., Pinhao, T., Torres, M., David, J. R., Vilela, C. R. and Carvalho, A. B. (2018). An investigation of Y chromosome incorporations in 400 species of Drosophila and related genera. PLoS Genet 14(11): e1007770. PubMed ID: 30388103

    Duxbury, E. M. L., Rostant, W. G. and Chapman, T. (2017). Manipulation of feeding regime alters sexual dimorphism for lifespan and reduces sexual conflict in Drosophila melanogaster. Proc Biol Sci 284(1854). PubMed ID: 28469030

    Early, A. M., Arguello, J. R., Cardoso-Moreira, M., Gottipati, S., Grenier, J. K. and Clark, A. G. (2016). Survey of global genetic diversity within the Drosophila immune system. Genetics [Epub ahead of print]. PubMed ID: 27815361

    Ebadi, H., Perry, M., Short, K., Klemm, K., Desplan, C., Stadler, P. F. and Mehta, A. (2018). Patterning the insect eye: From stochastic to deterministic mechanisms. PLoS Comput Biol 14(11): e1006363. PubMed ID: 30439954

    Elliott, K. H., Betini, G. S. and Norris, D. R. (2017). Fear creates an Allee effect: experimental evidence from seasonal populations. Proc Biol Sci 284(1857). PubMed ID: 28659452

    Elyashiv, E., Sattath, S., Hu, T. T., Strutsovsky, A., McVicker, G., Andolfatto, P., Coop, G. and Sella, G. (2016). A genomic map of the effects of linked selection in Drosophila. PLoS Genet 12: e1006130. PubMed ID: 27536991

    Faria, V.G., Martins, N.E., Magalhães, S., Paulo, T.F., Nolte, V., Schlötterer, C., Sucena, É and Teixeira, L. (2016). Drosophila adaptation to viral infection through defensive symbiont evolution. PLoS Genet 12: e1006297. PubMed ID: 27684942

    Filice, D. C. S. and Long, T. A. F. (2018). Genetic trade-offs between male reproductive traits in Drosophila melanogaster. Biol Lett 14(10). PubMed ID: 30333262

    Fraisse, C., Puixeu Sala, G. and Vicoso, B. (2018). Pleiotropy modulates the efficacy of selection in Drosophila melanogaster. Mol Biol Evol. PubMed ID: 30590559

    Fuller, Z. L., Haynes, G. D., Richards, S. and Schaeffer, S. W. (2016). Genomics of natural populations: How differentially expressed genes shape the evolution of chromosomal inversions in Drosophila pseudoobscura. Genetics [Epub ahead of print]. PubMed ID: 27401754

    Gerland, T.A., Sun, B., Smialowski, P., Lukacs, A., Thomae, A.W. and Imhof, A. (2017). The Drosophila speciation factor HMR localizes to genomic insulator sites. PLoS One 12: e0171798. PubMed ID: 28207793

    Garlovsky, M. D. and Snook, R. R. (2018). Persistent postmating, prezygotic reproductive isolation between populations. Ecol Evol 8(17): 9062-9073. PubMed ID: 30271566

    Gibert, J. M., Mouchel-Vielh, E., De Castro, S. and Peronnet, F. (2016). Phenotypic plasticity through transcriptional regulation of the evolutionary hotspot gene tan in Drosophila melanogaster. PLoS Genet 12: e1006218. PubMed ID: 27508387

    Gibert, J. M., Mouchel-Vielh, E. and Peronnet, F. (2018). Pigmentation pattern and developmental constraints: flight muscle attachment sites delimit the thoracic trident of Drosophila melanogaster. Sci Rep 8(1): 5328. PubMed ID: 29593305

    Goldman-Huertas, B., Mitchell, R. F., Lapoint, R. T., Faucher, C. P., Hildebrand, J. G. and Whiteman, N. K. (2015). Evolution of herbivory in Drosophilidae linked to loss of behaviors, antennal responses, odorant receptors, and ancestral diet. Proc Natl Acad Sci U S A 112(10): 3026-3031. PubMed ID: 25624509

    Griffin, P. C., Hangartner, S. B., Fournier-Level, A. and Hoffmann, A. A. (2016). Genomic trajectories to desiccation resistance: Convergence and divergence among replicate selected Drosophila lines. Genetics [Epub ahead of print]. PubMed ID: 28007884

    Graves, J. L., Jr., Hertweck, K. L., Phillips, M. A., Han, M. V., Cabral, L. G., Barter, T. T., Greer, L. F., Burke, M. K., Mueller, L. D. and Rose, M. R. (2017). Genomics of parallel experimental evolution in Drosophila. Mol Biol Evol [Epub ahead of print]. PubMed ID: 28087779

    Groth, B. R., Huang, Y., Monette, M. J. and Pool, J. E. (2018). Directional selection reduces developmental canalization against genetic and environmental perturbations in Drosophila wings. Evolution. PubMed ID: 29985527

    Gubala, A. M., Schmitz, J. F., Kearns, M. J., Vinh, T. T., Bornberg-Bauer, E., Wolfner, M. F. and Findlay, G. D. (2017). The goddard and saturn genes are essential for Drosophila male fertility and may have arisen de novo. Mol Biol Evol. PubMed ID: 28104747

    Gursky, V. V., Kozlov, K. N., Kulakovskiy, I. V., Zubair, A., Marjoram, P., Lawrie, D. S., Nuzhdin, S. V. and Samsonova, M. G. (2017). Translating natural genetic variation to gene expression in a computational model of the Drosophila gap gene regulatory network. PLoS One 12(9): e0184657. PubMed ID: 28898266

    Hackett, J. L., Wang, X., Smith, B. R. and Macdonald, S. J. (2016). Mapping QTL contributing to variation in posterior lobe morphology between strains of Drosophila melanogaster. PLoS One 11: e0162573. PubMed ID: 27606594

    Hanlon, S. L., Miller, D. E., Eche, S. and Hawley, R. S. (2018). Origin, composition, and structure of the supernumerary B chromosome of Drosophila melanogaster. Genetics. PubMed ID: 30249684

    Hart, M. L. I., Vu, B. L., Bolden, Q., Chen, K. T., Oakes, C. L., Zoronjic, L. and Meisel, R. P. (2018). Genes relocated between Drosophila chromosome arms evolve under relaxed selective constraints relative to non-relocated genes. J Mol Evol. PubMed ID: 29926120

    Harvanek, Z. M., Lyu, Y., Gendron, C. M., Johnson, J. C., Kondo, S., Promislow, D. E. L. and Pletcher, S. D. (2017). Perceptive costs of reproduction drive ageing and physiology in male Drosophila. Nat Ecol Evol 1(6): 152. PubMed ID: 28812624

    Hemmer, L. W. and Blumenstiel, J. P. (2016). Holding it together: rapid evolution and positive selection in the synaptonemal complex of Drosophila. BMC Evol Biol 16: 91. PubMed ID: 27150275

    Hickner, P. V., Rivaldi, C. L., Johnson, C. M., Siddappaji, M., Raster, G. J. and Syed, Z. (2016). The making of a pest: Insights from the evolution of chemosensory receptor families in a pestiferous and invasive fly, Drosophila suzukii. BMC Genomics 17: 648. PubMed ID: 27530109

    Highfill, C. A., Tran, J. H., Nguyen, S. K. T., Moldenhauer, T. R., Wang, X. and Macdonald, S. J. (2017). Naturally-segregating variation at Ugt86Dd contributes to nicotine resistance in Drosophila melanogaster. Genetics [Epub ahead of print]. PubMed ID: 28743761

    Hollis, B., Keller, L. and Kawecki, T. J. (2016). Sexual selection shapes development and maturation rates in Drosophila. Evolution [Epub ahead of print]. PubMed ID: 27883363

    Huang, Y. and Agrawal, A. F. (2016). Experimental evolution of gene expression and plasticity in alternative selective regimes. PLoS Genet 12: e1006336. PubMed ID: 27661078

    Izumitani, H. F., Kusaka, Y., Koshikawa, S., Toda, M. J. and Katoh, T. (2016). Phylogeography of the Subgenus Drosophila (Diptera: Drosophilidae): Evolutionary history of faunal divergence between the old and the new worlds. PLoS One 11: e0160051. PubMed ID: 27462734

    Jackson, B. C., Campos, J. L., Haddrill, P. R., Charlesworth, B. and Zeng, K. (2017). Variation in the intensity of selection on codon bias over time causes contrasting patterns of base composition evolution in Drosophila. Genome Biol Evol. PubMed ID: 28082609

    Jayaswal, V., Jimenez, J., Magie, R., Nguyen, K., Clifton, B., Yeh, S. and Ranz, J. M. (2018). A species-specific multigene family mediates differential sperm displacement in Drosophila melanogaster. Evolution 72(2): 399-403. PubMed ID: 29315521

    Jungreis, I., Chan, C. S., Waterhouse, R. M., Fields, G., Lin, M. F. and Kellis, M. (2016). Evolutionary dynamics of abundant stop codon readthrough. Mol Biol Evol [Epub ahead of print]. PubMed ID: 27604222

    Keais, G. L., Hanson, M. A., Gowen, B. E. and Perlman, S. J. (2017). X chromosome drive in a widespread Palearctic woodland fly, Drosophila testacea. J Evol Biol [Epub ahead of print]. PubMed ID: 28402000

    Khost, D. E., Eickbush, D. G. and Larracuente, A. M. (2017). Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster. Genome Res [Epub ahead of print]. PubMed ID: 28373483

    Kang, L., Garner, H. R., Price, D. K. and Michalak, P. (2017). Kang, L., Garner, H. R., Price, D. K. and Michalak, P. (2017). A test for gene flow among sympatric and allopatric Hawaiian picture-winged Drosophila. J Mol Evol [Epub ahead of print]. PubMed ID: 28492967

    Kursel, L. E. and Malik, H. S. (2017). Recurrent gene duplication leads to diverse repertoires of centromeric histones in Drosophila species. Mol Biol Evol [Epub ahead of print]. PubMed ID: 28333217

    Lack, J. B., Yassin, A., Sprengelmeyer, Q. D., Johanning, E. J., David, J. R. and Pool, J. E. (2016). Life history evolution and cellular mechanisms associated with increased size in high-altitude Drosophila. Ecol Evol 6: 5893-5906. PubMed ID: 27547363

    Lee, Y. C. G. and Karpen, G. H. (2017). Pervasive epigenetic effects of Drosophila euchromatic transposable elements impact their evolution. Elife 6. PubMed ID: 28695823

    Le Page, S., Sepil, I., Flintham, E., Pizzari, T., Carazo, P. and Wigby, S. (2017). Male relatedness and familiarity are required to modulate male-induced harm to females in Drosophila. Proc Biol Sci 284(1860). PubMed ID: 28794215

    Leung, W., et al. (2017). Retrotransposons are the major contributors to the expansion of the Drosophila ananassae Muller F element. G3 (Bethesda) 7(8):2439-2460. PubMed ID: 28667019

    Levin, T. C. and Malik, H. S. (2017). Rapidly evolving Toll-3/4 genes encode male-specific Toll-like receptors in Drosophila. Mol Biol Evol. PubMed ID: 28541576

    Lewis, S. H., Webster, C. L., Salmela, H. and Obbard, D. J. (2016). Repeated duplication of Argonaute2 is associated with strong selection and testis specialization in Drosophila. Genetics [Epub ahead of print]. PubMed ID: 27535930

    Librado, P. and Rozas, J. (2016). Weak polygenic selection drives the rapid adaptation of the chemosensory system: lessons from the upstream regions of the major gene families. Genome Biol Evol [Epub ahead of print]. PubMed ID: 27503297

    Lin, C. J., Hu, F., Dubruille, R., Vedanayagam, J., Wen, J., Smibert, P., Loppin, B. and Lai, E. C. (2018). The hpRNA/RNAi pathway is essential to resolve intragenomic conflict in the Drosophila male germline. Dev Cell 46(3): 316-326. PubMed ID: 30086302

    Lopez-Maestre, H., Carnelossi, E. A., Lacroix, V., Burlet, N., Mugat, B., Chambeyron, S., Carareto, C. M. and Vieira, C. (2017). Identification of misexpressed genetic elements in hybrids between Drosophila-related species. Sci Rep 7: 40618. PubMed ID: 28091568

    Mageeney, C. M., Kearse, M. G., Gershman, B. W., Pritchard, C. E., Colquhoun, J. M. and Ware, V. C. (2018). Functional interplay between ribosomal protein paralogues in the eRpL22 family in Drosophila melanogaster. Fly (Austin): 1-21. PubMed ID: 30465696

    Martinez, J., Tolosana, I., Ok, S., Smith, S., Snoeck, K., Day, J. P. and Jiggins, F. (2017). Symbiont strain is the main determinant of variation in Wolbachia-mediated protection against viruses across Drosophila species. Mol Ecol [Epub ahead of print]. PubMed ID: 28464440

    Mathur, V. and Schmidt, P.S. (2016). Adaptive patterns of phenotypic plasticity in laboratory and field environments in Drosophila melanogaster. Evolution [Epub ahead of print]. PubMed ID: 27925178

    Medina, P., Thornlow, B., Nielsen, R. and Corbett-Detig, R. (2018). Estimating the timing of multiple admixture pulses during local ancestry inference. Genetics. PubMed ID: 30206187

    Meiklejohn, C. D., Landeen, E. L., Gordon, K. E., Rzatkiewicz, T., Kingan, S. B., Geneva, A. J., Vedanayagam, J. P., Muirhead, C. A., Garrigan, D., Stern, D. L. and Presgraves, D. C. (2018). Gene flow mediates the role of sex chromosome meiotic drive during complex speciation. Elife 7. PubMed ID: 30543325

    Michalak, P., Kang, L., Sarup, P. M., Schou, M. F. and Loeschcke, V. (2017). Nucleotide diversity inflation as a genome-wide response to experimental lifespan extension in Drosophila melanogaster. BMC Genomics 18(1): 84. PubMed ID: 28088192

    Mittleman, B. E., Manzano-Winkler, B., Hall, J. B., Korunes, K. L. and Noor, M. A. (2017). The large X-effect on secondary sexual characters and the genetics of variation in sex comb tooth number in Drosophila subobscura. Ecol Evol 7(2): 533-540. PubMed ID: 28116050

    Morimoto, J., Ponton, F., Tychsen, I., Cassar, J. and Wigby, S. (2017). Interactions between the developmental and adult social environments mediate group dynamics and offspring traits in Drosophila melanogaster. Sci Rep 7(1): 3574. PubMed ID: 28620201

    Muto, L., Kamimura, Y., Tanaka, K. M. and Takahashi, A. (2018). An innovative ovipositor for niche exploitation impacts genital coevolution between sexes in a fruit-damaging Drosophila. Proc Biol Sci 285(1887). PubMed ID: 30257912

    Nourmohammad, A., Rambeau, J., Held, T., Kovacova, V., Berg, J. and Lassig, M. (2017). Adaptive evolution of gene expression in Drosophila. Cell Rep 20(6): 1385-1395. PubMed ID: 28793262

    Nozawa, M., Onizuka, K., Fujimi, M., Ikeo, K. and Gojobori, T. (2016). Accelerated pseudogenization on the neo-X chromosome in Drosophila miranda. Nat Commun 7: 13659. PubMed ID: 27897175

    O'Brien, E. K., Higgie, M., Reynolds, A., Hoffmann, A. A. and Bridle, J. R. (2017). Testing for local adaptation and evolutionary potential along altitudinal gradients in rainforest Drosophila: beyond laboratory estimates. Glob Chang Biol. PubMed ID: 28070978

    Osada, N., Miyagi, R. and Takahashi, A. (2017). Cis- and trans-regulatory effects on gene expression in a natural population of Drosophila melanogaster. Genetics [Epub ahead of print]. PubMed ID: 28615283

    Parker, D. J., Wiberg, R. A. W., Trivedi, U., Tyukmaeva, V. I., Gharbi, K., Butlin, R. K., Hoikkala, A., Kankare, M. and Ritchie, M. G. (2018). Inter- and intra-specific genomic divergence in Drosophila montana shows evidence for cold adaptation. Genome Biol Evol. PubMed ID: 30010752

    Pegoraro, M., Zonato, V., Tyler, E. R., Fedele, G., Kyriacou, C. P. and Tauber, E. (2017). Geographical analysis of diapause inducibility in European Drosophila melanogaster populations. J Insect Physiol 98: 238-244. PubMed ID: 28131702

    Peters, R. S., et al. (2017). Evolutionary History of the Hymenoptera. Curr Biol 27(7): 1013-1018. PubMed ID: 28343967

    Phillips, M. A., Long, A. D., Greenspan, Z. S., Greer, L. F., Burke, M. K., Villeponteau, B., Matsagas, K. C., Rizza, C. L., Mueller, L. D. and Rose, M. R. (2016). Genome-wide analysis of long-term evolutionary domestication in Drosophila melanogaster. Sci Rep 6: 39281. PubMed ID: 28004838

    Phillips, M. A., Rutledge, G. A., Kezos, J. N., Greenspan, Z. S., Talbott, A., Matty, S., Arain, H., Mueller, L. D., Rose, M. R. and Shahrestani, P. (2018). Effects of evolutionary history on genome wide and phenotypic convergence in Drosophila populations. BMC Genomics 19(1): 743. PubMed ID: 30305018

    Pischedda, A. and Chippindale, A. K. (2017). Direct benefits of choosing a high fitness mate can offset the indirect costs associated with intralocus sexual conflict. J Evolution [Epub ahead of print]. PubMed ID: 28369895

    Porcelli, D., Westram, A. M., Pascual, M., Gaston, K. J., Butlin, R. K. and Snook, R. R. (2016). Gene expression clines reveal local adaptation and associated trade-offs at a continental scale. Sci Rep 6: 32975. PubMed ID: 27599812

    Puerma, E., Orengo, D. J. and Aguade, M. (2016). The origin of chromosomal inversions as a source of segmental duplications in the Sophophora subgenus of Drosophila. Sci Rep 6: 30715. PubMed ID: 27470196

    Rajpurohit, S., Richardson, R., Dean, J., Vazquez, R., Wong, G. and Schmidt, P. S. (2016). Pigmentation and fitness trade-offs through the lens of artificial selection. Biol Lett 12(10). PubMed ID: 28120808

    Rice, G., Barmina, O., Hu, K. and Kopp, A. (2018). Evolving doublesex expression correlates with the origin and diversification of male sexual ornaments in the Drosophila immigrans species group. Evol Dev 20(2):78-88. PubMed ID: 29372584

    Rogers, R. L., Shao, L. and Thornton, K. R. (2017). Tandem duplications lead to novel expression patterns through exon shuffling in Drosophila yakuba. PLoS Genet 13(5): e1006795. PubMed ID: 28531189

    Rostant, W. G., Bowyer, J., Coupland, J., Facey, J., Hosken, D. J. and Wedell, N. (2017). Pleiotropic effects of DDT resistance on male size and behaviour. Behav Genet. PubMed ID: 28466236

    Roux, C., Fraisse, C., Romiguier, J., Anciaux, Y., Galtier, N. and Bierne, N. (2016). Shedding light on the grey zone of speciation along a continuum of genomic divergence. PLoS Biol 14(12): e2000234. PubMed ID: 28027292

    Said, I., Byrne, A., Serrano, V., Cardeno, C., Vollmers, C. and Corbett-Detig, R. (2018). Linked genetic variation and not genome structure causes widespread differential expression associated with chromosomal inversions. Proc Natl Acad Sci U S A 115(21): 5492-5497. Pubmed ID: 29735663

    Saint-Leandre, B., Clavereau, I., Hua-Van, A. and Capy, P. (2017). Transcriptional polymorphism of piRNA regulatory genes underlies the mariner activity in D. simulans testes. Mol Ecol [Epub ahead of print]. PubMed ID: 28401606

    Salazar-Jaramillo, L., Jalvingh, K. M., de Haan, A., Kraaijeveld, K., Buermans, H. and Wertheim, B. (2017). Inter- and intra-species variation in genome-wide gene expression of Drosophila in response to parasitoid wasp attack. BMC Genomics 18(1): 331. PubMed ID: 28449654

    Schmidt, J. M., Battlay, P., Gledhill-Smith, R. S., Good, R. T., Lumb, C., Fournier-Level, A. and Robin, C. (2017). Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel. Genetics [Epub ahead of print]. PubMed ID: 28935691

    Schor, I. E., Degner, J. F., Harnett, D., Cannavo, E., Casale, F. P., Shim, H., Garfield, D. A., Birney, E., Stephens, M., Stegle, O. and Furlong, E. E. (2017). Promoter shape varies across populations and affects promoter evolution and expression noise. Nat Genet. PubMed ID: 28191888

    Seeholzer, L. F., Seppo, M., Stern, D. L. and Ruta, V. (2018). Evolution of a central neural circuit underlies Drosophila mate preferences. Nature 559(7715):564-569. PubMed ID: 29995860

    Shahandeh, M. P., Pischedda, A. and Turner, T. L. (2017). Male mate choice via cuticular hydrocarbon pheromones drives reproductive isolation between Drosophila species. Evolution 72(1):123-135. PubMed ID: 29098691

    Shindey, R., Varma, V., Nikhil, K. L. and Sharma, V. K. (2017). Evolution of circadian rhythms in Drosophila melanogaster populations reared in constant light and dark regimes for over 330 generations. Chronobiol Int: 1-14. PubMed ID: 28156168

    Shropshire, J. D., On, J., Layton, E. M., Zhou, H. and Bordenstein, S. R. (2018). One prophage WO gene rescues cytoplasmic incompatibility in Drosophila melanogaster. Proc Natl Acad Sci U S A. PubMed ID: 29686091

    Siddiq, M. A., Loehlin, D. W., Montooth, K. L. and Thornton, J. W. (2017). Experimental test and refutation of a classic case of molecular adaptation in Drosophila melanogaster. Nat Ecol Evol 1(2): 25. PubMed ID: 28812605

    Signor, S.A., Liu, Y., Rebeiz, M. and Kopp, A. (2016). Genetic convergence in the evolution of male-specific color patterns in Drosophila. Curr Biol [Epub ahead of print]. PubMed ID: 27546578

    Simoes, P., Fragata, I., Seabra, S. G., Faria, G. S., Santos, M. A., Rose, M. R., Santos, M. and Matos, M. (2017). Predictable phenotypic, but not karyotypic, evolution of populations with contrasting initial history. Sci Rep 7(1): 913. PubMed ID: 28424494

    Sulak, M., Fong, L., Mika, K., Chigurupati, S., Yon, L., Mongan, N. P., Emes, R. D. and Lynch, V. J. (2016). TP53 copy number expansion is associated with the evolution of increased body size and an enhanced DNA damage response in elephants. Elife 5. PubMed ID: 27642012

    Syed, Z. A., Chatterjee, M., Samant, M. A. and Prasad, N. G. (2017). Reproductive isolation through experimental manipulation of sexually antagonistic coevolution in Drosophila melanogaster. Sci Rep 7(1): 3330. PubMed ID: 28611437

    Tauscher, P. M., Gui, J. and Shimmi, O. (2016). Adaptive protein divergence of BMP ligands takes place under developmental and evolutionary constraints. Development [Epub ahead of print]. PubMed ID: 27578781

    Tejeda, M. T., Arredondo, J., Liedo, P., Perez-Staples, D., Ramos-Morales, P. and Diaz-Fleischer, F. (2016). Reasons for success: rapid evolution for desiccation resistance and life-history changes in the polyphagous fly Anastrepha ludens. Evolution [Epub ahead of print]. PubMed ID: 27641541

    Tobler, R., Nolte, V. and Schlotterer, C. (2017). High rate of translocation-based gene birth on the Drosophila Y chromosome. Proc Natl Acad Sci U S A 114(44): 11721-11726. PubMed ID: 29078298

    Turissini, D. A., Comeault, A. A., Liu, G., Lee, Y. C. and Matute, D. R. (2017a). The ability of Drosophila hybrids to locate food declines with parental divergence. Evolution [Epub ahead of print]. PubMed ID: 28085186

    Turissini, D. A. and Matute, D. R. (2017b). Fine scale mapping of genomic introgressions within the Drosophila yakuba clade. PLoS Genet 13(9): e1006971. PubMed ID: 28873409

    Turelli, M., et al. (2018). Rapid Global Spread of wRi-like Wolbachia across Multiple Drosophila. Curr Biol 28(6): 963-971.e968. PubMed ID: 29526588

    Wei, K. H., Reddy, H. M., Rathnam, C., Lee, J., Lin, D., Ji, S., Mason, J. M., Clark, A. G. and Barbash, D. A. (2017). A Pooled Sequencing Approach Identifies a Candidate Meiotic Driver in Drosophila. Genetics [Epub ahead of print]. PubMed ID: 28258181

    Wensing, K. U., Koppik, M. and Fricke, C. (2017). Precopulatory but not postcopulatory male reproductive traits diverge in response to mating system manipulation in Drosophila melanogaster. Ecol Evol 7(23): 10361-10378. PubMed ID: 29238561

    Williams, C. M., McCue, M. D., Sunny, N. E., Szejner-Sigal, A., Morgan, T. J., Allison, D. B. and Hahn, D. A. (2016). Cold adaptation increases rates of nutrient flow and metabolic plasticity during cold exposure in Drosophila melanogaster. Proc Biol Sci 283. PubMed ID: 27605506

    Wolff, J. N., Gemmell, N. J., Tompkins, D. M. and Dowling, D. K. (2017). Introduction of a male-harming mitochondrial haplotype via 'Trojan Females' achieves population suppression in fruit flies. Elife 6. PubMed ID: 28467301

    Yablonovitch, A. L., Fu, J., Li, K., Mahato, S., Kang, L., Rashkovetsky, E., Korol, A. B., Tang, H., Michalak, P., Zelhof, A. C., Nevo, E. and Li, J. B. (2017). Regulation of gene expression and RNA editing in Drosophila adapting to divergent microclimates. Nat Commun 8(1): 1570. PubMed ID: 29146998

    Yang, B. and Wittkopp, P. J. (2017). Structure of the Transcriptional Regulatory Network Correlates with Regulatory Divergence in Drosophila. Mol Biol Evol [Epub ahead of print]. PubMed ID: 28333240

    Yassin, A., Delaney, E. K., Reddiex, A. J., Seher, T. D., Bastide, H., Appleton, N. C., Lack, J. B., David, J. R., Chenoweth, S. F., Pool, J. E. and Kopp, A. (2016). The pdm3 locus is a hotspot for recurrent evolution of female-limited color dimorphism in Drosophila. Curr Biol. PubMed ID: 27546577

    Yu, Y., Zhou, H., Kong, Y., Pan, B., Chen, L., Wang, H., Hao, P. and Li, X. (2016). The landscape of A-to-I RNA editome is shaped by both positive and purifying selection. PLoS Genet 12: e1006191. PubMed ID: 27467689

    Zhao, L. and Begun, D. J. (2017). Genomics of parallel adaptation at two timescales in Drosophila. PLoS Genet 13(10): e1007016. PubMed ID: 28968391

    Home page: The Interactive Fly © 1995, 1996 Thomas B. Brody, Ph.D.

    The Interactive Fly resides on the
    Society for Developmental Biology's Web server.