Stage-Specific Expression Profiling of Spermatogenesis Suggests that Meiotic Sex Chromosome Inactivation Drives Genomic Relocation of Testis-Expressed Genes
In Drosophila, genes expressed in males tend to accumulate on autosomes and are underrepresented on the X chromosome. In particular, genes expressed in testis have been observed to frequently relocate from the X chromosome to the autosomes. The inactivation of X-linked genes during male meiosis (i.e., meiotic sex chromosome inactivation—MSCI) was first proposed to explain male sterility caused by X-autosomal translocation in Drosophila, and more recently it was suggested that MSCI might provide the conditions under which selection would favor the accumulation of testis-expressed genes on autosomes. In order to investigate the impact of MSCI on Drosophila testis-expressed genes, we performed a global gene expression analysis of the three major phases of D. melanogaster spermatogenesis: mitosis, meiosis, and post-meiosis. First, we found evidence supporting the existence of MSCI by comparing the expression levels of X- and autosome-linked genes, finding the former to be significantly reduced in meiosis. Second, we observed that the paucity of X-linked testis-expressed genes was restricted to those genes highly expressed in meiosis. Third, we found that autosomal genes relocated through retroposition from the X chromosome were more often highly expressed in meiosis in contrast to their X-linked parents. These results suggest MSCI as a general mechanism affecting the evolution of some testis-expressed genes.
Published in the journal:
. PLoS Genet 5(11): e32767. doi:10.1371/journal.pgen.1000731
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1000731
Summary
In Drosophila, genes expressed in males tend to accumulate on autosomes and are underrepresented on the X chromosome. In particular, genes expressed in testis have been observed to frequently relocate from the X chromosome to the autosomes. The inactivation of X-linked genes during male meiosis (i.e., meiotic sex chromosome inactivation—MSCI) was first proposed to explain male sterility caused by X-autosomal translocation in Drosophila, and more recently it was suggested that MSCI might provide the conditions under which selection would favor the accumulation of testis-expressed genes on autosomes. In order to investigate the impact of MSCI on Drosophila testis-expressed genes, we performed a global gene expression analysis of the three major phases of D. melanogaster spermatogenesis: mitosis, meiosis, and post-meiosis. First, we found evidence supporting the existence of MSCI by comparing the expression levels of X- and autosome-linked genes, finding the former to be significantly reduced in meiosis. Second, we observed that the paucity of X-linked testis-expressed genes was restricted to those genes highly expressed in meiosis. Third, we found that autosomal genes relocated through retroposition from the X chromosome were more often highly expressed in meiosis in contrast to their X-linked parents. These results suggest MSCI as a general mechanism affecting the evolution of some testis-expressed genes.
Introduction
Sex chromosomes evolve differently than autosomes due to their distinct characteristics such as a lack of recombination between the homologues (X and Y), a different number of chromosome copies between sexes, and the proportion of heterochromatin [1]. Genomic and proteomic studies in Drosophila, mammals and worms reveal that male-biased genes, i.e. those that are more highly expressed in males than in females, are underrepresented in the X chromosome [2]–[7]. One possible mechanism contributing to this biased chromosomal distribution is the directional movement of testis-expressed genes out of the X chromosome. In Drosophila and mammals, there is a significant excess of genes retroposed from the X chromosome to the autosomes, and these genes are more likely to be expressed in testis [8]–[11]. Recently, in Drosophila, this observation was found to also apply for DNA-based gene duplication [12],[13].
Two evolutionary hypotheses have been proposed to explain the underrepresentation of genes expressed in males on the X chromosome. First, the meiotic sex chromosome inactivation (MSCI) hypothesis, can explain both the excess of retroposed genes from the X chromosome to the autosomes as well as the paucity of X-linked testis-expressed genes [9],[14]. According to this hypothesis, the inactivation of X-linked genes during male meiosis favors the accumulation of testis-expressed genes in autosomes where such genes can be expressed during the meiotic stage [14].
A second hypothesis, based on a model proposed by Rice [15] and by Charlesworth and co-authors [16], states that sex biased expression is driven by sexually antagonistic forces, i.e., the case of opposing selection pressure on the two sexes. In this model developed through the investigation of polymorphic equilibrium, dominant mutations with beneficial fitness effects in females, but detrimental effects in males, will have a higher probability of being fixed on the X chromosome [15], while recessive sexually antagonistic mutations will have a higher probability of being fixed on the autosomes [16]. The opposite pattern is expected for male-beneficial, female-detrimental mutations: dominant and recessive alleles will have higher probability of being fixed on the autosomes and on the X chromosome, respectively [15]. Recently, however, Patten and Haig [17] showed that male-beneficial mutations could also be fixed on the X chromosome even for some cases of dominant alleles. In addition, is important to note that Gibson and colleagues' theoretical and empirical findings [18] also suggest the X chromosome to be enriched in polymorphism for sexually antagonistic alleles.
In order to relate the sexual antagonism hypothesis to the male-biased gene chromosomal distribution, one must assume a model under which sexually antagonistic fitness variation contributes to sex-biased expression. In the model proposed by Rice [15], the paucity of X-linked male-biased genes would be a consequence of restrictive conditions for sexually antagonistic mutations. For instance, the model requires that: i) the majority of X-linked sexually antagonistic mutations are female-beneficial/male-detrimental; or ii) most male-beneficial/female-detrimental alleles are dominant [for a review, see [19]]. Although none of these requirements have been tested yet, the sexual antagonism hypothesis could explain the underrepresentation of male-biased genes on the Drosophila X chromosome [3],[20]. Most important, the observation of similar X chromosome underrepresentation found for somatic cells revealed that the chromosomal distribution of male-biased genes cannot be exclusively explained by MSCI [3],[20].
MSCI appears to be phylogenetically plastic: it has been observed in the X chromosome of male mammals and nematodes [21],[22], and in the Z chromosome of female birds [23]. Meiotic X chromosome inactivation in Drosophila males was first suggested as the cause of male sterility in mutants carrying X-autosomal translocations [14]. Recently, MSCI was investigated by inserting transgenic constructs expressing a testis-specific promoter, into different regions of the genome [24]. A reduction in the expression of reporter gene insertions into the X chromosome was observed, consistent with a MSCI model for Drosophila. In contrast to this result, however, a microarray-based study of male germline expression using spermatogenic arrest mutants did not detect a significant signal of MSCI [20].
In order to investigate the impact of MSCI on Drosophila testis-expressed genes, we conducted a global gene expression analysis of spermatogenesis in a developmental context. We compared X and autosome expression in three distinct regions of the testis, (the apical, mid- and posterior testis), that are correspondingly enriched with mitotic, meiotic and post-meiotic cells. It is important to note that in mammals, different studies [6],[25] have shown variability in gene expression reduction as a result of MSCI (e.g., a 60% to 80% X-linked expression reduction was observed in microarray-based experiments, [25]). In order to avoid an absolute definition of completely silenced expression, here we defined MSCI as a significant expression reduction in the male germline X chromosome.
Our results clearly show a small, but significant, reduction of X-linked gene expression in meiosis, consistent with MSCI. We also investigated the impact of MSCI on the chromosomal distribution of testis-expressed genes by analyzing their expression profile during spermatogenesis. MSCI was positively correlated with both the retrogene chromosomal movement and the underrepresentation of X-linked testis-biased genes. Our findings suggest MSCI as a general mechanism affecting the evolution of some testis-expressed genes.
Results/Discussion
Global gene expression analysis during spermatogenesis
To measure MSCI and investigate the impact of X inactivation on testis-expressed genes, we conducted a genome-wide expression analysis of wild-type males using three cell populations isolated from mitotic, meiotic and post-meiotic phases of spermatogenesis in Drosophila melanogaster, respectively (Materials and Methods, Table S1). Figure 1 details the precise regions of the testis that were isolated in this study. According to MSCI in mammals, we expected to find in Drosophila differential X-linked expression specific to spermatogenic meiotic cells [14],[20],[21],[24]. Therefore, our method allowed us to effectively measure the relative expression levels among the three phases and confirm MSCI.
Purification of Drosophila spermatogenic cells is not as straightforward as it is in mammals, where purified samples of specific germline phases can be isolated [6]. Therefore, our approach was to directly isolate testis regions enriched with RNAs from each of the three specific germline phases (Figure 1, Materials and Methods). While the cell types present at various stages of spermatogenesis are generally located in a gradient along the proximal-distal axis of the testis, most are not exclusively restricted to any one geographic region. Therefore, we used the cellular morphology of the three phases (Figure 1) as a guide to isolate cells enriched for the following cell populations: i) apical (mitotic) – enriched for apical hub stem cells, spermatogonial cysts with reduced levels of primary spermatocyte cysts; ii) mid-region (meiotic) – enriched for primary and secondary spermatocyte cysts; reduced levels of early spermatids; iii) distal (post-meiotic) – enriched for elongating spermatid cysts, individualized spermatid bundles and coiled spermatozoa, reduced levels of spermatocyte cyst cells. Careful dissection of these landmark regions produced exceptional reproducibility of the data as confirmed by the high correlation within replicates from the same spermatogenic phases (r 0.993, Figure S2).
In order to evaluate the cell content of mitotic, meiotic and post-meiotic samples, we examined our dataset for individual genes known to be expressed in specific spermatogenic stages. For instance, genes known to be involved in early spermatogenic processes are highly expressed in the mitotic phase, followed by decreased expression in meiosis and post-meiosis (Figure 2). For example, bag-of-marbles and benign gonial cell neoplasm genes act early in the Drosophila germ cells and their mutations prevent progression through spermatogenesis and cell differentiation [26]. Both genes were over-expressed in our mitotic sample (Figure 2), which, in overall, showed enrichment for genes involved in early spermatogenic processes (Fisher exact test, P 0.014). On the other hand, genes expressed in spermatocytes, such as K81, showed a peak of expression in meiosis (Figure 2). K81 is a paternal effect gene showing a post-fertilization phenotype but known to be exclusively expressed in primary spermatocytes [27]. Additionally, genes such as schuy and CG31858 showed low levels of mitotic expression followed by an increase in expression during meiosis, but peaking in the post-meiotic phase (Figure 2). Using in situ hybridization, a recent study showed that these genes are transcribed in post-meiosis [28]. These genes have been categorized and named as comet and cup genes according to their transcript localization patterns, since both are found in the end of elongated spermatids [28].
In summary, our samples were enriched with stage-specific cells and gene expression profiles can be accessed via the public database SpermPress (http://pondside.uchicago.edu/~longlab/spermpress). While our samples do not perfectly separate the three phases of spermatogenesis, our cell enrichment technique and the high degree of reproducibility in the replicate datasets should be enough to detect the reduction of X-linked expression during meiosis, as expected under MSCI.
MSCI in D. melanogaster
The existence of MSCI in Drosophila was recently investigated using reporter gene constructs whose expression was driven by a testis specific promoter [24]. In this study, the authors observed a significant reduction in the expression of constructs inserted into the X chromosome when compared to those inserted into the autosomes, and proposed this was consistent with the existence of MSCI. In contrast, a recent microarray survey of male germline expression using spermatogenic arrest mutants found no evidence for MSCI [20]. However, the mutants used in this study were completely depleted of later meiotic cells, compromising the overall investigation of meiotic MSCI [20].
We decided to directly confirm meiotic effects of MSCI in D. melanogaster by comparing the expression profile of the X chromosome and autosomes during the three stages of spermatogenesis. MSCI leads to a clear expression prediction: there should be a significant decrease in the expression of X-linked genes when compared to autosomal-linked genes exclusively during the meiotic phase.
Conventional microarray analyses based on average expression intensities, such as the one previously used to investigate the existence of X inactivation [20],[29] lose important information (such as variation and distribution of the expression signals), potentially compromising statistical hypothesis testing. Such limitations may be particularly important regarding the X chromosome expression reduction expected during MSCI [20]. There are two main statistical concerns we aimed to address: 1). how to deal with the variation unavoidably introduced by factors such as expression fluctuation among repeated experiments; and 2). how to deal with the small proportion of non-germline cells mixed in our samples. To approach these two issues, we introduced Bayesian statistical models for analyzing our data and testing the hypothesis of MSCI. This approach permitted the use of the data in its entirety to simultaneously estimate proportions of under-, over- and equally expressed genes in meiosis, while comparing X chromosome and autosomal distributions (Model A in Statistical Analyses, Methods in Text S1; Figure 3) [30]. The test avoided the specification of arbitrary thresholds commonly employed such as fold change methods [3],[31]. It also eliminated errors introduced by multiple hypothesis testing [32].
Our Bayesian approach modeled the proportions of genes from the X chromosome and autosomes that were over-, under- and equally expressed in meiosis relative to other spermatogenic phases. X inactivation could thus be observed as an excessive number of X-linked genes under-expressed in meiosis, i.e., with expression reduction in meiosis. On the other hand, if the X chromosome is active, we would expect to observe no difference between autosomal and X chromosome proportions of over-, under- and equally expressed genes.
We restricted our analyses to those genes that had been previously shown to be expressed in whole testis microarray experiments (Flyatlas presence call, Table S1). In our genome-wide phase analysis, we detected a significant under-expression of X-linked genes compared to autosome-linked genes in meiosis relative to mitosis (Figure 3A and 3B). In the meiotic phase, the proportion of over-expressed X-linked genes was significantly reduced when compared to the autosome-linked genes (Figure 3A; Bayesian P≤0.001), whereas the proportion of under-expressed X-linked genes was increased, as indicated by the distribution of X-linked and autosome-linked genes (Bayesian P≤0.001) (Figure 3B). The difference between the proportion of autosome-linked and X –linked genes under-expressed in meiosis (Figure 3B) is in agreement with the expression reduction of the X chromosome in meiosis expected by our previous definition of MSCI (Figure 3C). Similar results were obtained when using the complete dataset of genes, instead of only those already known to be expressed in testis (Figure S3C).
The clear separation of chromosome proportions for each class detected by our Bayesian model contrasts with the results obtained using an arbitrary twofold change threshold in the same data (Figure S3D, Methods in Text S1). This is a result of the increased power of the Bayesian approach for detecting genome-wide expression differences and of the negligible within-gene variability compared to between gene-variability produced by our experimental replicates.
Although we observed a significant reduction in meiotic expression on the X chromosome, it is clear that such reduction does not meet the textbook definition of “X inactivation”, that is, the complete silencing of expression. Indeed, we observed a fairly small (10%) average X-linked gene meiotic expression reduction when compared to genes in the autosomes, which does not indicate whole-scale inactivation of the X (94% vs. 84% for autosomal- and X-linked genes, respectively).
Moreover, the reduction observed in our experiments is considerably different (shows lower reduction magnitude) from the one obtained by Hense and co-authors using transgenic construct insertions [24]. The expression of transgenes inserted into the X observed by Hense et al., however, was universally reduced (between 3.4 to 10-fold) compared to transgenes inserted into the autosomes [24]. There are differences in the experimental properties that could account for the discrepancies obtained by the two experiments. First, Hense and co-authors analyzed construct insertions for the promoter of a gene highly and specifically expressed in testis [24], whereas our analysis includes the expression profile of all D. melanogaster gene products. This difference per se could generate discrepant observations of expression reduction. Indeed, Hense and co-authors [24] commented that their system may be especially sensitive for detecting X-inactivation because they used a minimal promoter from a gene that appears to be expressed relatively late during spermatogenesis. Second, the two studies used different methods to measure RNA intensity. Hense et al. measured RNA level by -Galactosidase activity and qRT-PCR, whereas we used a microarray technique [24]. These methods can present differences in resolution and expression variation. For instance, Hense and co-authors observed a lot of variation in expression reduction of the X chromosome, depending on the method used: -Galactosidase activity showed a 10-fold difference, whereas qRT-PCR showed 3.4-fold difference [24]. In our experiments, genes found to be under-expressed in meiosis showed an average two-fold expression reduction in meiosis compared to mitosis (Figure S4). However, individual gene expression reduction varied from 1.068 to 17.70 fold. Third, the heterogeneous cell makeup in our experiments could have obscured the signal on the microarray. Our samples were enriched with cells from different stages and therefore do not represent pure cell types. Other possible discrepancies also could account for our different results (e.g. transgene insertions can be affected more strongly by different chromatin structures, or X chromosome expression reduction can arise during later phases of meiosis).
Our global gene expression analysis, however, clearly showed that not all X-linked genes were affected in the same way, i.e., not all X-linked genes were under-expressed in meiosis. The fact that our samples were enriched with cells from different spermatogenesis stages rather than corresponding to pure samples could explain this result. On the other hand, MSCI in Drosophila might not be globally homogeneous, and therefore might be similar to what occurs in mammals [33].
Despite the discrepancies between our results and those of Hense et al. [24], it is important to note the similarities between them. Both showed a significant X chromosome expression reduction in male meiosis. The two studies used different yet complementary approaches and methods but reached similar conclusions. Our results form an independent and direct genomic test for the hypothesis of MSCI, which was indirectly supported by Hense and co-authors [24]. Moreover, the large-scale gene expression data of male meiosis provided by our study thus allowed us to investigate the evolutionary effects of meiotic expression reduction on X-linked testis-expressed genes, such as retrogene chromosomal relocations (see below).
The X chromosome has a paucity of testis male-biased genes (most likely male-biased germline cells) [3]–[5], and could be consequently enriched with non male-biased testis somatic genes. Therefore, one could argue that a lower proportion of somatic tissue in our meiotic sample (for example, absence of testis sheath), when compared to the mitotic sample, could account for the lower expression of X-linked genes. There are at least two arguments against this interpretation. First, most of the transcription in the testis occurs in germline cells and therefore somatic expression is unlikely to make a significant contribution [34]. Second, the paucity of male-biased X-linked genes was also observed for somatic tissues (in gonadectomized males) and therefore somatic testis genes probably show the same expression pattern [3]. Therefore, in all likelihood, differences in the amount of somatic and germline cells in our samples, did not make a major contribution to chromosomal expression differences.
Nonetheless, we decided to further investigate if the difference we observed in chromosomal expression was due to a somatic bias in our samples. We compared X chromosomal and autosomal expression exclusively for testis-biased genes selected from previous gonad comparison experiments [3],[4]. Using this approach, we were able to exclude the contribution of somatic genes. In complete agreement with our previous observations, we detected an excess of X-linked genes under-expressed in meiosis, thus providing supporting evidence for MSCI (Figure 4).
An alternative hypothesis to MSCI for our expression results is that the lowered expression of X chromosome genes in the meiotic phase is due to a loss of dosage compensation, which in Drosophila is achieved by hypertranscription of the X chromosome in males [35]. Previous analyses comparing ovary and testis expression have already established the existence of dosage compensation in the male germline [35]. Therefore, taking into account our results, it is unclear why dosage compensation should be limited to a period of the entire spermatogenesis process. The lower level of X chromosome expression restricted to the particular meiotic stage that we observed, however, is consistent with male germline X inactivation. Another piece of evidence against loss of dosage compensation related to lower expression in X-linked testis expressed genes was obtained by MSCI experiments done by Hense and co-authors [24]. In their work, the expression levels of transgenic constructs driven by a testis promoter inserted into the X were reduced compared to the same constructs inserted in only one of the homologous autosomal chromosomes (heterozygous insertion). If the phenomenon was due to loss of dosage compensation instead of MSCI, X-linked transgenic constructs should have been expressed at least at the same level as autosomal insertions, which was not observed [24].
At first glance, our global expression experiments and results are not able to empirically differentiate between MSCI and an ad hoc hypothesis of stage-specific loss of dosage compensation. If in fact there were a loss of dosage compensation occurring in meiosis, the results presented in the following sections would not change significantly. X chromosome down regulation in meiotic cells is associated with and impacts the general distribution of testis- biased genes and retrogenes.
Although we observed a significant reduction in expression for the X-linked loci in the meiotic sample when compared to the post-meiotic phase (Figure S5), we had hesitated to adopt an MSCI-based interpretation of this result. One reason for caution is that autoradiographic studies of Drosophila spermatogenesis have shown the absence of transcription in the post-meiotic phase [e.g., [34]]. The prevailing notion has been that most proteins required in the late stages of spermatogenesis are translated from mRNAs produced during early meiosis and stored in the cytoplasm [34],[36]. However, a recent study using in situ hybridization provided compelling evidence of post-meiotic transcription for 24 selected Drosophila genes [28]. The ambiguity associated between post-meiotic translation and transcription could confound any observations regarding the prevalence of X chromosome inactivation after meiosis.
MSCI and retrogenes movement out of the X chromosome in Drosophila
MSCI has been proposed to explain the accumulation of testis-expressed genes in autosomes [3],[9],[14]. Two specific predictions result from this hypothesis, one related to retrogenes, and the other related to testis-biased genes (see next section for the latter). The retrogene-based prediction refers to the disproportionate retroposition of X-linked genes to the autosomes, as a consequence of natural selection for expression of testis functions during MSCI [8],[9]. This prediction generates two testable expectations: (i) most autosomal retroposed genes that originated from X-linked parental genes should be over-expressed during the MSCI stage (i.e. meiosis); and (ii) retroposed genes should display complementary expression patterns during meiosis in relation to their X-linked parental genes. With our data, we have tested and confirmed both expectations.
We observed that approximately 63% of all retroposed genes were over-expressed during meiosis in comparison to mitosis, where X inactivation occurs. As observed for mammals retrogenes [37], both autosomal copies retroposed from the X and copies retroposed from other autosomes are more often over-expressed in meiosis than non-retrogene autosomal genes (32%, n = 7827). This indicates that any autosomal retrogene is more likely to be expressed in meiosis. One possible explanation is that the “hypertranscription” state of autosomal chromatin in meiotic cells might facilitate the expression of young retrocopies soon after their emergence [37]. However, confirming the first expectation, autosomal copies retroposed from the X chromosome were more often over-expressed in meiosis than genes retroposed between autosomes (gene proportions: 71% versus 58%, Bayesian P≤0.002). To test the second expectation, we measured the complementary expression of parental-retrogene pairs in meiosis, i.e. parental gene under-expression and retrogene over-expression in meiosis (Figure 5A and 5B). Specifically, we compared the expression of 27 X→A and 52 A→A parental-retrogene pairs [10] in the mitotic and meiotic phases (Figure 6, Table S2). An extension of the previous Bayesian analysis (Model B in Statistical Analyses, Methods in Text S1) revealed that, compared to A→A retrogenes, X→A retrogenes had a significantly higher proportion of complementary expression in meiosis, thus confirming the second expectation (Figure 5C, Bayesian P≤0.001).
The complementary expression pattern strongly suggested that MSCI occurs specifically during the meiotic phase, and revealed a spatial and temporal link between over-expression of the retrogene and inactivation of the parental gene on the X chromosome during meiosis. In addition, the observation of higher meiotic expression of genes retroposed from the X chromosome showed that the complementary expression pattern does not solely result from the inability of parental X-linked genes to evolve higher levels of meiotic expression. The most parsimonious scenario is that autosomal retrogenes are favored because they recover the function of an X-linked parental downregulated copy.
As mentioned in the previous section, our global transcriptome analysis used samples enriched for three phases of spermatogenesis, but also included accompanying somatic cell types, which could bias the expression profiles obtained. Additional evidence against a bias created by non-germline cells comes from the specific effects of MSCI on retrogenes. Reduced somatic contribution in our meiotic sample (e.g., absence of testis sheath) could potentially account for the excess of over-expressed autosomal genes. However, our analysis showed that two different groups of autosomal genes had significantly different expression patterns in meiosis. Autosomal genes retroposed from the X chromosome showed a significantly higher meiotic expression than those retroposed from the autosomes (71% versus 58%, Bayesian P 0.002). Such a difference would not be expected if the somatic contribution in our samples was a major factor or the expression profile is an artifact in our analysis. MSCI impacting testis-expressed genes is a more plausible explanation.
MSCI and the chromosomal distribution of Drosophila testis-biased genes
The existence of MSCI leads to another expectation: male-biased genes (those more highly expressed in males than in females) should accumulate in autosomes and, consequently, be underrepresented in the X chromosome [2]–[5]. There are no studies showing an effect of MSCI on the chromosomal distribution of male-biased genes. On the contrary, previous studies have shown that male-biased genes expressed not only in testis, but also in somatic cells are underrepresented in the X chromosome, which would suggest an explanation beyond MSCI (e.g., sexual antagonism or population genetic forces) [3],[20]. However, it is important to note that, in these studies, the proportion of somatic male-biased genes is an order of magnitude lower than the proportion of testis male-biased genes (<2% vs. <20%, Figure 2 in [3]), which would still allow an important role for MSCI.
In our expression comparison between meiotic and mitotic phases, approximately 60% of testis-biased genes (those more highly expressed in testes than in ovaries) were over-expressed in the meiotic phase (Figure 7C, Table S3). In order to test for the effect of MSCI on the chromosomal distribution of testis-biased genes, we analyzed the effects of mitosis and meiosis on the chromosomal distribution of testis-biased genes (Figure 7). The underrepresentation on the X chromosome of testis-biased genes (Figure 7A) is only observed for genes over-expressed in meiosis (Figure 7C) (Fisher exact test, P≤0.001). Testis-biased genes highly expressed in mitosis were not depleted in the X chromosome (Figure 7B). Our results show that both X chromosome inactivation and X-linked testis-biased gene underrepresentation occur only in meiosis. Therefore, it seems that the paucity of testis-biased genes on the X chromosome is affected by MSCI.
Our testis-biased gene analysis reveals a different chromosomal distribution of mitotic genes than the one observed using Drosophila spermatogenic arrest mutants [20]. In that study, genes expressed in mutant testes, enriched with mitotic cells, were found to be underrepresented in the X chromosome [20]. An important difference in our approach as opposed to the mutant analysis [20] was our ability to directly detect and compare gene expression levels (both higher and lower) in mitosis versus meiosis, instead of measuring gene expression levels only in mitotic cells in the mutant testis. Our analysis shows that there is not a paucity of X-linked genes that are more highly expressed in mitosis, therefore indicating the effects of MSCI. Another feature of our direct approach of global expression analysis of wild-type testis was that it obviated problems associated with the use of mutants, such as concomitant pleiotropic and physiological effects that may affect gene expression patterns.
Conclusion
In summary, our study coupled a global spermatogenic analysis with a Bayesian statistical method that overcame the limitations of conventional microarray comparison based on average expression intensities. We were able to detect significant expression reduction of X-linked genes during male meiotic phases. Although such reduction does not meet the absolute definition of gene silencing that may be unrealistic [25],[33], our results are in agreement with the detection of meiotic sex chromosome inactivation in males given the experimental methods currently available.
Our analysis also revealed a significant correlation between MSCI and retrogene chromosomal movement. Further, testis-biased genes are over-expressed during MSCI and the X-linked reduction is associated with the chromosomal distribution of those genes. These results suggest a critical impact of MSCI on the evolution of sex chromosomes in Drosophila. For example, the complementary nature of parental/retroposed gene pair expression during meiosis (Figure 5) illustrates the effects of MSCI on the origin, evolution and chromosomal localization of new testis-expressed retrogenes.
Previous studies using Drosophila ancestral X chromosome [9],[12],[13 and this work] and D. pseudoobscura neo-X chromosome [12],[13],[20] suggested that movement off of the X chromosome had significant role in reducing the proportion of male-biased genes in the X. Neo-X chromosome analysis also showed that other mechanisms such as gene gain and loss contributed to the paucity of X-linked male-biased genes [20]. The loss mechanism, for instance, could be obtained by intermediate pseudogenization of X-linked copies after gene duplication. In contrast, switch of expression profile between sexes (a male-biased gene turn to a female-biased gene) seems to be uncommon during the Drosophila gene evolution [20].
The impact of MSCI on the chromosomal distribution of male-biased genes in Drosophila can also be compared to the effect of MSCI in other organisms. In mammals, meiotic X inactivation affects the chromosomal distribution of testis-biased genes [37],[38]. A global gene expression analysis in different spermatogenic phases demonstrated that genes expressed during mitotic phases are enriched on the X chromosome, whereas those expressed later in spermatogenesis are enriched on autosomes [38]. Another study [37] demonstrated, through expression analyses along different spermatogenic stages that autosomal retrogenes specifically compensate for their X-linked parental gene that is silenced during meiosis.
Furthermore, although Drosophila and mammalian sex-chromosomes originated independently [39], there are similarities and differences among the forces shaping male gene evolution in the two systems [40]. In the case of sexual antagonism, the chromosomal distribution of male-biased genes in mammals could reflect the role of recessive alleles, since male genes expressed in mitosis tend to accumulate in the X chromosome [38]. In Drosophila, however, antagonistic forces acting on dominant mutations may contribute to demasculinization of the X chromosome, as revealed by somatic male-biased gene analyses [3],[20]. Note, however, that there many significant differences between flies and mammals (e.g., mechanisms of dosage compensation) [40]. As such, the statement that dominant/recessive sexual antagonistic alleles govern the chromosomal distribution in genes expressed in somatic cells could be an oversimplification. For genes expressed in meiotic cells, however, MSCI seems to play a similar role in both mammals and flies [37],[38 and this work] by driving the evolution of retrogenes and testis-biased genes. In summary, retrogene expression compensates parental X-linked inactivation and only late spermatogenic genes are depleted on the X. Therefore, taking mammalian and Drosophila data together, MSCI may be considered as a general mechanism and force impacting the evolution of testis-expressed genes.
Material and Methods
Isolation of tissue samples
All experiments used a wild-type strain of D. melanogaster originally collected near Tempe, Arizona (Wolbachia free strain [41]). Cells enriched for mitotic, meiotic and post-meiotic phases were obtained by dissection of apical, proximal and distal regions of the testis, respectively (Figure 1). Testes without seminal vesicles were dissected in PBS. Paired testes were separated using 0.25 mm diameter insect pins and only a single region was dissected from individual testis (Figure 1). This greatly increased the number of necessary dissections, but helped minimize contamination. Apical cells were obtained by separation of apical tips, whereas distal cells were obtained from the detachment of basal regions in order to avoid contamination with spermatocytes (Figure 1). Proximal cells were obtained as follows. First, individual testes placed in small drops had their distal region removed to allow the release of spermatid bundles, reducing the internal turgor in the testis. Then, the apical regions were removed and the spermatocytes were obtained by applying gradual pressure to the middle regions (proximal to the apical tip) in a posterior-anterior direction. During the procedure, contamination from distal regions was avoided by perpendicular positioning of the insect pin to prevent the exit of any remaining spermatid bundles as spermatocytes were teased out of the testis shealth. Cells from the different dissected regions were carefully pipetted to microcentrifuge tubes. For each of three replicate experiments, 250–500 testis dissections were used for RNA isolation. (For more detailed graphic protocols of dissecting testis expression data, see the database we prepared entitled “SpermPress: Drosophila spermatogenesis database” on the website (http:// pondside.uchicago. edu/∼longlab/spermpress).
Immunostaining techniques
Apical, proximal and distal regions of the testis were fixed in 95% ethanol followed by formaldehyde as previously described [42]. Indirect immunofluorescence staining was carried out using mouse α -tubulin (Sigma) and Cy3 donkey anti-mouse IgG (Jackson Laboratories) as primary and secondary antibodies, respectively. DNA was stained with Sytox Green (Amersham). Cells were incubated for 1 h at room temperature in solution containing 1%BSA, 1∶1500 and 1∶500 dilutions for primary and secondary antibodies, respectively. 2 mg/ml RNAse A was included in the primary antibody solution. Fixation and antibody stains were always followed by three wash steps in PBS-T. DNA stained in Sytox green (10 µM solution for 10 min at room temperature) was briefly washed before mounting on slides for visualization with a confocal microscope using Pro-long anti-fade media (Molecular Probes).
Arrays
Total RNA was extracted from the apical, proximal and distal germ cell populations using PicoPure™ RNA Isolation Kit (Arcturus). Three biological replicates were hybridized into Affymetrix Gene Chip Drosophila Genome 2.0 Arrays. cDNAs were synthesized according to the Invitrogen SuperScript RNA Amplification System, except for a modified in vitro transcription step where Biotin-Labeled cRNA was produced. Hybridization, scan and data processing were done using Affymetrix default protocols. Gene product expression was measured by hybridization intensity (log2) obtained using RMA background correction and quantile normalization (Bioconductor package in R). Full platform descriptions and data are available at the GEO under accession GSE18502. Individual gene expression profile can also be obtained at SpermPress database (http://pondside.uchicago.edu/~longlab/spermpress).
Selection of genes expressed in testis (FlyAtlas)
Confirmation of MSCI was done using only genes already known to be expressed in testis (Figure 3). We selected approximately nine thousand gene products present in at least 3 out of 4 Drosophila testis microarrays (Presence Call in Flyatlas [43], Table S1).
Testis-biased genes selection
Testis-biased genes were selected from a testis vs. ovary comparison contained in the Sebida database [44]. We selected 2608 genes that were classified as testis-biased genes in Drosophila gonad comparison datasets [3],[4] (Table S3). Testis-biased genes were assigned to 2268 and 331 gene product probe identifications located in the autosomes and in the X chromosome, respectively.
Parental-retrogene selection
Parental-retrogene pairs were selected from Bai and co-authors [10]. We were able to assign 91 pairs of probe IDs that corresponded to parental-retrogene pairs, excluding 3 pairs of duplication events after retrotransposition [10] (Table S2). We classified parental-retrogene pairs into retroposed “X→A” (n = 27), “A→A” (n = 52), “X→X” (n = 2) and “A→X” (n = 10) groups.
Statistical analyses
MSCI confirmation. In order to confirm X inactivation, we analyzed expression separately for X-linked and autosomal-linked genes. First, we compared meiotic gene expression to mitotic and to post-meiotic expressions, respectively. We could therefore classify genes as having higher (over), lower (under) or equal expression levels in meiosis compared to the other spermatogenic phases (mitosis or post-meiosis). Simultaneously, the proportions of genes in each class (over-, under- and equal expression) were estimated for X- and autosomal-linked genes. MSCI was detected by observing an excessive number of X-linked genes under-expressed in meiosis relative to any other phase. Excessive number in this case means a significantly higher proportion compared to autosomal-linked genes. In other words, we hypothesized that the X chromosome would possess more genes under-expressed in meiosis as compared to autosomes. Bayesian models were developed to estimate both chromosomal distributions as well as proportions of genes in expression classes [32],[45],[46]. More details about Bayesian estimation (Figure S1) and classification can be found in Model A in Statistical Analyses, Methods in Text S1.
Retrogene analysis
In order to verify if MSCI affects retrogene movement, we compared the spermatogenic expression of X->A and A->A parental-retrogene pairs. We expected that X->A pairs would have a higher proportion of complementary expression than A->A pairs. Complementary expression was defined as the under-expression of parental gene and the over-expression of the retrogene in meiosis relative to mitosis. Complementary expression was assessed by comparing mean expression intensities between mitosis and meiosis, which were jointly estimated [30],[45]. More specifically, we estimated the probability that a given pair showed meiotic retrogene over-expression and meiotic parental gene under-expression. All gene pairs in each group (X→A and A→A) were used simultaneously in our model estimation, therefore avoiding the need for multiple hypothesis testing corrections. More details about the Bayesian estimation of complementary expression [30],[45] can be found in Model B in Statistical Analyses, Methods in Text S1.
Testis-biased gene analyses
Two analyses were performed using the testis-biased gene dataset. First, we investigated MSCI by analyzing the proportions of testis-biased X-linked and autosome-linked genes under-expressed in meiosis (Figure 4). We used the under-expressed classes defined by Bayesian model A (Statistical Analyses, Methods in Text S1). Second, we assessed the testis-biased chromosomal proportion for mitotic and meiotic phases (Figure 7). Genes having higher expression in mitosis than in meiosis were considered as mitotic genes (Figure 7B), whereas those having higher expression in meiosis were considered as meiotic genes (Figure 7C). All testis-biased genes chromosomal proportions were parsed from the Sebida database [3],[4],[44]. Testis-biased gene chromosomal proportions were measured relative to the total number of genes by chromosome, where autosomes were pooled together.
Supporting Information
Zdroje
1. ViçosoB
CharlesworthB
2006 Evolution on the X chromosome: unusual patterns and processes. Nature Rev Genet 7 645 653
2. DorusS
BusbySA
GerikeU
ShabanowitzJ
HuntDF
2006 Genomic and functional evolution of the Drosophila melanogaster sperm proteome. Nature Genet 38 1440 1445
3. ParisiM
NuttallR
NaimanD
BouffardG
MalleyJ
2003 Paucity of genes on the Drosophila X chromosome showing male-biased expression. Science 299 697 700
4. ParisiM
NuttallR
EdwardsP
MinorJ
NaimanD
2004 A survey of ovary-, testis-, and soma-biased gene expression in Drosophila melanogaster adults. Genome Biol 5 R40
5. RanzJM
Castillo-DavisCI
MeiklejohnCD
HartlDL
2003 Sex-Dependent Gene Expression and Evolution of the Drosophila Transcriptome. Science 300 1742 1745
6. WangPJ
PageDC
McCarreyJR
2005 Differential expression of sex-linked and autosomal germ-cell-specific genes during spermatogenesis in the mouse. Hum Mol Genet 14 2911 2918
7. ReinkeV
GilIS
WardS
KazmerK
2004 Genome-wide germline-enriched and sex-biased expression profiles in Caenorhabditis elegans. Development 131 311 323
8. EmersonJJ
KaessmannH
BetránE
LongM
2004 Extensive gene traffic on the mammalian X chromosome. Science 303 537 540
9. BetránE
ThorntonK
LongM
2002 Retroposed new genes out of the X in Drosophila. Genome Res 12 1854 1859
10. BaiY
CasolaC
FeschotteC
BetránE
2006 Comparative genomics reveals a constant rate of origination and convergent acquisition of functional retrogenes in Drosophila. Genome Biol 8 R11
11. DaiH
YoshimatsuTF
LongM
2006 Retrogene movement within- and between-chromosomes in the evolution of Drosophila genomes. Gene 385 96 102
12. VibranovskiMD
ZhangY
LongM
2009 General gene movement off the X chromosome in the Drosophila genus. Genome Res 19 897 903
13. MeiselRP
HanMV
HahnMW
2009 A Complex Suite of Forces Drives Gene Traffic from Drosophila X Chromosomes. Genome Biol Evol 0 176 188
14. LifschytzE
LindsleyDL
1972 The role of X-chromosome inactivation during spermatogenesis. Proc Natl Acad Sci U S A 69 182 186
15. RiceWR
1984 Sex chromosomes and the evolution of sexual dimorphism. Evolution 38 735 742
16. CharlesworthB
CoyneAJ
BartonNH
1987 The relative rates of evolution of sex chromosomes and autosomes. Am Nat 130 113 146
17. PattenMM
HaigD
2009 Maintenance or loss of genetic variation under sexual and parental antagonism at a sex-linked locus. Evolution. In press
18. GibsonJR
ChippindaleAK
RiceWR
2002 The X chromosome is a hot spot for sexually antagonistic fitness variation. Proc Biol Sci 269 499 505
19. EllegrenH
ParschJ
2007 The evolution of sex-biased genes and sex-biased gene expression. Nat Rev Genet 8 689 698
20. SturgillD
ZhangY
ParisiM
OliverB
2007 Demasculinization of X chromosomes in the Drosophila genus. Nature 8 238 241
21. RichlerC
SoreqH
WahrmanJ
1992 X inactivation in mammalian testis is correlated with inactive X-specific transcription. Nature Genet 2 192 195
22. KellyWG
SchanerCE
DernburgAF
LeeMH
KimSK
2002 X-chromosome silencing in the germline of C. elegans. Development 129 479 492
23. SchoenmakersS
WassenaarE
HoogerbruggeJW
LavenJS
GrootegoedJA
2009 Female meiotic sex chromosome inactivation in chicken. PLoS Genet 5 e1000466 doi:10.1371/journal.pgen.1000466
24. HenseW
BainesJF
ParschJ
2007 X chromosome inactivation during Drosophila spermatogenesis. PLoS Biol 5 e273 doi:10.1371/journal.pbio.0050273
25. NguyenDK
DistecheCM
2006 Dosage compensation of the active X chromosome in mammals. Nat Genet 38 47 53
26. GönczyP
MatunisE
DiNardoS
1997 bag-of-marbles and benign gonial cell neoplasm act in the germline to restrict proliferation during Drosophila spermatogenesis. Development 124 4361 4371
27. LoppinB
LepetitD
DorusS
CoubleP
KarrTL
2005 Origin and neofunctionalization of a Drosophila paternal effect gene essential for zygoteviability. Curr Biol 15 87 93
28. BarreauC
BensonE
GudmannsdottirE
NewtonF
White-CooperH
2008 Post-meiotic transcription in Drosophila testes. Development 135 1897 1902
29. StraubT
BeckerPB
2007 Dosage compensation: the beginning and end of generalization. Nature Rev Genet 8 47 57
30. LopesHF
MüllerP
RavishankerN
2007 Bayesian computational methods in biomedical research.
KhattreeR
NaikDN
Computational Methods in Biomedical Research Boca Raton Chapman & Hall/CRC 211 259
31. JohnstonR
WangB
NuttallR
DoctoleroM
EdwardsP
2004 FlyGEM, a full transcriptome array platform for the Drosophila community. Genome Biol 5 R19
32. MüllerP
ParmigianiG
RiceK
2006 FDR and Bayesian multiple comparisons rules.
BernardoJM
BayarriMJ
BergerJO
DawidAP
HeckermanD
SmithAFM
WestM
Bayesian Statistics 8 (with discussion) Oxford Oxford University Press 349 370
33. CarrelL
WillardHF
2005 X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature 434 400 404
34. OlivieriG
OlivieriA
1965 Autoradiographic study of nucleic acid synthesis during spermatogenesis in Drosophila Melanogaster. Mutation Research 2 366 380
35. GuptaV
ParisiM
SturgillD
NuttallR
DoctoleroM
2006 Global analysis of X-chromosome dosage compensation. J Biol 5 3
36. SchäferM
NayerniaK
EngelW
SchäferU
1995 Translational control in spermatogenesis. Dev Biol 172 344 352
37. PotrzebowskiL
VinckenboschN
MarquesAC
ChalmelF
JégouB
2008 Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes. PLoS Biol 6 e80 doi:10.1371/journal.pbio.0060080
38. KhilPP
SmirnovaNA
RomanienkoPJ
Camerini-OteroRD
2004 The mouse X chromosome is enriched for sex-biased genes not subject to selection by meiotic sex chromosome inactivation. Nature Genet 36 642 646
39. CharlesworthB
1991 The evolution of sex chromosomes. Science 251 1030 1033
40. BachtrogD
2006 A dynamic view of sex chromosome evolution. Curr Opin Genet Dev 16 578 585
41. ClarkM
HeathBD
AndersonC
KarrTL
2006 Induced paternal effects mimic cytoplasmic incompatibility in Drosophila. Genetics 173 727 734
42. Ghosh-RoyA
KulkarniM
KumarV
ShirolikarS
RayK
2004 Cytoplasmic dynein-dynactin complex is required for spermatid growth but not axoneme assembly in Drosophila. Mol Biol Cell 15 2470 2483
43. ChintapalliVR
WangJ
DowJAT
2007 Using FlyAtlas to identify better Drosophila models of human disease. Nat Genet 39 715 720
44. GnadF
ParschJ
2006 Sebida: a database for the functional and evolutionary analysis of genes with sex-biased expression. Bioinformatics 22 2577 2579
45. DoK-A
MüllerP
VannucciM
2006 Bayesian Inference for Gene Expression and Proteomics. Cambridge Cambridge University Press
46. GamermanD
LopesHF
2006 Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference (2nd Edition). Boca Raton Chapman & Hall/CRC
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2009 Číslo 11
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- 10 Reasons to be Tantalized by the B73 Maize Genome
- X Chromosome: Expression and Escape
- Stage-Specific Expression Profiling of Spermatogenesis Suggests that Meiotic Sex Chromosome Inactivation Drives Genomic Relocation of Testis-Expressed Genes
- RNA-Silencing Enzymes Pol IV and Pol V in Maize: More than one Flavor?