Phylogenetic analysis of the repleta species group of the genus Drosophila using multiple sources of characters

Dietary Utilization Drives the Differentiation of Gut Bacterial Communities between Specialist and Generalist Drosophilid Flies

Gut bacteria play vital roles in the dietary detoxification, digestion, and nutrient supplementation of hosts during dietary specialization. The roles of gut bacteria in the host can be unveiled by comparing communities of specialist and generalist bacterial species. However, these species usually have a long evolutionary history, making it difficult to determine whether bacterial community differentiation is due to host dietary adaptation or phylogenetic divergence. In this regard, we investigated the bacterial communities from two Araceae-feeding Colocasiomyia species and further performed a meta-analysis by incorporating the published data from Drosophila bacterial community studies. The compositional and functional differentiation of bacterial communities was uncovered by comparing three (Araceae-feeding, mycophagous, and cactophilic) specialists with generalist flies. The compositional differentiation showed that Bacteroidetes and Firmicutes inhabited specialists, while more Proteobacteria lived in generalists. The functional prediction based on the bacterial community compositions suggested that amino acid metabolism and energy metabolism are overrepresented pathways in specialists and generalists, respectively. The differences were mainly associated with the higher utilization of structural complex carbohydrates, protein utilization, vitamin B₁₂ acquisition, and demand for detoxification in specialists than in generalists. The complementary roles of bacteria reveal a connection between gut bacterial communities and fly dietary specialization.

IMPORTANCE Gut bacteria may play roles in the dietary utilization of hosts, especially in specialist animals, during long-term host-microbe interaction. By comparing the gut bacterial communities between specialist and generalist drosophilid flies, we found that specialists harbor more bacteria linked to complex carbohydrate degradation, amino acid metabolism, vitamin B₁₂ formation, and detoxification than do generalists. This study reveals the roles of gut bacteria in drosophilid species in dietary utilization.

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.

Evolution of male genitalia in the Drosophila repleta species group (Diptera: Drosophilidae)

The Drosophila repleta group comprises more than one hundred species that inhabit several environments in the Neotropics and use different hosts as rearing and feeding resources. Rather homogeneous in their external morphology, they are generally distinguished by the male genitalia, seemingly their fastest evolving morphological trait, constituting an excellent model to study patterns of genital evolution in the context of a continental adaptive radiation. Although much is known about the evolution of animal genitalia at population level, surveys on macroevolutionary scale of this phenomenon are scarce. This study used a suite of phylogenetic comparative methods to elucidate the macroevolutionary patterns of genital evolution through deep time and large continental scales. Our results indicate that male genital size and some aspects of shape have been evolving by speciational evolution, probably due to the microevolutionary processes involved in species mate recognition. In contrast, several features of the aedeagus shape seemed to have evolved in a gradual fashion, with heterogeneous evolutionary phenotypic rates among clades. In general, the tempo of the evolution of aedeagus morphology was constant from the origin of the group until the Pliocene, when it accelerated in some clades that diversified mainly in this period. The incidence of novel ecological conditions in the tempo of aedeagus evolution and the relationship between species mate recognition and speciation in the Drosophila repleta group are discussed.

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.

Genomic changes associated with adaptation to arid environments in cactophilic Drosophila species

BACKGROUND

Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. Here we investigate the genetic changes associated with adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species.

RESULTS

We find disproportionately high rates of gene gains in internal branches in the species' phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species' differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one.

CONCLUSION

Adaptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change.

An Improved Genome Assembly for Drosophila navojoa, the Basal Species in the mojavensis Cluster

Three North American cactophilic Drosophila species, D. mojavensis, D. arizonae, and D. navojoa, are of considerable evolutionary interest owing to the shift from breeding in Opuntia cacti to columnar species. The 3 species form the "mojavensis cluster" of Drosophila. The genome of D. mojavensis was sequenced in 2007 and the genomes of D. navojoa and D. arizonae were sequenced together in 2016 using the same technology (Illumina) and assembly software (AllPaths-LG). Yet, unfortunately, the D. navojoa genome was considerably more fragmented and incomplete than its sister species, rendering it less useful for evolutionary genetic studies. The D. navojoa read dataset does not fully meet the strict insert size required by the assembler used (AllPaths-LG) and this incompatibility might explain its assembly problems. Accordingly, when we re-assembled the genome of D. navojoa with the SPAdes assembler, which does not have the strict AllPaths-LG requirements, we obtained a substantial improvement in all quality indicators such as N50 (from 84 kb to 389 kb) and BUSCO coverage (from 77% to 97%). Here we share a new, improved reference assembly for D. navojoa genome, along with a RNAseq transcriptome. Given the basal relationship of the Opuntia breeding D. navojoa to the columnar breeding D. arizonae and D. mojavensis, the improved assembly and annotation will allow researchers to address a range of questions associated with the genomics of host shifts, chromosomal rearrangements and speciation in this group.

Pupal emergence pattern in cactophilic Drosophila and the effect of host plants

Drosophila buzzatii and D. koepferae are sibling cactophilic species. The former breeds primarily on prickly pears (genus Opuntia) whereas the latter breeds on columnar cacti of the genera Cereus and Trichocereus, although with certain degree of niche overlapping. We examined the interspecific differences in diurnal temporal patterns of adult emergence from puparia and evaluated whether this behavior is affected by rearing in the different cactus hosts available in nature. We detected important host-dependent genetic variation for this trait differentially affecting the emergence schedule of these species. Diurnal pattern of emergence time was directly correlated with developmental time and negatively correlated with adult wing size, suggesting that early emergences are at least indirectly correlated with increased fitness. We discussed our results in terms of their putative effects on fitness and the genetic-metabolic pathways that would be presumably affected by host's nutritional-chemical differences.

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.

Genome-wide tests for introgression between cactophilic Drosophila implicate a role of inversions during speciation

Models of speciation‐with‐gene‐flow have shown that the reduction in recombination between alternative chromosome arrangements can facilitate the fixation of locally adaptive genes in the face of gene flow and contribute to speciation. However, it has proven frustratingly difficult to show empirically that inversions have reduced gene flow and arose during or shortly after the onset of species divergence rather than represent ancestral polymorphisms. Here, we present an analysis of whole genome data from a pair of cactophilic fruit flies, Drosophila mojavensis and D. arizonae, which are reproductively isolated in the wild and differ by several large inversions on three chromosomes. We found an increase in divergence at rearranged compared to colinear chromosomes. Using the density of divergent sites in short sequence blocks we fit a series of explicit models of species divergence in which gene flow is restricted to an initial period after divergence and may differ between colinear and rearranged parts of the genome. These analyses show that D. mojavensis and D. arizonae have experienced postdivergence gene flow that ceased around 270 KY ago and was significantly reduced in chromosomes with fixed inversions. Moreover, we show that these inversions most likely originated around the time of species divergence which is compatible with theoretical models that posit a role of inversions in speciation with gene flow.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.

Variations on a theme: diversification of cuticular hydrocarbons in a clade of cactophilic Drosophila

Background

We characterized variation and chemical composition of epicuticular hydrocarbons (CHCs) in the seven species of the Drosophila buzzatii cluster with gas chromatography/mass spectrometry. Despite the critical role of CHCs in providing resistance to desiccation and involvement in communication, such as courtship behavior, mating, and aggregation, few studies have investigated how CHC profiles evolve within and between species in a phylogenetic context. We analyzed quantitative differences in CHC profiles in populations of the D. buzzatii species cluster in order to assess the concordance of CHC differentiation with species divergence.

Results

Thirty-six CHC components were scored in single fly extracts with carbon chain lengths ranging from C₂₉ to C₃₉, including methyl-branched alkanes, n-alkenes, and alkadienes. Multivariate analysis of variance revealed that CHC amounts were significantly different among all species and canonical discriminant function (CDF) analysis resolved all species into distinct, non-overlapping groups. Significant intraspecific variation was found in different populations of D. serido suggesting that this taxon is comprised of at least two species. We summarized CHC variation using CDF analysis and mapped the first five CHC canonical variates (CVs) onto an independently derived period (per) gene + chromosome inversion + mtDNA COI gene for each sex. We found that the COI sequences were not phylogenetically informative due to introgression between some species, so only per + inversion data were used. Positive phylogenetic signal was observed mainly for CV1 when parsimony methods and the test for serial independence (TFSI) were used. These results changed when no outgroup species were included in the analysis and phylogenetic signal was then observed for female CV3 and/or CV4 and male CV4 and CV5. Finally, removal of divergent populations of D. serido significantly increased the amount of phylogenetic signal as up to four out of five CVs then displayed positive phylogenetic signal.

Conclusions

CHCs were conserved among species while quantitative differences in CHC profiles between populations and species were statistically significant. Most CHCs were species-, population-, and sex-specific. Mapping CHCs onto an independently derived phylogeny revealed that a significant portion of CHC variation was explained by species' systematic affinities indicating phylogenetic conservatism in the evolution of these hydrocarbon arrays, presumptive waterproofing compounds and courtship signals as in many other drosophilid species.

On the phylogeny of Mustelidae subfamilies: analysis of seventeen nuclear non-coding loci and mitochondrial complete genomes

BACKGROUND

Mustelidae, as the largest and most-diverse family of order Carnivora, comprises eight subfamilies. Phylogenetic relationships among these Mustelidae subfamilies remain argumentative subjects in recent years. One of the main reasons is that the mustelids represent a typical example of rapid evolutionary radiation and recent speciation event. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence and nuclear protein-coding data, herein we employ 17 nuclear non-coding loci (>15 kb), in conjunction with mt complete genome data (>16 kb), to clarify these enigmatic problems.

RESULTS

The combined nuclear intron and mt genome analyses both robustly support that Taxidiinae diverged first, followed by Melinae. Lutrinae and Mustelinae are grouped together in all analyses with strong supports. The position of Helictidinae, however, is enigmatic because the mt genome analysis places it to the clade uniting Lutrinae and Mustelinae, whereas the nuclear intron analysis favors a novel view supporting a closer relationship of Helictidinae to Martinae. This finding emphasizes a need to add more data and include more taxa to resolve this problem. In addition, the molecular dating provides insights into the time scale of the origin and diversification of the Mustelidae subfamilies. Finally, the phylogenetic performances and limits of nuclear introns and mt genes are discussed in the context of Mustelidae phylogeny.

CONCLUSION

Our study not only brings new perspectives on the previously obscured phylogenetic relationships among Mustelidae subfamilies, but also provides another example demonstrating the effectiveness of nuclear non-coding loci for reconstructing evolutionary histories in a group that has undergone rapid bursts of speciation.

Coevolution of male and female reproductive structures in Drosophila

The morphology of male genitalia whilst stable within species, exhibits huge interspecific variation. This variation is likely to be as a result of sexual selection due to the direct involvement of these reproductive structures in mating and sperm transfer. In contrast, internal soft tissue components of the genitalia are generally poorly investigated as they are not directly involved in physical and mechanical adequacy during sperm transfer. However, these soft tissue structures may also drive differential male-female interactions, particularly in internally fertilising organisms where females have the ability to store sperm and bias male reproductive success. In this paper we use the drosophila model to investigate the role of male and female reproductive elements in sexual selection. Our meta-analysis supplemented with additional new data clearly shows that within species, sperm length versus testis length, and sperm length versus seminal receptacle length, are highly correlated. Thus, independent of the phylogenetic relationship among species, gamete evolution is likely to result in sexual selection interactions that drive the evolution of internal reproductive components in both sexes. Our results and discussion of the literature highlight the importance of considering internal soft structures that may influence fertilisation, when investigating selective forces acting on the evolution of reproductive traits.

Cactophilic Drosophila in South America: a model for evolutionary studies

The Drosophila buzzatii cluster is composed of seven cactophilic species and their known geographical distribution encompasses the open vegetation diagonal, which includes the morphoclimatic Domains of the Caatinga, Chaco and Cerrado, which are situated between the Amazon and the Atlantic forests. Besides these areas, these cactophilic species are also found in a narrow strip along the Atlantic coast from northeastern Brazil to the southern tip of the country. The hypothesis of vicariant events, defining the core areas of each species, is proposed to explain the historical diversification for the cluster. The intraspecific analysis for the cluster shows a population structure with gene flow restricted by distance, range expansion with secondary contact resulting in introgression and simpatry, especially in the limits of the species distribution, polytypic populations and assortative mating in inter population experiments. There is a variation related to these events that depends on the species and geographic origin of the population analyzed. These events are, hypothetically, described as the results of expansion and retraction of the population ranges, as a consequence of their association with cacti, which theoretically follow the expansion and retraction of dry areas during the paleoclimatic oscillations in South America, as that promoted by the glacial cycles of the Quaternary. The Drosophila buzzatii cluster is divided into two groups. The first one is composed of D. buzzatii, a species that has a broad geographic distribution and no significant differentiation between its populations. The second is the Drosophila serido sibling set, which encompasses the others species and is characterized by a significant potential for differentiation.

Orthology, function and evolution of accessory gland proteins in the Drosophila repleta group

The accessory gland proteins (Acps) of Drosophila have become a model for the study of reproductive protein evolution. A major step in the study of Acps is to identify biological causes and consequences of the observed patterns of molecular evolution by comparing species groups with different biology. Here we characterize the Acp complement of Drosophila mayaguana, a repleta group representative. Species of this group show important differences in ecology and reproduction as compared to other Drosophila. Our results show that the extremely high rates of Acp evolution previously found are likely to be ubiquitous among species of the repleta group. These evolutionary rates are considerably higher than the ones observed in other Drosophila groups' Acps. This disparity, however, is not accompanied by major differences in the estimated number of Acps or in the functional categories represented as previously suggested. Among the genes expressed in accessory glands of D. mayaguana almost half are likely products of recent duplications. This allowed us to test predictions of the neofunctionalization model for gene duplication and paralog evolution in a more or less constrained timescale. We found that positive selection is a strong force in the early divergence of these gene pairs.

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.

Evidence of adaptive evolution of accessory gland proteins in closely related species of the Drosophila repleta group

Accessory gland proteins (Acps) are part of the seminal fluid of Drosophila species. These proteins have important reproductive functions, being responsible for the proper functioning of several steps of the fertilization process. Acps also contribute indirectly for the reproductive success of males by modulating female behavior. Evidence that Acps participate in sperm competition and sexual conflict includes findings that, on average, Acps have fast evolutionary rates, suggestive of adaptive evolution. This is especially true in species of the Drosophila repleta group. Nevertheless, only in a few occasions have robust statistical tests been used to determine whether observed evolutionary rates are in fact due to positive selection on amino acid substitutions between related species. Here we apply maximum likelihood tests for positive selection on 14 Acps of the D. repleta group. To increase statistical robustness, we use at least 8 sequences, all belonging to species of the Drosophila mulleri complex, for each gene analyzed. We found significant evidence of adaptive evolution for 10 of the tested genes. Among these, the ones with a conserved protein domain had positively selected sites within the functional region of the sequence. We also detected one instance of lineage-specific adaptive evolution in a clade formed by 2 sister species.

Molecular evolution and population genetics of two Drosophila mettleri cytochrome P450 genes involved in host plant utilization

Understanding the genetic basis of adaptation is one of the primary goals of evolutionary biology. The evolution of xenobiotic resistance in insects has proven to be an especially suitable arena for studying the genetics of adaptation, and resistant phenotypes are known to result from both coding and regulatory changes. In this study, we examine the evolutionary history and population genetics of two Drosophila mettleri cytochrome P450 genes that are putatively involved in the detoxification of alkaloids present in two of its cactus hosts: saguaro (Carnegiea gigantea) and senita (Lophocereus schottii). Previous studies demonstrated that Cyp28A1 was highly up-regulated following exposure to rotting senita tissue while Cyp4D10 was highly up-regulated following exposure to rotting saguaro tissue. Here, we show that a subset of sites in Cyp28A1 experienced adaptive evolution specifically in the D. mettleri lineage. Moreover, neutrality tests in several populations were also consistent with a history of selection on Cyp28A1. In contrast, we did not find evidence for positive selection on Cyp4D10, although this certainly does not preclude its involvement in host plant use. A surprising result that emerged from our population genetic analyses was the presence of significant genetic differentiation between flies collected from different host plant species (saguaro and senita) at Organ Pipe National Monument, Arizona, USA. This preliminary evidence suggests that D. mettleri may have evolved into distinctive host races that specialize on different hosts, a possibility that warrants further investigation.

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.

Grafting the molecular phylogenetic tree with morphological branches to reconstruct the evolutionary history of the genus Zaprionus (Diptera: Drosophilidae)

A molecular phylogeny for the drosophilid genus Zaprionus was inferred using a mitochondrial (CO-II) and a nuclear (Amyrel) gene using 22 available species. The combined molecular tree does not support the current classification, dubbed phylogenetic, based entirely upon a morphocline of forefemoral ornamentation. For species for which DNA was not available, phylogenetic positioning was only assigned using morphological characters. In order to avoid conflict between DNA and morphology in the combined analyses (supermatrix method), we developed a new method in which few morphological characters were sampled according to an a priori homoplasy assessment on the consensus molecular tree. At each internal node of the tree, a number of synapomorphies was determined, and species with no molecular sequences were grafted thereon. Analogously to tree vocabulary, we called our method 'morphological grafting'. New species groups and complexes were then defined in the light of our findings. Further, divergence times were estimated under a relaxed molecular clock, and historical biogeography was reconstructed under a maximum likelihood model. Zaprionus appears to be of recent origin in the Oriental region during the Late Miocene ( approximately 10 MYA), and colonization of Africa started shortly after ( approximately 7 MYA) via the maritime route of the Indian Ocean Islands. Most of the morphological and ecological diversification took place, later, in Western Africa during the Quaternary cyclic climatic changes. Furthermore, some species became recent invaders, with one, Zaprionus indianus, has successfully invaded South and North America during the last decade.

Evolutionary relationships of Drosophila mojavensis geographic host races and their sister species Drosophila arizonae

The cactophilic Drosophila mojavensis species group living in the deserts and dry tropical forests of the southwestern United States and Mexico provides a valuable system for studies in diversification and speciation. Rigorous studies of the relationships between host races of D. mojavensis and the relationships among the members of the species group (D. mojavensis, Drosophila arizona, and Drosophila navojoa) are lacking. We used mitochondrial CO1 sequence data to address the phylogenetics and population genetics of this species group. In this study we have found that the sister species D. mojavensis and D. arizonae share no mitochondrial haplotypes and thus show no evidence for recent introgression. We estimate the divergence time between D. mojavensis and D. arizonae to be between 1.91 and 2.97 million years ago. D. arizonae shows little structure in our population genetic analyses but there is phylogenetic differentiation between southeastern and northern populations of D. arizonae. Drosophila mojavensis shows significant population and phylogenetic structure across the four geographic regions of its distribution. The mitochondrial data support an origin of D. mojavensis on the mainland with early differentiation into the populations now found in the Mojave Desert and the Mainland Sonoran Desert and later colonization of the Baja Peninsula, in contrast to previous models. Also, the sister clade to D. mojavensis/D. arizonae includes D. navojoa and Drosophila huaylasi. By defining the genetic relationships among these populations, we provide a foundation for more sophisticated hypothesis testing regarding the timing of early speciation events and host switches in this species group.

Localization and characterization of X chromosome inversion breakpoints separating Drosophila mojavensis and Drosophila arizonae

Ectopic exchange between transposable elements or other repetitive sequences along a chromosome can produce chromosomal inversions. As a result, genome sequence studies typically find sequence similarity between corresponding inversion breakpoint regions. Here, we identify and investigate the breakpoint regions of the X chromosome inversion distinguishing Drosophila mojavensis and Drosophila arizonae. We localize one inversion breakpoint to 13.7 kb and localize the other to a 1-Mb interval. Using this localization and assuming microsynteny between Drosophila melanogaster and D. arizonae, we pinpoint likely positions of the inversion breakpoints to windows of less than 3000 bp. These breakpoints define the size of the inversion to approximately 11 Mb. However, in contrast to many other studies, we fail to find significant sequence similarity between the 2 breakpoint regions. The localization of these inversion breakpoints will facilitate future genetic and molecular evolutionary studies in this species group, an emerging model system for ecological genetics.

Molecular population genetics of the alpha-esterase5 gene locus in original and colonized populations of Drosophila buzzatii and its sibling Drosophila koepferae

Several studies have suggested that esterase-2 (EST-2) may be the target of natural selection in the cactophilic fly Drosophila buzzatii. In this work, we analyzed nucleotide variation in a fragment of alpha-esterase5 (alphaE5), the gene encoding EST-2, in original (Argentinian) and colonized (Australian) populations of D. buzzatii and in its sibling D. koepferae. Estimates of nucleotide heterozygosity in D. buzzatii were similar in Australia and Argentina, although we detected a loss of singletons in colonized populations, suggesting a moderate founder effect. Interspecific comparisons revealed that D. buzzatii was more polymorphic for nonsynonymous variation, whereas D. koepferae was more variable for synonymous and noncoding sites. The two major chromosomal arrangements (2st and 2j) in D. buzzatii displayed similar levels of nucleotide variation, whereas 2jz3 was monomorphic. The sequenced region allowed the discrimination of a greater number of EST-2 protein variants in the Australian sample than in the Argentinean sample. In D. koepferae, nucleotide variation in alphaE5 does not depart from neutral expectations, although tests of population structure were significant for silent variation. In contrast, D. buzzatii has probably undergone a recent population expansion in its South American range. In addition, the McDonald and Kreitman test revealed an excess of nonsynonymous polymorphism in both original and colonized populations of this species.

Testing chromosomal phylogenies and inversion breakpoint reuse in Drosophila

A combination of cytogenetic and bioinformatic procedures was used to test the chromosomal phylogeny relating Drosophila buzzatii with D. repleta. Chromosomes X and 2, harboring most of the inversions fixed between these two species, were analyzed. First, chromosomal segments conserved during the divergence of the two species were identified by comparative in situ hybridization to the D. repleta chromosomes of 180 BAC clones from a BAC-based physical map of the D. buzzatii genome. These conserved segments were precisely delimited with the aid of clones containing inversion breakpoints. Then GRIMM software was used to estimate the minimum number of rearrangements necessary to transform one genome into the other and identify all possible rearrangement scenarios. Finally, the most plausible inversion trajectory was tested by hybridizing 12 breakpoint-bearing BAC clones to the chromosomes of seven other species in the repleta group. The results show that chromosomes X and 2 of D. buzzatii and D. repleta differ by 12 paracentric inversions. Nine of them are fixed in chromosome 2 and entail two breakpoint reuses. Our results also show that the cytological relationship between D. repleta and D. mercatorum is closer than that between D. repleta and D. peninsularis, and we propose that the phylogenetic relationships in this lineage of the repleta group be reconsidered. We also estimated the rate of rearrangement between D. repleta and D. buzzatii and conclude that rates within the genus Drosophila vary substantially between lineages, even within a single species group.

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.

The mode of evolution of aggregation pheromones in Drosophila species

Aggregation pheromones are used by fruit flies of the genus Drosophila to assemble on breeding substrates, where they feed, mate and oviposit communally. These pheromones consist of species-specific blends of chemicals. Here, using a phylogenetic framework, we examine how differences among species in these pheromone blends have evolved. Theoretical predictions, genetic evidence, and previous empirical analysis of bark beetle species, suggest that aggregation pheromones do not evolve gradually, but via major, saltational shifts in chemical composition. Using pheromone data for 28 species of Drosophila we show that, unlike with bark beetles, the distribution of chemical components among species is highly congruent with their phylogeny, with closely related species being more similar in their pheromone blends than are distantly related species. This pattern is also strong within the melanogaster species group, but less so within the virilis species group. Our analysis strongly suggests that the aggregation pheromones of Drosophila exhibit a gradual, not saltational, mode of evolution. We propose that these findings reflect the function of the pheromones in the ecology of Drosophila, which does not hinge on species specificity of aggregation pheromones as signals.

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.

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.

Parallel evolution of chimeric fusion genes

To understand how novel functions arise, we must identify common patterns and mechanisms shaping the evolution of new genes. Here, we take advantage of data from three Drosophila genes, jingwei, Adh-Finnegan, and Adh-Twain, to find evolutionary patterns and mechanisms governing the evolution of new genes. All three of these genes are independently derived from Adh, which enabled us to use the extensive literature on Adh in Drosophila to guide our analyses. We discovered a fundamental similarity in the temporal, spatial, and types of amino acid changes that occurred. All three genes underwent rapid adaptive amino acid evolution shortly after they were formed, followed by later quiescence and functional constraint. These genes also show striking parallels in which amino acids change in the Adh region. We showed that these early changes tend to occur at amino acid residues that seldom, if ever, evolve in Drosophila Adh. Changes at these slowly evolving sites are usually associated with loss of function or hypomorphic mutations in Drosophila melanogaster. Our data indicate that shifting away from ancestral functions may be a critical step early in the evolution of chimeric fusion genes. We suggest that the patterns we observed are both general and predictive.

Evolutionary turnover of two pBuM satellite DNA subfamilies in the Drosophila buzzatii species cluster (repleta group): From alpha to alpha/beta arrays

The pBuM satellite DNA family was studied in seven Drosophila species from the buzzatii cluster (within the large Drosophila repleta group). The pBuM repeats are slightly AT-rich and show high levels of intraspecific sequence homogeneity. The pBuM family can be divided into two subfamilies. The pBuM-1 subfamily consists of tandemly arranged repetition units of approximately 190 bp, termed alpha. Alpha repeats were found in a high copy number in the genome of D. buzzatii, D. serido and D. antonietae. The pBuM-2 subfamily consists of tandemly arranged repetition units of 370 bp. Its origin is explained by an insertion of an approximately 180 bp foreign sequence (termed beta) in an alpha basic repeat unit, with subsequent homogenization/amplification events increasing its frequency. Alpha/beta repeats were found in a high copy number in the genome of D. serido, D. antonietae, D. seriema and D. gouveai. pBuM sequences were not detected in D. koepferae and D. borborema by hybridization experiments. The nucleotide analysis of 74 pBuM repeats revealed that apart from the beta insertion event, the evolution of the pBuM family has proceeded in a gradual fashion, mainly through accumulation and horizontal spread of nucleotide substitutions. Moreover, the data also indicate a faster evolutionary rate for the pBuM-2 subfamily than the pBuM-1 subfamily. Members of both subfamilies display a greater intraspecific than interspecific homogeneity, indicating a concerted mode of pBuM evolution. A scenario to explain the evolution of both satDNA subfamilies in the seven Drosophila species from the buzzatii cluster is proposed.

The origin and radiation of Macaronesian beetles breeding in Euphorbia: the relative importance of multiple data partitions and population sampling

Species-level phylogenies derived from many independent character sources and wide geographical sampling provide a powerful tool in assessing the importance of various factors associated with cladogenesis. In this study, we explore the relative importance of insular isolation and host plant switching in the diversification of a group of bark beetles (Curculionidae: Scolytinae) feeding and breeding in woody Euphor biaspurges. All species in the genus Aphanarthrumare each associated with only one species group of Euphorbia(succulents or one of three different arborescent groups), and the majority of species are endemic to one or several of the Macaronesian Islands. Hence, putative mechanisms of speciation could be assessed by identifying pairs of sister species in a phylogenetic analysis. We used DNA sequences from two nuclear and two mitochondrial genes, and morphological characters, to reconstruct the genealogical relationships among 92 individuals of 25 species and subspecies of Aphanarthrumand related genera. A stable tree topology was highly dependent on multiple character sources, but much less so on wide population sampling. However, multiple samples per species demonstrated one case of species paraphyly, as well as deep coalescence among three putative subspecies pairs. The phylogenetic analyses consistently placed the arborescent breeding and West African--Lanzarote-distributed species A. armatumin the most basal position in Aphanarthrum, rendering this genus paraphyletic with respect to Coleobothrus. Two major radiations followed, one predominantly African lineage of succulent feeding species, and one island radiation associated with arborescent host plants. Sister comparisons showed that most recent divergences occurred in allopatry on closely related hosts, with subsequent expansions obscuring more ancient events. Only 6 out of 24 cladogenetic events were associated with host switching, rendering geographical factors more important in recent diversification.

Phylogeny of the bears (Ursidae) based on nuclear and mitochondrial genes

The taxomic classification and phylogenetic relationships within the bear family remain argumentative subjects in recent years. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence data, herein we employ two nuclear single-copy gene segments, the partial exon 1 from gene encoding interphotoreceptor retinoid binding protein (IRBP) and the complete intron 1 from transthyretin (TTR) gene, in conjunction with previously published mt data, to clarify these enigmatic problems. The combined analyses of nuclear IRBP and TTR datasets not only corroborated prior hypotheses, positioning the spectacled bear most basally and grouping the brown and polar bear together but also provided new insights into the bear phylogeny, suggesting the sister-taxa association of sloth bear and sun bear with strong support. Analyses based on combination of nuclear and mt genes differed from nuclear analysis in recognizing the sloth bears as the earliest diverging species among the subfamily ursine representatives while the exact placement of the sun bear did not resolved. Asiatic and American black bears clustered as sister group in all analyses with moderate levels of bootstrap support and high posterior probabilities. Comparisons between the nuclear and mtDNA findings suggested that our combined nuclear dataset have the resolving power comparable to mtDNA dataset for the phylogenetic interpretation of the bear family. As can be seen from present study, the unanimous phylogeny for this recently derived family was still not produced and additional independent genetic markers were in need.

Hierarchical structure in the Drosophila mojavensis cluster (Diptera: Drosophilidae)

Sequences from three gene regions from the nuclear and mitochondrial genomes were used to examine the extent and nature of hierarchical structure in the Drosophilamojavensis cluster (Drosophila arizonae, D. mojavensis and D. navojoa) of the D. repleta species group. To determine the genetic divergence of these three species, sequence data were analyzed using maximum parsimony and population aggregation analysis. Individual and combined gene genealogies indicate that D. arizonae and D. mojavensis are neither diagnosable nor monophyletic with respect to one another. Although D. navojoa has differentiated from D. arizonae and D. mojavensis, as diagnosed by nuclear gene sequences, it may have undergone a reticulation event with D. arizonae. Our results suggest that either these taxa are still undergoing differentiation at the molecular level or have experienced gene flow in the recent past.

Evolution of water conservation mechanisms in Drosophila

Flies of the genus Drosophila inhabit a wide range of habitats, from the tropics to deserts to boreal forests. The primary physiological mechanism allowing Drosophila and other insects to survive in arid habitats is a reduction in rates of water loss. To understand mechanisms of water retention in greater detail, we investigated the three main routes by which Drosophila lose water: excretion, cuticular transpiration and respiratory loss through the spiracles. Excretory losses comprised <6% of total water flux and did not differ between xeric (cactophilic) and mesic species. No consistent relationship was observed between water-loss rates and the composition, physical properties or amounts of cuticular hydrocarbons, suggesting that cuticular transpiration did not differ among species from different habitats. Metabolic rates and water-loss rates were highly correlated. Cactophilic Drosophila were less active, and female cactophiles had lower metabolic rates than female mesic species of the same size. They were also more likely to exhibit a pattern of cyclic CO(2) release that may help to conserve water. We conclude that lower overall rates of water loss are achieved primarily by reduction of respiratory losses.