Phylogenetic relationships among species groups of the virilis-repleta radiation of Drosophila

Whole genome phylogenomics helps to resolve the phylogenetic position of the Zygothrica genus group (Diptera, Drosophilidae) and the causes of previous incongruences

Incomplete Lineage Sorting (ILS) and introgression are among the two main factors causing incongruence between gene and species trees. Advances in phylogenomic studies have allowed us to overcome most of these issues, providing reliable phylogenetic hypotheses while revealing the underlying evolutionary scenario. Across the last century, many incongruent phylogenetic reconstructions were recovered for Drosophilidae, employing a limited sampling of genetic markers or species. In these studies, the monophyly and the phylogenetic positioning of the Zygothrica genus group stood out as one of the most controversial questions. Thus, here, we addressed these issues using a phylogenomic approach, while accessing the influence of ILS and introgressions on the diversification of these species and addressing the spatio-temporal scenario associated with their evolution. For this task, the genomes of nine specimens from six Neotropical species belonging to the Zygothrica genus group were sequenced and evaluated in a phylogenetic framework encompassing other 39 species of Drosophilidae. Nucleotide and amino acid sequences recovered for a set of 2,534 single-copy genes by BUSCO were employed to reconstruct maximum likelihood (ML) concatenated and multi-species coalescent (MSC) trees. Likelihood mapping, quartet sampling, and reticulation tests were employed to infer the level and causes of incongruence. Lastly, a penalized-likelihood molecular clock strategy with fossil calibrations was performed to infer divergence times. Taken together, our results recovered the subdivision of Drosophila into six different lineages, one of which clusters species of the Zygothrica genus group (except for H. duncani). The divergence of this lineage was dated to Oligocene ∼ 31 Mya and seems to have occurred in the same timeframe as other key diversification within Drosophila. According to the concatenated and MSC strategies, this lineage is sister to the clade joining Drosophila (Siphlodora) with the Hawaiian Drosophila and Scaptomyza. Likelihood mapping, quartet sampling, reticulation reconstructions as well as introgression tests revealed that this lineage was the target of several hybridization events involving the ancestors of different Drosophila lineages. Thus, our results generally show introgression as a major source of previous incongruence. Nevertheless, the similar diversification times recovered for several of the Neotropical Drosophila lineages also support the scenario of multiple and simultaneous diversifications taking place at the base of Drosophilidae phylogeny, at least in the Neotropics.

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

Phylogeny of the Genus Drosophila

Understanding phylogenetic relationships among taxa is key to designing and implementing comparative analyses. The genus Drosophila, which contains over 1600 species, is one of the most important model systems in the biological sciences. For over a century, one species in this group, Drosophila melanogaster, has been key to studies of animal development and genetics, genome organization and evolution, and human disease. As whole-genome sequencing becomes more cost-effective, there is increasing interest in other members of this morphologically, ecologically, and behaviorally diverse genus. Phylogenetic relationships within Drosophila are complicated, and the goal of this paper is to provide a review of the recent taxonomic changes and phylogenetic relationships in this genus to aid in further comparative studies.

The Calibrated Phylogeny of the Drosophila fasciola Subgroup (D. repleta Group Wasserman) Indicates Neogene Diversification of Its Internal Branches

The species of the Drosophila fasciola subgroup Wasserman represent the dominant section of the Drosophila repleta group Wasserman in the American rainforests and have a broad geographical distribution in the New World. However, despite of its wide range, the D. fasciola subgroup is one of the most overlooked D. repleta subgroups. Here, we report a molecular phylogenetic analysis focused on the D. fasciola subgroup using two mitochondrial [cytochrome oxidase subunit I (COI), cytochrome oxidase subunit II (COII)] and two nuclear [elongation factor-1alpha F1 (EF-alphaF1) and transformer (tra)] genes. Overall, we found that this subgroup is a monophyletic taxon, subdivided into two main internal branches: named Fas1 and Fas2 clades. The diversification of these clades is estimated to have begun in the middle Miocene, around 12 Ma [95% high posterior density (HPD) 9.0-15 Ma], and might be associated with the colonization of South America by Central America populations after the closure of Isthmus of Panama due to the temporal congruence between these events. The terminal branches had their origins estimated to be in the Pliocene or the Plio-Pleistocene transition. For the later estimates, both the geomorphological influences and the climatic oscillations of the Pleistocene may have played a role in shaping the diversification of the D. fasciola group.

Phylogeography of the Subgenus Drosophila (Diptera: Drosophilidae): Evolutionary History of Faunal Divergence between the Old and the New Worlds

The current subgenus Drosophila (the traditional immigrans-tripunctata radiation) includes major elements of temperate drosophilid faunas in the northern hemisphere. Despite previous molecular phylogenetic analyses, the phylogeny of the subgenus Drosophila has not fully been resolved: the resulting trees have more or less varied in topology. One possible factor for such ambiguous results is taxon-sampling that has been biased towards New World species in previous studies. In this study, taxon sampling was balanced between Old and New World species, and phylogenetic relationships among 45 ingroup species selected from ten core species groups of the subgenus Drosophila were analyzed using nucleotide sequences of three nuclear and two mitochondrial genes. Based on the resulting phylogenetic tree, ancestral distributions and divergence times were estimated for each clade to test Throckmorton’s hypothesis that there was a primary, early-Oligocene disjunction of tropical faunas and a subsequent mid-Miocene disjunction of temperate faunas between the Old and the New Worlds that occurred in parallel in separate lineages of the Drosophilidae. Our results 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, our 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.

Phylogenetic patterns of geographical and ecological diversification in the subgenus Drosophila

Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. Members of the subgenus Drosophila are distributed across the globe and show a large diversity of ecological niches. Furthermore, taxonomic classification of Drosophila includes the rank radiation, which refers to closely related species groups. Nevertheless, it has never been tested if these taxonomic radiations correspond to evolutionary radiations. Here we present a study of the patterns of diversification of Drosophila to test for increased diversification rates in relation to the geographic and ecological diversification processes. For this, we have estimated and dated a phylogeny of 218 species belonging to the major species groups of the subgenus. The obtained phylogenies are largely consistent with previous studies and indicate that the major groups appeared during the Oligocene/Miocene transition or early Miocene, characterized by a trend of climate warming with brief periods of glaciation. Ancestral reconstruction of geographic ranges and ecological resource use suggest at least two dispersals to the Neotropics from the ancestral Asiatic tropical disribution, and several transitions to specialized ecological resource use (mycophagous and cactophilic). Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. However, diversification analyses show no significant support for adaptive radiations as a result of geographic dispersal or ecological resource shift. Also, cactophily has not resulted in an increase in the diversification rate of the repleta and related groups. It is thus concluded that the taxonomic radiations do not correspond to adaptive radiations.

A supermatrix-based molecular phylogeny of the family Drosophilidae

The genus Drosophila is diverse and heterogeneous and contains a large number of easy-to-rear species, so it is an attractive subject for comparative studies. The ability to perform such studies is currently compromised by the lack of a comprehensive phylogeny for Drosophila and related genera. The genus Drosophila as currently defined is known to be paraphyletic with respect to several other genera, but considerable uncertainty remains about other aspects of the phylogeny. Here, we estimate a phylogeny for 176 drosophilid (12 genera) and four non-drosophilid species, using gene sequences for up to 13 different genes per species (average: 4333 bp, five genes per species). This is the most extensive set of molecular data on drosophilids yet analysed. Phylogenetic analyses were conducted with maximum-likelihood (ML) and Bayesian approaches. Our analysis confirms that the genus Drosophila is paraphyletic with 100% support in the Bayesian analysis and 90% bootstrap support in the ML analysis. The subgenus Sophophora, which includes Drosophila melanogaster, is the sister clade of all the other subgenera as well as of most species of six other genera. This sister clade contains two large, well-supported subclades. The first subclade contains the Hawaiian Drosophila, the genus Scaptomyza, and the virilis-repleta radiation. The second contains the immigrans-tripunctata radiation as well as the genera Hirtodrosophila (except Hirtodrosophila duncani), Mycodrosophila, Zaprionus and Liodrosophila. We argue that these results support a taxonomic revision of the genus Drosophila.

Phylogeny of the Drosophila mesophragmatica group (Diptera, Drosophilidae): an example of Andean evolution

The mesophragmatica group of Drosophila belongs to the virilis-repleta radiation of the Drosophila subgenus. This group comprises 13 Neotropical species that are endemic to the South-American continent and seem to be fundamentally Andean in their distribution. The mesophragmatica-group phylogeny has been inferred previously by other authors based on morphological, cytological, and isozyme analyses. However, the relationships within the group have not yet been completely resolved, although its monophyletic origin has already been confirmed by molecular data. This work attempts to enhance the molecular approach to the relationships among the species of the mesophragmatica group, using both nuclear and mitochondrial markers. Phylogenetic analyses were performed using fragments of the nuclear alcohol dehydrogenase (Adh; 631 bp), alpha-methyldopa (Amd; 1211 bp), dopa-decarboxylase (Ddc; 1105 bp), and hunchback (Hb; 687 bp) genes and the mitochondrial cytochrome oxidase subunit II (COII; 672 bp) gene, and included a total of 4306 bp. The sequences obtained for eight representatives of the mesophragmatica group were analyzed both individually and in combination by distance methods, maximum parsimony, and maximum likelihood. Our results support subdivision of the mesophragmatica group into three main lineages: the first is composed of D. viracochi; the second comprises a clade grouping the sibling species D. pavani and D. gaucha; and the third encompasses D. gasici, D. brncici, and D. mesophragmatica. The best supported scenario suggests that D. viracochi is an early offshoot in the mesophragmatica group, with this and other early branchings occuring in the Pliocene/Pleistocene Epochs, possibly associated with Andean glacial refuges. Also based on the phylogenies obtained, we present a genealogical view of the evolution of previously described characters within the group.

Adaptation to long sperm in Drosophila: correlated development of the sperm roller and sperm packaging

Sperm are generally small and produced in huge numbers, but some species combine exaggerated sperm length with extremely limited numbers of sperm, an evolutionary trend that deviates from the theory of anisogamy. Sperm gigantism has arisen recurrently in various species, but insects exhibit the longest sperm, with some species of the Drosophilidae family producing sperm up to 6 cm in length. The anatomical, cytological, and physiological requirements for males to cope with these giant sperm were hitherto poorly understood. In this paper, we investigate the internal morphology of the male reproductive tract, and highlight specific features that may be linked to this increase in sperm size. We focus on species in the repleta group, within which sperm length varies by a factor of 35. An associated development of the sperm roller, a special twisting device inserted between the testis and the seminal vesicle, is demonstrated. Its length and the number of coils involved increase with sperm size, and it allows individual sperm to swell and roll into a spermatic pellet before reaching the seminal vesicle. This process occurs independently of and in addition to the sperm bundle coiling that takes place at the base of the testis. It is suggested that the emergence and development of the sperm roller may be a male adaptation to sperm gigantism.

Independent origins of new sex-linked chromosomes in the melanica and robusta species groups of Drosophila

Background

Recent translocations of autosomal regions to the sex chromosomes represent important systems for identifying the evolutionary forces affecting convergent patterns of sex-chromosome heteromorphism. Additions to the sex chromosomes have been reported in the melanica and robusta species groups, two sister clades of Drosophila. The close relationship between these two species groups and the similarity of their rearranged karyotypes motivates this test of alternative hypotheses; the rearranged sex chromosomes in both groups are derived through a common origin, or the rearrangements are derived through at least two independent origins. Here we examine chromosomal arrangement in representatives of the melanica and the robusta species groups and test these alternative hypotheses using a phylogenetic approach.

Results

Two mitochondrial and two nuclear gene sequences were used to reconstruct phylogenetic relationships of a set of nine ingroup species having fused and unfused sex chromosomes and representing a broad sample of both species groups. Different methods of phylogenetic inference, coupled with concurrent cytogenetic analysis, indicate that the hypothesis of independent origins of rearranged sex chromosomes within each species group is significantly more likely than the alternative hypothesis of a single common origin. An estimate tightly constrained around 8 My was obtained for the age of the rearranged sex chromosomes in the melanica group; however, a more loosely constrained estimate of 10–15 My was obtained for the age of the rearrangement in the robusta group.

Conclusion

Independent acquisition of new chromosomal arms by the sex chromosomes in the melanica and robusta species groups represents a case of striking convergence at the karyotypic level. Our findings indicate that the parallel divergence experienced by newly sex-linked genomic regions in these groups represents an excellent system for studying the tempo of sex chromosome evolution.

Molecular phylogeny of the Drosophila virilis section (Diptera: Drosophilidae) based on mitochondrial and nuclear sequences

Regardless of the well-documented virilis species group, most groups of the Drosophila virilis section have not been completely studied at molecular level since it was suggested. Therefore, phylogenetic relationships among and within species groups of the virilis section are generally unknown. In present paper, the complete mitochondrial ND2 gene and fragment of COI gene in combination with a nuclear gene, Adh coding region, were used to derive the most extensive molecular phylogeny to date for the Drosophila virilis section. A total of 111 individuals covering 61 species were sampled in this study. Novel phylogenetic findings included (1) support for the paraphyly of the melanica and robusta species group and at least two subgroups of the robusta species group, the lacertosa and okadai subgroups, were distinguished as paraphyletic taxa. In addition, (2) present results revealed the sister relationship between D. moriwakii and the robusta subgroup, conflicting with current taxonomy regarding D. moriwakii, which was shifted from the robusta species group to the melanica group. (3) In contrast to the robusta and melanica species groups, monophyly of the polychaeta species group, the angor group and the virilis group was confirmed, respectively. However, the monophyletic quadrisetata species group was resolved with uncertainty. (4) Our analyses of combined data set suggested close relationship between the quadrisetata species group and the unpublished clefta group, and the okadai subgroup is sister to the clade comprising of the quadrisetata and clefta species groups. Within the virilis section, D. fluvialis and three tropical species groups, the polychaeta group, the angor group and the repleta group, are found to branch off earlier than other ingroup taxa. This suggests that the virilis section might have originated in the Old World tropics. Besides, the derived status of the close affinities of the quadrisetata group, the clefta group, and the melanica and robusta groups is probably the result of their adaptation to forests between subtropical and cool-temperate climate. Based on the consideration of the phylogenetic placement of the species of the virilis section, we suggest that at least five independent migrations occurred from the Old World to the New World.

Evolutionary genetics of reproductive behavior in Drosophila: connecting the dots

Species of the genus Drosophila exhibit enormous variation in all of their reproductive behaviors: resource use and specialization, courtship signaling, sperm utilization, and female remating. The genetic bases of this variability and its evolution are poorly understood. At the same time, Drosophila comparative genomics now has developed to a point at which approaches previously only possible with D. melanogaster can be exploited to address these questions. We have taken advantage of the known phylogenetic relationships of this group of flies not only to place these behaviors in an evolutionary framework, but to provide a roadmap for future genetic studies.

Dopa decarboxylase: a model gene-enzyme system for studying development, behavior, and systematics

Throughout its long evolutionary history, the Dopa decarboxylase gene (Ddc) has acquired a variety of functions in insects. The enzyme (DDC) catalyzes the production of the neural transmitters dopamine and serotonin. Not surprisingly, evidence of the enzyme's involvement in the behavior of insects is beginning to accumulate. In addition, DDC plays a role in wound healing, parasite defense, pigmentation, and cuticle hardening. A high degree of sequence conservation has allowed comparisons of the Ddc-coding regions from various insects, facilitating a number of recent studies on insect systematics. This review outlines the diverse functions of Ddc and illustrates how studies of this model system address many questions on insect neurobiology, developmental biology, and systematics.

Phylogeny and molecular evolution of the Drosophila hydei subgroup (Drosophila repleta group) inferred from the Xanthine dehydrogenase gene

The hydei subgroup (Drosophila repleta group) consists of seven species divided into two complexes: bifurca and hydei, whose phylogenetic relationships are not well understood. To evaluate the molecular phylogeny of this subgroup, we analyzed 2085 bp of coding sequence of the Xanthine dehydrogenase gene in six available species of the hydei subgroup, with Drosophila buzzatii and Drosophila mulleri as an outgroup. For phylogenetic reconstruction we adopted a maximum-likelihood framework, based on the adjustment of descriptive models of nucleotide substitution to real data. We employed distance-based and weighted parsimony methods to construct candidate phylogenies. In all cases, we obtained only one completely resolved tree with strong statistical support for each node, that shows a phylogeny that is partially discordant with the proposed systematics of the subgroup. This tree suggests that the two species complexes are paraphyletic, as opposed to classic phylogenies using morphologic and cytologic traits. This discordance is discussed in relation to its implication for the evolutionary history of the hydei subgroup.

Molecular phylogeny of the subgenus Drosophila (Diptera, Drosophilidae) with an emphasis on Neotropical species and groups: A nuclear versus mitochondrial gene approach

The genus Drosophila has played an essential role in many biological studies during the last 100 years but much controversy and many incompletely addressed issues still remain to be elucidated regarding the phylogeny of this genus. Because information on the Neotropical species contained in the subgenus Drosophila is particularly incomplete, with this taxonomic group being underrepresented in many studies, we designed a study to answer some evolutionary questions related to these species. We subjected at least 41 Drosophilidae taxa to a phylogenetic analysis using a 516-base pair (bp) fragment of the alpha-methyldopa (Amd) nuclear gene and a 672 bp fragment of the mitochondrial cytochrome oxidase subunit II (COII) gene both individually and in combination. We found that the subgenus Drosophila is paraphyletic and subdivided into two main clusters: the first containing species traditionally placed in the virilis-repleta radiation and the second assembling species of the immigrans-Hirtodrosophila radiation. Inside the first of these clusters we could detect the monophyly of both the flavopilosa (the sister-clade of the annulimana group) and the mesophragmatica (closely related to the repleta group) species groups. Concerning the immigrans-Hirtodrosophila lineage, Zaprionus, Liodrosophila, Samoaia, and Hirtodrosophila were the early offshoots, followed by the immigrans, quinaria, testacea, and funebris species groups. The tripunctata radiation appears to be a derived clade, composed of a paraphyletic tripunctata group, intimately interposed with members of the cardini, guarani, and guaramunu species groups. Overall, the COII gene yielded a poor phylogenetic performance when compared to the Amd gene, the evolutionary hypothesis of which agreed with the total evidence tree. This phenomenon can be explained by the fast saturation of transitional substitutions in COII, due to strong biases in both base composition and substitution patterns, as also by its great among-site rate variation heterogeneity.

Automated measurement of Drosophila wings

Background

Many studies in evolutionary biology and genetics are limited by the rate at which phenotypic information can be acquired. The wings of Drosophila species are a favorable target for automated analysis because of the many interesting questions in evolution and development that can be addressed with them, and because of their simple structure.

Results

We have developed an automated image analysis system (WINGMACHINE) that measures the positions of all the veins and the edges of the wing blade of Drosophilid flies. A video image is obtained with the aid of a simple suction device that immobilizes the wing of a live fly. Low-level processing is used to find the major intersections of the veins. High-level processing then optimizes the fit of an a priori B-spline model of wing shape. WINGMACHINE allows the measurement of 1 wing per minute, including handling, imaging, analysis, and data editing. The repeatabilities of 12 vein intersections averaged 86% in a sample of flies of the same species and sex.

Comparison of 2400 wings of 25 Drosophilid species shows that wing shape is quite conservative within the group, but that almost all taxa are diagnosably different from one another. Wing shape retains some phylogenetic structure, although some species have shapes very different from closely related species. The WINGMACHINE system facilitates artificial selection experiments on complex aspects of wing shape. We selected on an index which is a function of 14 separate measurements of each wing. After 14 generations, we achieved a 15 S.D. difference between up and down-selected treatments.

Conclusion

WINGMACHINE enables rapid, highly repeatable measurements of wings in the family Drosophilidae. Our approach to image analysis may be applicable to a variety of biological objects that can be represented as a framework of connected lines.