Unraveling the Genetic Basis for the Rapid Diversification of Male Genitalia between Drosophila Species.

In the last 240,000 years, males of the Drosophila simulans species clade have evolved striking differences in the morphology of their epandrial posterior lobes and claspers (surstyli). These appendages are used for grasping the female during mating and so their divergence is most likely driven by sexual selection. Mapping studies indicate a highly polygenic and generally additive genetic basis for these morphological differences. However, we have limited understanding of the gene regulatory networks that control the development of genital structures and how they evolved to result in this rapid phenotypic diversification. Here, we used new D. simulans/D. mauritiana introgression lines on chromosome arm 3L to generate higher resolution maps of posterior lobe and clasper differences between these species. We then carried out RNA-seq on the developing genitalia of both species to identify the expressed genes and those that are differentially expressed between the two species. This allowed us to test the function of expressed positional candidates during genital development in D. melanogaster. We identified several new genes involved in the development and possibly the evolution of these genital structures, including the transcription factors Hairy and Grunge. Furthermore, we discovered that during clasper development Hairy negatively regulates tartan (trn), a gene known to contribute to divergence in clasper morphology. Taken together, our results provide new insights into the regulation of genital development and how this has evolved between species.

tartan underlies the evolution of Drosophila male genital morphology.

Male genital structures are among the most rapidly evolving morphological traits and are often the only features that can distinguish closely related species. This process is thought to be driven by sexual selection and may reinforce species separation. However, while the genetic bases of many phenotypic differences have been identified, we still lack knowledge about the genes underlying evolutionary differences in male genital organs and organ size more generally. The claspers (surstyli) are periphallic structures that play an important role in copulation in insects. Here, we show that divergence in clasper size and bristle number between Drosophila mauritiana and Drosophila simulans is caused by evolutionary changes in tartan (trn), which encodes a transmembrane leucine-rich repeat domain protein that mediates cell-cell interactions and affinity. There are no fixed amino acid differences in trn between D. mauritiana and D. simulans, but differences in the expression of this gene in developing genitalia suggest that cis-regulatory changes in trn underlie the evolution of clasper morphology in these species. Finally, analyses of reciprocal hemizygotes that are genetically identical, except for the species from which the functional allele of trn originates, determined that the trn allele of D. mauritiana specifies larger claspers with more bristles than the allele of D. simulans Therefore, we have identified a gene underlying evolutionary change in the size of a male genital organ, which will help to better understand not only the rapid diversification of these structures, but also the regulation and evolution of organ size more broadly.

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles.

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

Gene regulatory network architecture in different developmental contexts influences the genetic basis of morphological evolution.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.

Human eye conditions: insights from the fly eye.

The fruit fly Drosophila melanogaster has served as an excellent model to study and understand the genetics of many human diseases from cancer to neurodegeneration. Studying the regulation of growth, determination and differentiation of the compound eyes of this fly, in particular, have provided key insights into a wide range of diseases. Here we review the regulation of the development of fly eyes in light of shared aspects with human eye development. We also show how understanding conserved regulatory pathways in eye development together with the application of tools for genetic screening and functional analyses makes Drosophila a powerful model to diagnose and characterize the genetics underlying many human eye conditions, such as aniridia and retinitis pigmentosa. This further emphasizes the importance and vast potential of basic research to underpin applied research including identifying and treating the genetic basis of human diseases.

Sexual dimorphism and natural variation within and among species in the Drosophila retinal mosaic.

BACKGROUND: Insect compound eyes are composed of ommatidia, which contain photoreceptor cells that are sensitive to different wavelengths of light defined by the specific rhodopsin proteins that they express. The fruit fly Drosophila melanogaster has several different ommatidium types that can be localised to specific retinal regions, such as the dorsal rim area (DRA), or distributed stochastically in a mosaic across the retina, like the 'pale' and 'yellow' types. Variation in these ommatidia patterns very likely has important implications for the vision of insects and could underlie behavioural and environmental adaptations. However, despite the detailed understanding of ommatidia specification in D. melanogaster, the extent to which the frequency and distribution of the different ommatidium types vary between sexes, strains and species of Drosophila is not known. RESULTS: We investigated the frequency and distribution of ommatidium types based on rhodopsin protein expression, and the expression levels of rhodopsin transcripts in the eyes of both sexes of different strains of D. melanogaster, D. simulans and D. mauritiana. We found that while the number of DRA ommatidia was invariant, Rh3 expressing ommatidia were more frequent in the larger eyes of females compared to the males of all species analysed. The frequency and distribution of ommatidium types also differed between strains and species. The D. simulans strain ZOM4 has the highest frequency of Rh3 expressing ommatidia, which is associated with a non-stochastic patch of pale and odd-coupled ommatidia in the dorsal-posterior of their eyes. CONCLUSIONS: Our results show that there is striking variation in the frequency and distribution of ommatidium types between sexes, strains and species of Drosophila. This suggests that evolutionary changes in the underlying regulatory mechanisms can alter the distribution of ommatidium types to promote or restrict their expression in specific regions of the eye within and between species, and that this could cause differences in vision among these flies.

Sox21b underlies the rapid diversification of a novel male genital structure between Drosophila species.

The emergence and diversification of morphological novelties is a major feature of animal evolution.1,2,3,4,5,6,7,8,9 However, relatively little is known about the genetic basis of the evolution of novel structures and the mechanisms underlying their diversification. The epandrial posterior lobes of male genitalia are a novelty of particular Drosophila species.10,11,12,13 The lobes grasp the female ovipositor and insert between her abdominal tergites and, therefore, are important for copulation and species recognition.10,11,12,14,15,16,17 The posterior lobes likely evolved from co-option of a Hox-regulated gene network from the posterior spiracles10 and have since diversified in morphology in the D. simulans clade, in particular, over the last 240,000 years, driven by sexual selection.18,19,20,21 The genetic basis of this diversification is polygenic but, to the best of our knowledge, none of the causative genes have been identified.22,23,24,25,26,27,28,29,30 Identifying the genes underlying the diversification of these secondary sexual structures is essential to understanding the evolutionary impact on copulation and species recognition. Here, we show that Sox21b negatively regulates posterior lobe size. This is consistent with expanded Sox21b expression in D. mauritiana, which develops smaller posterior lobes than D. simulans. We tested this by generating reciprocal hemizygotes and confirmed that changes in Sox21b underlie posterior lobe evolution between these species. Furthermore, we found that posterior lobe size differences caused by the species-specific allele of Sox21b significantly affect copulation duration. Taken together, our study reveals the genetic basis for the sexual-selection-driven diversification of a novel morphological structure and its functional impact on copulatory behavior.

Genetic and developmental analysis of differences in eye and face morphology between Drosophila simulans and Drosophila mauritiana.

Eye and head morphology vary considerably among insects and even between closely related species of Drosophila. Species of the D. melanogaster subgroup, and other Drosophila species, exhibit a negative correlation between eye size and face width (FW); for example, D. mauritiana generally has bigger eyes composed of larger ommatidia and conversely a narrower face than its sibling species. To better understand the evolution of eye and head morphology, we investigated the genetic and developmental basis of differences in eye size and FW between male D. mauritiana and D. simulans. QTL mapping of eye size and FW showed that the major loci responsible for the interspecific variation in these traits are localized to different genomic regions. Introgression of the largest effect QTL underlying the difference in eye size resulted in flies with larger eyes but no significant difference in FW. Moreover,introgression of a QTL region on the third chromosome that contributes to the FW difference between these species affected FW, but not eye size. We also observed that this difference in FW is detectable earlier in the development of the eye­antennal disc than the difference in the size of the retinal field. Our results suggest that different loci that act at different developmental stages underlie changes in eye size and FW. Therefore, while there is a negative correlation between these traits in Drosophila, we show genetically that they also have the potential to evolve independently and this may help to explain the evolution of these traits in other insects.

Extensive paternal mtDNA leakage in natural populations of Drosophila melanogaster.

Strict maternal inheritance is considered a hallmark of animal mtDNA. Although recent reports suggest that paternal leakage occurs in a broad range of species, it is still considered an exceptionally rare event. To evaluate the impact of paternal leakage on the evolution of mtDNA, it is essential to reliably estimate the frequency of paternal leakage in natural populations. Using allele-specific real-time quantitative PCR (RT-qPCR), we show that heteroplasmy is common in natural populations with at least 14% of the individuals carrying multiple mitochondrial haplotypes. However, the average frequency of the minor mtDNA haplotype is low (0.8%), which suggests that this pervasive heteroplasmy has not been noticed before due to a lack of power in sequencing surveys. Based on the distribution of mtDNA haplotypes in the offspring of heteroplasmic mothers, we found no evidence for strong selection against one of the haplotypes. We estimated that the rate of paternal leakage is 6% and that at least 100 generations are required for complete sorting of mtDNA haplotypes. Despite the high proportion of heteroplasmic individuals in natural populations, we found no evidence for recombination between mtDNA molecules, suggesting that either recombination is rare or recombinant haplotypes are counter-selected. Our results indicate that evolutionary studies using mtDNA as a marker might be biased by paternal leakage in this species.

DrosoPhyla: Resources for Drosophilid Phylogeny and Systematics.

The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4,400 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the genera, subgenera, and species groups in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred a species tree with highly supported groups for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.

Multiple hybridization events between Drosophila simulans and Drosophila mauritiana are supported by mtDNA introgression.

The study of speciation has advanced considerably in the last decades because of the increased application of molecular tools. In particular, the quantification of gene flow between recently diverged species could be addressed. Drosophila simulans and Drosophila mauritiana diverged, probably allopatrically, from a common ancestor approximately 250,000 years ago. However, these species share one mitochondrial DNA (mtDNA) haplotype indicative of a recent episode of introgression. To study the extent of gene flow between these species, we took advantage of a large sample of D. mauritiana and employed a range of different markers, i.e. nuclear and mitochondrial sequences, and microsatellites. This allowed us to detect two new mtDNA haplotypes (MAU3 and MAU4). These haplotypes diverged quite recently from haplotypes of the siII group present in cosmopolitan populations of D. simulans. The mean divergence time of the most diverged haplotype (MAU4) is approximately 127,000 years, which is more than 100,000 years before the assumed speciation time. Interestingly, we also found some evidence for gene flow at the nuclear level because an excess of putatively neutral loci shows significantly reduced differentiation between D. simulans and D. mauritiana. Our results suggest that these species are exchanging genes more frequently than previously thought.

Characterization of the Genetic Architecture Underlying Eye Size Variation Within Drosophila melanogaster and Drosophila simulans.

The compound eyes of insects exhibit striking variation in size, reflecting adaptation to different lifestyles and habitats. However, the genetic and developmental bases of variation in insect eye size is poorly understood, which limits our understanding of how these important morphological differences evolve. To address this, we further explored natural variation in eye size within and between four species of the Drosophila melanogaster species subgroup. We found extensive variation in eye size among these species, and flies with larger eyes generally had a shorter inter-ocular distance and vice versa We then carried out quantitative trait loci (QTL) mapping of intra-specific variation in eye size and inter-ocular distance in both D. melanogaster and D. simulans This revealed that different genomic regions underlie variation in eye size and inter-ocular distance in both species, which we corroborated by introgression mapping in D. simulans This suggests that although there is a trade-off between eye size and inter-ocular distance, variation in these two traits is likely to be caused by different genes and so can be genetically decoupled. Finally, although we detected QTL for intra-specific variation in eye size at similar positions in D. melanogaster and D. simulans, we observed differences in eye fate commitment between strains of these two species. This indicates that different developmental mechanisms and therefore, most likely, different genes contribute to eye size variation in these species. Taken together with the results of previous studies, our findings suggest that the gene regulatory network that specifies eye size has evolved at multiple genetic nodes to give rise to natural variation in this trait within and among species.

Nonrandom Wolbachia infection status of Drosophila melanogaster strains with different mtDNA haplotypes.

Wolbachia are maternally inherited bacteria, which typically spread in the host population by inducing cytoplasmic incompatibility (CI). In Drosophila melanogaster, Wolbachia is quite common but CI is variable, with most of the studies reporting low levels of CI. Surveying mitochondrial DNA (mtDNA) variation and infection status in a worldwide D. melanogaster collection, we found that the Wolbachia infection was not randomly distributed among flies with different mtDNA haplotypes. This preferential infection of some mtDNA haplotypes could be caused by a recent spread of mtDNA haplotypes associated with the infection. The comparison of contemporary D. melanogaster samples with lines collected more than 50 years ago shows that indeed one haplotype with a high incidence of Wolbachia infection has increased in frequency. Consistent with this observation, we found that the acquisition of a Wolbachia infection in a population from Crete was accompanied with an almost complete mtDNA replacement, with the Wolbachia-associated haplotype becoming abundant. Although it is difficult to identify the evolutionary forces causing the global increase of wMel, the parallel sweep of Wolbachia and an mtDNA haplotype suggests a fitness advantage of the Wolbachia infection.

Contrasting patterns of natural variation in global Drosophila melanogaster populations.

Despite the popularity of Drosophila melanogaster in functional and evolutionary genetics, the global pattern of natural variation has not yet been comprehensively described in this species. For the first time, we report a combined survey using neutral microsatellites and mitochondrial sequence variation jointly. Thirty-five populations originating from five continents were compared. In agreement with previous microsatellite studies, sub-Saharan African populations were the most variable ones. Consistent with previous reports of a single 'out of Africa' habitat expansion, we found that non-African populations contained a subset of the African alleles. The pattern of variation detected for the mitochondrial sequences differed substantially. The most divergent haplotypes were detected in the Mediterranean region while Africa harboured most haplotypes, which were all closely related. In the light of the well-established African origin of D. melanogaster, our results cast severe doubts about the suitability of mtDNA for biogeographic inference in this model organism.

Genetic architecture and functional characterization of genes underlying the rapid diversification of male external genitalia between Drosophila simulans and Drosophila mauritiana.

Male sexual characters are often among the first traits to diverge between closely related species and identifying the genetic basis of such changes can contribute to our understanding of their evolutionary history. However, little is known about the genetic architecture or the specific genes underlying the evolution of male genitalia. The morphology of the claspers, posterior lobes, and anal plates exhibit striking differences between Drosophila mauritiana and D. simulans. Using QTL and introgression-based high-resolution mapping, we identified several small regions on chromosome arms 3L and 3R that contribute to differences in these traits. However, we found that the loci underlying the evolution of clasper differences between these two species are independent from those that contribute to posterior lobe and anal plate divergence. Furthermore, while most of the loci affect each trait in the same direction and act additively, we also found evidence for epistasis between loci for clasper bristle number. In addition, we conducted an RNAi screen in D. melanogaster to investigate if positional and expression candidate genes located on chromosome 3L, are also involved in genital development. We found that six of these genes, including components of Wnt signaling and male-specific lethal 3 (msl3), regulate the development of genital traits consistent with the effects of the introgressed regions where they are located and that thus represent promising candidate genes for the evolution these traits.

A perspective on micro-evo-devo: progress and potential.

The term "micro-evo-devo" refers to the combined study of the genetic and developmental bases of natural variation in populations and the evolutionary forces that have shaped this variation. It thus represents a synthesis of the fields of evolutionary developmental biology and population genetics. As has been pointed out by several others, this synthesis can provide insights into the evolution of organismal form and function that have not been possible within these individual disciplines separately. Despite a number of important successes in micro-evo-devo, however, it appears that evo devo and population genetics remain largely separate spheres of research, limiting their ability to address evolutionary questions. This also risks pushing contemporary evo devo to the fringes of evolutionary biology because it does not describe the causative molecular changes underlying evolution or the evolutionary forces involved. Here we reemphasize the theoretical and practical importance of micro-evo-devo as a strategy for understanding phenotypic evolution, review the key recent insights that it has provided, and present a perspective on both the potential and the remaining challenges of this exciting interdisciplinary field.

Evolution of mir-92a underlies natural morphological variation in Drosophila melanogaster.

Identifying the genetic mechanisms underlying phenotypic change is essential to understanding how gene regulatory networks and ultimately the genotype-to-phenotype map evolve. It is recognized that microRNAs (miRNAs) have the potential to facilitate evolutionary change [1-3]; however, there are no known examples of natural morphological variation caused by evolutionary changes in miRNA expression. Therefore, the contribution of miRNAs to evolutionary change remains unknown [1, 4]. Drosophila melanogaster subgroup species display a portion of trichome-free cuticle on the femur of the second leg called the "naked valley." It was previously shown that Ultrabithorax (Ubx) is involved in naked valley variation between D. melanogaster and D. simulans [5, 6]. However, naked valley size also varies among populations of D. melanogaster, ranging from 1,000 up to 30,000 µm(2). We investigated the genetic basis of intraspecific differences in the naked valley in D. melanogaster and found that neither Ubx nor shavenbaby (svb) [7, 8] contributes to this morphological difference. Instead, we show that changes in mir-92a expression underlie the evolution of naked valley size in D. melanogaster through repression of shavenoid (sha) [9]. Therefore, our results reveal a novel mechanism for morphological evolution and suggest that modulation of the expression of miRNAs potentially plays a prominent role in generating organismal diversity.

Evolution of eye morphology and rhodopsin expression in the Drosophila melanogaster species subgroup.

A striking diversity of compound eye size and shape has evolved among insects. The number of ommatidia and their size are major determinants of the visual sensitivity and acuity of the compound eye. Each ommatidium is composed of eight photoreceptor cells that facilitate the discrimination of different colours via the expression of various light sensitive Rhodopsin proteins. It follows that variation in eye size, shape, and opsin composition is likely to directly influence vision. We analyzed variation in these three traits in D. melanogaster, D. simulans and D. mauritiana. We show that D. mauritiana generally has larger eyes than its sibling species, which is due to a combination of larger ommatidia and more ommatidia. In addition, intra- and inter-specific differences in eye size among D. simulans and D. melanogaster strains are mainly caused by variation in ommatidia number. By applying a geometric morphometrics approach to assess whether the formation of larger eyes influences other parts of the head capsule, we found that an increase in eye size is associated with a reduction in the adjacent face cuticle. Our shape analysis also demonstrates that D. mauritiana eyes are specifically enlarged in the dorsal region. Intriguingly, this dorsal enlargement is associated with enhanced expression of rhodopsin 3 in D. mauritiana. In summary, our data suggests that the morphology and functional properties of the compound eyes vary considerably within and among these closely related Drosophila species and may be part of coordinated morphological changes affecting the head capsule.

Variability patterns and positively selected sites at the gametophytic self-incompatibility pollen SFB gene in a wild self-incompatible Prunus spinosa (Rosaceae) population.

Current models for the generation of new gametophytic self-incompatibility specificities require that neutral variability segregates within specificity classes. Furthermore, one of the models predicts greater ratios of nonsynonymous to synonymous substitutions in pollen than in pistil specificity genes. All models assume that new specificities arise by mutation only. To test these models, 21 SFB (the pollen S-locus) alleles from a wild Prunus spinosa (Rosaceae) population were obtained. For seven of these, the corresponding S-haplotype was also characterized. The SFB data set was also used to identify positively selected sites. Those sites are likely to be the ones responsible for defining pollen specificities. Of the 23 sites identified as being positively selected, 21 are located in the variable (including a new region described here) and hypervariable regions. Little variability is found within specificity classes. There is no evidence for selective sweeps being more frequent in pollen than in pistil specificity genes. The S-RNase and the SFB genes have only partially correlated evolutionary histories. None of the models is compatible with the variability patterns found in the SFB and the S-haplotype data.