A Major Locus Controls a Genital Shape Difference Involved in Reproductive Isolation Between Drosophila yakuba and Drosophila santomea

A developmental atlas of male terminalia across twelve species of Drosophila

How complex morphologies evolve is one of the central questions in evolutionary biology. Observing the morphogenetic events that occur during development provides a unique perspective on the origins and diversification of morphological novelty. One can trace the tissue of origin, emergence, and even regression of structures to resolve murky homology relationships between species. Here, we trace the developmental events that shape some of the most diverse organs in the animal kingdom-the male terminalia (genitalia and analia) of Drosophilids. Male genitalia are known for their rapid evolution with closely related species of the Drosophila genus demonstrating vast variation in their reproductive morphology. We used confocal microscopy to monitor terminalia development during metamorphosis in twelve related species of Drosophila. From this comprehensive dataset, we propose a new staging scheme for pupal terminalia development based on shared developmental landmarks, which allows one to align developmental time points between species. We were able to trace the origin of different substructures, find new morphologies and suggest possible homology of certain substructures. Additionally, we demonstrate that posterior lobe is likely originated prior to the split between the Drosophila melanogaster and the Drosophila yakuba clade. Our dataset opens up many new directions of research and provides an entry point for future studies of the Drosophila male terminalia evolution and development.

Resolving between novelty and homology in the rapidly evolving phallus of Drosophila

The genitalia present some of the most rapidly evolving anatomical structures in the animal kingdom, possessing a variety of parts that can distinguish recently diverged species. In the Drosophila melanogaster group, the phallus is adorned with several processes, pointed outgrowths, that are similar in size and shape between species. However, the complex three-dimensional nature of the phallus can obscure the exact connection points of each process. Previous descriptions based upon adult morphology have primarily assigned phallic processes by their approximate positions in the phallus and have remained largely agnostic regarding their homology relationships. In the absence of clearly identified homology, it can be challenging to model when each structure first evolved. Here, we employ a comparative developmental analysis of these processes in eight members of the melanogaster species group to precisely identify the tissue from which each process forms. Our results indicate that adult phallic processes arise from three pupal primordia in all species. We found that in some cases the same primordia generate homologous structures whereas in other cases, different primordia produce phenotypically similar but remarkably non-homologous structures. This suggests that the same gene regulatory network may have been redeployed to different primordia to induce phenotypically similar traits. Our results highlight how traits diversify and can be redeployed, even at short evolutionary scales.

[Small scale evolution]

Small-scale evolution or microevolution concerns evolution at the intra-specific level or between closely related species. At the intra-specific level, it allows the analysis of the evolutionary forces at work: mutation, genetic drift, migration and selection. Moreover, because of the short evolutionary time, it is easier to identify the genetic basis of observed phenotypic differences. Most studies focus on current populations but more and more analyses are performed on ancient DNA. This provides important information for tracing the history of populations and also allows the reconstruction of phenotypes of individuals that disappeared several thousand years ago. In this short review, I present studies showing how pre-zygotic or post-zygotic barriers involved in species formation are set up using the example of the geographical barrier due to the formation of the Isthmus of Panama and that of the heterochromatin divergence in Drosophilidae. I also describe the different approaches that have been used to identify the genetic basis of well known phenotypic variations: candidate gene approach (about melanism in felines), QTL mapping (variation in the number of lateral bone plates in sticklebacks), association study (pigmentation in the Asian ladybird). Finally, I illustrate the key impact of natural selection with the iconic example of the evolution of the beak of Galapagos finches, and the role of certain developmental genes in its morphological diversification.

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome-wide dynamics of speciation-through-introgression.

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.

A morphological trait involved in reproductive isolation between Drosophila sister species is sensitive to temperature

Male genitalia are usually extremely divergent between closely related species, but relatively constant within one species. Here we examine the effect of temperature on the shape of the ventral branches, a male genital structure involved in reproductive isolation, in the sister species Drosophila santomea and Drosophila yakuba. We designed a semi-automatic measurement machine learning pipeline that can reliably identify curvatures and landmarks based on manually digitized contours of the ventral branches. With this method, we observed that temperature does not affect ventral branches in D. yakuba but that in D. santomea ventral branches tend to morph into a D. yakuba-like shape at lower temperature. We found that male genitalia structures involved in reproductive isolation can be relatively variable within one species and can resemble the shape of closely related species' genitalia through plasticity to temperature. Our results suggest that reproductive isolation mechanisms can be dependent on the environmental context.

Evolutionary expansion of apical extracellular matrix is required for the elongation of cells in a novel structure

One of the fundamental gaps in our knowledge of how novel anatomical structures evolve is understanding the origins of the morphogenetic processes that form these features. Here, we traced the cellular development of a recently evolved morphological novelty, the posterior lobe of D. melanogaster. We found that this genital outgrowth forms through extreme increases in epithelial cell height. By examining the apical extracellular matrix (aECM), we also uncovered a vast matrix associated with the developing genitalia of lobed and non-lobed species. Expression of the aECM protein Dumpy is spatially expanded in lobe-forming species, connecting the posterior lobe to the ancestrally derived aECM network. Further analysis demonstrated that Dumpy attachments are necessary for cell height increases during posterior lobe development. We propose that the aECM presents a rich reservoir for generating morphological novelty and highlights a yet unseen role for aECM in regulating extreme cell height.

An Atlas of Transcription Factors Expressed in Male Pupal Terminalia of Drosophila melanogaster

During development, transcription factors and signaling molecules govern gene regulatory networks to direct the formation of unique morphologies. As changes in gene regulatory networks are often implicated in morphological evolution, mapping transcription factor landscapes is important, especially in tissues that undergo rapid evolutionary change. The terminalia (genital and anal structures) of Drosophila melanogaster and its close relatives exhibit dramatic changes in morphology between species. While previous studies have identified network components important for patterning the larval genital disc, the networks governing adult structures during pupal development have remained uncharted. Here, we performed RNA-seq in whole Drosophila melanogaster male terminalia followed by in situ hybridization for 100 highly expressed transcription factors during pupal development. We find that the male terminalia are highly patterned during pupal stages and that specific transcription factors mark separate structures and substructures. Our results are housed online in a searchable database (https://flyterminalia.pitt.edu/) as a resource for the community. This work lays a foundation for future investigations into the gene regulatory networks governing the development and evolution of Drosophila terminalia.

Quantifying the extent of morphological homoplasy: A phylogenetic analysis of 490 characters in Drosophila

Homoplasy is a fundamental phenomenon in evolutionary biology but an appraisal of its extent at the morphological level is still lacking. Here, we analyzed the evolution of 490 morphological characters conceptualized among 56 drosophilid species. We found that two thirds of morphological changes were homoplastic and that the level of homoplasy depended on the stage of development and the type of the organ, with the adult terminalia being the least homoplastic. In spite of its predominance at the character change level, homoplasy accounts for only ∼13% of between species similarities in pairwise comparisons. These results provide empirical insights on the limits of morphological changes and the frequency of recurrent evolution.

Genetic basis of species-specific genitalia reveals role in species diversification

The diversity of genital morphology among closely related animals with internal fertilization is well known, but the genetic backgrounds are unclear. Here, we show that, in Carabus (Ohomopterus) beetles showing correlated evolution of male and female genital parts, only a few major quantitative trait loci (QTLs) determine differences in genital dimensions between sister species, and sequence divergence is pronounced in the genomic regions containing genital QTLs. The major QTLs for male and female genital dimensions reside in different locations within the same linkage group, implying that coevolution between the sexes is only loosely constrained and can respond to sexually antagonistic selection. The same genomic regions containing the major QTLs show elevated divergence between three pairs of parapatric species with marked differences in genital parts. Our study demonstrates that species diversification can follow coevolution of genitalia between the sexes, even without tight linkage of loci affecting male and female genital dimensions.

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

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

A potential pitfall in studies of biological shape: Does size matter?

The number of published studies using geometric morphometrics (GM) for analysing biological shape has increased steadily since the beginning of the 1990s, covering multiple research areas such as ecology, evolution, development, taxonomy and palaeontology. Unfortunately, we have observed that many published studies using GM do not evaluate the potential allometric effects of size on shape, which normally require consideration or assessment. This might lead to misinterpretations and flawed conclusions in certain cases, especially when size effects explain a large part of the shape variation. We assessed, for the first time and in a systematic manner, how often published studies that have applied GM consider the potential effects of allometry on shape. We reviewed the 300 most recent published papers that used GM for studying biological shape. We also estimated how much of the shape variation was explained by allometric effects in the reviewed papers. More than one-third (38%) of the reviewed studies did not consider the allometric component of shape variation. In studies where the allometric component was taken into account, it was significant in 88% of the cases, explaining up to 87.3% of total shape variation. We believe that one reason that may cause the observed results is a misunderstanding of the process that superimposes landmark configurations, i.e. the Generalized Procrustes Analysis, which removes isometric effects of size on shape, but not allometric effects. Allometry can be a crucial component of shape variation. We urge authors to address, and report, size effects in studies of biological shape. However, we do not propose to always remove size effects, but rather to evaluate the research question with and without the allometric component of shape variation. This approach can certainly provide a thorough understanding of how much size contributes to the observed shaped variation.

Genetic and Transgenic Reagents for Drosophila simulans, D. mauritiana, D. yakuba, D. santomea, and D. virilis

Species of the Drosophila melanogaster species subgroup, including the species D. simulans, D. mauritiana, D. yakuba, and D. santomea, have long served as model systems for studying evolution. However, studies in these species have been limited by a paucity of genetic and transgenic reagents. Here, we describe a collection of transgenic and genetic strains generated to facilitate genetic studies within and between these species. We have generated many strains of each species containing mapped piggyBac transposons including an enhanced yellow fluorescent protein (EYFP) gene expressed in the eyes and a ϕC31 attP site-specific integration site. We have tested a subset of these lines for integration efficiency and reporter gene expression levels. We have also generated a smaller collection of other lines expressing other genetically encoded fluorescent molecules in the eyes and a number of other transgenic reagents that will be useful for functional studies in these species. In addition, we have mapped the insertion locations of 58 transposable elements in D. virilis that will be useful for genetic mapping studies.

Drosophila yakuba mayottensis, a new model for the study of incipient ecological speciation

A full understanding of how ecological factors drive the fixation of genetic changes during speciation is obscured by the lack of appropriate models with clear natural history and powerful genetic toolkits. In a recent study, we described an early stage of ecological speciation in a population of the generalist species Drosophila yakuba (melanogaster subgroup) on the island of Mayotte (Indian Ocean). On this island, flies are strongly associated with the toxic fruits of noni (Morinda citrifolia) and show a partial degree of pre-zygotic reproductive isolation. Here, I mine the nuclear and mitochondrial genomes and provide a full morphological description of this population. Only 29 nuclear sites (< 4 × 10^-7 of the genome) are fixed in this population and absent from 3 mainland populations and the closest relative D. santomea, but no mitochondrial or morphological character distinguish Mayotte flies from the mainland. This result indicates that physiological and behavioral traits may evolve faster than morphology at the early stages of speciation. Based on these differences, the Mayotte population is designated as a new subspecies, Drosophila yakuba mayottensis subsp. nov., and its strong potential in understanding the genetics of speciation and plant-insect interactions is discussed.

Emerging developmental genetic model systems in holometabolous insects

The number of insect species that are amenable to functional genetic studies is growing rapidly and provides many new research opportunities in developmental and evolutionary biology. The holometabolous insects represent a disproportionate percentage of animal diversity and are thus well positioned to provide model species for a wide variety of developmental processes. Here we discuss emerging holometabolous models, and review some recent breakthroughs. For example, flies and midges were found to use structurally unrelated long-range pattern organizers, butterflies and moths revealed extensive pattern formation during oogenesis, new imaging possibilities in the flour beetle Tribolium castaneum showed how embryos break free of their extraembryonic membranes, and the complex genetics governing interspecies difference in head shape were revealed in Nasonia wasps.

The Genetic Basis of Baculum Size and Shape Variation in Mice

The rapid divergence of male genitalia is a preeminent evolutionary pattern. This rapid divergence is especially striking in the baculum, a bone that occurs in the penis of many mammalian species. Closely related species often display diverse baculum morphology where no other morphological differences can be discerned. While this fundamental pattern of evolution has been appreciated at the level of gross morphology, nearly nothing is known about the genetic basis of size and shape divergence. Quantifying the genetic basis of baculum size and shape variation has been difficult because these structures generally lack obvious landmarks, so comparing them in three dimensions is not straightforward. Here, we develop a novel morphometric approach to quantify size and shape variation from three-dimensional micro-CT scans taken from 369 bacula, representing 75 distinct strains of the BXD family of mice. We identify two quantitative trait loci (QTL) that explain ∼50% of the variance in baculum size, and a third QTL that explains more than 20% of the variance in shape. Together, our study demonstrates that baculum morphology may diverge relatively easily, with mutations at a few loci of large effect that independently modulate size and shape. Based on a combination of bioinformatic investigations and new data on RNA expression, we prioritized these QTL to 16 candidate genes, which have hypothesized roles in bone morphogenesis and may enable future genetic manipulation of baculum morphology.

High throughput preparation of fly genomic DNA in 96-well format using a paint-shaker

Sample homogenization is an essential step for genomic DNA extraction, with multiple downstream applications in Molecular Biology. Genotyping hundreds or thousands of samples requires an automation of this homogenization step, and high throughput homogenizer equipment currently costs 7000 euros or more. We present an apparatus for homogenization of individual Drosophila adult flies in 96-well micro-titer dishes, which was built from a small portable paint-shaker (F5 portable paint-shaker, Ushake). Single flies are disrupted in each well that contains extraction buffer and a 4-mm metal ball. Our apparatus can hold up to five 96-well micro-titer plates. Construction of the homogenizer apparatus takes about 3–4 days, and all equipment can be obtained from a home improvement store. The total material cost is approximately 700 euros including the paint-shaker. We tested the performance of our apparatus using the ZR-96 Quick-gDNA™ kit (Zymo Research) homogenization buffer and achieved nearly complete tissue homogenization after 15 minutes of shaking. PCR tests did not detect any cross contamination between samples of neighboring wells. We obtained on average 138 ng of genomic DNA per fly, and DNA quality was adequate for standard PCR applications. In principle, our tissue homogenizer can be used for isolation of DNA suitable for library production and high throughput genotyping by Multiplexed Shotgun Genotyping (MSG), as well as RNA isolation from single flies. The sample adapter can also hold and shake other items, such as centrifuge tubes (15–50 mL) or small bottles.