Faster-X Evolution of Gene Expression in Drosophila

DNA sequences on X chromosomes often have a faster rate of evolution when compared to similar loci on the autosomes, and well articulated models provide reasons why the X-linked mode of inheritance may be responsible for the faster evolution of X-linked genes. We analyzed microarray and RNA–seq data collected from females and males of six Drosophila species and found that the expression levels of X-linked genes also diverge faster than autosomal gene expression, similar to the “faster-X” effect often observed in DNA sequence evolution. Faster-X evolution of gene expression was recently described in mammals, but it was limited to the evolutionary lineages shortly following the creation of the therian X chromosome. In contrast, we detect a faster-X effect along both deep lineages and those on the tips of the Drosophila phylogeny. In Drosophila males, the dosage compensation complex (DCC) binds the X chromosome, creating a unique chromatin environment that promotes the hyper-expression of X-linked genes. We find that DCC binding, chromatin environment, and breadth of expression are all predictive of the rate of gene expression evolution. In addition, estimates of the intraspecific genetic polymorphism underlying gene expression variation suggest that X-linked expression levels are not under relaxed selective constraints. We therefore hypothesize that the faster-X evolution of gene expression is the result of the adaptive fixation of beneficial mutations at X-linked loci that change expression level in cis. This adaptive faster-X evolution of gene expression is limited to genes that are narrowly expressed in a single tissue, suggesting that relaxed pleiotropic constraints permit a faster response to selection. Finally, we present a conceptional framework to explain faster-X expression evolution, and we use this framework to examine differences in the faster-X effect between Drosophila and mammals.     

In Vivo RNAi Rescue in Drosophila melanogaster with Genomic Transgenes from Drosophila pseudoobscura

Background

Systematic, large-scale RNA interference (RNAi) approaches are very valuable to systematically investigate biological processes in cell culture or in tissues of organisms such as Drosophila. A notorious pitfall of all RNAi technologies are potential false positives caused by unspecific knock-down of genes other than the intended target gene. The ultimate proof for RNAi specificity is a rescue by a construct immune to RNAi, typically originating from a related species.

Methodology/Principal Findings

We show that primary sequence divergence in areas targeted by Drosophila melanogaster RNAi hairpins in five non-melanogaster species is sufficient to identify orthologs for 81% of the genes that are predicted to be RNAi refractory. We use clones from a genomic fosmid library of Drosophila pseudoobscura to demonstrate the rescue of RNAi phenotypes in Drosophila melanogaster muscles. Four out of five fosmid clones we tested harbour cross-species functionality for the gene assayed, and three out of the four rescue a RNAi phenotype in Drosophila melanogaster.

Conclusions/Significance

The Drosophila pseudoobscura fosmid library is designed for seamless cross-species transgenesis and can be readily used to demonstrate specificity of RNAi phenotypes in a systematic manner.     

Sex chromosome evolution: life,death and repetitive DNA

Dimorphic sex chromosomes create problems. Males of many species, including Drosophila, are heterogametic, with dissimilar X and Y chromosomes. The essential process of dosage compensation modulates the expression of X-linked genes in one sex to maintain a constant ratio of X to autosomal expression. This involves the regulation of hundreds of dissimilar genes whose only shared property is chromosomal address. Drosophila males dosage compensate by up regulating X-linked genes 2 fold. This is achieved by the Male Specific Lethal (MSL) complex, which is recruited to genes on the X chromosome and modifies chromatin to increase expression. How the MSL complex is restricted to X-linked genes remains unknown. Recent studies of sex chromosome evolution have identified a central role for 2 types of repetitive elements in X recognition. Helitrons carrying sites that recruit the MSL complex have expanded across the X chromosome in at least one Drosophila species.¹ Our laboratory found that siRNA from an X-linked satellite repeat promotes X recognition by a yet unknown mechanism.² The recurring adoption of repetitive elements as X-identify elements suggests that the large and mysterious fraction of the genome called “junk” DNA is actually instrumental in the evolution of sex chromosomes.     

Dosage Compensation of X-Linked Muller Element F Genes but Not X-Linked Transgenes in the Australian Sheep Blowfly

In most animals that have X and Y sex chromosomes, chromosome-wide mechanisms are used to balance X-linked gene expression in males and females. In the fly Drosophila melanogaster, the dosage compensation mechanism also generally extends to X-linked transgenes. Over 70 transgenic lines of the Australian sheep blowfly Lucilia cuprina have been made as part of an effort to develop male-only strains for a genetic control program of this major pest of sheep. All lines carry a constitutively expressed fluorescent protein marker gene. In all 12 X-linked lines, female larvae show brighter fluorescence than male larvae, suggesting the marker gene is not dosage compensated. This has been confirmed by quantitative RT-PCR for selected lines. To determine if endogenous X-linked genes are dosage compensated, we isolated 8 genes that are orthologs of genes that are on the fourth chromosome in D. melanogaster. Recent evidence suggests that the D. melanogaster fourth chromosome, or Muller element F, is the ancestral X chromosome in Diptera that has reverted to an autosome in Drosophila species. We show by quantitative PCR of male and female DNA that 6 of the 8 linkage group F genes reside on the X chromosome in L. cuprina. The other two Muller element F genes were found to be autosomal in L. cuprina, whereas two Muller element B genes were found on the same region of the X chromosome as the L. cuprina orthologs of the D. melanogaster Ephrin and gawky genes. We find that the L. cuprina X chromosome genes are equally expressed in males and females (i.e., fully dosage compensated). Thus, unlike in Drosophila, it appears that the Lucilia dosage compensation system is specific for genes endogenous to the X chromosome and cannot be co-opted by recently arrived transgenes.     

A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

Out of the testis,into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome‐scale sequence and multi‐tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well‐known “out of the testis” origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary‐expressed genes. Therefore, I propose a novel “into the ovary” hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined “out of the testis, into the ovary” evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex‐specific adaptation in Drosophila.     

The non-dosage compensated Lsp1α gene of Drosophila melanogaster escapes acetylation by MOF in larval fat body nuclei, but is flanked by two dosage compensated genes

Background

In Drosophila melanogaster dosage compensation of most X-linked genes is mediated by the male-specific lethal (MSL) complex, which includes MOF. MOF acetylates histone H4 at lysine 16 (H4K16ac). The X-linked Larval serum protein one α (Lsp1α) gene has long been known to be not dosage compensated. Here we have examined possible explanations for why the Lsp1α gene is not dosage compensated.

Results

Quantitative RNase protection analysis showed that the genes flanking Lsp1α are expressed equally in males and females and confirmed that Lsp1α is not dosage compensated. Unlike control X-linked genes, Lsp1α was not enriched for H4K16ac in the third instar larval fat body, the tissue in which the gene is actively expressed. X-linked Lsp1α promoter-lacZ reporter transgenes are enriched for H4K16ac in third instar larval fat body. An X-linked reporter gene bracketed by Lsp1α flanking regions was dosage compensated. One of the genes flanking Lsp1α is expressed in the same tissue. This gene shows a modest enrichment for H4K16ac but only at the part of the gene most distant from Lsp1α. Phylogenetic analyses of the sequences of the genomes of 12 Drosophila species shows that Lsp1α is only present within the melanogaster subgroup of species.

Conclusion

Lsp1α is not modified by the MSL complex but is in a region of the X chromosome that is regulated by the MSL complex. The high activity or tissue-specificity of the Lsp1α promoter does not prevent regulation by the MSL complex. The regions flanking Lsp1α do not appear to block access by the MSL complex. Lsp1α appears to have recently evolved within the melanogaster subgroup of Drosophila species. The most likely explanation for why Lsp1α is not dosage compensated is that the gene has not evolved a mechanism to independently recruit the MSL complex, possibly because of its recent evolutionary origin, and because there appears to be a low level of bound MSL complex in a nearby gene that is active in the same tissue.     

Adaptive Evolution and the Birth of CTCF Binding Sites in the Drosophila Genome

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.     

Comparative Genomic Hybridization (CGH) Reveals a Neo-X Chromosome and Biased Gene Movement in Stalk-Eyed Flies (Genus Teleopsis)

Chromosomal location has a significant effect on the evolutionary dynamics of genes involved in sexual dimorphism, impacting both the pattern of sex-specific gene expression and the rate of duplication and protein evolution for these genes. For nearly all non-model organisms, however, knowledge of chromosomal gene content is minimal and difficult to obtain on a genomic scale. In this study, we utilized Comparative Genomic Hybridization (CGH), using probes designed from EST sequence, to identify genes located on the X chromosome of four species in the stalk-eyed fly genus Teleopsis. Analysis of log₂ ratio values of female-to-male hybridization intensities from the CGH microarrays for over 3,400 genes reveals a strongly bimodal distribution that clearly differentiates autosomal from X-linked genes for all four species. Genotyping of 33 and linkage mapping of 28 of these genes in Teleopsis dalmanni indicate the CGH results correctly identified chromosomal location in all cases. Syntenic comparison with Drosophila indicates that 90% of the X-linked genes in Teleopsis are homologous to genes located on chromosome 2L in Drosophila melanogaster, suggesting the formation of a nearly complete neo-X chromosome from Muller element B in the dipteran lineage leading to Teleopsis. Analysis of gene movement both relative to Drosophila and within Teleopsis indicates that gene movement is significantly associated with 1) rates of protein evolution, 2) the pattern of gene duplication, and 3) the evolution of eyespan sexual dimorphism. Overall, this study reveals that diopsids are a critical group for understanding the evolution of sex chromosomes within Diptera. In addition, we demonstrate that CGH is a useful technique for identifying chromosomal sex-linkage and should be applicable to other organisms with EST or partial genomic information.     

Comparison of the oxidative phosphorylation (OXPHOS) nuclear genes in the genomes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae

Background

In eukaryotic cells, oxidative phosphorylation (OXPHOS) uses the products of both nuclear and mitochondrial genes to generate cellular ATP. Interspecies comparative analysis of these genes, which appear to be under strong functional constraints, may shed light on the evolutionary mechanisms that act on a set of genes correlated by function and subcellular localization of their products.

Results

We have identified and annotated the Drosophila melanogaster, D. pseudoobscura and Anopheles gambiae orthologs of 78 nuclear genes encoding mitochondrial proteins involved in oxidative phosphorylation by a comparative analysis of their genomic sequences and organization. We have also identified 47 genes in these three dipteran species each of which shares significant sequence homology with one of the above-mentioned OXPHOS orthologs, and which are likely to have originated by duplication during evolution. Gene structure and intron length are essentially conserved in the three species, although gain or loss of introns is common in A. gambiae. In most tissues of D. melanogaster and A. gambiae the expression level of the duplicate gene is much lower than that of the original gene, and in D. melanogaster at least, its expression is almost always strongly testis-biased, in contrast to the soma-biased expression of the parent gene.

Conclusions

Quickly achieving an expression pattern different from the parent genes may be required for new OXPHOS gene duplicates to be maintained in the genome. This may be a general evolutionary mechanism for originating phenotypic changes that could lead to species differentiation.     

Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

Assessing the impact of comparative genomic sequence data on the functional annotation of the Drosophila genome

Background

It is widely accepted that comparative sequence data can aid the functional annotation of genome sequences; however, the most informative species and features of genome evolution for comparison remain to be determined.

Results

We analyzed conservation in eight genomic regions (apterous, even-skipped, fushi tarazu, twist, and Rhodopsins 1, 2, 3 and 4) from four Drosophila species (D. erecta, D. pseudoobscura, D. willistoni, and D. littoralis) covering more than 500 kb of the D. melanogaster genome. All D. melanogaster genes (and 78-82% of coding exons) identified in divergent species such as D. pseudoobscura show evidence of functional constraint. Addition of a third species can reveal functional constraint in otherwise non-significant pairwise exon comparisons. Microsynteny is largely conserved, with rearrangement breakpoints, novel transposable element insertions, and gene transpositions occurring in similar numbers. Rates of amino-acid substitution are higher in uncharacterized genes relative to genes that have previously been studied. Conserved non-coding sequences (CNCSs) tend to be spatially clustered with conserved spacing between CNCSs, and clusters of CNCSs can be used to predict enhancer sequences.

Conclusions

Our results provide the basis for choosing species whose genome sequences would be most useful in aiding the functional annotation of coding and cis-regulatory sequences in Drosophila. Furthermore, this work shows how decoding the spatial organization of conserved sequences, such as the clustering of CNCSs, can complement efforts to annotate eukaryotic genomes on the basis of sequence conservation alone.     

Zinc Finger Binding Motifs Do Not Explain Recombination Rate Variation within or between Species of Drosophila

In humans and mice, the Cys₂His₂ zinc finger protein PRDM9 binds to a DNA sequence motif enriched in hotspots of recombination, possibly modifying nucleosomes, and recruiting recombination machinery to initiate Double Strand Breaks (DSBs). However, since its discovery, some researchers have suggested that the recombinational effect of PRDM9 is lineage or species specific. To test for a conserved role of PRDM9-like proteins across taxa, we use the Drosophila pseudoobscura species group in an attempt to identify recombination associated zinc finger proteins and motifs. We leveraged the conserved amino acid motifs in Cys₂His₂ zinc fingers to predict nucleotide binding motifs for all Cys₂His₂ zinc finger proteins in Drosophila pseudoobscura and identified associations with empirical measures of recombination rate. Additionally, we utilized recombination maps from D. pseudoobscura and D. miranda to explore whether changes in the binding motifs between species can account for changes in the recombination landscape, analogous to the effect observed in PRDM9 among human populations. We identified a handful of potential recombination-associated sequence motifs, but the associations are generally tenuous and their biological relevance remains uncertain. Furthermore, we found no evidence that changes in zinc finger DNA binding explains variation in recombination rate between species. We therefore conclude that there is no protein with a DNA sequence specific human-PRDM9-like function in Drosophila. We suggest these findings could be explained by the existence of a different recombination initiation system in Drosophila.     

Yeast Srp1, a nuclear protein related toDrosophila and mouse pendulin,is required for normal migration,division, and integrity of nuclei during mitosis

This paper describes genes from yeast and mouse with significant sequence similarities to aDrosophila gene that encodes the blood cell tumor suppressor pendulin. The protein encoded by the yeast gene, Srp1p, and mouse pendulin share 42% and 51% amino acid identity withDrosophila pendulin, respectively. All three proteins consist of 10.5 degenerate tandem repeats of ～ 42 amino acids each. Similar repeats occur in a superfamily of proteins that includes theDrosophila Armadillo protein. All three proteins contain a consensus sequence for a bipartite nuclear localization signal (NLS) in the N-terminal domain, which is not part of the repeat structure. Confocal microscopic analysis of yeast cells stained with antibodies against Srp1p reveals that this protein is intranuclear throughout the cell cycle. Targeted gene disruption shows thatSRP1 is an essential gene. Despite their sequence similarities,Drosophila and mouse pendulin are unable to rescue the lethality of anSRP1 disruption. We demonstrate that yeast cells depleted of Srp1p arrest in mitosis with a G2 content of DNA. Arrested cells display abnormal structures and orientations of the mitotic spindles, aberrant segregation of the chromatin and the nuclei, and threads of chromatin emanating from the bulk of nuclear DNA. This phenotype suggests that Srplp is required for the normal function of microtubules and the spindle pole bodies, as well as for nuclear integrity. We suggest that Srp1p interacts with multiple components of the cell nucleus that are required for mitosis and discuss its functional similarities to, and differences fromDrosophila pendulin.     

An evolutionary model for the duplication and divergence of esterase genes in Drosophila

Summary The esterase 5 (Est-5 = gene, EST 5 = protein) enzyme in Drosophila pseudoobscura is encoded by one of three paralogous genes, Est-5A, Est5B, and Est-5C, that are tightly clustered on the right arm of the X chromosome. The homologous Est-6 locus in Drosophila melanogaster has only one paralogous neighbor, Est-P. Comparisons of coding and flanking DNA sequences among the three D. pseudoobscura and two D. melanogaster genes suggest that two paralogous genes were present before the divergence of D. pseudoobscura from D. melanogaster and that, later, a second duplication occurred in D. pseudoobscura. Nucleotide sequences of the coding regions of the three D. pseudoobscura genes showed 78–85% similarity in pairwise comparisons, whereas the relatedness between Est-6 and Est-P was only 67%. The higher degree of conservation in D. pseudoobscura likely results from the comparatively recent divergence of Est-5B and Est-5C and from possible gene conversion events between Est-5A and Est-5B. Analyses of silent and replacement site differences in the two exons of the paralogous and orthologous genes in each species indicate that common selective forces are acting on all five loci. Further evidence for common purifying selective constraints comes from the conservation of hydropathy profiles and proposed catalytic residues. However, different levels of amino acid substitution between the paralogous genes in D. melanogaster relative to those in D. pseudoobscura suggest that interspecific differences in selection also exist.Offprint requests to: R.C. Richmond