Rapid evolutionary change, constraints and the maintenance of polymorphism in natural populations of Drosophila melanogaster

Allele frequencies can shift rapidly within natural populations. Under certain conditions, repeated rapid allele frequency shifts can lead to the long-term maintenance of polymorphism. In recent years, studies of the model insect Drosophila melanogaster have suggested that this phenomenon is more common than previously believed and is often driven by some form of balancing selection, such as temporally fluctuating or sexually antagonistic selection. Here we discuss some of the general insights into rapid evolutionary change revealed by large-scale population genomic studies, as well as the functional and mechanistic causes of rapid adaptation uncovered by single-gene studies. As an example of the latter, we consider a regulatory polymorphism of the D. melanogaster fezzik gene. Polymorphism at this site has been maintained at intermediate frequency over an extended period of time. Regular observations from a single population over a period of 7 years revealed significant differences in the frequency of the derived allele and its variance across collections between the sexes. These patterns are highly unlikely to arise from genetic drift alone or from the action of sexually antagonistic or temporally fluctuating selection individually. Instead, the joint action of sexually antagonistic and temporally fluctuating selection can best explain the observed rapid and repeated allele frequency shifts. Temporal studies such as those reviewed here further our understanding of how rapid changes in selection can lead to the long-term maintenance of polymorphism as well as improve our knowledge of the forces driving and limiting adaptation in nature.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

Dynamics and stage-specificity of between-population gene expression divergence in the Drosophila melanogaster larval fat body

Gene expression variation is pervasive across all levels of organismal organization, including development. Few studies, however, have examined variation in developmental transcriptional dynamics among populations, or how it contributes to phenotypic divergence. Indeed, the evolution of gene expression dynamics when both the evolutionary and temporal timescale are comparatively short remains relatively uncharacterized. Here, we examined coding and non-coding gene expression in the fat body of an ancestral African and a derived European Drosophila melanogaster population across three developmental stages spanning ten hours of larval development. Between populations, expression divergence was largely stage-specific. We detected higher expression variation during the late wandering stage, which may be a general feature of this stage. During this stage, we also detected higher and more extensive lncRNA expression in Europe, suggesting that lncRNA expression may be more important in derived populations. Interestingly, the temporal breadth of protein-coding and lncRNA expression became more restricted in the derived population. Taken together with the signatures of potential local adaptation that we detected at the sequence level in 9-25% of candidate genes (those showing evidence of expression divergence between populations), this finding suggests that gene expression becomes more developmental stage-specific during adaptation to new environments. We further used RNAi to identify several candidate genes that likely contribute to known phenotypic divergence between these populations. Our results shed light on the evolution and dynamics of expression variation over short developmental and evolutionary timescales, and how this variation contributes to population and phenotypic divergence.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Natural variation in the transcriptional response of Drosophila melanogaster to oxidative stress

Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species’ ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3′ untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3′ untranslated region.

Sexual antagonism, temporally fluctuating selection, and variable dominance affect a regulatory polymorphism in Drosophila melanogaster

Understanding how genetic variation is maintained within species is a major goal of evolutionary genetics that can shed light on the preservation of biodiversity. Here, we examined the maintenance of a regulatory single-nucleotide polymorphism (SNP) of the X-linked Drosophila melanogaster gene fezzik. The derived variant at this site is at intermediate frequency in many worldwide populations but absent in populations from the ancestral species range in sub-Saharan Africa. We collected and genotyped wild-caught individuals from a single European population biannually over a period of 5 years, which revealed an overall difference in allele frequency between the sexes and a consistent change in allele frequency across seasons in females but not in males. Modeling based on the observed allele and genotype frequencies suggested that both sexually antagonistic and temporally fluctuating selection may help maintain variation at this site. The derived variant is predicted to be female-beneficial and mostly recessive; however, there was uncertainty surrounding our dominance estimates and long-term modeling projections suggest that it is more likely to be dominant. By examining gene expression phenotypes, we found that phenotypic dominance was variable and dependent upon developmental stage and genetic background, suggesting that dominance may be variable at this locus. We further determined that fezzik expression and genotype are associated with starvation resistance in a sex-dependent manner, suggesting a potential phenotypic target of selection. By characterizing the mechanisms of selection acting on this SNP, our results improve our understanding of how selection maintains genetic and phenotypic variation in natural populations.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.

The Influence of Chromosomal Environment on X-Linked Gene Expression in Drosophila melanogaster

Sex chromosomes often differ from autosomes with respect to their gene expression and regulation. In Drosophila melanogaster, X-linked genes are dosage compensated by having their expression upregulated in the male soma, a process mediated by the X-chromosome-specific binding of the dosage compensation complex (DCC). Previous studies of X-linked gene expression found a negative correlation between a gene’s male-to-female expression ratio and its distance to the nearest DCC binding site in somatic tissues, including head and brain, which suggests that dosage compensation influences sex-biased gene expression. A limitation of the previous studies, however, was that they focused on endogenous X-linked genes and, thus, could not disentangle the effects of chromosomal position from those of gene-specific regulation. To overcome this limitation, we examined the expression of an exogenous reporter gene inserted at many locations spanning the X chromosome. We observed a negative correlation between the male-to-female expression ratio of the reporter gene and its distance to the nearest DCC binding site in somatic tissues, but not in gonads. A reporter gene’s location relative to a DCC binding site had greater influence on its expression than the local regulatory elements of neighboring endogenous genes, suggesting that intra-chromosomal variation in the strength of dosage compensation is a major determinant of sex-biased gene expression. Average levels of sex-biased expression did not differ between head and brain, but there was greater positional effect variation in the brain, which may explain the observed excess of endogenous sex-biased genes located on the X chromosome in this tissue.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Population Genetic and Functional Analysis of a cis-Regulatory Polymorphism in the DrosophilamelanogasterMetallothionein A gene

Although gene expression can vary extensively within and among populations, the genetic basis of this variation and the evolutionary forces that maintain it are largely unknown. In Drosophilamelanogaster, a 49-bp insertion/deletion (indel) polymorphism in the Metallothionein A (MtnA) gene is associated with variation in MtnA expression and oxidative stress tolerance. To better understand the functional and evolutionary significance of this polymorphism, we investigated it in several worldwide populations. In a German population, the deletion was present at a high and stable frequency over multiple seasons and years, and was associated with increased MtnA expression. There was, however, no evidence that the polymorphism was maintained by overdominant, seasonally fluctuating, or sexually antagonistic selection. The deletion was rare in a population from the species' ancestral range in sub-Saharan Africa and is likely the result of non-African admixture, suggesting that it spread to high frequency following the species' out-of-Africa expansion. Using data from a North American population, we found that the deletion was associated with MtnA expression and tolerance to oxidative stress induced by menadione sodium bisulfite. Our results are consistent with the deletion being selectively favored in temperate populations due to the increased MtnA expression and oxidative stress tolerance that it confers.

Regulation of the X Chromosome in the Germline and Soma of Drosophila melanogaster Males

During the evolution of heteromorphic sex chromosomes, the sex-specific Y chromosome degenerates, while the X chromosome evolves new mechanisms of regulation. Using bioinformatic and experimental approaches, we investigate the expression of the X chromosome in Drosophila melanogaster. We observe nearly complete X chromosome dosage compensation in male somatic tissues, but not in testis. The X chromosome contains disproportionately fewer genes with high expression in testis than the autosomes, even after accounting for the lack of dosage compensation, which suggests that another mechanism suppresses their expression in the male germline. This is consistent with studies of reporter genes and transposed genes, which find that the same gene has higher expression when autosomal than when X-linked. Using a new reporter gene that is expressed in both testis and somatic tissues, we find that the suppression of X-linked gene expression is limited to genes with high expression in testis and that the extent of the suppression is positively correlated with expression level.

Functional characterization of adaptive variation within a cis-regulatory element influencing Drosophila melanogaster growth

Gene expression variation is a major contributor to phenotypic diversity within species and is thought to play an important role in adaptation. However, examples of adaptive regulatory polymorphism are rare, especially those that have been characterized at both the molecular genetic level and the organismal level. In this study, we perform a functional analysis of the Drosophila melanogaster CG9509 enhancer, a cis-regulatory element that shows evidence of adaptive evolution in populations outside the species’ ancestral range in sub-Saharan Africa. Using site-directed mutagenesis and transgenic reporter gene assays, we determined that 3 single nucleotide polymorphisms are responsible for the difference in CG9509 expression that is observed between sub-Saharan African and cosmopolitan populations. Interestingly, while 2 of these variants appear to have been the targets of a selective sweep outside of sub-Saharan Africa, the variant with the largest effect on expression remains polymorphic in cosmopolitan populations, suggesting it may be subject to a different mode of selection. To elucidate the function of CG9509, we performed a series of functional and tolerance assays on flies in which CG9509 expression was disrupted. We found that CG9509 plays a role in larval growth and influences adult body and wing size, as well as wing loading. Furthermore, variation in several of these traits was associated with variation within the CG9509 enhancer. The effect on growth appears to result from a modulation of active ecdysone levels and expression of growth factors. Taken together, our findings suggest that selection acted on 3 sites within the CG9509 enhancer to increase CG9509 expression and, as a result, reduce wing loading as D. melanogaster expanded out of sub-Saharan Africa.

Much of the phenotypic variation that is observed within species is thought to be caused by variation in gene expression. Variants within cis-regulatory elements, which affect the expression of nearby genes within the same DNA strand, are thought to be an abundant resource upon which natural selection can act. Understanding the functional consequences of adaptive cis-regulatory changes is important, as it can help elucidate the mechanisms underlying phenotypic evolution in general and provide insight into the development and maintenance of biodiversity. However, functional analyses of these types of changes remain rare. Here we present a functional analysis of an adaptively evolving enhancer element of a D. melanogaster gene called CG9509, of previously unknown function. We show that 3 single nucleotide polymorphisms located within the enhancer of this gene are responsible for an increase in CG9509 expression in cosmopolitan populations (outside of south and central Africa) relative to sub-Saharan populations, which include ancestral populations. We further show that CG9509 is involved in the regulation of growth rate and body size determination and propose that the CG9509 enhancer underwent positive selection to reduce wing loading as the species expanded out of sub-Saharan Africa.

X-linkage is not a general inhibitor of tissue-specific gene expression in Drosophila melanogaster

As a consequence of its difference in copy number between males and females, the X chromosome is subject to unique evolutionary forces and gene regulatory mechanisms. Previous studies of Drosophila melanogaster have shown that the expression of X-linked, testis-specific reporter genes is suppressed in the male germline. However, it is not known whether this phenomenon is restricted to testis-expressed genes or if it is a more general property of genes with tissue-specific expression, which are also underrepresented on the X chromosome. To test this, we compared the expression of three tissue-specific reporter genes (ovary, accessory gland and Malpighian tubule) inserted at various autosomal and X-chromosomal locations. In contrast to testis-specific reporter genes, we found no reduction of X-linked expression in any of the other tissues. In accessory gland and Malpighian tubule, we detected higher expression of the X-linked reporter genes, which suggests that they are at least partially dosage compensated. We found no difference in the tissue-specificity of X-linked and autosomal reporter genes. These findings indicate that, in general, the X chromosome is not a detrimental environment for tissue-specific gene expression and that the suppression of X-linked expression is limited to the male germline.

De Novo Transcriptome Assembly and Sex-Biased Gene Expression in the Cyclical Parthenogenetic Daphnia galeata

Daphnia species have become models for ecological genomics and exhibit interesting features, such as high phenotypic plasticity and a densely packed genome with many lineage-specific genes. They are also cyclic parthenogenetic, with alternating asexual and sexual cycles and environmental sex determination. Here, we present a de novo transcriptome assembly of over 32,000 D. galeata genes and use it to investigate gene expression in females and spontaneously produced males of two clonal lines derived from lakes in Germany and the Czech Republic. We find that only a low percentage (18%) of genes shows sex-biased expression and that there are many more female-biased gene (FBG) than male-biased gene (MBG). Furthermore, FBGs tend to be more conserved between species than MBGs in both sequence and expression. These patterns may be a consequence of cyclic parthenogenesis leading to a relaxation of purifying selection on MBGs. The two clonal lines show considerable differences in both number and identity of sex-biased genes, suggesting that they may have reproductive strategies differing in their investment in sexual reproduction. Orthologs of key genes in the sex determination and juvenile hormone pathways, which are thought to be important for the transition from asexual to sexual reproduction, are present in D. galeata and highly conserved among Daphnia species.

An Indel Polymorphism in the MtnA 3' Untranslated Region Is Associated with Gene Expression Variation and Local Adaptation in Drosophila melanogaster

Insertions and deletions (indels) are a major source of genetic variation within species and may result in functional changes to coding or regulatory sequences. In this study we report that an indel polymorphism in the 3’ untranslated region (UTR) of the metallothionein gene MtnA is associated with gene expression variation in natural populations of Drosophila melanogaster. A derived allele of MtnA with a 49-bp deletion in the 3' UTR segregates at high frequency in populations outside of sub-Saharan Africa. The frequency of the deletion increases with latitude across multiple continents and approaches 100% in northern Europe. Flies with the deletion have more than 4-fold higher MtnA expression than flies with the ancestral sequence. Using reporter gene constructs in transgenic flies, we show that the 3' UTR deletion significantly contributes to the observed expression difference. Population genetic analyses uncovered signatures of a selective sweep in the MtnA region within populations from northern Europe. We also find that the 3’ UTR deletion is associated with increased oxidative stress tolerance. These results suggest that the 3' UTR deletion has been a target of selection for its ability to confer increased levels of MtnA expression in northern European populations, likely due to a local adaptive advantage of increased oxidative stress tolerance.

Although molecular variation is abundant in natural populations, understanding how this variation affects organismal phenotypes that are subject to natural selection remains a major challenge in the field of evolutionary genetics. Here we show that a deletion mutation in a noncoding region of the Drosophila melanogaster Metallothionein A gene leads to a significant increase in gene expression and increases survival under oxidative stress. The deletion is in high frequency in three distinct geographic regions: in northern European populations, in northern populations along the east coast of North America, and in southern populations along the east coast of Australia. In northern European populations the deletion shows population genetic signatures of recent positive selection. Thus, we provide evidence for a regulatory polymorphism that underlies local adaptation in natural populations.

Variation in the X:Autosome Distribution of Male-Biased Genes among Drosophila melanogaster Tissues and Its Relationship with Dosage Compensation

Genes that are expressed differently between males and females (sex-biased genes) often show a nonrandom distribution in their genomic location, particularly with respect to the autosomes and the X chromosome. Previous studies of Drosophila melanogaster found a general paucity of male-biased genes on the X chromosome, although this is mainly limited to comparisons of whole flies or body segments containing the reproductive organs. To better understand the chromosomal distribution of sex-biased genes in various tissues, we used a common analysis framework to analyze microarray and RNA sequence data comparing male and female gene expression in individual tissues (brain, Malpighian tubule, and gonads), composite structures (head and gonadectomized carcass), and whole flies. Although there are relatively few sex-biased genes in the brain, there is a strong and highly significant enrichment of male-biased genes on the X chromosome. A weaker enrichment of X-linked male-biased genes is seen in the head, suggesting that most of this signal comes from the brain. In all other tissues, there is either no departure from the random expectation or a significant paucity of male-biased genes on the X chromosome. The brain and head also differ from other tissues in that their male-biased genes are significantly closer to binding sites of the dosage compensation complex. We propose that the interplay of dosage compensation and sex-specific regulation can explain the observed differences between tissues and reconcile disparate results reported in previous studies.

Adaptive divergence of a transcriptional enhancer between populations of Drosophila melanogaster

As species colonize new habitats they must adapt to the local environment. Much of this adaptation is thought to occur at the regulatory level; however, the relationships among genetic polymorphism, expression variation and adaptation are poorly understood. Drosophila melanogaster, which expanded from an ancestral range in sub-Saharan Africa around 15 000 years ago, represents an excellent model system for studying regulatory evolution. Here, we focus on the gene CG9509, which differs in expression between an African and a European population of D. melanogaster. The expression difference is caused by variation within a transcriptional enhancer adjacent to the CG9509 coding sequence. Patterns of sequence variation indicate that this enhancer was the target of recent positive selection, suggesting that the expression difference is adaptive. Analysis of the CG9509 enhancer in new population samples from Europe, Asia, northern Africa and sub-Saharan Africa revealed that sequence polymorphism is greatly reduced outside the ancestral range. A derived haplotype absent in sub-Saharan Africa is at high frequency in all other populations. These observations are consistent with a selective sweep accompanying the range expansion of the species. The new data help identify the sequence changes responsible for the difference in enhancer activity.

Population and sex differences in Drosophila melanogaster brain gene expression

Background

Changes in gene regulation are thought to be crucial for the adaptation of organisms to their environment. Transcriptome analyses can be used to identify candidate genes for ecological adaptation, but can be complicated by variation in gene expression between tissues, sexes, or individuals. Here we use high-throughput RNA sequencing of a single Drosophila melanogaster tissue to detect brain-specific differences in gene expression between the sexes and between two populations, one from the ancestral species range in sub-Saharan Africa and one from the recently colonized species range in Europe.

Results

Relatively few genes (<100) displayed sexually dimorphic expression in the brain, but there was an enrichment of sex-biased genes, especially male-biased genes, on the X chromosome. Over 340 genes differed in brain expression between flies from the African and European populations, with the inter-population divergence being highly correlated between males and females. The differentially expressed genes included those involved in stress response, olfaction, and detoxification. Expression differences were associated with transposable element insertions at two genes implicated in insecticide resistance (Cyp6g1 and CHKov1).

Conclusions

Analysis of the brain transcriptome revealed many genes differing in expression between populations that were not detected in previous studies using whole flies. There was little evidence for sex-specific regulatory adaptation in the brain, as most expression differences between populations were observed in both males and females. The enrichment of genes with sexually dimorphic expression on the X chromosome is consistent with dosage compensation mechanisms affecting sex-biased expression in somatic tissues.

Newly evolved genes: moving from comparative genomics to functional studies in model systems. How important is genetic novelty for species adaptation and diversification?

Genes are gained and lost over the course of evolution. A recent study found that over 1,800 new genes have appeared during primate evolution and that an unexpectedly high proportion of these genes are expressed in the human brain. But what are the molecular functions of newly evolved genes and what is their impact on an organism's fitness? The acquisition of new genes may provide a rich source of genetic diversity that fuels evolutionary innovation. Although gene manipulation experiments are not feasible in humans, studies in model organisms, such as Drosophila melanogaster, have shown that new genes can quickly become integrated into genetic networks and become essential for survival or fertility. Future studies of new genes, especially chimeric genes, and their functions will help determine the role of genetic novelty in the adaptation and diversification of species.

Rate of amino acid substitution is influenced by the degree and conservation of male-biased transcription over 50 myr of Drosophila evolution

Sex-biased gene expression (i.e., the differential expression of genes between males and females) is common among sexually reproducing species. However, genes often differ in their sex-bias classification or degree of sex bias between species. There is also an unequal distribution of sex-biased genes (especially male-biased genes) between the X chromosome and the autosomes. We used whole-genome expression data and evolutionary rate estimates for two different Drosophilid lineages, melanogaster and obscura, spanning an evolutionary time scale of around 50 Myr to investigate the influence of sex-biased gene expression and chromosomal location on the rate of molecular evolution. In both lineages, the rate of protein evolution correlated positively with the male/female expression ratio. Genes with highly male-biased expression, genes expressed specifically in male reproductive tissues, and genes with conserved male-biased expression over long evolutionary time scales showed the fastest rates of evolution. An analysis of sex-biased gene evolution in both lineages revealed evidence for a "fast-X" effect in which the rate of evolution was greater for X-linked than for autosomal genes. This pattern was particularly pronounced for male-biased genes. Genes located on the obscura "neo-X" chromosome, which originated from a recent X-autosome fusion, showed rates of evolution that were intermediate between genes located on the ancestral X-chromosome and the autosomes. This suggests that the shift to X-linkage led to an increase in the rate of molecular evolution.

Dissecting gene expression in mosquito

Gene expression is known to vary extensively among tissues and between sexes. However, detailed descriptions of tissue- and sex-specific gene expression are available for only a few model organisms. A new study published in BMC Genomics presents such a data set for the mosquito, Anopheles gambiae, which is the vector of human malaria. In addition to providing a valuable resource for the community of mosquito researchers, the study allows comparative transcriptomic studies of dipteran insects to be extended over 250 million years of evolution, since the divergence of A. gambiae and Drosophila melanogaster.

Population transcriptomics of Drosophila melanogaster females

BACKGROUND

Variation at the level of gene expression is abundant in natural populations and is thought to contribute to the adaptive divergence of populations and species. Gene expression also differs considerably between males and females. Here we report a microarray analysis of gene expression variation among females of 16 Drosophila melanogaster strains derived from natural populations, including eight strains from the putative ancestral range in sub-Saharan Africa and eight strains from Europe. Gene expression variation among males of the same strains was reported previously.

RESULTS

We detected relatively low levels of expression polymorphism within populations, but much higher expression divergence between populations. A total of 569 genes showed a significant expression difference between the African and European populations at a false discovery rate of 5%. Genes with significant over-expression in Europe included the insecticide resistance gene Cyp6g1, as well as genes involved in proteolysis and olfaction. Genes with functions in carbohydrate metabolism and vision were significantly over-expressed in the African population. There was little overlap between genes expressed differently between populations in females and males.

CONCLUSIONS

Our results suggest that adaptive changes in gene expression have accompanied the out-of-Africa migration of D. melanogaster. Comparison of female and male expression data indicates that the vast majority of genes differing in expression between populations do so in only one sex and suggests that most regulatory adaptation has been sex-specific.

On the utility of short intron sequences as a reference for the detection of positive and negative selection in Drosophila

The detection of selection, both positive and negative, acting on a DNA sequence or class of nucleotide sites requires comparison with a reference sequence that is unaffected by selection. In Drosophila, recent findings of widespread selective constraint, as well as adaptive evolution, in both coding and noncoding regions highlight the difficulties in choosing such a reference sequence. Here, we investigate the utility of short intron sequences as a reference for the detection of selection. For a set of 119 Drosophila melanogaster genes containing 195 short introns (<or=120 bp), we analyzed polymorphism and divergence at 1) 4-fold synonymous sites, 2) all sites of introns <or=120 bp, 3) all sites of introns <or=65 bp, 4) bases 8-30 of introns <or=120 bp, and 5) bases 8-30 of introns <or=65 bp. The last class of sites shows the highest levels of both interspecific divergence and intraspecific polymorphism, suggesting that these sites are under the least selective constraint. Bases 8-30 of introns <or=65 bp also have the lowest ratio of divergence to polymorphism, which may indicate that a small proportion of substitutions in the other classes of sites are the result of adaptive evolution. Although there is little signal of selection on the primary sequence of short introns, patterns of insertion-deletion polymorphism and divergence suggest that both positive and negative selection act to maintain an optimal intron length.

Molecular evolution of sex-biased genes in the Drosophila ananassae subgroup

Background

Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura.

Results

Sex-biased expression is generally conserved between D. melanogaster and D. ananassae, with the majority of genes exhibiting the same bias in the two species. However, about one-third of the genes have either gained or lost sex-biased expression in one of the species and a small proportion of genes (~4%) have changed bias from one sex to the other. The male-biased genes of D. ananassae show evidence of positive selection acting at the protein level. However, the signal of adaptive protein evolution for male-biased genes is not as strong in D. ananassae as it is in D. melanogaster and is limited to genes with conserved male-biased expression in both species. Within D. ananassae, a significant signal of adaptive evolution is also detected for female-biased and unbiased genes.

Conclusions

Our findings extend previous observations of widespread adaptive protein evolution to an independent Drosophila lineage, the D. ananassae subgroup. However, the rate of adaptive evolution is not greater for male-biased genes than for female-biased or unbiased genes, which suggests that there are differences in sex-biased gene evolution between the two lineages.

The transcriptional response of Drosophila melanogaster to infection with the sigma virus (Rhabdoviridae)

Background

Bacterial and fungal infections induce a potent immune response in Drosophila melanogaster, but it is unclear whether viral infections induce an antiviral immune response. Using microarrays, we examined the changes in gene expression in Drosophila that occur in response to infection with the sigma virus, a negative-stranded RNA virus (Rhabdoviridae) that occurs in wild populations of D. melanogaster.

Principal Findings

We detected many changes in gene expression in infected flies, but found no evidence for the activation of the Toll, IMD or Jak-STAT pathways, which control immune responses against bacteria and fungi. We identified a number of functional categories of genes, including serine proteases, ribosomal proteins and chorion proteins that were overrepresented among the differentially expressed genes. We also found that the sigma virus alters the expression of many more genes in males than in females.

Conclusions

These data suggest that either Drosophila do not mount an immune response against the sigma virus, or that the immune response is not controlled by known immune pathways. If the latter is true, the genes that we identified as differentially expressed after infection are promising candidates for controlling the host's response to the sigma virus.

The influence of demography and weak selection on the McDonald-Kreitman test: an empirical study in Drosophila

The McDonald-Kreitman (MK) test, which compares the ratio of polymorphism to divergence at nonsynonymous and synonymous sites, is frequently used to detect adaptive evolution in protein-coding sequences. Because the two classes of sites share a common evolutionary history, the MK test is thought to be robust to most demographic factors. However, weak selection on nonsynonymous sites can bias the MK test, especially when a species' effective population size has not been constant. Here, we present an empirical analysis of the influence of demography on the MK test by comparing test results for a common set of 136 genes, including a set of sex-biased genes that shows a strong signal of adaptive evolution, in two Drosophila melanogaster populations: an ancestral population from Africa and a derived population from Europe. The latter has undergone a relatively recent bottleneck, which has reduced its effective population size. We find that the MK test has less power to detect positive selection in the European population for two reasons. First, the overall reduced level of standing variation decreases the statistical power of the test. Second, the segregation of slightly deleterious nonsynonymous mutations biases the MK test away from detecting positive selection. The latter effect is stronger for X-linked genes, which have experienced the greatest reduction in effective population size outside of Africa, and also leads to the underestimation of rates of adaptive protein evolution by multilocus implementations of the MK test. Interestingly, a subset of autosomal female-biased genes shows an increased signal of adaptive evolution in the European population. This is inconsistent with currently accepted demographic scenarios and may reflect female-specific changes in selective constraint following the colonization of non-African habitats.

Effects of X-linkage and sex-biased gene expression on the rate of adaptive protein evolution in Drosophila

Patterns of polymorphism and divergence in Drosophila protein-coding genes suggest that a considerable fraction of amino acid differences between species can be attributed to positive selection and that genes with sex-biased expression, that is, those expressed predominantly in one sex, have especially high rates of adaptive evolution. Previous studies, however, have been restricted to autosomal sex-biased genes and, thus, do not provide a complete picture of the evolutionary forces acting on sex-biased genes across the genome. To determine the effects of X-linkage on sex-biased gene evolution, we surveyed DNA sequence polymorphism and divergence in 45 X-linked genes, including 17 with male-biased expression, 13 with female-biased expression, and 15 with equal expression in the 2 sexes. Using both single- and multilocus tests for selection, we found evidence for adaptive evolution in both groups of sex-biased genes. The signal of adaptive evolution was particularly strong for X-linked male-biased genes. A comparison with data from 91 autosomal genes revealed a "fast-X" effect, in which the rate of adaptive evolution was greater for X-linked than for autosomal genes. This effect was strongest for male-biased genes but could be seen in the other groups as well. A genome-wide analysis of coding sequence divergence that accounted for sex-biased expression also uncovered a fast-X effect for male-biased and unbiased genes, suggesting that recessive beneficial mutations play an important role in adaptation.

Gene expression variation in African and European populations of Drosophila melanogaster

BACKGROUND

Differences in levels of gene expression among individuals are an important source of phenotypic variation within populations. Recent microarray studies have revealed that expression variation is abundant in many species, including Drosophila melanogaster. However, previous expression surveys in this species generally focused on a small number of laboratory strains established from derived populations. Thus, these studies were not ideal for population genetic analyses.

RESULTS

We surveyed gene expression variation in adult males of 16 D. melanogaster strains from two natural populations, including an ancestral African population and a derived European population. Levels of expression polymorphism were nearly equal in the two populations, but a higher number of differences was detected when comparing strains between populations. Expression variation was greatest for genes associated with few molecular functions or biological processes, as well as those expressed predominantly in males. Our analysis also identified genes that differed in expression level between the European and African populations, which may be candidates for adaptive regulatory evolution. Genes involved in flight musculature and fatty acid metabolism were over-represented in the list of candidates.

CONCLUSION

Overall, stabilizing selection appears to be the major force governing gene expression variation within populations. However, positive selection may be responsible for much of the between-population expression divergence. The nature of the genes identified to differ in expression between populations may reveal which traits were important for local adaptation to the European and African environments.

Properties of the Kinesin-1 motor DdKif3 from Dictyostelium discoideum

The amoeba Dictyostelium discoideum possesses genes for 13 different kinesins. Here we characterize DdKif3, a member of the Kinesin-1 family. Kinesin-1 motors form homodimers that can move micrometer-long distances on microtubules using the energy derived from ATP hydrolysis. We expressed recombinant motors in Escherichia coli and tested them in different in vitro assays. Full-length and truncated Kif3 motors were active in gliding and ATPase assays. They showed a strong dependence on ionic strength. Like the full-length motor, the truncated DdKif3-592 motor (aa 1-592; comprising motor domain, neck, and partial stalk) reached its maximum speed of around 2.0micrcom s(-1) at a potassium acetate concentration of 200mM. The shortened DdKif3-342 motor (aa 1-342; comprising motor domain, partial neck) showed a high ATP turnover, comparable to that of the fungal Kinesin-1, Nkin. Results from the duty cycle calculations and gliding assays indicate that DdKif3 is a processive motor. A GFP-fusion protein revealed a mainly cytoplasmic localization of DdKif3. Immunofluorescence staining makes an association with the endoplasmic reticulum or mitochondria unlikely. Despite a similar phylogenetic distance to both metazoa and fungi, in terms of its biochemical properties DdKif3 revealed a closer similarity to fungal than animal kinesins.

X chromosome inactivation during Drosophila spermatogenesis

Genes with male- and testis-enriched expression are under-represented on the Drosophila melanogaster X chromosome. There is also an excess of retrotransposed genes, many of which are expressed in testis, that have “escaped” the X chromosome and moved to the autosomes. It has been proposed that inactivation of the X chromosome during spermatogenesis contributes to these patterns: genes with a beneficial function late in spermatogenesis should be selectively favored to be autosomal in order to avoid inactivation. However, conclusive evidence for X inactivation in the male germline has been lacking. To test for such inactivation, we used a transgenic construct in which expression of a lacZ reporter gene was driven by the promoter sequence of the autosomal, testis-specific ocnus gene. Autosomal insertions of this transgene showed the expected pattern of male- and testis-specific expression. X-linked insertions, in contrast, showed only very low levels of reporter gene expression. Thus, we find that X linkage inhibits the activity of a testis-specific promoter. We obtained the same result using a vector in which the transgene was flanked by chromosomal insulator sequences. These results are consistent with global inactivation of the X chromosome in the male germline and support a selective explanation for X chromosome avoidance of genes with beneficial effects late in spermatogenesis.

During spermatogenesis, the X chromosome is inactivated in the male germline (sperm cells), thereby silencing, or inactivating, genes residing on the X chromosome. X chromosome silencing is thought to be common among species with XY sex determination and has important implications for genome evolution. For example, genes with increased expression in the male tend to be under-represented on the X chromosome, and many testes-specific genes have been “retrotransposed,” or moved, from the sex to autosomal chromosomes. However, compelling evidence for X chromosome inactivation in the fruit fly Drosophila has been lacking. Here, we used a transgenic technique to test for male germline X inactivation in this important model organism. We randomly inserted a “reporter gene” whose expression requires a regulatory element for an autosomal testis-specific gene into multiple autosomal and X-chromosomal locations. We found that autosomal insertions of the reporter gene have significantly higher expression in the male germline than X-linked insertions. This pattern holds for two different transgenes with nearly 50 independent insertions, providing strong evidence for X chromosome inactivation during spermatogenesis. The silencing of X-linked gene expression in the male germline may contribute to the observed paucity of male-expressed genes on the X chromosome and the excess of retrotransposed genes that have moved from the X chromosome to the autosomes.

Genes expressed in the Drosophila testes are under-represented on the X chromosome to avoid the global inactivation of the X chromosome that occurs in the male germ line.

The evolution of sex-biased genes and sex-biased gene expression

Differences between males and females in the optimal phenotype that is favoured by selection can be resolved by the evolution of differential gene expression in the two sexes. Microarray experiments have shown that such sex-biased gene expression is widespread across organisms and genomes. Sex-biased genes show unusually rapid sequence evolution, are often labile in their pattern of expression, and are non-randomly distributed in the genome. Here we discuss the characteristics and expression of sex-biased genes, and the selective forces that shape this previously unappreciated source of phenotypic diversity. Sex-biased gene expression has implications beyond just evolutionary biology, including for medical genetics.

Prevalence of positive selection among nearly neutral amino acid replacements in Drosophila

We have estimated the selective effects of amino acid replacements in natural populations by comparing levels of polymorphism in 91 genes in African populations of Drosophila melanogaster with their divergence from Drosophila simulans. The genes include about equal numbers whose level of expression in adults is greater in males, greater in females, or approximately equal in the sexes. Markov chain Monte Carlo methods were used to sample key parameters in the stationary distribution of polymorphism and divergence in a model in which the selective effect of each nonsynonymous mutation is regarded as a random sample from some underlying normal distribution whose mean may differ from one gene to the next. Our analysis suggests that approximately 95% of all nonsynonymous mutations that could contribute to polymorphism or divergence are deleterious, and that the average proportion of deleterious amino acid polymorphisms in samples is approximately 70%. On the other hand, approximately 95% of fixed differences between species are positively selected, although the scaled selection coefficient (N(e)s) is very small. We estimate that approximately 46% of amino acid replacements have N(e)s < 2, approximately 84% have N(e)s < 4, and approximately 99% have N(e)s < 7. Although positive selection among amino acid differences between species seems pervasive, most of the selective effects could be regarded as nearly neutral. There are significant differences in selection between sex-biased and unbiased genes, which relate primarily to the mean of the distributions of mutational effects and the fraction of slightly deleterious and weakly beneficial mutations that are fixed.

Widespread adaptive evolution of Drosophila genes with sex-biased expression

Many genes in higher eukaryotes show sexually dimorphic expression, and these genes tend to be among the most divergent between species. In most cases, however, it is not known whether this rapid divergence is caused by positive selection or if it is due to a relaxation of selective constraint. To distinguish between these two possibilities, we surveyed DNA sequence polymorphism in 91 Drosophila melanogaster genes with male-, female-, or nonsex-biased expression and determined their divergence from the sister species D. simulans. Using several single- and multilocus statistical tests, we estimated the type and strength of selection influencing the evolution of the proteins encoded by genes of each expression class. Adaptive evolution, as indicated by a relative excess of nonsynonymous divergence between species, was common among the sex-biased genes (both male and female). Male-biased genes, in particular, showed a strong and consistent signal of positive selection, while female-biased genes showed more variation in the type of selection they experience. Genes expressed equally in the two sexes, in contrast, showed no evidence for adaptive evolution between D. melanogaster and D. simulans. This suggests that sexual selection and intersexual coevolution are the major forces driving genetic differentiation between species.

Sebida: a database for the functional and evolutionary analysis of genes with sex-biased expression

We describe Sebida, a database of genes with sex-biased expression. The database integrates results from multiple, independent microarray studies comparing male and female gene expression in Drosophila melanogaster, Drosophila simulans and Anopheles gambiae. Sebida uses standard nomenclature, which allows individual genes to be compared across different microarray platforms and to be queried by gene name, symbol, or annotation number. In addition to ratios of male/female expression for each gene, Sebida also contains information useful for evolutionary studies, such as local recombination rate, degree of codon bias and interspecific divergence at synonymous and non-synonymous sites.

AVAILABILITY

Sebida can be accessed at http://www.sebida.de

Positive correlation between evolutionary rate and recombination rate in Drosophila genes with male-biased expression

Previous studies have shown that genes that are expressed predominantly or exclusively in males tend to evolve rapidly in comparison to other genes. In most cases, however, it is unknown whether this rapid evolution is the result of increased positive (or sexual) selection on male-expressed traits or if it is due to a relaxation of selective constraints. To distinguish between these two possibilities, we analyzed the relationship between the nonsynonymous substitution rate (dN) and local recombination rate for 343 Drosophila genes that were classified as male, female, or nonsex biased in their expression. For the male-biased genes, a positive correlation between dN and recombination rate was observed. This can be explained by an increased rate of adaptive evolution in regions of higher recombination due to a reduction of Hill-Robertson interference. In contrast, the correlation between dN and recombination rate was negative for both female- and nonsex-biased genes, suggesting that these genes are primarily subject to purifying selection, which is expected to be less effective in regions of reduced recombination.

Patterns of synonymous codon usage in Drosophila melanogaster genes with sex-biased expression

The nonrandom use of synonymous codons (codon bias) is a well-established phenomenon in Drosophila. Recent reports suggest that levels of codon bias differ among genes that are differentially expressed between the sexes, with male-expressed genes showing less codon bias than female-expressed genes. To examine the relationship between sex-biased gene expression and level of codon bias on a genomic scale, we surveyed synonymous codon usage in 7276 D. melanogaster genes that were classified as male-, female-, or non-sex-biased in their expression in microarray experiments. We found that male-biased genes have significantly less codon bias than both female- and non-sex-biased genes. This pattern holds for both germline and somatically expressed genes. Furthermore, we find a significantly negative correlation between level of codon bias and degree of sex-biased expression for male-biased genes. In contrast, female-biased genes do not differ from non-sex-biased genes in their level of codon bias and show a significantly positive correlation between codon bias and degree of sex-biased expression. These observations cannot be explained by differences in chromosomal distribution, mutational processes, recombinational environment, gene length, or absolute expression level among genes of the different expression classes. We propose that the observed codon bias differences result from differences in selection at synonymous and/or linked nonsynonymous sites between genes with male- and female-biased expression.

Identification of a locus under complex positive selection in Drosophila simulans by haplotype mapping and composite-likelihood estimation

The recent action of positive selection is expected to influence patterns of intraspecific DNA sequence variation in chromosomal regions linked to the selected locus. These effects include decreased polymorphism, increased linkage disequilibrium, and an increased frequency of derived variants. These effects are all expected to dissipate with distance from the selected locus due to recombination. Therefore, in regions of high recombination, it should be possible to localize a target of selection to a relatively small interval. Previously described patterns of intraspecific variation in three tandemly arranged, testes-expressed genes (janusA, janusB, and ocnus) in Drosophila simulans included all three of these features. Here we expand the original sample and also survey nucleotide polymorphism at three neighboring loci. On the basis of recombination events between derived and ancestral alleles, we localize the target of selection to a 1.5-kb region surrounding janusB. A composite-likelihood-ratio test based on the spatial distribution and frequency of derived polymorphic variants corroborates this result and provides an estimate of the strength of selection. However, the data are difficult to reconcile with the simplest model of positive selection, whereas a new composite-likelihood method suggests that the data are better described by a model in which the selected allele has not yet gone to fixation.

Functional analysis of Drosophila melanogaster gene regulatory sequences by transgene coplacement

The function of putative regulatory sequences identified by comparative genomics can be elucidated only through experimentation. Here the effectiveness of using heterologous gene constructs and transgene coplacement to characterize regulatory sequence function is demonstrated. This method shows that a sequence in the Adh 3'-untranslated region negatively regulates expression, independent of gene or chromosomal context.

Pleiotropic effect of disrupting a conserved sequence involved in a long-range compensatory interaction in the Drosophila Adh gene

Recent advances in experimental analyses of the evolution of RNA secondary structures suggest a more complex scenario than that typically considered by Kimura's classical model of compensatory evolution. In this study, we examine one such case in more detail. Previous experimental analysis of long-range compensatory interactions between the two ends of Drosophila Adh mRNA failed to fit the classical model of compensatory evolution. To further investigate and verify long-range pairing in Drosophila Adh with respect to models of compensatory evolution and its potential functional role, we introduced site-directed mutations in the Drosophila melanogaster Adh gene. We explore two alternative hypotheses for why previous analysis of long-range compensatory interactions failed to fit the classical model. Specifically, we investigate whether the disruption of a conserved short-range pairing within Adh exon 2 has an effect on Adh expression or if there is a dual functional role of a conserved sequence in the 3'-UTR in both long-range pairing and the negative regulation of Adh expression. We find that a classical result was not observed due to the pleiotropic effect of changing a nucleotide involved in both long-range base pairing and the negative regulation of gene expression.

Molecular evolution of sex-biased genes in Drosophila

Studies of morphology, interspecific hybridization, protein/DNA sequences, and levels of gene expression have suggested that sex-related characters (particularly those involved in male reproduction) evolve rapidly relative to non-sex-related characters. Here we report a general comparison of evolutionary rates of sex-biased genes using data from cDNA microarray experiments and comparative genomic studies of Drosophila. Comparisons of nonsynonymous/synonymous substitution rates (d(N)/d(S)) between species of the D. melanogaster subgroup revealed that genes with male-biased expression had significantly faster rates of evolution than genes with female-biased or unbiased expression. The difference was caused primarily by a higher d(N) in the male-biased genes. The same pattern was observed for comparisons among more distantly related species. In comparisons between D. melanogaster and D. pseudoobscura, genes with highly biased male expression were significantly more divergent than genes with highly biased female expression. In many cases, orthologs of D. melanogaster male-biased genes could not be identified in D. pseudoobscura through a Blast search. In contrast to the male-biased genes, there was no clear evidence for accelerated rates of evolution in female-biased genes, and most comparisons indicated a reduced rate of evolution in female-biased genes relative to unbiased genes. Male-biased genes did not show an increased ratio of nonsynonymous/synonymous polymorphism within D. melanogaster, and comparisons of polymorphism/divergence ratios suggest that the rapid evolution of male-biased genes is caused by positive selection.

Selective Constraints on Intron Evolution in Drosophila

Intron sizes show an asymmetrical distribution in a number of organisms, with a large number of “short” introns clustered around a minimal intron length and a much broader distribution of longer introns. In Drosophila melanogaster, the short intron class is centered around 61 bp. The narrow length distribution suggests that natural selection may play a role in maintaining intron size. A comparison of 15 orthologous introns among species of the D. melanogaster subgroup indicates that, in general, short introns are not under greater DNA sequence or length constraints than long introns. There is a bias toward deletions in all introns (deletion/insertion ratio is 1.66), and the vast majority of indels are of short length (&lt;10 bp). Indels occurring on the internal branches of the phylogenetic tree are significantly longer than those occurring on the terminal branches. These results are consistent with a compensatory model of intron length evolution in which slightly deleterious short deletions are frequently fixed within species by genetic drift, and relatively rare larger insertions that restore intron length are fixed by positive selection. A comparison of paralogous introns shared among duplicated genes suggests that length constraints differ between introns within the same gene. The janusA, janusB, and ocnus genes share two short introns derived from a common ancestor. The first of these introns shows significantly fewer indels than the second intron, although the two introns show a comparable number of substitutions. This indicates that intron-specific selective constraints have been maintained following gene duplication, which preceded the divergence of the D. melanogaster species subgroup.

Rapid evolution of male-biased gene expression in Drosophila

A number of genes associated with sexual traits and reproduction evolve at the sequence level faster than the majority of genes coding for non-sex-related traits. Whole genome analyses allow this observation to be extended beyond the limited set of genes that have been studied thus far. We use cDNA microarrays to demonstrate that this pattern holds in Drosophila for the phenotype of gene expression as well, but in one sex only. Genes that are male-biased in their expression show more variation in relative expression levels between conspecific populations and two closely related species than do female-biased genes or genes with sexually monomorphic expression patterns. Additionally, elevated ratios of interspecific expression divergence to intraspecific expression variation among male-biased genes suggest that differences in rates of evolution may be due in part to natural selection. This finding has implications for our understanding of the importance of sexual dimorphism for speciation and rates of phenotypic evolution.

Patterns of DNA sequence variation suggest the recent action of positive selection in the janus-ocnus region of Drosophila simulans

Levels of nucleotide polymorphism in three paralogous Drosophila simulans genes, janusA (janA), janusB (janB), and ocnus (ocn), were surveyed by DNA sequencing. The three genes lie in tandem within a 2.5-kb region of chromosome arm 3R. In a sample of eight alleles from a worldwide distribution we found a significant departure from neutrality by several statistical tests. The most striking feature of this sample was that in a 1.7-kb region containing the janA and janB genes, 30 out of 31 segregating sites contained variants present only once in the sample, and 29 of these unique variants were found in the same allele. A restriction survey of an additional 28 lines of D. simulans revealed strong linkage disequilibrium over the janA-janB region and identified six more alleles matching the rare haplotype. Among the rare alleles, the level of DNA sequence variation was typical for D. simulans autosomal genes and showed no departure from neutrality. In addition, the rare haplotype was more similar to the D. melanogaster sequence, indicating that it was the ancestral form. These results suggest that the derived haplotype has risen to high worldwide frequency relatively recently, most likely as a result of natural selection.

Molecular evolution of the ocnus and janus genes in the Drosophila melanogaster species subgroup

Genes involved in male fertility are potential targets for sexual selection, and their evolution may play a role in reproductive isolation and speciation. Here we describe a new Drosophila melanogaster gene, ocnus (ocn), that encodes a protein abundant in testes nuclear extracts. RT-PCR indicates that ocn transcription is limited to males and is specific to testes. ocn shares homology with another testis-specific gene, janusB (janB), and is located just distal to janB on chromosome 3. The two genes also share homology with the adjacent janusA (janA) gene, suggesting that multiple duplication events have occurred within this region of the genome. We cloned and sequenced these three genes from species of the D. melanogaster species subgroup. Phylogenetic analysis based on protein-encoding sequences predicts a duplication pattern of janA --> janA janB --> janA janB ocn, with the latter event occurring after the divergence of the D. melanogaster and Drosophila obscura species groups. We found significant heterogeneity in the rates of evolution among the three genes within the D. melanogaster species subgroup as measured by the ratio of nonsynonymous to synonymous substitutions, suggesting that diversification of gene function followed each duplication event and that each gene evolved under different selective constraints. All three genes showed faster rates of evolution than genes encoding proteins with metabolic function. These results are consistent with previous studies that have detected an increased rate of evolution in genes with reproductive function.

Stacking and stability of RNA duplexes containing fluorobenzene and fluorobenzimidazole nucleosides

Six different fluorobenzene or fluorobenzimidazole ribonucleosides and one abasic site were incorporated in oligoribonucleotides. Individual contributions of base stacking and solvation of the modified nucleosides could be determined. In fluorobenzene.fluorobenzimidazole-modified base pairs a duplex stabilizing force was found that points to a weak F...H hydrogen bond. The lipophilicity of the unprotected nucleosides were investigated by determination of 1-octanol water partition coefficients.