Developmental bias predicts 60 million years of wing shape evolution.

The degree to which developmental biases affect trait evolution is subject to much debate. Here, we first quantify fluctuating asymmetry as a measure of developmental variability, i.e., the propensity of developmental systems to create some phenotypic variants more often than others, and show that it predicts phenotypic and standing genetic variation as well as deep macroevolutionary divergence in wing shape in sepsid flies. Comparing our data to the findings of a previous study demonstrates that developmental variability in the sepsid fly Sepsis punctum strongly aligns with mutational, standing genetic, and macroevolutionary variation in the Drosophilidae--a group that diverged from the sepsid lineage ca. 64 My ago. We also find that developmental bias in S. punctum wing shape aligns with the effects of allometry, but less so with putatively adaptive thermal plasticity and population differentiation along latitude. Our findings demonstrate that developmental bias in fly wings predicts evolvability and macroevolutionary trajectories on a much greater scale than previously appreciated but also suggest that causal explanations for such alignments may go beyond simple constraint hypotheses.

Skeleton phylogeny reconstructed with transcriptomes for the tribe Drosophilini (Diptera: Drosophilidae).

The family Drosophilidae is one of the most important model systems in evolutionary biology. Thanks to advances in high-throughput sequencing technology, a number of molecular phylogenetic analyses have been undertaken by using large data sets of many genes and many species sampled across this family. Especially, recent analyses using genome sequences have depicted the family-wide skeleton phylogeny with high confidence. However, the taxon sampling is still insufficient for minor lineages and non-Drosophila genera. In this study, we carried out phylogenetic analyses using a large number of transcriptome-based nucleotide sequences, focusing on the largest, core tribe Drosophilini in the Drosophilidae. In our analyses, some noise factors against phylogenetic reconstruction were taken into account by removing putative paralogy from the datasets and examining the effects of missing data, i.e. gene occupancy and site coverage, and incomplete lineage sorting. The inferred phylogeny has newly resolved the following phylogenetic positions/relationships at the genomic scale: (i) the monophyly of the subgenus Siphlodora including Zaprionus flavofasciatus to be transferred therein; (ii) the paraphyly of the robusta and melanica species groups within a clade comprised of the robusta, melanica and quadrisetata groups and Z. flavofasciatus; (iii) Drosophila curviceps (representing the curviceps group), D. annulipes (the quadrilineata subgroup of the immigrans group) and D. maculinotata clustered into a clade sister to the Idiomyia + Scaptomyza clade, forming together the expanded Hawaiian drosophilid lineage; (iv) Dichaetophora tenuicauda (representing the lineage comprised of the Zygothrica genus group and Dichaetophora) placed as the sister to the clade of the expanded Hawaiian drosophilid lineage and Siphlodora; and (v) relationships of the subgenus Drosophila and the genus Zaprionus as follows: (Zaprionus, (the quadrilineata subgroup, ((D. sternopleuralis, the immigrans group proper), (the quinaria radiation, the tripunctata radiation)))). These results are to be incorporated into the so-far published phylogenomic tree as a backbone (constraint) tree for grafting much more species based on sequences of a limited number of genes. Such a comprehensive, highly confident phylogenetic tree with extensive and dense taxon sampling will provide an essential framework for comparative studies of the Drosophilidae.

Competitive history shapes rapid evolution in a seasonal climate.

Eco-evolutionary dynamics will play a critical role in determining species' fates as climatic conditions change. Unfortunately, we have little understanding of how rapid evolutionary responses to climate play out when species are embedded in the competitive communities that they inhabit in nature. We tested the effects of rapid evolution in response to interspecific competition on subsequent ecological and evolutionary trajectories in a seasonally changing climate using a field-based evolution experiment with Drosophila melanogaster Populations of D. melanogaster were either exposed, or not exposed, to interspecific competition with an invasive competitor, Zaprionus indianus, over the summer. We then quantified these populations' ecological trajectories (abundances) and evolutionary trajectories (heritable phenotypic change) when exposed to a cooling fall climate. We found that competition with Z. indianus in the summer affected the subsequent evolutionary trajectory of D. melanogaster populations in the fall, after all interspecific competition had ceased. Specifically, flies with a history of interspecific competition evolved under fall conditions to be larger and have lower cold fecundity and faster development than flies without a history of interspecific competition. Surprisingly, this divergent fall evolutionary trajectory occurred in the absence of any detectible effect of the summer competitive environment on phenotypic evolution over the summer or population dynamics in the fall. This study demonstrates that competitive interactions can leave a legacy that shapes evolutionary responses to climate even after competition has ceased, and more broadly, that evolution in response to one selective pressure can fundamentally alter evolution in response to subsequent agents of selection.

Phylotranscriptomics Reveals Discordance in the Phylogeny of Hawaiian Drosophila and Scaptomyza (Diptera: Drosophilidae).

Island radiations present natural laboratories for studying the evolutionary process. The Hawaiian Drosophilidae are one such radiation, with nearly 600 described species and substantial morphological and ecological diversification. These species are largely divided into a few major clades, but the relationship between clades remains uncertain. Here, we present new assembled transcriptomes from 12 species across these clades, and use these transcriptomes to resolve the base of the evolutionary radiation. We recover a new hypothesis for the relationship between clades, and demonstrate its support over previously published hypotheses. We then use the evolutionary radiation to explore dynamics of concordance in phylogenetic support, by analyzing the gene and site concordance factors for every possible topological combination of major groups. We show that high bootstrap values mask low evolutionary concordance, and we demonstrate that the most likely topology is distinct from the topology with the highest support across gene trees and from the topology with highest support across sites. We then combine all previously published genetic data for the group to estimate a time-calibrated tree for over 300 species of drosophilids. Finally, we digitize dozens of published Hawaiian Drosophilidae descriptions, and use this to pinpoint probable evolutionary shifts in reproductive ecology as well as body, wing, and egg size. We show that by examining the entire landscape of tree and trait space, we can gain a more complete understanding of how evolutionary dynamics play out across an island radiation.

Evolution of Olfactory Receptors Tuned to Mustard Oils in Herbivorous Drosophilidae.

The diversity of herbivorous insects is attributed to their propensity to specialize on toxic plants. In an evolutionary twist, toxins betray the identity of their bearers when herbivores coopt them as cues for host-plant finding, but the evolutionary mechanisms underlying this phenomenon are poorly understood. We focused on Scaptomyza flava, an herbivorous drosophilid specialized on isothiocyanate (ITC)-producing (Brassicales) plants, and identified Or67b paralogs that were triplicated as mustard-specific herbivory evolved. Using in vivo heterologous systems for the expression of olfactory receptors, we found that S. flava Or67bs, but not the homologs from microbe-feeding relatives, responded selectively to ITCs, each paralog detecting different ITC subsets. Consistent with this, S. flava was attracted to ITCs, as was Drosophila melanogaster expressing S. flava Or67b3 in the homologous Or67b olfactory circuit. ITCs were likely coopted as olfactory attractants through gene duplication and functional specialization (neofunctionalization and subfunctionalization) in S. flava, a recently derived herbivore.

Temporal flexibility of gene regulatory network underlies a novel wing pattern in flies.

Organisms have evolved endless morphological, physiological, and behavioral novel traits during the course of evolution. Novel traits were proposed to evolve mainly by orchestration of preexisting genes. Over the past two decades, biologists have shown that cooption of gene regulatory networks (GRNs) indeed underlies numerous evolutionary novelties. However, very little is known about the actual GRN properties that allow such redeployment. Here we have investigated the generation and evolution of the complex wing pattern of the fly Samoaia leonensis We show that the transcription factor Engrailed is recruited independently from the other players of the anterior-posterior specification network to generate a new wing pattern. We argue that partial cooption is made possible because 1) the anterior-posterior specification GRN is flexible over time in the developing wing and 2) this flexibility results from the fact that every single gene of the GRN possesses its own functional time window. We propose that the temporal flexibility of a GRN is a general prerequisite for its possible cooption during the course of evolution.

A small survey of introduced Zaprionus indianus (Diptera: Drosophilidae) in orchards of the eastern United States.

The African fig fly, Zaprionus indianus (Gupta), is a generalist fruit fly that typically breeds in decaying fruits from over 70 plant species. The species has spread globally from its native range in tropical Africa, becoming an invasive pest on ripening figs in Brazil. First reported in the United States in 2005 in Florida, Z. indianus has since been documented as far north as Canada and is hypothesized to recolonize northwards from southern refugia each year. We sampled drosophilid communities over the growing season at 2 orchards in Virginia from 2020 to 2022 and 11 orchards along the East Coast during the fall of 2022 to quantify the abundance of Z. indianus relative to other drosophilids across locations, seasons, and fruit crops. Massachusetts had the northernmost population, with no Z. indianus detected in Maine and no correlation between latitude and relative abundance. Variation in Z. indianus relative abundance was high between nearby orchards and abundance was higher on peaches relative to apples within orchards. Comparisons of seasonal abundance curves between 2 Virginia orchards showed similar dynamics across years with individuals first detected around July and becoming absent around December, with peaks in late summer and mid-fall. The variation in seasonal and latitudinal abundance shown here highlights a need for broader sampling to accurately characterize the range, spread, and environmental tolerances of Z. indianus in North America.

Evolution and development of male-specific leg brushes in Drosophilidae.

The origin, diversification, and secondary loss of sexually dimorphic characters are common in animal evolution. In some cases, structurally and functionally similar traits have evolved independently in multiple lineages. Prominent examples of such traits include the male-specific grasping structures that develop on the front legs of many dipteran insects. In this report, we describe the evolution and development of one of these structures, the male-specific "sex brush." The sex brush is composed of densely packed, irregularly arranged modified bristles and is found in several distantly related lineages in the family Drosophilidae. Phylogenetic analysis using 250 genes from over 200 species provides modest support for a single origin of the sex brush followed by many secondary losses; however, independent origins of the sex brush cannot be ruled out completely. We show that sex brushes develop in very similar ways in all brush-bearing lineages. The dense packing of brush hairs is explained by the specification of bristle precursor cells at a near-maximum density permitted by the lateral inhibition mechanism, as well as by the reduced size of the surrounding epithelial cells. In contrast to the female and the ancestral male condition, where bristles are arranged in stereotypical, precisely spaced rows, cell migration does not contribute appreciably to the formation of the sex brush. The complex phylogenetic history of the sex brush can make it a valuable model for investigating coevolution of sex-specific morphology and mating behavior.

Evolution and genomic basis of the plant-penetrating ovipositor: a key morphological trait in herbivorous Drosophilidae.

Herbivorous insects are extraordinarily diverse, yet are found in only one-third of insect orders. This skew may result from barriers to plant colonization, coupled with phylogenetic constraint on plant-colonizing adaptations. The plant-penetrating ovipositor, however, is one trait that surmounts host plant physical defences and may be evolutionarily labile. Ovipositors densely lined with hard bristles have evolved repeatedly in herbivorous lineages, including within the Drosophilidae. However, the evolution and genetic basis of this innovation has not been well studied. Here, we focused on the evolution of this trait in Scaptomyza, a genus sister to Hawaiian Drosophila, that contains a herbivorous clade. Our phylogenetic approach revealed that ovipositor bristle number increased as herbivory evolved in the Scaptomyza lineage. Through a genome-wide association study, we then dissected the genomic architecture of variation in ovipositor bristle number within S. flava. Top-associated variants were enriched for transcriptional repressors, and the strongest associations included genes contributing to peripheral nervous system development. Individual genotyping supported the association at a variant upstream of Gαi, a neural development gene, contributing to a gain of 0.58 bristles/major allele. These results suggest that regulatory variation involving conserved developmental genes contributes to this key morphological trait involved in plant colonization.

Insect cross-tolerance to freezing and drought stress: role of metabolic rearrangement.

The accumulation of trehalose has been suggested as a mechanism underlying insect cross-tolerance to cold/freezing and drought. Here we show that exposing diapausing larvae of the drosophilid fly, Chymomyza costata to dry conditions significantly stimulates their freeze tolerance. It does not, however, improve their tolerance to desiccation, nor does it significantly affect trehalose concentrations. Next, we use metabolomics to compare the complex alterations to intermediary metabolism pathways in response to three environmental factors with different ecological meanings: environmental drought (an environmental stressor causing mortality), decreasing ambient temperatures (an acclimation stimulus for improvement of cold hardiness), and short days (an environmental signal inducing diapause). We show that all three factors trigger qualitatively similar metabolic rearrangement and a similar phenotypic outcome-improved larval freeze tolerance. The similarities in metabolic response include (but are not restricted to) the accumulation of typical compatible solutes and the accumulation of energy-rich molecules (phosphagens). Based on these results, we suggest that transition to metabolic suppression (a state in which chemical energy demand is relatively low but need for stabilization of macromolecules is high) represents a common axis of metabolic pathway reorganization towards accumulation of non-toxic cytoprotective compounds, which in turn stimulates larval freeze tolerance.

Genetic Diversity and Thermal Performance in Invasive and Native Populations of African Fig Flies.

During biological invasions, invasive populations can suffer losses of genetic diversity that are predicted to negatively impact their fitness/performance. Despite examples of invasive populations harboring lower diversity than conspecific populations in their native range, few studies have linked this lower diversity to a decrease in fitness. Using genome sequences, we show that invasive populations of the African fig fly, Zaprionus indianus, have less genetic diversity than conspecific populations in their native range and that diversity is proportionally lower in regions of the genome experiencing low recombination rates. This result suggests that selection may have played a role in lowering diversity in the invasive populations. We next use interspecific comparisons to show that genetic diversity remains relatively high in invasive populations of Z. indianus when compared with other closely related species. By comparing genetic diversity in orthologous gene regions, we also show that the genome-wide landscape of genetic diversity differs between invasive and native populations of Z. indianus indicating that invasion not only affects amounts of genetic diversity but also how that diversity is distributed across the genome. Finally, we use parameter estimates from thermal performance curves for 13 species of Zaprionus to show that Z. indianus has the broadest thermal niche of measured species, and that performance does not differ between invasive and native populations. These results illustrate how aspects of genetic diversity in invasive species can be decoupled from measures of fitness, and that a broad thermal niche may have helped facilitate Z. indianus's range expansion.

Genomic signatures of admixture and selection are shared among populations of Zaprionus indianus across the western hemisphere.

Introduced species have become an increasingly common component of biological communities around the world. A central goal in invasion biology is therefore to identify the demographic and evolutionary factors that underlie successful introductions. Here we use whole genome sequences, collected from populations in the native and introduced range of the African fig fly, Zaprionus indianus, to quantify genetic relationships among them, identify potential sources of the introductions, and test for selection at different spatial scales. We find that geographically widespread populations in the western hemisphere are genetically more similar to each other than to lineages sampled across Africa, and that these populations share a mixture of alleles derived from differentiated African lineages. Using patterns of allele-sharing and demographic modelling we show that Z. indinaus have undergone a single expansion across the western hemisphere with admixture between African lineages predating this expansion. We also find support for selection that is shared across populations in the western hemisphere, and in some cases, with a subset of African populations. This suggests either that parallel selection has acted across a large part of Z. indianus's introduced range; or, more parsimoniously, that Z. indianus has experienced selection early on during (or prior-to) its expansion into the western hemisphere. We suggest that the range expansion of Z. indianus has been facilitated by admixture and selection, and that management of this invasion could focus on minimizing future admixture by controlling the movement of individuals within this region rather than between the western and eastern hemisphere.

Comparative mitogenomics of Drosophilidae and the evolution of the Zygothrica genus group (Diptera, Drosophilidae).

The Zygothrica genus group of Drosophilidae encompasses more than 437 species and five genera. Although knowledge regarding its diversity has increased, uncertainties about its monophyly and position within Drosophilidae remain. Genomic approaches have been widely used to address different phylogenetic questions and analyses involving the mitogenome have revealed a cost-efficient tool to these studies. Thus, this work aims to characterize mitogenomes of three species of the Zygothrica genus group (from the Hirtodrosophila, Paraliodrosophila and Zygothrica genera), while comparing them with orthologous sequences from other 23 Drosophilidae species and addressing their phylogenetic position. General content concerning gene order and overlap, nucleotide composition, start and stop codon, codon usage and tRNA structures were compared, and phylogenetic trees were constructed under different datasets. The complete mitogenomes characterized for H. subflavohalterata affinis H002 and P. antennta present the PanCrustacea gene order with 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, 13 protein coding genes and an A+T rich region with two T-stretched elements. Some peculiarities such as the almost complete overlap of genes tRNAH/ND4, tRNAF/ND5 and tRNAS2/ND1 are reported for different Drosophilidae species. Non-canonical secondary structures were encountered for tRNAS1 and tRNAY, revealing patterns that apply at different phylogenetic scales. According to the best depiction of the mitogenomes evolutionary history, the three Neotropical species of the Zygothrica genus group encompass a monophyletic lineage sister to Zaprionus, composing with this genus a clade that is sister to the Drosophila subgenus.

Viral susceptibility across host species is largely independent of dietary protein to carbohydrate ratios.

The likelihood of a successful host shift of a parasite to a novel host species can be influenced by environmental factors that can act on both the host and parasite. Changes in nutritional resource availability have been shown to alter pathogen susceptibility and the outcome of infection in a range of systems. Here, we examined how dietary protein to carbohydrate altered susceptibility in a large cross-infection experiment. We infected 27 species of Drosophilidae with an RNA virus on three food types of differing protein to carbohydrate ratios. We then measured how viral load and mortality across species was affected by changes in diet. We found that changes in the protein:carbohydrate in the diet did not alter the outcomes of infection, with strong positive inter-species correlations in both viral load and mortality across diets, suggesting no species-by-diet interaction. Mortality and viral load were strongly positively correlated, and this association was consistent across diets. This suggests changes in diet may give consistent outcomes across host species, and may not be universally important in determining host susceptibility to pathogens.

Karyotype Description of Two Andean Species of the guarani Group of Drosophila (Díptera: Drosophilidae) and Cytological Notes.

The guarani group of Drosophila genus (Diptera: Drosophilidae) is formed by 24 species however the relationship of these species is not clear. In the present study are described the karyotypes of Drosophila sachapuyu Peñafiel and Rafael, 2018 and Drosophila zamorana Peñafiel and Rafael, 2018, two Andean species members of the guarani group. Mitotic chromosomes from cerebral ganglia of third stand larval were obtained by thermal shock and cell suspension techniques. The karyotype of D. sachapuyu, presents 2n = 10 (4R, 1V; X = R, Y = R) while D. zamorana exhibits karyotype 2n = 12 (5R, 1V; X = V, Y = R).

A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly.

BACKGROUND: The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics. RESULTS: Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS). CONCLUSIONS: Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

Insect mitochondria as targets of freezing-induced injury.

Many insects survive internal freezing, but the great complexity of freezing stress hinders progress in understanding the ultimate nature of freezing-induced injury. Here, we use larvae of the drosophilid fly, Chymomyza costata to assess the role of mitochondrial responses to freezing stress. Respiration analysis revealed that fat body mitochondria of the freeze-sensitive (non-diapause) phenotype significantly decrease oxygen consumption upon lethal freezing stress, while mitochondria of the freeze-tolerant (diapausing, cold-acclimated) phenotype do not lose respiratory capacity upon the same stress. Using transmission electron microscopy, we show that fat body and hindgut mitochondria swell, and occasionally burst, upon exposure of the freeze-sensitive phenotype to lethal freezing stress. By contrast, mitochondrial swelling is not observed in the freeze-tolerant phenotype exposed to the same stress. We hypothesize that mitochondrial swelling results from permeability transition of the inner mitochondrial membrane and loss of its barrier function, which causes osmotic influx of cytosolic water into the matrix. We therefore suggest that the phenotypic transition to diapause and cold acclimation could be associated with adaptive changes that include the protection of the inner mitochondrial membrane against permeability transition and subsequent mitochondrial swelling. Accumulation of high concentrations of proline and other cryoprotective substances might be a part of such adaptive changes as we have shown that freezing-induced mitochondrial swelling was abolished by feeding the freeze-sensitive phenotype larvae on a proline-augmented diet.

Changes in temperature alter the potential outcomes of virus host shifts.

Host shifts-where a pathogen jumps between different host species-are an important source of emerging infectious disease. With on-going climate change there is an increasing need to understand the effect changes in temperature may have on emerging infectious disease. We investigated whether species' susceptibilities change with temperature and ask if susceptibility is greatest at different temperatures in different species. We infected 45 species of Drosophilidae with an RNA virus and measured how viral load changes with temperature. We found the host phylogeny explained a large proportion of the variation in viral load at each temperature, with strong phylogenetic correlations between viral loads across temperature. The variance in viral load increased with temperature, while the mean viral load did not. This suggests that as temperature increases the most susceptible species become more susceptible, and the least susceptible less so. We found no significant relationship between a species' susceptibility across temperatures, and proxies for thermal optima (critical thermal maximum and minimum or basal metabolic rate). These results suggest that whilst the rank order of species susceptibilities may remain the same with changes in temperature, some species may become more susceptible to a novel pathogen, and others less so.

Host shifts result in parallel genetic changes when viruses evolve in closely related species.

Host shifts, where a pathogen invades and establishes in a new host species, are a major source of emerging infectious diseases. They frequently occur between related host species and often rely on the pathogen evolving adaptations that increase their fitness in the novel host species. To investigate genetic changes in novel hosts, we experimentally evolved replicate lineages of an RNA virus (Drosophila C Virus) in 19 different species of Drosophilidae and deep sequenced the viral genomes. We found a strong pattern of parallel evolution, where viral lineages from the same host were genetically more similar to each other than to lineages from other host species. When we compared viruses that had evolved in different host species, we found that parallel genetic changes were more likely to occur if the two host species were closely related. This suggests that when a virus adapts to one host it might also become better adapted to closely related host species. This may explain in part why host shifts tend to occur between related species, and may mean that when a new pathogen appears in a given species, closely related species may become vulnerable to the new disease.

Phylogenetic relationships between fungus-associated Neotropical species of the genera Hirtodrosophila, Mycodrosophila and Zygothrica (Diptera, Drosophilidae), with insights into the evolution of breeding sites usage.

The Neotropical region harbors an astonishing diversity of species, but still encompasses the least studied biogeographic region of the world. These properties apply for different taxonomic groups, and can be exemplified by drosophilids. In fact, high levels of cryptic diversity have recently been discovered for Neotropical species of the Zygothrica genus group, but relationships among these species, or them and other Drosophilidae species still remains to be addressed. Therefore, the aim of this study was to evaluate the phylogenetic relationships between fungus-associated Neotropical species of the genera Hirtodrosophila, Mycodrosophila and Zygothrica, which together with Paramycodrosophila and Paraliodrosophila compose the Zygothrica genus group. For this, fragments of the mitochondrial cytochrome oxidase subunits I (COI) and II (COII) genes, and the nuclear alpha methyldopa (Amd) and dopa decarboxylase (Ddc) genes were newly characterized for 43 Neotropical specimens of fungus-associated drosophilids, and analyzed in the context of 51 additional Drosophilinae sequences plus one Steganinae outgroup. Based on the resulting phylogeny, the evolution of breeding sites usage was also evaluated through ancestral character reconstructions. Our results revealed the Zygothrica genus group as a monophyletic lineage of Drosophila that branches after the subgenera Sophophora and Drosophila. Within this lineage, Mycodrosophila species seem to encompass the early offshoot, followed by a grade of Hirtodrosophila species, with derived branches mostly occupied by representatives of Zygothrica. This genus, in particular, was subdivided into five major clades, two of which include species of Hirtodrosophila, whose generic status needs to be reevatuated. According to our results, the use of fungi as breeding sites encompasses a symplesiomorphy for the Zygothrica genus group, since one of the recovered clades is currently specialized in using flowers as breeding sites whereas a sole species presents a reversal to the use of fruits of a plant of Gentianales. So, in general, this study supports the paraphyly of Drosophila in relation to fungus-associated Neotropical species of Drosophilidae, providing the first molecular insights into the phylogenetic patterns related to the evolution of this diverse group of species and some of its characteristic traits.