Positive selection at the binding sites of the male-specific lethal complex involved in dosage compensation in Drosophila

Intermolecular interactions drive protein adaptive and co-adaptive evolution at both species and population levels

Proteins are the building blocks for almost all the functions in cells. Understanding the molecular evolution of proteins and the forces that shape protein evolution is essential in understanding the basis of function and evolution. Previous studies have shown that adaptation frequently occurs at the protein surface, such as in genes involved in host–pathogen interactions. However, it remains unclear whether adaptive sites are distributed randomly or at regions associated with particular structural or functional characteristics across the genome, since many proteins lack structural or functional annotations. Here, we seek to tackle this question by combining large-scale bioinformatic prediction, structural analysis, phylogenetic inference, and population genomic analysis of Drosophila protein-coding genes. We found that protein sequence adaptation is more relevant to function-related rather than structure-related properties. Interestingly, intermolecular interactions contribute significantly to protein adaptation. We further showed that intermolecular interactions, such as physical interactions, may play a role in the coadaptation of fast-adaptive proteins. We found that strongly differentiated amino acids across geographic regions in protein-coding genes are mostly adaptive, which may contribute to the long-term adaptive evolution. This strongly indicates that a number of adaptive sites tend to be repeatedly mutated and selected throughout evolution in the past, present, and maybe future. Our results highlight the important roles of intermolecular interactions and coadaptation in the adaptive evolution of proteins both at the species and population levels.

Rapid Evolution of Autosomal Binding Sites of the Dosage Compensation Complex in Drosophila melanogaster and Its Association With Transcription Divergence

How pleiotropy influences evolution of protein sequence remains unclear. The male-specific lethal (MSL) complex in Drosophila mediates dosage compensation by 2-fold upregulation of the X chromosome in males. Nevertheless, several MSL proteins also bind autosomes and likely perform functions not related to dosage compensation. Here, we study the evolution of MOF, MSL1, and MSL2 biding sites in Drosophila melanogaster and its close relative Drosophila simulans. We found pervasive expansion of the MSL binding sites in D. melanogaster, particularly on autosomes. The majority of these newly-bound regions are unlikely to function in dosage compensation and associated with an increase in expression divergence between D. melanogaster and D. simulans. While dosage-compensation related sites show clear signatures of adaptive evolution, these signatures are even more marked among autosomal regions. Our study points to an intriguing avenue of investigation of pleiotropy as a mechanism promoting rapid protein sequence evolution.

The Effect of Hybridization on Dosage Compensation in Member Species of the Anopheles gambiae Species Complex

Dosage compensation has evolved in concert with Y-chromosome degeneration in many taxa that exhibit heterogametic sex chromosomes. Dosage compensation overcomes the biological challenge of a “half dose” of X chromosome gene transcripts in the heterogametic sex. The need to equalize gene expression of a hemizygous X with that of autosomes arises from the fact that the X chromosomes retain hundreds of functional genes that are actively transcribed in both sexes and interact with genes expressed on the autosomes. Sex determination and heterogametic sex chromosomes have evolved multiple times in Diptera, and in each case the genetic control of dosage compensation is tightly linked to sex determination. In the Anopheles gambiae species complex (Culicidae), maleness is conferred by the Y-chromosome gene Yob, which despite its conserved role between species is polymorphic in its copy number between them. Previous work demonstrated that male An. gambiae s.s. males exhibit complete dosage compensation in pupal and adult stages. In the present study, we have extended this analysis to three sister species in the An. gambiae complex: An. coluzzii, An. arabiensis, and An. quadriannulatus. In addition, we analyzed dosage compensation in bi-directional F1 hybrids between these species to determine if hybridization results in the mis-regulation and disruption of dosage compensation. Our results confirm that dosage compensation operates in the An. gambiae species complex through the hypertranscription of the male X chromosome. Additionally, dosage compensation in hybrid males does not differ from parental males, indicating that hybridization does not result in the mis-regulation of dosage compensation.

GFZF, a Glutathione S-Transferase Protein Implicated in Cell Cycle Regulation and Hybrid Inviability, Is a Transcriptional Coactivator

The core promoters of protein-encoding genes play a central role in regulating transcription. M1BP is a transcriptional activator that associates with a core promoter element known as Motif 1 that resides at thousands of genes in Drosophila To gain insight into how M1BP functions, we identified an interacting protein called GFZF. GFZF had been previously identified in genetic screens for factors involved in maintenance of hybrid inviability, the G₂-M DNA damage checkpoint, and RAS/mitogen-activated protein kinase (MAPK) signaling, but its contribution to these processes was unknown. Here, we show that GFZF resides in the nucleus and functions as a transcriptional coactivator. In addition, we show that GFZF is a glutathione S-transferase (GST). Thus, GFZF is the first transcriptional coactivator with intrinsic GST activity, and its identification as a transcriptional coactivator provides an explanation for its role in numerous biological processes.

Evolution of dosage compensation under sexual selection differs between X and Z chromosomes

Complete sex chromosome dosage compensation has more often been observed in XY than ZW species. In this study, using a population genetic model and the chicken transcriptome, we assess whether sexual conflict can account for this difference. Sexual conflict over expression is inevitable when mutation effects are correlated across the sexes, as compensatory mutations in the heterogametic sex lead to hyperexpression in the homogametic sex. Coupled with stronger selection and greater reproductive variance in males, this results in slower and less complete evolution of Z compared with X dosage compensation. Using expression variance as a measure of selection strength, we find that, as predicted by the model, dosage compensation in the chicken is most pronounced in genes that are under strong selection biased towards females. Our study explains the pattern of weak dosage compensation in ZW systems, and suggests that sexual selection plays a major role in shaping sex chromosome dosage compensation.

Complete sex chromosome dosage compensation is largely limited to male heterogametic species, with the majority of female heterogametic species displaying incomplete dosage compensation. Here, the authors show that sexual conflict over gene expression combined with sexual selection in males can explain this pattern.

The evolution of sex chromosomes in organisms with separate haploid sexes

The evolution of dimorphic sex chromosomes is driven largely by the evolution of reduced recombination and the subsequent accumulation of deleterious mutations. Although these processes are increasingly well understood in diploid organisms, the evolution of dimorphic sex chromosomes in haploid organisms (U/V) has been virtually unstudied theoretically. We analyze a model to investigate the evolution of linkage between fitness loci and the sex-determining region in U/V species. In a second step, we test how prone nonrecombining regions are to degeneration due to accumulation of deleterious mutations. Our modeling predicts that the decay of recombination on the sex chromosomes and the addition of strata via fusions will be just as much a part of the evolution of haploid sex chromosomes as in diploid sex chromosome systems. Reduced recombination is broadly favored, as long as there is some fitness difference between haploid males and females. The degeneration of the sex-determining region due to the accumulation of deleterious mutations is expected to be slower in haploid organisms because of the absence of masking. Nevertheless, balancing selection often drives greater differentiation between the U/V sex chromosomes than in X/Y and Z/W systems. We summarize empirical evidence for haploid sex chromosome evolution and discuss our predictions in light of these findings.

Recurrent turnover of chromosome-specific satellites in Drosophila

Repetitive DNA are DNA sequences that are repeated multiple times in the genome and normally considered nonfunctional. Several studies predict that the rapid evolution of chromosome-specific satellites led to hybrid incompatibilities and speciation. Interestingly, in Drosophila, the X and dot chromosomes share a unique and noteworthy property: They are identified by chromosome-specific binding proteins and they are particularly involved in genetic incompatibilities between closely related species. Here, I show that the X and dot chromosomes are overpopulated by certain repetitive elements that undergo recurrent turnover in Drosophila species. The portion of the X and dot chromosomes covered by such satellites is up to 52 times and 44 times higher than in other chromosomes, respectively. In addition, the newly evolved X chromosome in D. pseudoobscura (the chromosomal arm XR) has been invaded by the same satellite that colonized the ancestral X chromosome (chromosomal arm XL), whereas the autosomal homologs in other species remain mostly devoid of satellites. Contrarily, the Müller element F in D. ananassae, homolog to the dot chromosome in D. melanogaster, has no overrepresented DNA sequences compared with any other chromosome. The biology and evolutionary patterns of the characterized satellites suggest that they provide both chromosomes with some kind of structural identity and are exposed to natural selection. The rapid satellite turnover fits some speciation models and may explain why these two chromosomes are typically involved in hybrid incompatibilities.

Conservation and de novo acquisition of dosage compensation on newly evolved sex chromosomes in Drosophila

Dosage compensation has arisen in response to the evolution of distinct male (XY) and female (XX) karyotypes. In Drosophila melanogaster, the MSL complex increases male X transcription approximately twofold. X-specific targeting is thought to occur through sequence-dependent binding to chromatin entry sites (CESs), followed by spreading in cis to active genes. We tested this model by asking how newly evolving sex chromosome arms in Drosophila miranda acquired dosage compensation. We found evidence for the creation of new CESs, with the analogous sequence and spacing as in D. melanogaster, providing strong support for the spreading model in the establishment of dosage compensation.

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.

The MOF chromobarrel domain controls genome-wide H4K16 acetylation and spreading of the MSL complex

The histone H4 lysine 16 (H4K16)-specific acetyltransferase MOF is part of two distinct complexes involved in X chromosome dosage compensation and autosomal transcription regulation. Here we show that the MOF chromobarrel domain is essential for H4K16 acetylation throughout the Drosophila genome and is required for spreading of the male-specific lethal (MSL) complex on the X chromosome. The MOF chromobarrel domain directly interacts with nucleic acids and potentiates MOF's enzymatic activity after chromatin binding, making it a unique example of a chromo-like domain directly controlling acetylation activity in vivo. We also show that the Drosophila-specific N terminus of MOF has evolved to perform sex-specific functions. It modulates nucleosome binding and HAT activity and controls MSL complex assembly, thus regulating MOF function in dosage compensation. We propose that MOF has been especially tailored to achieve tight regulation of its enzymatic activity and enable its dual role on X and autosomes.

msl2 mRNA is bound by free nuclear MSL complex in Drosophila melanogaster

In Drosophila, the global increase in transcription from the male X chromosome to compensate for its monosomy is mediated by the male-specific lethal (MSL) complex consisting of five proteins and two non-coding RNAs, roX1 and roX2. After an initial sequence-dependent recognition by the MSL complex of 150-300 high affinity sites, the spread to the majority of the X-linked genes depends on local MSL-complex concentration and active transcription. We have explored whether any additional RNA species are associated with the MSL complex. No additional roX RNA species were found, but a strong association was found between a spliced and poly-adenylated msl2 RNA and the MSL complex. Based on our results, we propose a model in which a non-chromatin-associated partial or complete MSL-complex titrates newly transcribed msl2 mRNA and thus regulates the amount of available MSL complex by feedback. This represents a novel mechanism in chromatin structure regulation.

Does positive selection drive transcription factor binding site turnover? A test with Drosophila cis-regulatory modules

Transcription factor binding site(s) (TFBS) gain and loss (i.e., turnover) is a well-documented feature of cis-regulatory module (CRM) evolution, yet little attention has been paid to the evolutionary force(s) driving this turnover process. The predominant view, motivated by its widespread occurrence, emphasizes the importance of compensatory mutation and genetic drift. Positive selection, in contrast, although it has been invoked in specific instances of adaptive gene expression evolution, has not been considered as a general alternative to neutral compensatory evolution. In this study we evaluate the two hypotheses by analyzing patterns of single nucleotide polymorphism in the TFBS of well-characterized CRM in two closely related Drosophila species, Drosophila melanogaster and Drosophila simulans. An important feature of the analysis is classification of TFBS mutations according to the direction of their predicted effect on binding affinity, which allows gains and losses to be evaluated independently along the two phylogenetic lineages. The observed patterns of polymorphism and divergence are not compatible with neutral evolution for either class of mutations. Instead, multiple lines of evidence are consistent with contributions of positive selection to TFBS gain and loss as well as purifying selection in its maintenance. In discussion, we propose a model to reconcile the finding of selection driving TFBS turnover with constrained CRM function over long evolutionary time.

Local effects of limited recombination: historical perspective and consequences for population estimates of adaptive evolution

Recent years have witnessed the integration of theoretical advances in population genetics with large-scale analyses of complete genomes, with a growing number of studies suggesting pervasive natural selection that includes frequent deleterious as well as adaptive mutations. In finite populations, however, mutations under selection alter the fate of genetically linked mutations (the so-called Hill-Robertson effect). Here we review the evolutionary consequences of selection at linked sites (linked selection) focusing on its effects on nearby nucleotides in genomic regions with nonreduced recombination. We argue that these local effects of linkage may account for differences in selection intensity among genes. We also show that even high levels of recombination are unlikely to remove all effects of linked selection, causing a reduction in the polymorphism to divergence ratio (r(pd)) at neutral sites. Because a number of methods employed to estimate the magnitude and frequency of adaptive mutations take reduced r(pd) as evidence of positive selection, ignoring local linkage effects may lead to misleading estimates of the proportion of adaptive substitutions and estimates of positive selection. These biases are caused by employing methods that do not account for local variation in the relative effective population size (N(e)) caused by linked selection.

Third chromosome candidate genes for conspecific sperm precedence between D. simulans and D. mauritiana

BACKGROUND

Male - female incompatibilities can be critical in keeping species as separate and discrete units. Premating incompatibilities and postzygotic hybrid sterility/inviability have been widely studied as isolating barriers between species. In recent years, a number of studies have brought attention to postmating prezygotic barriers arising from male - male competition and male - female interactions. Yet little is known about the genetic basis of postmating prezygotic isolation barriers between species.

RESULTS

Using D. simulans lines with mapped introgressions of D. mauritiana into their third chromosome, we find at least two D. mauritiana introgressions causing male breakdown in competitive paternity success. Eighty one genes within the mapped introgressed regions were identified as broad-sense candidates on the basis of male reproductive tract expression and male-related function. The list of candidates was narrowed down to five genes based on differences in male reproductive tract expression between D. simulans and D. mauritiana. Another ten genes were confirmed as candidates using evidence of adaptive gene coding sequence diversification in the D. simulans and/or D. mauritiana lineage. Our results show a complex genetic basis for conspecific sperm precedence, with evidence of gene interactions between at least two third chromosome loci. Pleiotropy is also evident from correlation between conspecific sperm precedence and female induced fecundity and the identification of candidate genes that might exert an effect through genetic conflict and immunity.

CONCLUSIONS

We identified at least two loci responsible for conspecific sperm precedence. A third of candidate genes within these two loci are located in the 89B cytogenetic position, highlighting a possible major role for this chromosome position during the evolution of species specific adaptations to postmating prezygotic reproductive challenges.

A sequence motif enriched in regions bound by the Drosophila dosage compensation complex

Background

In Drosophila melanogaster, dosage compensation is mediated by the action of the dosage compensation complex (DCC). How the DCC recognizes the fly X chromosome is still poorly understood. Characteristic sequence signatures at all DCC binding sites have not hitherto been found.

Results

In this study, we compare the known binding sites of the DCC with oligonucleotide profiles that measure the specificity of the sequences of the D. melanogaster X chromosome. We show that the X chromosome regions bound by the DCC are enriched for a particular type of short, repetitive sequences. Their distribution suggests that these sequences contribute to chromosome recognition, the generation of DCC binding sites and/or the local spreading of the complex. Comparative data indicate that the same sequences may be involved in dosage compensation in other Drosophila species.

Conclusions

These results offer an explanation for the wild-type binding of the DCC along the Drosophila X chromosome, contribute to delineate the forces leading to the establishment of dosage compensation and suggest new experimental approaches to understand the precise biochemical features of the dosage compensation system.

Progress and prospects toward our understanding of the evolution of dosage compensation

In many eukaryotic organisms, gender is determined by a pair of heteromorphic sex chromosomes. Degeneration of the non-recombining Y chromosome is a general facet of sex chromosome evolution. Selective pressure to restore expression levels of X-linked genes relative to autosomes accompanies Y-chromosome degeneration, thus driving the evolution of dosage compensation mechanisms. This review focuses on evolutionary aspects of dosage compensation, in light of recent advances in comparative and functional genomics that have substantially increased our understanding of the molecular mechanisms of dosage compensation and how it evolved. We review processes involved in sex chromosome evolution, and discuss the dynamic interaction between Y degeneration and the acquisition of dosage compensation. We compare mechanisms of dosage compensation and the origin of dosage compensation genes between different taxa and comment on sex chromosomes that apparently lack compensation mechanisms. Finally, we discuss how dosage compensation systems can also influence the evolution of well-established sex chromosomes.

Genetic testing of the hypothesis that hybrid male lethality results from a failure in dosage compensation

Several recent studies have suggested that F(1) hybrid male lethality in crosses between Drosophila melanogaster and D. simulans is due to a failure in dosage compensation, caused by incompatibilities between D. simulans dosage compensation proteins and the D. melanogaster X chromosome. Contrary to the predictions of this hypothesis, mutations in four essential D. melanogaster dosage compensation genes are shown here to moderately increase rather than decrease hybrid male viability.

A surrogate approach to study the evolution of noncoding DNA elements that organize eukaryotic genomes

Comparative genomics provides a facile way to address issues of evolutionary constraint acting on different elements of the genome. However, several important DNA elements have not reaped the benefits of this new approach. Some have proved intractable to current day sequencing technology. These include centromeric and heterochromatic DNA, which are essential for chromosome segregation as well as gene regulation, but the highly repetitive nature of the DNA sequences in these regions make them difficult to assemble into longer contigs. Other sequences, like dosage compensation X chromosomal sites, origins of DNA replication, or heterochromatic sequences that encode piwi-associated RNAs, have proved difficult to study because they do not have recognizable DNA features that allow them to be described functionally or computationally. We have employed an alternate approach to the direct study of these DNA elements. By using proteins that specifically bind these noncoding DNAs as surrogates, we can indirectly assay the evolutionary constraints acting on these important DNA elements. We review the impact that such "surrogate strategies" have had on our understanding of the evolutionary constraints shaping centromeres, origins of DNA replication, and dosage compensation X chromosomal sites. These have begun to reveal that in contrast to the view that such structural DNA elements are either highly constrained (under purifying selection) or free to drift (under neutral evolution), some of them may instead be shaped by adaptive evolution and genetic conflicts (these are not mutually exclusive). These insights also help to explain why the same elements (e.g., centromeres and replication origins), which are so complex in some eukaryotic genomes, can be simple and well defined in other where similar conflicts do not exist.

Drosophila dosage compensation: a complex voyage to the X chromosome

Dosage compensation is the crucial process that equalizes gene expression from the X chromosome between males (XY) and females (XX). In Drosophila, the male-specific lethal (MSL) ribonucleoprotein complex mediates dosage compensation by upregulating transcription from the single male X chromosome approximately twofold. A key challenge is to understand how the MSL complex distinguishes the X chromosome from autosomes. Recent studies suggest that this occurs through a multi-step targeting mechanism that involves DNA sequence elements and epigenetic marks associated with transcription. This review will discuss the relative contributions of sequence elements and transcriptional marks to the complete pattern of MSL complex binding.

Evolution of the Drosophila nuclear pore complex results in multiple hybrid incompatibilities

Speciation often involves the evolution of incompatible gene interactions that cause sterility or lethality in hybrids between populations. These so-called hybrid incompatibilities occur between two or more functionally divergent loci. We show that the nucleoporin 160kDa (Nup160) gene of the fruitfly Drosophila simulans is incompatible with one or more factors on the D. melanogaster X chromosome, causing hybrid lethality. Nup160 encodes a nuclear pore complex protein and shows evidence of adaptive evolution. Furthermore, the protein encoded by Nup160 directly interacts with that of another hybrid lethality gene, Nup96, indicating that at least two lethal hybrid incompatibility genes have evolved as byproducts of divergent coevolution among interacting components of the Drosophila nuclear pore complex.

Comparative analysis of five immunity-related genes reveals different levels of adaptive evolution in the virilis and melanogaster groups of Drosophila

Five immunity-related genes previously reported to be evolving under positive selection in Drosophila melanogaster and D. simulans have been analysed across the Drosophila genus using two types of approaches, random-site and branch-site likelihood models as well as the proportion of synonymous and non-synonymous variation within and between species. Different selective pressures have been detected in the sample of genes, one showing evidence for adaptive evolution across the phylogeny of Drosophila and two showing lineage-specific positive selection. Furthermore, amino-acid sites identified as being under positive selection in the melanogaster and the virilis groups are different, suggesting that the evolution of the proteins in these two divergent groups may have been shaped by different pathogens.