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 profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. This study shows that a 24 h age difference of Drosophila simulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. These results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework (Hsu, 2019).
Divergent populations across different environments are exposed to critical sensory information related to locating a host or mate, as well as avoiding predators and pathogens. These sensory signals generate evolutionary changes in neuroanatomy and behavior; however, few studies have investigated patterns of neural architecture that occur between sensory systems, or that occur within large groups of closely-related organisms. This study examine 62 species within the genus Drosophila and describes an inverse resource allocation between vision and olfaction, which was consistently observe at the periphery, within the brain, as well as during larval development. This sensory variation was noted across the entire genus and appears to represent repeated, independent evolutionary events, where one sensory modality is consistently selected for at the expense of the other. Moreover, evidence is provided of a developmental genetic constraint through the sharing of a single larval structure, the eye-antennal imaginal disc. In addition, the ecological implications of visual or olfactory bias were examined, including the potential impact on host-navigation and courtship (Keesey, 2019).
Thermal adaptation is typically detected by examining the tolerance of a few populations to extreme temperatures within a single life stage. However, the extent to which adaptation occurs among many different populations might depend on the tolerance of multiple life stages and the average temperature range that the population experiences. This study examined local adaptation to native temperature conditions in eleven populations of the well-known cosmopolitan fruit fly, Drosophila melanogaster. These populations were sampled from across the global range of D. melanogaster. Traits related to fitness were measured during each life stage to determine if certain stages are more sensitive to changes in temperature than others. D. melanogaster appeared to show local adaptation to native temperatures during the egg, larval, and adult life stages, but not the pupal stage. This suggests that across the entire distribution of D. melanogaster, certain life stages might be locally adapted to native temperatures, while other stages might use phenotypic plasticity or tolerance to a wide range of temperatures experienced in the native environment of this species (Austin, 2019).
During evolution, organisms have acquired variable feeding habits. Some species are nutritional generalists that adapt to various food resources, while others are specialists, feeding on specific resources. However, much remains to be discovered about how generalists adapt to diversified diets. Larvae of the generalists Drosophila melanogaster and D. simulans develop on three diets with different nutrient balances, whereas specialists D. sechellia and D. elegans cannot develop on carbohydrate-rich diets. The generalist D. melanogaster downregulates the expression of diverse metabolic genes systemically by transforming growth factor beta (TGF-beta)/Activin signaling, maintains metabolic homeostasis, and successfully adapts to the diets. In contrast, the specialist D. sechellia expresses those metabolic genes at higher levels and accumulates various metabolites on the carbohydrate-rich diet, culminating in reduced adaptation. Phenotypic similarities and differences strongly suggest that the robust carbohydrate-responsive regulatory systems are evolutionarily retained through genome-environment interactions in the generalists and contribute to their nutritional adaptabilities (Watanabe, 2019).
Climate warming is threatening biodiversity worldwide. Climate specialists such as alpine species are especially likely to be vulnerable. Adaptation by rapid evolution is the only long-term option for survival of many species, but the adaptive evolutionary potential of heat resistance has not been assessed in an alpine invertebrate. This study show that the alpine fly Drosophila nigrosparsa cannot readily adapt to heat stress. Heat-exposed flies from a regime with increased ambient temperature and a regime with increased temperature plus artificial selection for heat tolerance were less heat tolerant than the control group. Increased ambient temperature affected negatively both fitness and competitiveness. Ecological niche models predicted the loss of three quarters of the climatically habitable areas of this fly by the end of this century. These findings suggest that, alongside with other climate specialists, species from mountainous regions are highly vulnerable to climate warming and unlikely to adapt through evolutionary genetic changes (Kinzner, 2019).
Relationships between an organism and its environment can be fundamental in the understanding how populations change over time and species arise. Local ecological conditions can shape variation at multiple levels, among these are the evolutionary history and trajectories of coding genes. This study examines the rate of molecular evolution at protein-coding genes throughout the genome in response to host adaptation in the cactophilic Drosophila mojavensis. These insects are intimately associated with cactus necroses, developing as larvae and feeding as adults in these necrotic tissues. Drosophila mojavensis is composed of four isolated populations across the deserts of western North America and each population has adapted to utilize different cacti that are chemically, nutritionally, and structurally distinct. High coverage Illumina sequencing was performed on three previously unsequenced populations of D. mojavensis. Genomes were assembled using the previously sequenced genome of D. mojavensis from Santa Catalina Island (USA) as a template. Protein coding genes were aligned across all four populations and rates of protein evolution were determined for all loci using a several approaches. Loci that exhibited elevated rates of molecular evolution tend to be shorter, have fewer exons, low expression, be transcriptionally responsive to cactus host use and have fixed expression differences across the four cactus host populations. Fast evolving genes were involved with metabolism, detoxification, chemosensory reception, reproduction and behavior. Results of this study give insight into the process and the genomic consequences of local ecological adaptation (Allan, 2019).
Population genomic data has revealed patterns of genetic variation associated with adaptation in many taxa. Yet understanding the adaptive process that drives such patterns is challenging; it requires disentangling the ecological agents of selection, determining the relevant timescales over which evolution occurs, and elucidating the genetic architecture of adaptation. Doing so for the adaptation of hosts to their microbiome is of particular interest with growing recognition of the importance and complexity of host-microbe interactions. This study tracked the pace and genomic architecture of adaptation to an experimental microbiome manipulation in replicate populations of Drosophila melanogaster in field mesocosms. Shifts in microbiome composition altered population dynamics and led to divergence between treatments in allele frequencies, with regions showing strong divergence found on all chromosomes. Moreover, at divergent loci previously associated with adaptation across natural populations, the more common allele was found in fly populations experimentally enriched for a certain microbial group was also more common in natural populations with high relative abundance of that microbial group. These results suggest that microbiomes may be an agent of selection that shapes the pattern and process of adaptation and, more broadly, that variation in a single ecological factor within a complex environment can drive rapid, polygenic adaptation over short timescales (Rudman, 2019).
Drosophila melanogaster recently spread from its tropical origin in Africa and became a cosmopolitan species that has adapted to a wide range of different thermal environments, including temperate climates. An important limiting factor of temperate climates has probably been their low and varying temperatures. The transcriptional output of genes can vary across temperatures, which might have been detrimental while settling in temperate environments. The reduction of temperature-sensitive expression of functionally important genes to ensure consistent levels of gene expression might have been relevant while adapting to such environments. This study focuses on the gene vestigial (vg) whose product is a key factor in wing development. Evidence is provided that temperature-sensitivity of vg has been buffered in populations from temperate climates. Temperature-sensitivity of vg gene expression was investigated in six natural populations, including four temperate populations (three from Europe and one from high-altitude Africa), and two tropical populations from the ancestral species range. All temperate populations exhibited a lower degree of temperature-induced expression plasticity than the tropical populations (Voigt, 2019).
Drosophila melanogaster is thought to originate from sub-Saharan Africa. This study documents the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets were found. Combining these genomes with existing data, this study probed the history of this species' geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began approximately 10,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,400 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. These findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments (Sprengelmeyer, 2019).
In animals, the most common type of RNA editing is the deamination of adenosines (A) into inosines (I). Because inosines base-pair with cytosines (C), they are interpreted as guanosines (G) by the cellular machinery and genomically encoded G alleles at edited sites mimic the function of edited RNAs. The contribution of this hardwiring effect on genome evolution remains obscure. This study looked for population genomics signatures of adaptive evolution associated with A-to-I RNA edited sites in humans and Drosophila melanogaster. Single nucleotide polymorphisms at edited sites occur 3 (humans) to 15 times (Drosophila) more often than at unedited sites, the nucleotide G is virtually the unique alternative allele at edited sites and G alleles segregate at higher frequency at edited sites than at unedited sites. This study study reveals that a significant fraction of coding synonymous and nonsynonymous as well as silent and intergenic A-to-I RNA editing sites are likely adaptive in the distantly related human and Drosophila lineages (Popitsch, 2020).
The pervasive occurrence of sexual dimorphism demonstrates different adaptive strategies of males and females. While different reproductive strategies of the two sexes are well-characterized, very little is known about differential functional requirements of males and females in their natural habitats. The impact environmental change on the selection response was studied in both sexes. Exposing replicated Drosophila populations to a novel temperature regime, sex-specific changes were demonstrated in gene expression, metabolic and behavioral phenotypes in less than 100 generations. This indicates not only different functional requirements of both sexes in the new environment but also rapid sex-specific adaptation. Supported by computer simulations it is proposed that altered sex-biased gene regulation from standing genetic variation, rather than new mutations, is the driver of rapid sex-specific adaptation. This discovery of environmentally driven divergent functional requirements of males and females has important implications-possibly even for gender aware medical treatments (Hsu, 2020).
Neuronal activity is temperature-sensitive and affects behavioral traits important for individual fitness, such as locomotion and courtship. Yet not enough is known about the evolutionary response of neuronal phenotypes in new temperature environments. This study used long-term experimental evolution of Drosophila simulans populations exposed to novel temperature regimes. A direct relationship was demonstrated between thermal selective pressure and the evolution of neuronally expressed molecular and behavioral phenotypes. Several essential neuronal genes evolve lower expression at high temperatures and higher expression at low temperatures, with dopaminergic neurons standing out by displaying the most consistent expression change across independent replicates. The link between evolved gene expression and behavioral changes was functionally validated by pharmacological intervention in the experimentally evolved D. simulans populations as well as by genetically triggered expression changes of key genes in D. melanogaster. Since natural temperature clines confirm these results for Drosophila and Anopheles populations, it is concluded that neuronal dopamine evolution is a key factor for temperature adaptation (Jaksic, 2020).
Many organisms enter a dormant state in their life cycle to deal with predictable changes in environments over the course of a year. The timing of dormancy is therefore a key seasonal adaptation, and it evolves rapidly with changing environments. The hypothesis that differences in the timing of seasonal activity are driven by differences in the rate of development during diapause was tested in Rhagoletis pomonella, a fly specialized to feed on fruits of seasonally limited host plants. Transcriptomes from the central nervous system across a time series during diapause show consistent and progressive changes in transcripts participating in diverse developmental processes, despite a lack of gross morphological change. Moreover, population genomic analyses suggested that many genes of small effect enriched in developmental functional categories underlie variation in dormancy timing and overlap with gene sets associated with development rate in Drosophila melanogaster. These transcriptional data also suggested that a recent evolutionary shift from a seasonally late to a seasonally early host plant drove more rapid development during diapause in the early fly population. Moreover, genetic variants that diverged during the evolutionary shift were also enriched in putative cis regulatory regions of genes differentially expressed during diapause development. Overall, these data suggest polygenic variation in the rate of developmental progression during diapause contributes to the evolution of seasonality in R. pomonella. Patterns that suggest hourglass-like developmental divergence early and late in diapause development are discussed along with an important role for hub genes in the evolution of transcriptional divergence (Dowle, 2020).
One hypothesis for the function of sleep is that it serves as a mechanism to conserve energy. Recent studies have suggested that increased sleep can be an adaptive mechanism to improve survival under food deprivation in Drosophila melanogaster. To test the generality of this hypothesis, Sleep and its plastic response to starvation was compared in a temperate and tropical population of Drosophila melanogaster. Flies from the temperate population were found to be more starvation resistant, and it was hypothesized that they would engage in behaviors that are considered to conserve energy, including increased sleep and reduced movement. Surprisingly, temperate flies slept less and moved more when they were awake compared to tropical flies, both under fed and starved conditions, therefore sleep did not correlate with population-level differences in starvation resistance. In contrast, total sleep and percent change in sleep when starved were strongly positively correlated with starvation resistance within the tropical population, but not within the temperate population. Thus, unexpectedly complex relationships between starvation and sleep were observed that vary both within and across populations. These observations falsify the simple hypothesis of a straightforward relationship between sleep and energy conservation. The hypothesis that starvation is correlated with metabolic phenotypes was tested by investigating stored lipid and carbohydrate levels; stored metabolites were found to partially contributed towards variation starvation resistance. These findings demonstrate that the function of sleep under starvation can rapidly evolve on short timescales and raise new questions about the physiological correlates of sleep and the extent to which variation in sleep is shaped by natural selection (Sarikaya, 2020).
Ecological adaptation is frequently inferred by the comparison of natural populations from different environments. Nevertheless, the inference of the selective forces suffers the challenge that many environmental factors covary. With well-controlled environmental conditions, experimental evolution provides a powerful approach to complement the analysis of natural populations. On the other hand, it is apparent that laboratory conditions differ in many ways from natural environments, which raises the question to what extent selection responses in experimental evolution studies can inform about adaptation processes in the wild. This study compared the expression profiles of replicated Drosophila melanogaster populations which have been exposed to two distinct temperature regimes (18/28 °C and 10/20 °C) in the laboratory for more than 80 generations. Using gene-wise differential expression analysis and co-expression network analysis, 541 genes and three co-regulated gene modules were identified that evolved in the same direction in both temperature regimes, and most of these changes probably reflect an adaptation to the space constrain or diurnal temperature fluctuation that is common in both selection regimes. 203 genes and seven modules evolved temperature-specific expression changes. Remarkably, a significant overlap was detected of these temperature-adaptive genes/modules from experimental evolution with temperature-adaptive genes inferred from natural Drosophila populations covering two different temperature clines. It is concluded that well-designed experimental evolution studies are a powerful tool to dissect evolutionary responses (Hsu, 2020).
Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, this study examined the genetic basis and evolutionary history of ovarian diapause. Outbred females were exposed to different temperatures and day lengths, ovarian development was characterized for over 2800 flies, and their complete, phased genomes were reconstructed. Diapause, scored at two different developmental cutoffs, was found to be modest heritability, and hundreds of SNPs associated with each of the two phenotypes were identified. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal variation in diapause seen in North America. Alleles associated with diapause are segregating in Zambia, suggesting that variation in diapause relies on ancestral polymorphisms, and both pro- and anti-diapause alleles have experienced selection in North America. Finally, outdoor mesocosms were used to track diapause under natural conditions. Hybrid swarms reared outdoors were found to evolve increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. These results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space (Erickson, 2020).
The repeatability or predictability of evolution is a central question in evolutionary biology and most often addressed in experimental evolution studies. This study inferred how genetically heterogeneous natural systems acquire the same molecular changes to address how genomic background affects adaptation in natural populations. In particular, advantage was taken of independently formed neo-sex chromosomes in Drosophila species that have evolved dosage compensation by co-opting the dosage-compensation male-specific lethal (MSL) complex to study the mutational paths that have led to the acquisition of hundreds of novel binding sites for the MSL complex in different species. This complex recognizes a conserved 21-bp GA-rich sequence motif that is enriched on the X chromosome, and newly formed X chromosomes recruit the MSL complex by de novo acquisition of this binding motif. Recently formed sex chromosomes were identified in the D. melanica and D. robusta species groups by genome sequencing and generate genomic occupancy maps of the MSL complex to infer the location of novel binding sites. Diverse mutational paths were utilized in each species to evolve hundreds of de novo binding motifs along the neo-X, including expansions of microsatellites and transposable element (TE) insertions. However, the propensity to utilize a particular mutational path differs between independently formed X chromosomes and appears to be contingent on genomic properties of that species, such as simple repeat or TE density. This establishes the 'genomic environment' as an important determinant in predicting the outcome of evolutionary adaptations (Ellison, 2019).
This study took advantage of naturally occurring variation in sex chromosome karyotype in Drosophila species to study independent replicates of solving the same evolutionary challenge: to dosage compensate newly formed neo-X chromosomes by acquiring hundreds of MSL-binding sites in response to Y degeneration (Ellison, 2019).
The independent acquisition of dosage compensation in Drosophila allows several important questions in evolutionary biology and gene regulation to be addressed: first, how repeatable is evolution? Evolutionary biologists have long debated the predictability of the evolutionary process. At one extreme, evolution could be highly idiosyncratic and unpredictable, since the survival of the fittest could occur along a great number of forking paths. Alternatively, constraints on evolution may force independent lineages to frequently converge on the same genetic solutions for the same evolutionary challenge. Second, how do regulatory networks evolve? And what is the contribution of TEs to regulatory evolution? Evolutionary innovations and adaptations often require rapid and concerted changes in regulation of gene expression at many loci. TEs constitute the most dynamic part of eukaryotic genomes, and the dispersal of TEs that contain a regulatory element may allow for the same regulatory motif to be recruited at many genomic locations, thereby drawing multiple genes into the same regulatory network. Third, what makes a binding motif functional? The genomes of complex organisms encompass megabases of DNA, and regulatory molecules must distinguish specific targets within this vast landscape. Regulatory factors typically identify their targets through sequence-specific interactions with the underlying DNA, but they typically bind only a fraction of the candidate genomic regions containing their specific target sequence motif. An unresolved mystery in regulatory evolution is what drives the specificity of binding to a subset of genomic regions that all appear to have a sequence that matches the consensus binding motif (Ellison, 2019).
Several features make dosage compensation in Drosophila a promising system to tackle these questions. The genetic architecture for most adaptations -- especially those involving regulatory changes -- as well as the timing and exact selective forces driving them is generally little understood. In contrast, detailed knowledge is available of the molecular mechanism of dosage compensation in Drosophila. The cis- and trans-acting components of this regulatory network and the regulatory motif for targeting the MSL complex to the X are known. Clear expectations are available of which genomic regions should acquire dosage compensation and about the timing and the evolutionary forces that drive wiring of hundreds of genes into the dosage-compensation network on newly evolved X chromosomes. Specifically, Y degeneration is a general facet of sex chromosome evolution, creating selective pressures to up-regulate X-linked genes in males. Dosage compensation should thus only evolve on neo-X chromosomes whose neo-Y homologs have started to degenerate and should evolve simultaneously or shortly after substantial gene loss has occurred on the neo-Y. Indeed, comparative data in Drosophila support this model of dosage-compensation evolution. Drosophila species with partially eroded neo-Y chromosomes exist that have not yet evolved MSL-mediated dosage compensation, including D. busckii and D. albomicans, lending empirical support to the notion that dosage compensation evolves in response to Y degeneration and not the other way round. Thus, a refined understanding of how, when, why, and where dosage compensation in Drosophila evolves makes this an ideal model system to study the repeatability of evolution and the evolution of regulatory networks (Ellison, 2019).
Results from evolution experiments indicate that although evolution is not identical in replicate populations, there is an important degree of predictability. Experimentally evolved populations under controlled, identical conditions consistently show parallelism in which mutations in certain genes are repeatedly selected. However, organisms adapting to similar environments are not genetically identical, but their genome instead carries the legacy of their unique evolutionary trajectory, raising the question of how genomic differences affect genetic parallelism (Ellison, 2019).
Sex chromosome-autosome fusions have independently created neo-sex chromosomes in different Drosophila lineages. This provides everal independent replicates to study how, on the molecular level, evolution has solved the same adaptive challenge: acquiring hundreds of binding sites to recruit the MSL complex to newly formed X chromosomes. This allows quantification of how much variation there is, both within and between species, in the underlying mutational paths to acquire hundreds of MSL-binding sites on neo-X chromosomes and identify genomic contingencies that will influence the repeatability of evolutionary trajectories. Importantly, neo-sex chromosomes of Drosophila are evolutionarily young (between 0.1-15 MY old), which allows, in many cases, inferring of the causative mutations that have resulted in the gain of a regulatory element and decipher the evolutionary processes at work to draw hundreds of genes into a new regulatory network (Ellison, 2019).
The results suggest that the evolution of MSL-binding sites is highly opportunistic but contingent on genomic background. In particular, each independently evolved neo-X chromosome was found to use a diverse set of mutational pathways to acquire MSL-binding sites on a new neo-X chromosome, ranging from microsatellite expansions to the utilization of presites to TE insertions. However, different lineages differ with regards to the frequency of which mutational paths are most often followed to acquire novel binding sites, and this propensity may depend on the genomic background. In particular, species found with the higher density of simple repeats are more prone to utilize expansions in GA microsatellites to gain a novel MSL-binding site. In contrast, D. robusta has an elevated TE density compared to its sibling species, and it was found that the dispersal of a TE has played an important role in the acquisition of MSL-binding sites on its neo-X chromosome. Thus, this suggests that the genomic background of a species predisposes it to evolve along a particular path, yet the evolutionary process is random and resourceful with regards to utilizing a variety of mutations to create novel MSL-binding sites. However, different phenotypes show drastic differences in their underlying genetic architecture, and the importance of genomic background likely differs among traits and (Ellison, 2019).
Evolutionary innovations and adaptations often require rapid and concerted changes in regulation of gene expression at many loci. It has been suggested that TEs play a key role in rewiring regulatory networks, since the dispersal of TEs that contain a regulatory element may allow for the same regulatory motif to be recruited at many genomic locations. A handful of recent studies have implicated TEs as drivers of key evolutionary innovations, including placentation in mammals or rewiring the core regulatory network of human embryonic stem cells. While these studies demonstrate that TEs can, in principle, contribute to the creation or rewiring of regulatory networks, they do not address the question of how often regulatory elements evolve by TE insertions versus by other mutations. That is, the importance of TEs in contributing to regulatory evolution is not known. Quantification of the role of TEs would require a priori knowledge of how and when regulatory networks evolve and a detailed molecular understanding of which genes are being drawn into a regulatory network and how. As discussed above, these parameters are well understood for dosage compensation in flies (Ellison, 2019).
Previous work in D. miranda has shown that a helitron TE was recruited into the dosage-compensation network at two independent time points. The younger 1.5-MY-old neo-X chromosome of D. miranda is in the process of evolving dosage compensation, and dozens of new CESs on this chromosome were created by insertions of the ISX element. The domesticated ISX TE gained a novel MRE motif by a 10-bp deletion in the ISY element, which is a highly abundant TE in the D. miranda genome. The remnants of a related (but different) TE at CES was found on the older neo-X of this species (which formed roughly 13-15 MY ago), but the TE was too eroded to reconstruct its evolutionary history. This study, identified another domesticated TE that was utilized to deliver MSL-binding sites to a newly formed neo-X chromosome, but no significant TE contribution was found for MSL-binding site evolution in two independent neo-X chromosomes (Ellison, 2019).
The data shed light on the question of when TEs are expected to be important in regulatory evolution. For TEs to contribute to regulatory rewiring, two conditions have to be met: a regulatory element (or a progenitor sequence that can easily mutate into the required binding motif) needs to be present in the TE, and TE needs to be active in the genome (and not yet silenced by the host machinery). TEs undergo a characteristic life cycle in which they invade a new species (or escape the genome defense by mutation) and transpose until they are silenced by the host genome. Once a TE is robustly repressed, it no longer can serve as a vehicle to disperse regulatory elements, so the time window when a particular TE family can be domesticated is probably short and needs to coincide with a necessity to disperse regulatory motifs. A high TE burden does increase that chance, but at a cost: maintaining active TEs in the genome allows a rapid response to evolutionary challenges but also creates a major source of genomic mutation, illegitimate recombination, genomic rearrangements, and genome size inflation (Ellison, 2019).
The current findings support this view of a TE tradeoff. The ISY element in D. miranda is the most highly abundant transposon in the D. miranda genome and is massively contributing to the degeneration of the neo-Y in this specie. Indeed, the genomic analysis has revealed >20,000 novel insertions of the ISY element on the neo-Y, often within genes. Yet, it contained a sequence that was only one mutational step away from a functional MSL-binding site (that is, a single 10-bp deletion), and domestication of this element allowed for the rapid dispersal of functional binding sites for the MSL complex along the neo-X. The domestication of the TE in D. robusta occurred too long ago for to reconstruct its exact evolutionary history and the potential damage its mobilization may have caused while it was active. However, consistent with a tradeoff that the host genome faces, it was found that D. robusta has a higher TE density than its sister species and also a considerably larger genome size, yet a TE contributed to wiring hundreds of genes into the dosage-compensation network on its neo-X (Ellison, 2019).
Perhaps surprisingly, in many instances, it was not possible to detect specific mutations that would generate a novel MSL binding motif. Instead, it was found that functional MSL-binding sites are derived from presites containing the GA-rich motif that was already present in an ancestor in which the neo-X is autosomal and in which these sequences do not recruit the MSL complex. The MSL binding motif is only modestly enriched on the X chromosome compared to the autosomes (only approximately 2-fold), and only a small fraction of putative binding sites are actually bound by the MSL complex. The dosage-compensation machinery shares this characteristic with many other sequence-specific binding factors whose predicted target motifs are often in vast excess to the sites actually utilized. It has been speculated that other genomic aspects, such as chromatin context or the 3D organization of the genome, could help to distinguish between utilized and nonutilized copies of a motif. The finding that a large number of sites can acquire the ability to recruit the MSL complex, without any apparent associated changes at the DNA level, supports the view that epigenetic modifications or changes to the 3D architecture of the genome help to ultimately determine which putative binding sites in the genome are actually utilized. In D. melanogaster, the X chromosome has a unique satellite DNA composition, and it was suggested that these repeats play a primary role in determining X identity during dosage compensation. Furthermore, localization of the MSL complex to MREs is dependent on an additional cofactor, the CLAMP protein. CLAMP binds directly to GA-rich MRE sequences and targets MSL to the X chromosome but also binds to GA-rich sequence elements throughout the genome. Recent work has shown that variability in sequence composition within similar GA-rich motifs drive specificity for CLAMP binding, and variation within seemingly similar cis elements may also drive context-specific targeting of the MSL complex. Future investigations of changes in the chromatin level, the repeat content, and the genomic architecture of these newly formed sex chromosomes will help to resolve this outstanding question (Ellison, 2019).
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).
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).
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).
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).
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).
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).
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).
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).
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).
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)
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).
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)
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).
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).
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).
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).
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).
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).
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).
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).
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).
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).
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).
Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. This study investigated the genetic changes associated with
adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species. Disproportionately high rates of gene gains were found in internal branches in the species' phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species' differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one. It is concluded that daptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change (Rane, 2019).
Phenotypic plasticity can allow organisms to respond to environmental changes by producing better matching phenotypes without any genetic change. Because of this, plasticity is predicted to be a major mechanism by which a population can survive the initial stage of colonizing a novel environment. This prediction was tested by challenging wild Drosophila melanogaster with increasingly extreme larval environments and then examining expression of alcohol dehydrogenase (ADH) and its relationship to larval survival in the first generation of encountering a novel environment. Most families responded in the adaptive direction of increased ADH activity in higher alcohol environments and families with higher plasticity were also more likely to survive in the highest alcohol environment. Thus, plasticity of ADH activity was positively selected in the most extreme environment and was a key trait influencing fitness. Furthermore, there was significant heritability of ADH plasticity that can allow plasticity to evolve in subsequent generations after initial colonization. The adaptive value of plasticity, however, was only evident in the most extreme environment and had little impact on fitness in less extreme environments. The results provide one of the first direct tests of the adaptive role of phenotypic plasticity in colonizing a novel environment (Wang, 2019).
The genetic architecture of adaptive traits is of key importance to predict evolutionary responses. Most adaptive traits are polygenic-i.e., result from selection on a large number of genetic loci-but most molecularly characterized traits have a simple genetic basis. This discrepancy is best explained by the difficulty in detecting small allele frequency changes (AFCs) across many contributing loci. To resolve this, laboratory natural selection was used to detect signatures for selective sweeps and polygenic adaptation. This study exposed 10 replicates of a Drosophila simulans population to a new temperature regime and uncovered a polygenic architecture of an adaptive trait with high genetic redundancy among beneficial alleles. Convergent responses were observed for several phenotypes-e.g., fitness, metabolic rate, and fat content-and a strong polygenic response (99 selected alleles; mean s = 0.059). However, each of these selected alleles increased in frequency only in a subset of the evolving replicates. Different evolutionary paradigms were discerned based on the heterogeneous genomic patterns among replicates. Redundancy and quantitative trait (QT) paradigms fitted the experimental data better than simulations assuming independent selective sweeps. These results show that natural D. simulans populations harbor a vast reservoir of adaptive variation facilitating rapid evolutionary responses using multiple alternative genetic pathways converging at a new phenotypic optimum. This key property of beneficial alleles requires the modification of testing strategies in natural populations beyond the search for convergence on the molecular level (Barghi, 2019).
The cohabitation of Drosophila melanogaster with humans is nearly ubiquitous. Though it has been well-established that this fly species originated in sub-Saharan Africa, and only recently has spread globally, many details of its swift expansion remain unclear. Elucidating the demographic history of D. melanogaster provides a unique opportunity to investigate how human movement might have impacted patterns of genetic diversity in a commensal species, as well as providing neutral null models for studies aimed at identifying genomic signatures of local adaptation. This study used whole-genome data from five populations (Africa, North America, Europe, Central Asia, and the South Pacific) to carry out demographic inferences, with particular attention to the inclusion of migration and admixture. The importance of these parameters for model fitting is demonstrated, and how previous estimates of divergence times are likely to be significantly underestimated as a result of not including them is shown. Finally, how human movement along early shipping routes might have shaped the present-day population structure of D. melanogaster is discussed (Arguello, 2019).
Gene expression profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. This study shows that a 24 h age difference of Drosophila simulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. These results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework (Hsu, 2019).
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).
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).
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).
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).
Genes involved in immune defense against pathogens provide some of the most well-known examples of both directional and balancing selection. Antimicrobial peptides (AMPs) are innate immune effector genes, playing a key role in pathogen clearance in many species, including Drosophila. Conflicting lines of evidence have suggested AMPs may be under directional, balancing or purifying selection. This study used both a linear model and control gene-based approach to show that balancing selection is an important force shaping AMP diversity in Drosophila. In D. melanogaster, this is most clearly observed in ancestral African populations. Furthermore, the signature of balancing selection is even more striking once background selection has been accounted for. Balancing selection also acts on AMPs in D. mauritiana, an isolated island endemic separated from D. melanogaster by about 4 million years of evolution. This suggests that balancing selection may be broadly acting to maintain adaptive diversity in Drosophila AMPs, as has been found in other taxa (Chapman, 2019).
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).
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).
Orphan genes, lacking detectable homologs in outgroup species, typically represent 10-30% of eukaryotic genomes. This study investigated the roots of orphan gene emergence in the Drosophila genus. Across the annotated proteomes of twelve species, 6297 orphan genes were found within 4953 taxon-specific clusters of orthologs. By inferring the ancestral DNA as non-coding for between 550 and 2467 (8.7-39.2%) of these genes, this study describes for the first time how de novo emergence contributes to the abundance of clade-specific Drosophila genes. In support of them having functional roles, it was shown that de novo genes have robust expression and translational support. However, the distinct nucleotide sequences of de novo genes, which have characteristics intermediate between intergenic regions and conserved genes, reflect their recent birth from non-coding DNA. It was found that de novo genes encode more disordered proteins than both older genes and intergenic regions. Together, these results suggest that gene emergence from non-coding DNA provides an abundant source of material for the evolution of new proteins. Following gene birth, gradual evolution over large evolutionary timescales moulds sequence properties towards those of conserved genes, resulting in a continuum of properties whose starting points depend on the nucleotide sequences of an initial pool of novel genes (Heames, 2020).
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).
Insect odorant receptor (Or) genes determine the responses of sensory neurons that mediate critical behaviours. The Drosophila melanogaster Or22 locus represents an interesting example of molecular evolution, with high levels of sequence divergence and copy number variation between D. melanogaster and other Drosophila species, and a corresponding high level of variability in the responses of the neuron it controls, ab3A. However, the link between Or22 molecular and functional diversity has not been established. This study shows that several naturally occurring Or22 variants generate major shifts in neuronal response properties. The molecular changes were determined that underpin these response shifts, one of which represents a chimaeric gene variant previously suggested to be under natural selection. In addition it was shown that several alternative molecular genetic mechanisms have evolved for ensuring that where there is more than one gene copy at this locus, only one functional receptor is generated. These data thus provide a causal link between the striking levels of phenotypic neuronal response variation found in natural populations of D. melanogaster and genetic variation at the Or22 locus. Since neuronal responses govern animal behaviour, it is predicted that Or22 may be a key player in underlying one or more olfactory-driven behaviours of significant adaptive importance (Shaw, 2019).
Drosophila suzukii Matsumura (Diptera: Drosophilidae) is a significant invasive pest in soft-skin fruits and berries in Asia, Europe, and North and South America. Many herbivorous insects use multiple cues for host selection, particularly olfactory and visual stimuli. The visual system of closely-related Drosophila melanogaster is well-documented, expressing strong sensitivity to short-wavelength colors (ultraviolet to green) and only limited sensitivity to long-wavelength colors (red to infrared). The results suggest that D. suzukii have limited ability to distinguish red consistent with visual sensitivity range within the melanogaster subgroup. It is proposed that color contrast rather than color appearance may be of greater importance in orientation and attraction. It is proposed that differences in reflectance between light wavelengths important for color opponency are key to color discrimination to provide color contrast between foreground and background, as occurs between fruit and foliage, during host-finding (Little, 2019).
The behaviors of closely related species can be remarkably different, and these differences have important ecological and evolutionary consequences. While the recent boom in genotype-phenotype studies has led to a greater understanding of the genetic architecture and evolution of a variety of traits, studies identifying the genetic basis of behaviors are, comparatively, still lacking. This is likely because they are complex and environmentally sensitive phenotypes, making them difficult to measure reliably for association studies. The Drosophila species complex holds promise for addressing these challenges, as the behaviors of closely related species can be readily assayed in a common environment. This study investigated the genetic basis of an evolved behavioral difference, pupation site choice, between Drosophila melanogaster and D. simulans. In this study, A significant contribution was demonstrated of the X chromosome to the difference in pupation site choice behavior between these species. Using a panel of X-chromosome deficiencies, the majority of the X chromosome was screened for causal loci, and two regions were identified associated with this X-effect. Gene disruption and RNAi data were collecting supporting a single gene that affects pupation behavior within each region: Fas2 and tilB. Finally, differences in tilB expression were shown to correlate with the differences in pupation site choice behavior between species. This evidence associating two genes with differences in a complex, environmentally sensitive behavior represents the first step towards a functional and evolutionary understanding of this behavioral divergence (Pischedda, 2019).
Trade-off between vision and olfaction, the fact that investment in one correlates with decreased investment in the other, has been demonstrated by a wealth of comparative studies. However, there is still no empirical evidence suggesting how these two sensory systems coevolve, i.e. simultaneously or alternatively. The "Dark-flies" (Drosophila melanogaster) constitute a unique model to investigate such relation since they have been reared in the dark since 1954, approximately 60 years (~1500 generations). To observe how vision and olfaction evolve, populations of Dark-flies were reared in normal lighting conditions for 1 (DF1G) and 65 (DF65G) generations. The sizes of the visual (optic lobes, OLs) and olfactory (antennal lobes, ALs) primary centres, as well as the rest of the brain, and compared the results with the original and its genetically most similar strain (Oregon flies). Whereas the ALs decreased in size, the OLs (together with the brain) increased in size in the Dark-flies returned back to the light, both in the DF1G and DF65G. These results experimentally show that trade-off between vision and olfaction occurs simultaneously, and suggests that there are possible genetic and epigenetic processes regulating the size of both optic and antennal lobes. Furthermore, although the Dark-flies were able to mate and survive in the dark with a reduced neural investment, individuals being returned to the light seem to have been selected with reinvestment in visual capabilities despite a potential higher energetic cost (Ozer, 2020).
The evolution of animal behaviour is poorly understood. Despite numerous correlations between interspecific divergence in behaviour and nervous system structure and function, demonstrations of the genetic basis of these behavioural differences remain rare. This study develop a neurogenetic model, Drosophila sechellia, a species that displays marked differences in behaviour compared to its close cousin Drosophila melanogaster that are linked to its extreme specialization on noni fruit (Morinda citrifolia). Using calcium imaging, this study identified olfactory pathways in D. sechellia that detect volatiles emitted by the noni host. This mutational analysis indicates roles for different olfactory receptors in long- and short-range attraction to noni, and cross-species allele-transfer experiments demonstrate that the tuning of one of these receptors is important for species-specific host-seeking. The molecular determinants of this functional change were identified, and their evolutionary origin and behavioural importance were characterized. Circuit tracing was performed in the D. sechellia brain, and receptor adaptations were found to be accompanied by increased sensory pooling onto interneurons as well as species-specific central projection patterns. This work reveals an accumulation of molecular, physiological and anatomical traits that are linked to behavioural divergence between species, and defines a model for investigating speciation and the evolution of the nervous system (Auer, 2020).
Females of many animal species mate several times with different males (polyandry), whereas females of some species mate with a single male (monandry) only once. Little is known about the mechanisms by which these different mating systems evolve. Females of Drosophila prolongata mate serially, unlike Drosophila melanogaster females that refuse to remate for several days after their first mating (remating suppression [RS]). Nevertheless, interestingly, nonvirgin D. prolongata females refuse to remate with males that are prohibited from performing their species-specific courtship behavior, leg vibration (LV), suggesting that LV overrides RS making it cryptic in D. prolongata. This study examined how long this cryptic RS persists. Surprisingly, it was sustained for at least 2 weeks, showing that RS is substantially augmented in D. prolongata compared to that of D. melanogaster. The two most closely related species to D. prolongata, Drosophila rhopaloa and Drosophila carrolli, do not perform LV and showed augmented RS, supporting the idea that augmented RS could have evolved before LV was acquired. These results suggested that D. prolongata females are intrinsically monandrous, whereas the newly evolved courtship behavior makes them polyandrous. This is a rare case in which a proximate mechanism of polyandry evolution from monandry is demonstrated (Minekawa, 2020).
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). 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).
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).
It is unclear where in the nervous system evolutionary changes tend to occur. To localize the source of neural evolution that has generated divergent behaviors, a new approach was developed to label and functionally manipulate homologous neurons across Drosophila species. Homologous descending neurons were examined that drive courtship song in two species that sing divergent song types, and relevant evolutionary changes were localized in circuit function downstream of the intrinsic physiology of these descending neurons. This evolutionary change causes different species to produce divergent motor patterns in similar social contexts. Artificial stimulation of these descending neurons drives multiple song types, suggesting that multifunctional properties of song circuits may facilitate rapid evolution of song types (Ding, 2019).
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).
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).
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).
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).
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).
Germ lines are the cell lineages that give rise to the sperm and eggs in animals. The germ lines first arise from primordial germ cells (PGCs) during embryogenesis: these form from either a presumed derived mode of preformed germ plasm (inheritance) or from an ancestral mechanism of inductive cell-cell signalling (induction). Numerous genes involved in germ line specification and development have been identified and functionally studied. However, little is known about the molecular evolutionary dynamics of germ line genes in metazoan model systems. This study examined the molecular evolution of germ line genes within three metazoan model systems. These include the genus Drosophila (N=34 genes, inheritance), the fellow insect Apis (N=30, induction), and their more distant relative Caenorhabditis (N=23, inheritance). Using multiple species and established phylogenies in each genus, this study reports that germ line genes exhibited marked variation in the constraint on protein sequence divergence (dN/dS) and codon usage bias (CUB) within each genus. Importantly, it was found that de novo lineage-specific inheritance (LSI) genes in Drosophila (osk, pgc) and in Caenorhabditis (pie-1, pgl-1), which are essential to germ plasm functions under the derived inheritance mode, displayed rapid protein sequence divergence relative to the other germ line genes within each respective genus. This may reflect the evolution of specialized germ plasm functions and/or low pleiotropy of LSI genes, features not shared with other germ line genes. In addition, it was observed that the relative ranking of dN/dS and of CUB between genera were each more strongly correlated between Drosophila and Caenorhabditis, from different phyla, than between Drosophila and its insect relative Apis, suggesting taxonomic differences in how germ line genes have evolved. Taken together, the present results advance understanding of the evolution of animal germ line genes within three well-known metazoan models. Further, the findings provide insights to the molecular evolution of germ line genes with respect to LSI status, pleiotropy, adaptive evolution as well as PGC-specification mode (Whittle, 2019a).
At the very end of the larval stage Drosophila expectorate a glue secreted by their salivary glands to attach themselves to a substrate while pupariating. The glue is a mixture of apparently unrelated proteins, some of which are highly glycosylated and possess internal repeats. Because species adhere to distinct substrates (i.e. leaves, wood, rotten fruits), glue genes are expected to evolve rapidly. Available genome sequences and PCR-sequencing of regions of interest were used to investigate the glue genes in 20 Drosophila species. A new gene in addition was discovered to the seven glue genes annotated in D. melanogaster. A phase 1 intron was identified at a conserved position present in five of the eight glue genes of D. melanogaster, suggesting a common origin for those glue genes. A slightly significant rate of gene turnover was inferred. Both the number of repeats and the repeat sequence were found to diverge rapidly, even between closely related species. High repeat number variation was also detected at the intrapopulation level in D. melanogaster. It is concluded that most conspicuous signs of accelerated evolution are found in the repeat regions of several glue genes (Da Lage, 2019).
The subgenus Sophophora of Drosophila, which includes D. melanogaster, is an important model for the study of molecular evolution, comparative genomics, and evolutionary developmental biology. Numerous phylogenetic studies have examined species relationships in the well-known melanogaster, obscura, willistoni, and saltans species groups, as well as the relationships among these clades. In contrast, other species groups of Sophophora have been relatively neglected and have not been subjected to molecular phylogenetic analysis. This study focuses on the endemic African Drosophila fima and dentissima lineages. Both these clades fall within the broadly defined melanogaster species group, but are otherwise distantly related to each other. The new phylogeny supports pervasive divergent and convergent evolution of male-specific grasping structures (sex combs). The implications of these results are discussed for defining the boundaries of the melanogaster species group, and weigh the relative merits of "splitting" and "lumping" approaches to the taxonomy of this key model system (Kopp, 2019).
Evolutionary convergence is a core issue in the study of adaptive evolution, as well as a highly debated topic at present. Few studies have analyzed this issue using a 'real-time' or evolutionary trajectory approach. Do populations that are initially differentiated converge to a similar adaptive state when experiencing a common novel environment? Drosophila subobscura populations founded from different locations and years showed initial differences and variation in evolutionary rates in several traits during short-term (approximately 20 generations) laboratory adaptation. This study has extended that analysis to 40 more generations to analyze (1) how differences in evolutionary dynamics among populations change between shorter and longer time spans, and (2) whether evolutionary convergence occurs after 60 generations of evolution in a common environment. Substantial variation was found in longer term evolutionary trajectories and differences between short- and longer term evolutionary dynamics. Although pervasive patterns of convergence were observed toward the character values of long-established populations, populations still remain differentiated for several traits at the final generations analyzed. This pattern might involve transient divergence, as is reported in some cases, indicating that more generations should lead to final convergence. These findings highlight the importance of longer term studies for understanding convergent evolution (Simoes, 2019).
Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. This study focused on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, this study sought to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so population samples of D. melanogaster transcriptomes were generated from five tissues: accessory gland, testis, larval salivary gland, female head, and first instar larva. These data were combined with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster. A large number of gain or loss of expression phenotypes were observed, most of which were polymorphic within D. melanogaster. Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, many more novel expression phenotypes were found in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue specific expression. Finally, in addition to qualitatively novel expression phenotypes, genes exhibiting major quantitative expression divergence in the D. melanogaster lineage were identified (Cridland, 2020).
Synergistic epistasis for deleterious alleles is relevant to the mutation load paradox and the evolution of sex and recombination. Some studies have shown evidence of synergistic epistasis for spontaneous or induced deleterious mutations appearing in mutation-accumulation experiments. However, many newly arising mutations may not actually be segregating in natural populations because of the erasing action of natural selection. A demonstration of synergistic epistasis for naturally segregating alleles can be achieved by means of inbreeding depression studies, as deleterious recessive allelic effects are exposed in inbred lines. This paper reports the results of two independent inbreeding experiments carried out with two different populations of Drosophila melanogaster. The results show a consistent accelerated inbreeding depression for fitness, suggesting synergistic epistasis among deleterious alleles. Computer simulations were performed assuming different possible models of epistasis and mutational parameters for fitness, finding some of them to be compatible with the results observed. These results suggest that synergistic epistasis for deleterious mutations not only occurs among newly arisen spontaneous or induced mutations, but also among segregating alleles in natural populations (Dominguez-Garcia, 2019).
Predicting how species will respond to selection pressures requires understanding the factors that constrain their evolution. This study used genome engineering of Drosophila to investigate constraints on the repeated evolution of unrelated herbivorous insects to toxic cardiac glycosides, which primarily occurs via a small subset of possible functionally-relevant substitutions to Na(+),K(+)-ATPase. Surprisingly, frequently observed adaptive substitutions were found at at two sites, 111 and 122, are lethal when homozygous and adult heterozygotes exhibit dominant neural dysfunction. A phylogenetically correlated substitution, A119S, was found that partially ameliorates the deleterious effects of substitutions at 111 and 122. Despite contributing little to cardiac glycoside-insensitivity in vitro, A119S, like substitutions at 111 and 122, substantially increases adult survivorship upon cardiac glycoside exposure. These results demonstrate the importance of epistasis in constraining adaptive paths. Moreover, by revealing distinct effects of substitutions in vitro and in vivo, the results underscore the importance of evaluating the fitness of adaptive substitutions and their interactions in whole organisms (Taverner, 2019).
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).
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).
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).
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).
The impact of different reproductive barriers on species or population isolation may vary in different stages of speciation depending on evolutionary forces acting within species and through species' interactions. Genetic incompatibilities between interacting species are expected to reinforce prezygotic barriers in sympatric populations and lead to cascade reinforcement between conspecific populations living within and outside the areas of sympatry. These predictions were tested and this work studied whether and how the strength and target of reinforcement between Drosophila montana and Drosophila flavomontana vary between sympatric populations with different histories and species abundances. All barriers between D. montana females and D. flavomontana males were nearly complete, while in the reciprocal cross strong postzygotic isolation was accompanied by prezygotic barriers whose strength varied according to population composition. Sexual isolation between D. flavomontana females and D. montana males was increased in long-established sympatric populations, where D. flavomontana is abundant, while postmating prezygotic (PMPZ) barriers were stronger in populations where this species is a new invader and still rare and where female discrimination against heterospecific males was lower. Strengthening of sexual and PMPZ barriers in this cross also induced cascade reinforcement of respective barriers between D. flavomontana populations, which is a classic signature of reinforcement process (Poikela, 2019).
Lifetime reproductive capacity is a critical fitness component. In insects, female reproductive capacity is largely determined by the number of ovarioles, the egg-producing subunits of the ovary. Recent work has provided insights into ovariole number regulation in Drosophila melanogaster. However, whether mechanisms discovered under laboratory conditions explain evolutionary variation in natural populations is an outstanding question. This study investigated potential effects of ecology on the developmental processes underlying ovariole number evolution among Hawaiian Drosophila, a large adaptive radiation wherein the highest and lowest ovariole numbers of the family have evolved within 25 million years. Previous studies proposed that ovariole number correlated with oviposition substrate but sampled largely one clade of these flies and were limited by a provisional phylogeny and the available comparative methods. This hypothesis was tested by applying phylogenetic modeling to an expanded sampling of ovariole numbers and substrate types and shows support for these predictions across all major groups of Hawaiian Drosophila, wherein ovariole number variation is best explained by adaptation to specific substrates. Furthermore, oviposition substrate evolution is linked to changes in the allometric relationship between body size and ovariole number. Finally, evidence is provided that the major changes in ovarian cell number that regulate D. melanogaster ovariole number also regulate ovariole number in Hawaiian drosophilids. Thus, evidence is provided that this remarkable adaptive radiation is linked to evolutionary changes in a key reproductive trait regulated at least partly by variation in the same developmental parameters that operate in the model species D. melanogaster (Sarikaya, 2019).
When two species interbreed, the resulting hybrid offspring are often sterile, with the heterogametic (e.g. XY) hybrid usually being more severely affected. The prevailing theory for this pattern of sterility evokes divergent changes in separate lineages having maladaptive interactions when placed together in a hybrid individual, with recessive factors on the sex chromosome interacting with dominant factors on the autosomes. The effect of these interactions on gametogenesis should not be uniform across species pairs unless genetic divergence follows the same paths in different lineages or if a specific stage of gametogenesis is more susceptible to detrimental genetic interactions. A detailed cellular characterization of hybrid male sterility was performed across three recently diverged species pairs of Drosophila. Across all three pairs, sterile hybrid sperm are alive but exhibit rapid nuclear de-condensation with age, with active, but non-differentiated, mitochondria. Surprisingly, all three sets of interspecies hybrids produce half of the number of sperm per round of spermatogenesis, with each sperm cell containing two tails. Non-disjunction failures during meiosis I were identified as the likely cause. Thus, errors during meiosis I may be a general phenomenon underlying Drosophila male sterility, indicating either a heightened sensitivity of this spermatogenic stage to failure, or a basis to sterility other than the prevailing model (Kanippayoor, 2020).
The application of Wolbachia in insect pest and vector control requires the establishment of genotypically stable host associations. The cytoplasmic incompatibility (CI) inducing Wolbachia strain wCer2 naturally occurs in the cherry fruit fly Rhagoletis cerasi as co-infection with other strains and was transferred to other fruit flies by embryonic microinjections. wCer2 genome data was obtained from its native and three novel hosts, Drosophila simulans, Drosophila melanogaster and Ceratitis capitata and assessed its genome stability, characteristics, and CI factor (cif) genes. De novo assembly was successful from Wolbachia cell-enriched singly infected D. simulans embryos, with minimal host and other bacterial genome traces. The low yield of Wolbachia sequence reads from total genomic extracts of one multiply infected R. cerasi pupa and one singly infected C. capitata adult limited de novo assemblies but was sufficient for comparative analyses. Across hosts wCer2 was stable in genome synteny and content. Polymorphic nucleotide sites were found in wCer2 of each host, however only one nucleotide was different between R. cerasi and C. capitata, and none between replicated D. simulans lines. The wCer2 genome is highly similar to wAu (D. simulans), wMel (D. melanogaster) and wRec (Drosophila recens). In contrast to wMel and wRec (each with one cif gene pair) and wAu (without any cif genes), wCer2 has three pairs of Type I cif genes and one Type IV cifB gene without a cifA complement. This may explain previously reported CI patterns of wCer2, including incomplete rescue of its own CI modification in three novel host species (Morrow, 2020).
Recently diverged populations in the early stages of speciation offer an opportunity to understand mechanisms of isolation and their relative contributions. Drosophila willistoni is a tropical species with broad distribution from Argentina to the southern United States, including the Caribbean islands. A postzygotic barrier between northern populations (North America, Central America, and the northern Caribbean islands) and southern populations (South American and the southern Caribbean islands) has been recently documented and used to propose the existence of two different subspecies. This study has identified premating isolation between populations regardless of their subspecies status. No evidence of postmating prezygotic isolation was found, and this study proceeded to characterize hybrid male sterility between the subspecies. Sterile male hybrids transfer an ejaculate that is devoid of sperm but causes elongation and expansion of the female uterus. In sterile male hybrids, bulging of the seminal vesicle appears to impede the movement of the sperm toward the sperm pump, where sperm normally mixes with accessory gland products. The results highlight a unique form of hybrid male sterility in Drosophila that is driven by a mechanical impediment to transfer sperm rather than by an abnormality of the sperm itself. Interestingly, this form of sterility is reminiscent of a form of infertility (azoospermia) that is caused by lack of sperm in the semen due to blockages that impede the sperm from reaching the ejaculate (Davis, 2020).
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).
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).
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).
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).
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).
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).
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).
Cytological studies revealed that the number of chromosomes and their organization varies across species. The increasing availability of whole genome sequences of multiple species across specific phylogenies has confirmed and greatly extended these cytological observations. In the Drosophila genus, the ancestral karyotype consists of five rod-like acrocentric chromosomes (Muller elements A to E) and one dot-like chromosome (element F), each exhibiting a generally conserved gene content. Chromosomal fusions and paracentric inversions are thus the major contributors, respectively, to chromosome number variation among species and to gene order variation within chromosomal element. The subobscura cluster of Drosophila consists in three species that retain the genus ancestral karyotype and differ by a reduced number of fixed inversions. This study used cytological information and the D. guanche genome sequence to identify and molecularly characterize the breakpoints of inversions that became fixed since the D. guanche-D. subobscura split. The results have led to a proposal of a modified version of the D. guanche cytological map of its X chromosome, and to establish that (i) most inversions became fixed in the D. subobscura lineage and (ii) the order in which the four X chromosome overlapping inversions occurred and became fixed (Orengo. 2019).
It has been hypothesized that individually-rare hidden structural variants (SVs) could account for a significant fraction of variation in complex traits. This study identified more than 20,000 euchromatic SVs from 14 Drosophila melanogaster genome assemblies, of which ~40% are invisible to high specificity short-read genotyping approaches. SVs are common, with 31.5% of diploid individuals harboring a SV in genes larger than 5kb, and 24% harboring multiple SVs in genes larger than 10kb. SV minor allele frequencies are rarer than amino acid polymorphisms, suggesting that SVs are more deleterious. A number of functionally important genes were shown to harbor previously hidden structural variants likely to affect complex phenotypes. Furthermore, SVs are overrepresented in candidate genes associated with quantitative trait loci mapped using the Drosophila Synthetic Population Resource. It is concluded that SVs are ubiquitous, frequently constitute a heterogeneous allelic series, and can act as rare alleles of large effect (Chakraborty, 2019).
The Drosophila obscura species group shows dramatic variation in karyotype, including transitions among sex chromosomes. Members of the affinis and pseudoobscura subgroups contain a neo-X chromosome (a fusion of the X with an autosome), and it was shown that ancestral Y genes have become autosomal in species harboring the neo-X. Detailed analysis of species in the pseudoobscura subgroup revealed that ancestral Y genes became autosomal through a translocation to the small dot chromosome. This study shows that the Y-dot translocation is restricted to the pseudoobscura subgroup, and translocation of ancestral Y genes in the affinis subgroup likely followed a different route. Most ancestral Y genes appear to have translocated to unique autosomal or X-linked locations in different taxa of the affinis subgroup, and a dynamic model of sex chromosome formation and turnover in the obscura species group is proposed. The results suggest that Y genes can find unique paths to escape unfavorable genomic environments that form after sex chromosome-autosome fusions (Bracewell, 2020).
Chromosomal inversions are among the primary drivers of genome structure evolution in a wide range of natural populations. While there is an impressive array of theory and empirical analyses that have identified conditions under which inversions can be positively selected, comparatively little data is available on the fitness impacts of these genome structural rearrangements themselves. Because inversion breakpoints can disrupt functional elements and alter chromatin domains, the precise positioning of an inversion's breakpoints can strongly affect its fitness. This study compared the fine-scale distribution of low frequency inversion breakpoints with those of high frequency inversions and inversions that have gone to fixation between Drosophila species. A number of differences were identified among frequency classes that may influence inversion fitness. In particular, breakpoints that are proximal to insulator elements, generate large tandem duplications, and minimize impacts on gene coding spans are more prevalent in high frequency and fixed inversions than in rare inversions. The data suggest that natural selection acts to preserve both genes and larger cis-regulatory networks in the occurrence and spread of rearrangements. These factors may act to limit the availability of high fitness arrangements when suppressed recombination is favorable (McBroome, 2020).
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).
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).
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 It is common to find considerable genetic variation in susceptibility to infection in natural populations. This study has investigated whether natural selection increases this variation by testing whether host populations show more genetic variation in susceptibility to pathogens that they naturally encounter than novel pathogens. In a large cross-infection experiment involving four species of Drosophila and four host-specific viruses, greater genetic variation was always found in susceptibility to viruses that had coevolved with their host. The genetic architecture of resistance was examined in one host species, finding that there are more major-effect genetic variants in coevolved host-pathogen interactions. It is concluded that selection by pathogens has increased genetic variation in host susceptibility, and much of this effect is caused by the occurrence of major-effect resistance polymorphisms within populations (Duxbury, 2019).
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).