Sequence variation of alcohol dehydrogenase (Adh) paralogs in cactophilic Drosophila

Phylogenomics provides insights into the evolution of cactophily and host plant shifts in Drosophila

Cactophilic species of the Drosophila buzzatii cluster (repleta group) comprise an excellent model group to investigate genomic changes underlying adaptation to extreme climate conditions and host plants. In particular, these species form a tractable system to study the transition from chemically simpler breeding sites (like prickly pears of the genus Opuntia) to chemically more complex hosts (columnar cacti). Here, we report four highly contiguous genome assemblies of three species of the buzzatii cluster. Based on this genomic data and inferred phylogenetic relationships, we identified candidate taxonomically restricted genes (TRGs) likely involved in the evolution of cactophily and cactus host specialization. Functional enrichment analyses of TRGs within the buzzatii cluster identified genes involved in detoxification, water preservation, immune system response, anatomical structure development, and morphogenesis. In contrast, processes that regulate responses to stress, as well as the metabolism of nitrogen compounds, transport, and secretion were found in the set of species that are columnar cacti dwellers. These findings are in line with the hypothesis that those genomic changes brought about key mechanisms underlying the adaptation of the buzzatii cluster species to arid regions in South America.

A Bayesian Mutation-Selection Framework for Detecting Site-Specific Adaptive Evolution in Protein-Coding Genes

In recent years, codon substitution models based on the mutation–selection principle have been extended for the purpose of detecting signatures of adaptive evolution in protein-coding genes. However, the approaches used to date have either focused on detecting global signals of adaptive regimes—across the entire gene—or on contexts where experimentally derived, site-specific amino acid fitness profiles are available. Here, we present a Bayesian site-heterogeneous mutation–selection framework for site-specific detection of adaptive substitution regimes given a protein-coding DNA alignment. We offer implementations, briefly present simulation results, and apply the approach on a few real data sets. Our analyses suggest that the new approach shows greater sensitivity than traditional methods. However, more study is required to assess the impact of potential model violations on the method, and gain a greater empirical sense its behavior on a broader range of real data sets. We propose an outline of such a research program.

Computational Sequence Analysis of Inversion Breakpoint Regions in the Cactophilic Drosophila mojavensis Lineage

We investigated rates of chromosomal evolution in Drosophila mojavensis using whole-genome sequence information from D. mojavensis, Drosophila buzzatii, and Drosophila virilis. Drosophila mojavensis is a cactophilic species of the repleta group living under extreme ecological conditions in the deserts of the Southwestern United States and Northwestern México. The genome of D. buzzatii, another member of the repleta group, was recently sequenced and the largest scaffolds anchored to all chromosomes using diverse procedures. Chromosome organization between D. mojavensis and D. buzzatii was compared using MUMmer and GRIMM software. Our results corroborate previous cytological analyses that indicated chromosome 2 differed between these 2 species by 10 inversions, chromosomes X and 5 differed by one inversion each, and chromosome 4 was homosequential. In contrast, we found that chromosome 3 differed by 5 inversions instead of the expected 2 that were previously inferred by cytological analyses. Thirteen of these inversions occurred in the D. mojavensis lineage: 12 are fixed and one of them is a polymorphic inversion previously described in populations from Sonora and Baja California, México. We previously investigated the breakpoints of chromosome 2 inversions fixed in D. mojavensis. Here we characterized the breakpoint regions of the 5 inversions found in chromosome 3 in order to infer the molecular mechanism that generated each inversion and its putative functional consequences. Overall, our results reveal a number of gene alterations at the inversion breakpoints with putative adaptive consequences that point to natural selection as the cause for fast chromosomal evolution in D. mojavensis.

Genomic analysis of the four ecologically distinct cactus host populations of Drosophila mojavensis

BACKGROUND

Relationships between an organism and its environment can be fundamental in the understanding how populations change over time and species arise. Local ecological conditions can shape variation at multiple levels, among these are the evolutionary history and trajectories of coding genes. This study examines the rate of molecular evolution at protein-coding genes throughout the genome in response to host adaptation in the cactophilic Drosophila mojavensis. These insects are intimately associated with cactus necroses, developing as larvae and feeding as adults in these necrotic tissues. Drosophila mojavensis is composed of four isolated populations across the deserts of western North America and each population has adapted to utilize different cacti that are chemically, nutritionally, and structurally distinct.

RESULTS

High coverage Illumina sequencing was performed on three previously unsequenced populations of D. mojavensis. Genomes were assembled using the previously sequenced genome of D. mojavensis from Santa Catalina Island (USA) as a template. Protein coding genes were aligned across all four populations and rates of protein evolution were determined for all loci using a several approaches.

CONCLUSIONS

Loci that exhibited elevated rates of molecular evolution tend to be shorter, have fewer exons, low expression, be transcriptionally responsive to cactus host use and have fixed expression differences across the four cactus host populations. Fast evolving genes were involved with metabolism, detoxification, chemosensory reception, reproduction and behavior. Results of this study give insight into the process and the genomic consequences of local ecological adaptation.

Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

Background

Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. Here we aim to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae.

Results

We used a phylogenetic and population genetic analysis framework to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. We found strong evidence of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains.

Conclusions

The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1250-x) contains supplementary material, which is available to authorized users.

Behavioral evolution accompanying host shifts in cactophilic Drosophila larvae

For plant utilizing insects, the shift to a novel host is generally accompanied by a complex set of phenotypic adaptations. Many such adaptations arise in response to differences in plant chemistry, competitive environment, or abiotic conditions. One less well-understood factor in the evolution of phytophagous insects is the selective environment provided by plant shape and volume. Does the physical structure of a new plant host favor certain phenotypes? Here, we use cactophilic Drosophila, which have colonized the necrotic tissues of cacti with dramatically different shapes and volumes, to examine this question. Specifically, we analyzed two behavioral traits in larvae, pupation height, and activity that we predicted might be related to the ability to utilize variably shaped hosts. We found that populations of D. mojavensis living on lengthy columnar or barrel cactus hosts have greater activity and pupate higher in a laboratory environment than populations living on small and flat prickly pear cactus cladodes. Crosses between the most phenotypically extreme populations suggest that the genetic architectures of these behaviors are distinct. A comparison of activity in additional cactophilic species that are specialized on small and large cactus hosts shows a consistent trend. Thus, we suggest that greater motility and an associated tendency to pupate higher in the laboratory are potential larval adaptations for life on a large plant where space is more abundant and resources may be more sparsely distributed.

Identification of misexpressed genetic elements in hybrids between Drosophila-related species

Crosses between close species can lead to genomic disorders, often considered to be the cause of hybrid incompatibility, one of the initial steps in the speciation process. How these incompatibilities are established and what are their causes remain unclear. To understand the initiation of hybrid incompatibility, we performed reciprocal crosses between two species of Drosophila (D. mojavensis and D. arizonae) that diverged less than 1 Mya. We performed a genome-wide transcriptomic analysis on ovaries from parental lines and on hybrids from reciprocal crosses. Using an innovative procedure of co-assembling transcriptomes, we show that parental lines differ in the expression of their genes and transposable elements. Reciprocal hybrids presented specific gene categories and few transposable element families misexpressed relative to the parental lines. Because TEs are mainly silenced by piwi-interacting RNAs (piRNAs), we hypothesize that in hybrids the deregulation of specific TE families is due to the absence of such small RNAs. Small RNA sequencing confirmed our hypothesis and we therefore propose that TEs can indeed be major players of genome differentiation and be implicated in the first steps of genomic incompatibilities through small RNA regulation.

Genome Evolution in Three Species of Cactophilic Drosophila

We report genomes of two species of cactophilic Drosophila: Drosophila arizonae and D. navojoa These two are the closest relatives of D. mojavensis, forming the D. mojavensis cluster. D. mojavensis and D. arizonae diverged from D. navojoa ∼5.8 Mya, while the split between D. arizonae and D. mojavensis is more recent, at 1.5 Mya. Together the three genomes provide opportunities to examine genomic changes associated with speciation and host shifts in this ecologically defined group of flies. The three species are also separated by fixed inversion differences in three of their six chromosomes. While the levels of nucleotide divergence in the colinear chromosomes are significantly lower than in the inverted chromosomes, consistent with a past role of the inversions in preventing gene flow, the patterns differ among the inverted chromosomes when the locations of nucleotides inside or outside of the inversions are considered. For Muller element E, there is greater divergence external to the inversion breakpoints. For Muller A, the divergence is slightly higher inside the inversions, while for Muller B, the breakpoints and hence the difference in substitutions in relation to the inversions could not be determined. The differences among the inverted chromosomes, especially once the breakpoints are clearly established, could aid in dating the origins of the inversions.

Molecular evolution of candidate genes involved in post-mating-prezygotic reproductive isolation

Traits involved in post-copulatory interactions between the sexes may evolve rapidly as a result of sexual selection and/or sexual conflict, leading to post-mating-prezygotic (PMPZ) reproductive isolating barriers between diverging lineages. Although the importance of PMPZ isolation is recognized, the molecular basis of such incompatibilities is not well understood. Here, we investigate molecular evolution of a subset of Drosophila mojavensis and Drosophila arizonae reproductive tract genes. These include genes that are transcriptionally regulated by conspecific mating in females, many of which are misregulated in heterospecific crosses, and a set of male genes whose transcripts are transferred to females during mating. As a group, misregulated female genes are not more divergent and do not appear to evolve under different selection pressures than other female reproductive genes. Male transferred genes evolve at a higher rate than testis-expressed genes, and at a similar rate compared to accessory gland protein genes, which are known to evolve rapidly. Four of the individual male transferred genes show patterns of divergent positive selection between D. mojavensis and D. arizonae. Three of the four genes belong to the sperm-coating protein-like family, including an ortholog of antares, which influences female fertility and receptivity in Drosophila melanogaster. Synthesis of these molecular evolutionary analyses with transcriptomics and predicted functional information makes these genes candidates for involvement in PMPZ reproductive incompatibilities between D. mojavensis and D. arizonae.

Ecological genomics of host shifts in Drosophila mojavensis

Advances in next-generation sequencing technologies have liberated our dependency on model laboratory species for answering genomic and transcriptomic level questions. These new techniques have dramatically expanded our breadth of study organisms and have allowed the analysis of species from diverse ecological environments. One such species is the cactophilic Drosophila mojavensis that inhabits the deserts of western North America. These insects feed and develop in the necrotic cacti, feeding largely on the microflora of the necrotic plant tissues. Drosophila mojavensis is composed of four geographically and ecologically separated populations. Each population (Baja California peninsula, mainland Sonoran Desert, Mojave Desert and Santa Catalina Island) utilizes the necrotic tissues of distinct cactus species. The differences in the nutritional and chemical composition of the necroses include a set of toxic compounds to which resident population must adapt. These ecological differences have facilitated many of the life history, behavior, physiological and genetic differences between the cactus host populations. Genomic resources have allowed investigators to examine the genomic and transcriptional level changes associated with the local adaptation of the four D. mojavensis populations, thereby providing further understanding of the genetic mechanism of adaptation and its role in the divergence of ecologically distinct populations.

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.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

Divergence in olfactory host plant preference in D. mojavensis in response to cactus host use

Divergence in host adaptive traits has been well studied from an ecological and evolutionary perspective, but identification of the proximate mechanisms underlying such divergence is less well understood. Behavioral preferences for host plants are often mediated by olfaction and shifts in preference may be accompanied by changes in the olfactory system. In this study, we examine the evolution of host plant preferences in cactophilic Drosophila mojavensis that feeds and breeds on different cacti throughout its range. We show divergence in electrophysiological responses and olfactory behavior among populations with host plant shifts. Specifically, significant divergence was observed in the Mojave Desert population that specializes on barrel cactus. Differences were observed in electrophysiological responses of the olfactory organs and in behavioral responses to barrel cactus volatiles. Together our results suggest that the peripheral nervous system has changed in response to different ecological environments and that these changes likely contribute to divergence among D. mojavensis populations.

Transcriptome-wide expression variation associated with environmental plasticity and mating success in cactophilic Drosophila mojavensis

Ecological speciation occurs with the adaptation of populations to different environments and concurrent evolution of reproductive isolation. Phenotypic plasticity might influence both ecological adaptation and reproductive traits. We examined environment-specific gene expression and male mating success in cactophilic Drosophila mojavensis using transcriptome sequencing. This species exhibits cactus-dependent mating success across different species of host plants, with genotype-by-environment interactions for numerous traits. We cultured flies from egg to eclosion on two natural cactus hosts and surveyed gene expression in adult males that were either successful or unsuccessful in achieving copulation in courtship trials. We identified gene expression differences that included functions involved with metabolism, most likely related to chemical differences between host cactus species. Several epigenetic-related functions were identified that might play a role in modulating gene expression in adults due to host cactus effects on larvae, and mating success. Cactus-dependent mating success involved expression differences of genes implicated in translation, transcription, and nervous system development. This suggests a role of neurological function genes in the mating success of D. mojavensis males. Together, these results suggest that the influence of environmental variation on mating success via regulation of gene expression might be an important aspect of ecological speciation.

Evolutionary toxicogenomics: diversification of the Cyp12d1 and Cyp12d3 genes in Drosophila species

Gene duplication and divergence are overwhelmingly considered to be the primary mechanisms by which cytochrome P450 monooxygenases (P450s) have radiated into a large and diverse gene superfamily. To address how environmental stress drives the fixation and diversification of gene duplications, we have analyzed Cyp12d1 and Cyp12d3, a pair of duplicated genes found in the sequenced Drosophila genomes of the melanogaster group. The paralog Cyp12d3, which is not found in Drosophila melanogaster, is basal to the melanogaster group, after it split from the obscura group (ca. 50 mya), and has a significant signature of positive selection in two species (D. sechellia and D. ananassae). Examination of the Cyp12d1 region in D. melanogaster wildtype and isoline populations revealed variation both in copy number and sequence, including splice-site variations, which certainly alter gene function. Further investigations of several strains have identified three cases in which differences in the Cyp12d1 gene region are associated with the differences in transcript abundance and transcriptional responses to the environmental stresses that have not been seen for other detoxificative loci. Together, these data highlight the value of using both macro- and microevolutionary approaches in studying the duplication and divergence events associated with detoxification genes and lay important groundwork for future studies in the field of evolutionary toxicogenomics, which uses the principles of phylogenetic analysis to predict possible enzymatic functions.

Model-based comparisons of phylogeographic scenarios resolve the intraspecific divergence of cactophilic Drosophila mojavensis

The cactophilic fly Drosophila mojavensis exhibits considerable intraspecific genetic structure across allopatric geographic regions and shows associations with different host cactus species across its range. The divergence between these populations has been studied for more than 60years, yet their exact historical relationships have not been resolved. We analysed sequence data from 15 intronic X-linked loci across populations from Baja California, mainland Sonora-Arizona and Mojave Desert regions under an isolation-with-migration model to assess multiple scenarios of divergence. We also compared the results with a pre-existing sequence data set of eight autosomal loci. We derived a population tree with Baja California placed at its base and link their isolation to Pleistocene climatic oscillations. Our estimates suggest the Baja California population diverged from an ancestral Mojave Desert/mainland Sonora-Arizona group around 230,000-270,000years ago, while the split between the Mojave Desert and mainland Sonora-Arizona populations occurred one glacial cycle later, 117,000-135,000years ago. Although we found these three populations to be effectively allopatric, model ranking could not rule out the possibility of a low level of gene flow between two of them. Finally, the Mojave Desert population showed a small effective population size, consistent with a historical population bottleneck. We show that model-based inference from multiple loci can provide accurate information on the historical relationships of closely related groups allowing us to set into historical context a classic system of incipient ecological speciation.

Gene alterations at Drosophila inversion breakpoints provide prima facie evidence for natural selection as an explanation for rapid chromosomal evolution

Background

Chromosomal inversions have been pervasive during the evolution of the genus Drosophila, but there is significant variation between lineages in the rate of rearrangement fixation. D. mojavensis, an ecological specialist adapted to a cactophilic niche under extreme desert conditions, is a chromosomally derived species with ten fixed inversions, five of them not present in any other species.

Results

In order to explore the causes of the rapid chromosomal evolution in D. mojavensis, we identified and characterized all breakpoints of seven inversions fixed in chromosome 2, the most dynamic one. One of the inversions presents unequivocal evidence for its generation by ectopic recombination between transposon copies and another two harbor inverted duplications of non-repetitive DNA at the two breakpoints and were likely generated by staggered single-strand breaks and repair by non-homologous end joining. Four out of 14 breakpoints lay in the intergenic region between preexisting duplicated genes, suggesting an adaptive advantage of separating previously tightly linked duplicates. Four out of 14 breakpoints are associated with transposed genes, suggesting these breakpoints are fragile regions. Finally two inversions contain novel genes at their breakpoints and another three show alterations of genes at breakpoints with potential adaptive significance.

Conclusions

D. mojavensis chromosomal inversions were generated by multiple mechanisms, an observation that does not provide support for increased mutation rate as explanation for rapid chromosomal evolution. On the other hand, we have found a number of gene alterations at the breakpoints with putative adaptive consequences that directly point to natural selection as the cause of D. mojavensis rapid chromosomal evolution.

Sequence signatures of a recent chromosomal rearrangement in Drosophila mojavensis

The X-chromosome inversion, Xe, distinguishes Drosophila mojavensis and D. arizonae. Earlier work mapped the breakpoints of this inversion to large intervals and provided hypotheses for the locations of the breakpoints within 3000-bp intergenic regions on the D. mojavensis genome sequence assembly. Here, we sequenced these regions directly in the putatively ancestral D. arizonae X-chromosome. We find that the two inversion breakpoints are near an inverted gene duplication and a common repetitive element, respectively, and these features were likely present in the non-inverted ancestral chromosome on the D. mojavensis lineage. Contrary to an earlier hypothesis, the inverted gene duplication appears to predate the inversion. We find no sequence similarity between the breakpoint regions in the D. mojavensis ancestor, excluding an ectopic-exchange model of chromosome rearrangements. We also found no evidence that staggered single-strand breaks caused the inversion. We suggest these features may have contributed to the chromosomal breakages resulting in this inversion.

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.

The molecular basis of host adaptation in cactophilic Drosophila: molecular evolution of a glutathione S-transferase gene (GstD1) in Drosophila mojavensis

Drosophila mojavensis is a cactophilic fly endemic to the northwestern deserts of North America. This species includes four genetically isolated cactus host races each individually specializing on the necrotic tissues of a different cactus species. The necrosis of each cactus species provides the resident D. mojavensis populations with a distinct chemical environment. A previous investigation of the role of transcriptional variation in the adaptation of D. mojavensis to its hosts produced a set of candidate loci that are differentially expressed in response to host shifts, and among them was glutathione S-transferase D1 (GstD1). In both D. melanogaster and Anopheles gambiae, GstD1 has been implicated in the resistance of these species to the insecticide dichloro-diphenyl-trichloroethane (DDT). The pattern of sequence variation of the GstD1 locus from all four D. mojavensis populations, D. arizonae (sister species), and D. navojoa (outgroup) has been examined. The data suggest that in two populations of D. mojavensis GstD1 has gone through a period of adaptive amino acid evolution. Further analyses indicate that of the seven amino acid fixations that occurred in the D. mojavensis lineage, two of them occur in the active site pocket, potentially having a significant effect on substrate specificity and in the adaptation to alternative cactus hosts.

Multilocus nuclear sequences reveal intra- and interspecific relationships among chromosomally polymorphic species of cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae, a pair of sibling species endemic to North America, constitute an important model system to study ecological genetics and the evolution of reproductive isolation. This species pair can produce fertile hybrids in some crosses and are sympatric in a large part of their ranges. Despite the potential for hybridization in nature, however, evidence of introgression has not been rigorously sought. Further, the evolutionary relationships within and among the geographically distant populations of the two species have not been characterized in detail using high-resolution molecular studies. Both species have six chromosomes: five large acrocentrics and one 'dot' chromosome. Fixed inversion differences between the species exist in three chromosomes (X, 2 and 3) while three are colinear (4, 5 and 6), suggesting that were introgression to occur, it would be most likely in the colinear chromosomes. We utilized nucleotide sequence variation at multiple loci on five chromosomes to test for evidence of introgression, and to test various scenarios for the evolutionary relationships of these two species and their populations. While we do not find evidence of recent introgression, loci in the colinear chromosomes appear to have participated in exchange in the past. We also found considerable population structure within both species. The level of differentiation discovered among D. arizonae populations was unexpectedly high and suggests that its populations, as well as those of D. mojavensis, may be themselves undergoing incipient speciation and merit further attention.

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.

Functional genomics of cactus host shifts in Drosophila mojavensis

Understanding the genetic basis of adaptation to novel environments remains one of the major challenges confronting evolutionary biologists. While newly developed genomic approaches hold considerable promise for addressing this overall question, the relevant tools have not often been available in the most ecologically interesting organisms. Our study organism, Drosophila mojavensis, is a cactophilic Sonoran Desert endemic utilizing four different cactus hosts across its geographical range. Its well-known ecology makes it an attractive system in which to study the evolution of gene expression during adaptation. As a cactophile, D. mojavensis oviposits in the necrotic tissues of cacti, therefore exposing larvae and even adults to the varied and toxic compounds of rotting cacti. We have developed a cDNA microarray of D. mojavensis to examine gene expression associated with cactus host use. Using a population from the Baja California population we examined gene expression differences of third instar larvae when reared in two chemically distinct cactus hosts, agria (Stenocereus gummosus, native host) vs. organpipe (Stenocereus thurberi, alternative host). We have observed differential gene expression associated with cactus host use in genes involved in metabolism and detoxification.

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.

Microsatellite variation among diverging populations of Drosophila mojavensis

Divergence and speciation may occur by various means, depending on the particular history, selective environments, and genetic composition of populations. In Drosophila mojavensis, a good model of incipient speciation, understanding the population genetic structure within this group facilitates our ability to understand the context in which reproductive isolation among populations is developing. Here we report the genetic structure and relationships of D. mojavensis populations at nuclear loci. We surveyed 29 populations throughout the distribution of D. mojavensis for four microsatellite loci to differentiation among populations of this species. These loci reveal four distinct geographical regions of D. mojavensis populations in the south-western United States and north-western Mexico--(i) Baja California peninsula (Baja), (ii) Sonora, Mexico-southern Arizona, United States (Sonora), (iii) Mojave Desert and Grand Canyon (Mojave), and (iv) Santa Catalina Island (Catalina). While all regions show strong isolation, Mojave and Catalina are highly diverged from other regions. Within any region, populations are largely homogenous over broad geographical distances. Based on the population structure, we find clear geographical barriers to gene flow appear to have a strong effect in isolating populations across regions for this species.

Sympatry, allopatry and sexual isolation between Drosophila mojavensis and D. arizonae

Populations of the North American cactophilic fruitfly Drosophila mojavensis and its sibling species D. arizonae exist both in sympatry and in allopatry. Females of D. arizonae, regardless of their population of origin, are effectively completely isolated behaviorally from D. mojavensis males. On the other hand, females of D. mojavensis from the sympatric populations in Sonora, Mexico exhibit significantly stronger premating isolation from D. arizonae males than do D. mojavensis females from allopatric populations from the Baja California peninsula. Earlier studies interpreted these limited observations as support for reinforcement. Since the time of those studies, additional allopatric populations of D. mojavensis have been collected from southern California and from Santa Catalina Island, off the coast of southern California. Here, we tested the prediction that if sympatry is in fact associated with increased isolation in D. mojavensis, these additional allopatric populations also should show, relative to the sympatric ones, less isolation from D. arizonae. Our results are consistent with this prediction and suggest that isolation is in fact stronger in sympatry.

Multilocus test for introgression between the cactophilic species Drosophila mojavensis and Drosophila arizonae

Information obtained from laboratory studies regarding the efficacy of barriers to gene flow (reproductive isolation) between species is often incomplete or misleading, so detailed genetic analyses are needed to determine whether hybridization and introgression occur in nature. Previous laboratory studies of the cactophilic species Drosophila mojavensis and Drosophila arizonae suggest that reproductive isolation is incomplete and that gene flow may occur in sympatry. We sampled 18 nuclear and one mitochondrial loci from multiple populations of D. arizonae and D. mojavensis to test for the signature of recent or historic gene flow between these two species. We located chromosomal regions that were inverted between these species and analyzed those regions independently of others. Statistical tests for introgression using all loci or only collinear loci failed to reject expectations of an isolation model. Further tests using average nucleotide differences between species and phylogenetic analyses also failed to find support for introgression between D. mojavensis and D. arizonae. Additional ecological and behavioral studies of these species in their natural habitats are required to explain why the signature of gene flow was not detected at the DNA sequence level in populations when laboratory studies suggest such gene flow should be possible.

Ribosomal intergenic spacer (IGS) length variation across the Drosophilinae (Diptera: Drosophilidae)

Background

The intergenic spacer of the ribosomal genes in eukaryotes (IGS) contains duplications of the core transcription promoter. The number of these duplicated promoters, as measured by the IGS length, appears to be correlated with growth rate and development time in several distantly related taxa. In the present study, we examined IGS length variation across a number of species of Drosophila to determine the amount of variation in this trait across different evolutionary time scales. Furthermore, we compared the usefulness of two methods commonly used to determine IGS length: Southern Blot Hybridization (SB) and Polymerase Chain Reaction (PCR).

Results

Our results show broad variation in IGS length across the genus Drosophila, but closely related species had similar IGS lengths. Our results also suggest that PCR tends to underestimate the true IGS size when the size is greater than 5 kb, and that this degree of underestimation is greater as the IGS size increases.

Conclusion

Broad variation in IGS length occurs across large evolutionary divergences in the subfamily Drosophilinae. Although average IGS length has been shown to evolve rapidly under artificial selection, closely related taxa generally have similar average IGS lengths. Our comparison of methods suggests that without previous knowledge of the DNA sequence of the IGS and flanking regions, both methods be used to accurately measure IGS length.

Molecular population genetics of accessory gland protein genes and testis-expressed genes in Drosophila mojavensis and D. arizonae

Molecular population genetic investigation of Drosophila male reproductive genes has focused primarily on melanogaster subgroup accessory gland protein genes (Acp's). Consistent with observations from male reproductive genes of numerous taxa, Acp's evolve more rapidly than nonreproductive genes. However, within the Drosophila genus, large data sets from additional types of male reproductive genes and from different species groups are lacking. Here we report findings from a molecular population genetics analysis of male reproductive genes of the repleta group species, Drosophila arizonae and D. mojavensis. We find that Acp's have dramatically higher average pairwise Ka/Ks (0.93) than testis-enriched genes (0.19) and previously reported melanogaster subgroup Acp's (0.42). Overall, 10 of 19 Acp's have Ka/Ks > 1 either in nonpolarized analyses or in at least one lineage of polarized analyses. Of the nine Acp's for which outgroup data were available, average Ka/Ks was considerably higher in D. mojavensis (2.08) than in D. arizonae (0.87). Contrasts of polymorphism and divergence suggest that adaptive protein evolution at Acp's is more common in D. mojavensis than in D. arizonae.

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

BACKGROUND

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

RESULTS

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

CONCLUSION

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.

Inhibition of alcohol dehydrogenase after 2-propanol exposure in different geographic races of Drosophila mojavensis: lack of evidence for selection at the Adh-2 locus

High frequencies of the fast allele of alcohol dehydrogenase-2 (Adh-2F) are found in populations of Drosophila mojavensis that inhabit the Baja California peninsula (race BII) whereas the slow allele (Adh-2S) predominates at most other localities within the species' geographic range. Race BII flies utilize necrotic tissue of pitaya agria cactus (Stenocereus gummosus) which contains high levels of 2-propanol, whereas flies from most other localities utilize different cactus hosts in which 2-propanol levels are low. To test if 2-propanol acts as a selective force on Adh-2 genotype, or whether some other yet undetermined genetic factor is responsible, mature males of D. mojavensis lines derived from the Grand Canyon (race A) and Santa Catalina Island (race C), each with individuals homozygous for Adh-2F and Adh-2S, were exposed to 2-propanol for 24 h and ADH-2 specific activity was then determined on each genotype. Flies from five other localities homozygous for either the fast or slow allele also were examined. Results for all reported races of D. mojavensis were obtained. 2-propanol exposure inhibited ADH-2 specific activity in both genotypes from all localities, but inhibition was significantly less in two populations of race BII flies homozygous for Adh-2F. When F/F and S/S genotypes in flies from the same locality were compared, both genotypes showed high 2-propanol inhibition that was not statistically different, indicating that the F/F genotype alone does not provide a benefit against the inhibitory effects of 2-propanol. ADH-1 activity in female ovaries was inhibited less by 2-propanol than ADH-2. These results do not support the hypothesis that 2-propanol acts as a selective factor favoring the Adh-2F allele.

Activity variation in alcohol dehydrogenase paralogs is associated with adaptation to cactus host use in cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae are species of cactophilic flies that share a recent duplication of the alcohol dehydrogenase (Adh) locus. One paralog (Adh-2) is expressed in adult tissues and the other (Adh-1) in larvae and ovaries. Enzyme activity measurements of the ADH-2 amino acid polymorphism in D. mojavensis suggest that the Fast allozyme allele has a higher activity on 2-propanol than 1-propanol. The Fast allele was found at highest frequency in populations that utilize hosts with high proportions of 2-propanol, while the Slow allele is most frequent in populations that utilize hosts with high proportions of 1-propanol. This suggests that selection for ADH-2 allozyme alleles with higher activity on the most abundant alcohols is occurring in each D. mojavensis population. In the other paralog, ADH-1, significant differences between D. mojavensis and D. arizonae are associated with a previously shown pattern of adaptive protein evolution in D. mojavensis. Examination of protein sequences showed that a large number of amino acid fixations between the paralogs have occurred in catalytic residues. These changes are potentially responsible for the significant difference in substrate specificity between the paralogs. Both functional and sequence variation within and between paralogs suggests that Adh has played an important role in the adaptation of D. mojavensis and D. arizonae to their cactophilic life.

Comparisons of maize pericarp color1 alleles reveal paralogous gene recombination and an organ-specific enhancer region

The maize (Zea mays) p1 (for pericarp color1) gene encodes an R2R3 Myb-like transcription factor that regulates the flavonoid biosynthetic pathway in floral organs, most notably kernel pericarp and cob. Alleles of the p1 gene condition distinct tissue-specific pigmentation patterns; to elucidate the molecular basis of these allele-specific expression patterns, we characterized two novel P1-rw (for red pericarp/white cob) alleles, P1-rw1077 and P1-rw751Ac. Structural analysis of P1-rw1077 indicated that this allele was generated by recombination between p1 and the tightly linked paralogous gene, p2. In the resulting gene, the p1 coding sequence was replaced by the p2 coding sequence, whereas the flanking p1 regulatory sequences remained largely preserved. The red pericarp color specified by P1-rw1077 suggests that the p1- and p2-encoded proteins are functionally equivalent as regulatory factors in the flavonoid biosynthesis pathway. Sequence analysis shows that the P1-rw1077 allele lacks a 386-bp sequence in a distal enhancer region 5 kb upstream of the transcription start site. An independently derived P1-rw allele contains an Ac insertion into the same sequence, indicating that this site likely contains cob glume-specific regulatory elements.

Coalescent processes and relaxation of selective constraints leading to contrasting genetic diversity at paralogs AtHVA22d and AtHVA22e in Arabidopsis thaliana

Duplicate loci offer a very powerful system for understanding the complicated genome structure and adaptive evolution of a gene family. In this study, the genetic variation at paralogs AtHVA22d and AtHVA22e, members of an ABA- and stress-inducible gene family, is examined in the selfing Arabidopsis thaliana. Population genetic analysis indicates contrasting levels of nucleotide diversity at overall exon sequence and nonsynonymous sites between AtHVA22d (pi = 0.00337, pi(rep) = 0.00158) and AtHVA22e (pi = 0.00054, pi(rep) = 0.00023). The fact of Ka/Ks ratios significantly less than 1 in all sequences indicates that both genes are functional and subjected to purifying selection. In addition, rooted at barley HVA22, accelerated evolution is detected at replacement changes in the AtHVA22d locus, indicating relaxation of purifying selection after gene duplication. However, relative rate tests reveal no deviation from the neutrality at synonymous sites between the two paralogs. Based on clock-like evolution, the rate of synonymous substitution is estimated at 1.83 x 10(-9) substitutions per site per year; and the divergence of the two paralogs is traced to 90 MYA, coinciding with a period of the diversification of angiosperms. Given no codon usage bias in both genes, natural selection alone cannot account for the 6.4-fold differences in the nucleotide variation at synonymous sites between the two paralogs. Random processes resulting in different coalescence times, 3.65 MYA at AtHVA22d vs. 1.20 MYA at AtHVA22e, may have predominantly contributed to the evident differences of the genetic diversity. Partially nonoverlapping modes of expression between the two functional paralogs suggest a subfunctionalization hypothesis for explaining the fates of duplicate loci.

Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene

Evidence for the evolutionary maintenance of a hairpin structure possibly involved in intron processing had been found in intron 1 of the alcohol dehydrogenase gene (Adh) in diverse Drosophila species. In this study, the putative hairpin structure was evaluated systematically in Drosophila melanogaster by elimination of either side of the stem using site-directed mutagenesis. The effects of these mutations and the compensatory double mutant on intron splicing efficiency and ADH protein production were assayed in Drosophila melanogaster Schneider L2 cells and germ-line transformed adult flies. Mutations that disrupt the putative hairpin structure right upstream of the intron branch point were found to cause a significant reduction in both splicing efficiency and ADH protein production. In contrast, the compensatory double mutant that restores the putative hairpin structure was indistinguishable from the WT in both splicing efficiency and ADH level. It was also observed by mutational analysis that a more stable secondary structure (with a longer stem) in this intron decreases both splicing efficiency and ADH protein production. Implications for RNA secondary structure and intron evolution are discussed.