Document Detail

Incomplete sex chromosome dosage compensation in the Indian meal moth, Plodia interpunctella, based on de novo transcriptome assembly.
Jump to Full Text
MedLine Citation:
PMID:  23034217     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
Males and females experience differences in gene dose for loci in the nonrecombining region of heteromorphic sex chromosomes. If not compensated, this leads to expression imbalances, with the homogametic sex on average exhibiting greater expression due to the doubled gene dose. Many organisms with heteromorphic sex chromosomes display global dosage compensation mechanisms, which equalize gene expression levels between the sexes. However, birds and Schistosoma have been previously shown to lack chromosome-wide dosage compensation mechanisms, and the status in other female heterogametic taxa including Lepidoptera remains unresolved. To further our understanding of dosage compensation in female heterogametic taxa and to resolve its status in the lepidopterans, we assessed the Indian meal moth, Plodia interpunctella. As P. interpunctella lacks a complete reference genome, we conducted de novo transcriptome assembly combined with orthologous genomic location prediction from the related silkworm genome, Bombyx mori, to compare Z-linked and autosomal gene expression levels for each sex. We demonstrate that P. interpunctella lacks complete Z chromosome dosage compensation, female Z-linked genes having just over half the expression level of males and autosomal genes. This finding suggests that the Lepidoptera and possibly all female heterogametic taxa lack global dosage compensation, although more species will need to be sampled to confirm this assertion.
Authors:
Peter W Harrison; Judith E Mank; Nina Wedell
Related Documents :
6961457 - Quantitative karyotyping of human chromosomes by dual beam flow cytometry.
23845247 - Sex determination: how worms count to two.
86937 - Hoffman's violet and dahlia as specific stains for animal chromosomes.
1724417 - Anti-kinetochore staining for single laser, bivariate flow sorting of indian muntjac ch...
23412657 - Regulation and coordination of nuclear envelope and nuclear pore complex assembly.
18752197 - Comparison of maize similarity and dissimilarity genetic coefficients based on microsat...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't    
Journal Detail:
Title:  Genome biology and evolution     Volume:  4     ISSN:  1759-6653     ISO Abbreviation:  Genome Biol Evol     Publication Date:  2012  
Date Detail:
Created Date:  2012-11-07     Completed Date:  2013-04-12     Revised Date:  2013-10-24    
Medline Journal Info:
Nlm Unique ID:  101509707     Medline TA:  Genome Biol Evol     Country:  England    
Other Details:
Languages:  eng     Pagination:  1118-26     Citation Subset:  IM    
Affiliation:
Research Department of Genetics, Evolution and Environment, University College London, United Kingdom. p.w.harrison@ucl.ac.uk
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Animals
Dosage Compensation, Genetic / genetics*
Female
Gene Dosage
Gene Expression / physiology
Male
Moths / genetics*
Recombination, Genetic* / genetics,  physiology
Sex Characteristics
Sex Chromosomes / genetics*
Transcriptome / genetics*
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Genome Biol Evol
Journal ID (iso-abbrev): Genome Biol Evol
Journal ID (publisher-id): gbe
Journal ID (hwp): gbe
ISSN: 1759-6653
Publisher: Oxford University Press
Article Information
Download PDF
© The Author(s) 2012. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.
creative-commons:
Accepted Day: 26 Month: 9 Year: 2012
Print publication date: Year: 2012
Electronic publication date: Day: 3 Month: 10 Year: 2012
collection publication date: Year: 2012
pmc-release publication date: Day: 3 Month: 10 Year: 2012
Volume: 4 Issue: 11
First Page: 1118 Last Page: 1126
PubMed Id: 23034217
ID: 3514961
DOI: 10.1093/gbe/evs086
Publisher Id: evs086

Incomplete Sex Chromosome Dosage Compensation in the Indian Meal Moth, Plodia interpunctella, Based on De Novo Transcriptome Assembly
Peter W. Harrison1*
Judith E. Mank1
Nina Wedell2
1Research Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
2University of Exeter, Daphne du Maurier Centre for Ecology and Conservation, College of Life and Environmental Sciences, Cornwall, United Kingdom
Correspondence: *Corresponding author: E-mail: p.w.harrison@ucl.ac.uk.
Associate editor: Marta Wayne

Introduction

There is an inherent imbalance in gene dosage for sex-linked genes between the sexes (Bachtrog et al. 2011). In female heterogametic taxa, such as birds, snakes, and lepidopterans, ZW females have one fewer copy of all Z-linked genes compared with ZZ males. A converse situation exists in male heterogametic taxa, such as mammals and Drosophila, where XY males have half the gene dose for X-linked genes compared with XX females. Differences in gene dose result in differences in expression for many genes (Schlattl et al. 2011; Zhou et al. 2011), and this could produce differences in expression for many hundreds of genes in systems with highly differentiated sex chromosomes.

Although the relationship between dose and expression is complex (Papp et al. 2003; Birchler et al. 2005, 2007), it is clear that some genes are dosage sensitive (Redon et al. 2006; Estivill and Armengol 2007; Makino and McLysaght 2010), and for these genes, changes in gene dose can result in severe phenotypes. Ohno (1967) suggested that this imbalance in gene dosage would result in selection for a dosage compensation mechanism to restore expression of the single X or Z chromosome in the heterogametic sex to the diploid level that existed before degradation of the Y or W gene content. Originally envisaged to apply to the entirety of the sex chromosome, it is clear that animals actually exhibit a range of sex chromosome dosage compensation (Mank 2009; Mank et al. 2011). At one end of the spectrum are organisms with nearly complete chromosome-wide compensating mechanisms, such as Drosophila and Caenorhabditis elegans (Lucchesi 1973; McDonel et al. 2006; Ercan et al. 2007), where the average expression of the X chromosome between the sexes, and between the single X and the diploid autosomes in the heterogametic sex, is equal. At the other end of the spectrum are species with only partial compensation, where the majority of genes on the X or Z are uncompensated in the heterogametic sex. In these cases, only a subset of sex-linked genes shows balanced expression between the sexes, and the average expression of the single X or Z in the heterogametic sex is less than the autosomal average and also less than the X or Z chromosome average in the homogametic sex. Examples of incomplete dosage compensation have been documented in birds (Itoh et al. 2007; Naurin et al. 2011; Wolf and Bryk 2011), platypus (Deakin et al. 2008), and Schistosoma (Vicoso and Bachtog 2011). There are also species in the middle of the spectrum, such as Tribolium (Prince et al. 2010), and potentially the eutherian mammals (Julien et al. 2012, Pessia et al. 2012), though there is still an unresolved debate regarding the status of dosage compensation in the placental mammals, with some studies indicating complete dosage compensation and others showing only partial compensation (Xiong et al. 2010; Deng et al. 2011; Kharchenko et al. 2011; Lin et al. 2012).

This variation in sex chromosome dosage compensation begs interesting questions about why some clades evolve complex mechanisms of dosage compensation for whole chromosomes and others simply compensate dosage sensitive genes (Mank 2009; Mank and Ellegren 2009; Wright and Mank 2012). To address these questions, we require information about the status of dosage compensation in a range of animals, both male and female heterogametic. The Lepidoptera, comprising butterflies and moths, share an orthologous Z chromosome, although the W is not conserved (Sahara et al. 2012). The status of dosage compensation of the Lepidopteran Z chromosome is currently unclear. Initial evidence from the silkworm, Bombyx mori, indicates incomplete dosage compensation, similar to that found in birds (Arunkumar et al. 2009; Zha et al. 2009). However, a subsequent reanalysis of the same data set, which corrected for a number of microarray biases, called this result into question, suggesting that Z:A ratios are less than 1 in both females and males (Walters and Hardcastle 2011). Unfortunately, the data set was ultimately deemed too noisy to be conclusive.

As a large clade with ancient and highly distinct Z and W sex chromosomes, the Lepidoptera represent a valuable group to investigate the status of sex chromosome dosage compensation. To resolve whether chromosome-wide Z-linked gene dosage compensation occurs in Lepidoptera, we built a de novo transcriptome assembly from Illumina RNA-Seq data from the Indian meal moth, Plodia interpunctella, which lacks a complete draft reference genome. The three fully sequenced Lepidopteran genomes, Bombyx (Mita et al. 2004), Monarch (Zhan et al. 2011), and Heliconius (The Heliconius Genome Consortium 2012), show a remarkable conservation of Z chromosome gene content across both butterflies and moths. We were therefore able to use orthology from the fully sequenced B. mori genome to assign putative genomic location for assembled transcripts as either Z-linked or autosomal, and we then estimated male and female expression levels for all putative 1:1 orthologs. The results show that not only is this strategy effective for assessing genomic location and expression level in the absence of full genome sequences, but more importantly that Plodia lacks complete sex chromosome dosage compensation.


Materials and Methods
Sample Preparation and Sequencing

We collected whole-body samples from 20 male and 20 virgin female Indian meal moths (P. interpunctella) at 1 day of age. Larvae were reared at a standard density (10 larvae in 150 ml of medium) using standard protocols on a diet of bran midlings, yeast, honey, and glycerol and were maintained at 28°C with a 16:8 h light–dark cycle (Ingleby et al. 2010). These samples were divided into four nonoverlapping sex-specific pools, each consisting of 10 same-sex individuals, thus creating two replicate pools for each sex. Pooling is useful in that it reduces replicate variance, and therefore increases statistical power, although at the cost of the ability to track individuals. RNA was extracted using an RNeasy Qiagen kit and prepared for sequencing by the University of Exeter Biosciences Sequencing Service using standard Illumina protocols.

Whole transcriptome shotgun sequencing (RNA-Seq) was performed on an Illumina Genome Analyzer II, with a separate lane for each of the four pools. The resulting 76-bp paired-end reads were assessed for quality using FastQC (http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc, last accessed November 2, 2011) and filtered with FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html, last accessed August 5, 2011) to exclude reads where less than 90% of positions base calls had quality scores of Q20 or above. Additionally, any reads identified as containing residual adaptor sequence were removed before transcriptome assembly. Postfiltering, there were on average more than 24 million reads per pool.

De Novo Transcriptome Assembly

We built a de novo transcriptome assembly for P. interpunctella with the filtered RNA-Seq reads using Trinity (Grabherr et al. 2011). Trinity does not rely on orthology but rather builds the transcriptome de novo based on read information alone. Trinity clusters the sequence data into a series of de Brujin graphs, used to process similar reads that correspond to each gene, locus, or set of analogous genes within the transcriptome. These graphs are processed to determine full-length transcripts, identifying alternatively spliced isoforms and paralogous genes. We combined the four pools of filtered paired-end reads in a single de novo transcriptome assembly to enable orthology determination and comparisons of the generated contigs across the four samples. We used RSEM (Li and Dewey 2011) to map the filtered read sets back to the resulting transcriptome assembly to obtain read count estimates for each of the 58,979 contigs for each of the four pools.

Previous microarray studies typically filtered out all genes with expression below the background hybridization threshold, effectively removing loci with low or no expression. The digital nature of RNA-Seq data raises the potential for inappropriate filtering thresholds. In particular to this study, because gene silencing rates have been shown to be different between sex chromosomes and autosomes, and because this can affect estimates of dosage compensation (Kharchenko et al. 2011), we discarded all contigs that did not have at least four reads per million mappable reads in at least two samples. This avoids any potential problems with different proportions of contigs on the Z and autosomes without significant expression and also reduces the number of partial and erroneous contigs that typically have low read coverage. This resulted in 17,246 contigs. The absence of filtering in a previous study of sex chromosome dosage compensation (Xiong et al. 2010) was heavily criticized, and this called into question their conclusions based on their low-stringency filters (Castagné et al. 2011; Deng et al. 2011; Kharchenko et al. 2011; Lin et al. 2011). However, although filtering to remove genes with minimal expression is still debated (He et al. 2011), filtering to remove nonexpressed and inactivated genes is an important step when comparing sex chromosomes with autosomes. This controversy highlights the significance of the type of filtering that is performed and the effect this can have on the conclusions drawn from this type of expression study.

Orthology and Chromosome Location

We compared Plodia contigs with the genes from the annotated genome of the silkworm B. mori (Mita et al. 2004) obtained from KAIKObase version 3.2.1 (Shimomura et al. 2009). Putative orthology was determined using a reciprocal best-hit approach with a Basic Local Alignment Search Tool (BLAST) e-value cutoff of 1 × 1010 (Altschul et al. 1990). The 1:1 orthology prediction between Bombyx and Plodia allowed us to use the chromosome positions of Bombyx genes as a putative chromosome location for each of the Plodia contigs. Plodia contigs that reciprocally matched Bombyx Z genes were classed as Z-borne contigs and contigs that matched Bombyx autosomal genes were termed autosomal in Plodia. All subsequent analyses were based solely on inferred 1:1 orthologs, with the remaining contigs discarded. This allowed us to identify 350 Plodia contigs with putative Z-linkage and 7,997 putative autosomal contigs. Although contigs do not necessarily correspond directly to single genes, the methods of the Trinity assembly (Grabherr et al. 2011), the fact that we removed all contigs <200 bp and the exclusive use of 1:1 orthologs means that most Plodia contigs for the expression study represent single genes.

Expression Comparisons of Z versus Autosomes in Both Sexes

We calculated fragments per kilobase per million mappable reads (FPKM), a factor that corrects for contig length and variations in read depth across samples (Mortazavi et al. 2008) for 1:1 Plodia orthologs for each pool individually. Because the correlation between same-sex replicates was very high (male r2 = 0.830, female r2 = 0.869), we combined replicates and calculated mean male and female FPKM for each ortholog by averaging across same-sex replicate pools. Repeatability of FPKM estimates for same-sex replicate pools was also estimated with Spearman’s ρ, an estimate of rank order correlation, and replicates showed high similarity for both male (Spearman’s ρ = 0.906; P < 2.2 × 1016) and female (Spearman’s ρ = 0.928; P < 2.2 × 1016) pools.

The FPKM-corrected expression levels of Z-linked and autosomal contigs were compared within each sex. Mann–Whitney U (MWU) was employed to test for significant differences between average Z and autosome expression levels. A resampling test was also employed to determine the probability that a subset of the autosomal contigs could yield the significant differences observed between Z and autosome contigs in females. Random samples of 350 female autosomal contigs (corresponding to the number of Z-linked orthologs in our data set) were selected as pseudo-Z linked sets, to compare to the remaining autosomal contigs. Average expression for the two sets was calculated and MWU employed to test for significance between the groups. This additional evidence supports the predictions of genomic location and orthology, was replicated 1,000 times, and the number of significant results was recorded with exact permutation P values (Phipson and Smyth 2010).

In addition to Z to autosome ratios within each sex, we plotted the differences in expression profile between male and female Z-linked versus autosomal expression, first log2 transforming the data to aid visualization, then fitting linear regression lines. Z-linked contig expression was further analyzed by separating Z-linked contigs into four quartiles based on male expression levels. Male-to-female Z-linked expression was compared, and MWU was employed to test for significant differences between the sexes Z-linked expression within each quartile.

Gene Ontology

We tested whether there was any over- or underabundance of functional Gene Ontology (GO) terms (Ashburner et al. 2000) for the Z chromosome compared with the autosomes and for contigs exhibiting male- or female-limited expression using the Drosophila melanogaster GO. We used BLAST to compare each Plodia contig to the annotated D. melanogaster genome annotation, taking the highest ranked BLAST hit [release 5.44 (Adams et al. 2000)] for all BLASTs where the highest hit passed a threshold of P < 1 × 105. Z-chromosome, male-limited, and female-limited target lists were separately compared with all expressed contigs using GOrilla to identify enriched GO terms with a threshold of P < 1 × 105 with multiple hypothesis correction (Eden et al. 2007, 2009).


Results

Using RNA-Seq data from four nonoverlapping sex-specific pools, each containing 10 moth individuals, we constructed a de novo transcriptome assembly for the non-model Lepidopteran, the Indian meal moth, P. interpunctella, using Trinity (Grabherr et al. 2011). Synteny was inferred from the silk moth, B. mori, annotated genome sequence (Mita et al. 2004; Shimomura et al. 2009) for 1:1 orthologs, yielding 350 putative Z-linked contigs and 7,997 putative autosomal contigs. Correlations between FPKM-corrected expression levels in same-sex replicate pools were highly significant for both the male (r2 = 0.830; Spearman’s ρ = 0.906; P < 2.2 × 1016) and female (r2 = 0.869; Spearman’s ρ = 0.928; P < 2.2 × 1016) pools.

Replicate pools were combined to provide average autosomal and Z-linked expression (table 1 and fig. 1). Dosage compensation is a mechanism by which the heterogametic sex, in the case of P. interpunctella, the female, would hypertranscribe Z-linked loci to balance expression of the single Z chromosome with the diploid autosomes (Ohno 1967; Mank 2009; Mank et al. 2011). It is therefore most appropriate to test for dosage compensation by comparing average Z expression with average autosomal expression (Z:A) for males and females separately. Male average Z:A of 0.954 was not significantly different from one (MWU; P = ns), in stark contrast to female Z:A of 0.534 (MWU; P = 2.2 × 1016; table 1), as would be expected in the absence of dosage compensation mechanisms.

To support the orthology predictions of genomic location in Plodia, we used a resampling approach to determine the probability that a random subsample of 350 female autosomal contigs would yield at least as much expression difference as we observed between female Z-linked and autosomal contigs. We compared average expression between autosomal and 1,000 resampled pseudo-Z contig sets, and in these comparisons, none of the resamples showed as much expression difference as our inferred Z and autosomal data set (P = 1 × 104, 1,000 replicates). Comparing the expression of the pseudo-Z sets with the observed expression difference alone does not provide sufficient statistical rigor, as even if the resampling results are not as extreme as the real data, the differences between the pseudo-Z contig sets and the remaining female autosomal contigs could still be significant. Therefore, we also assessed statistical difference between our autosomal and resampled pseudo-Z replicates with MWU tests using a cutoff of P < 5 × 102. Few of the resampled pseudo-Z sets showed any significant difference to the autosomal data (P = 4.4 × 102, 1,000 replicates), again indicating that Z contigs on average had been inferred correctly and that our data were sufficient to assess dosage compensation differences between the Z chromosome and autosomes. The high level of synteny between two other Lepidopteran genomes, the Monarch (Zhan et al. 2011) and Heliconius (The Heliconius Genome Consortium 2012) butterflies with Bombyx, adds further support to our methodology of using orthology to predict Z-linkage. These combined lines of evidence add support for our genomic location predictions.

Although it is most appropriate to compare Z:A ratios within each sex separately, comparing male and female expression is useful in estimating the degree of dosage effect. Figure 2 illustrates the dosage effects of the Z chromosome, with female expression of Z-linked contigs averaging 0.58× that of males due to their halved gene dose. Figure 2 also suggests that the degree of dosage effect is different for lowly versus highly expressed contigs. We investigated this further by comparing male and female FPKM for Z-linked contigs based on average expression in males (fig. 3). This analysis, which separated Z-linked contigs into quartiles of increasing male expression, demonstrates that dosage effects are least manifest for lowly expressed contigs (MWU; P = ns) and increase for each quartile, with significant differences between males and females for the second (MWU; P = 4.41 × 103), third (MWU; P = 4.22 × 109), and fourth (MWU; P = 6.77 × 1011) quartiles. No functional difference could be detected using GO analyses between Z-linked and autosomal contigs or between sex-limited contigs and autosomes (supplementary material S1, Supplementary Material online). Interestingly, sex-limited autosomal contigs showed high expression levels, with male-limited autosomal contigs averaging 110.72 FPKM and female-limited autosomal contigs averaging 375.78 FPKM, compared with the autosomal average (58.79 FPKM in males and 60.98 FPKM in females).


Discussion

Our results indicate that Plodia lacks complete Z chromosome dosage compensation, with average female expression of Z-linked contigs slightly more than half that of male Z chromosome expression and slightly more than half of the average autosomal expression levels. This suggests that incomplete dosage compensation may be the norm in Lepidoptera, although data on additional species, and a resolution of the debate surrounding the Bombyx data (Zha et al. 2009; Walters and Hardcastle 2011), will be needed to confirm the generality of this pattern throughout the clade.

Plodia exhibits a similar pattern of partial, gene-by-gene compensation observed in both birds (Arnold et al. 2008) and Schistosoma (Vicoso and Bachtrog 2011), the female heterogametic animal taxa that have been previously assessed. This indicates that Z chromosomes may not generally be associated with global chromosome-wide mechanisms of dosage compensation that exist in some male heterogametic systems. The effects of Z chromosome dose in Plodia, the imbalance in Z to autosome expression levels, are somewhat greater (0.53) than the dosage effects observed in birds, which range from 0.6 to 0.75 (Mank and Ellegren 2009; Wolf and Bryk 2011). Our estimate of dosage effect in Plodia is also lower than the estimate from Bombyx, which ranges from 0.67 to 0.83 depending on the tissue sampled (Zha et al. 2009). However, this estimate may be confounded by data filtering and noise or difference in sampling strategy (Walters and Hardcastle 2011). It is worth noting that the majority of lepidopterans lack recombination in females (Marec and Traut 1993; Yoshido et al. 2005; Wang 2011), just as recombination is absent in heterogametic XY males in Drosophila (Morgan 1914).

Methodological Remarks

To better understand the evolution of sex chromosomes and dosage compensation mechanisms, we require a broader understanding of the status of dosage compensation in a diverse range of male and female heterogametic taxa. However, the desired diversity in organisms does not always match the availability of complete reference genomes with which to assess dosage compensation. Genome sequencing has focused on particular groups, and although more and more genomes are becoming available, there are notable gaps. Additionally, a large number of genome projects have either only sequenced the homogametic sex or have not sequenced in sufficient depth to accurately reconstruct the sex chromosomes in enough detail to enable mapping and quantification of RNA-Seq data. Our study has overcome these barriers by conducting de novo transcriptomic sequence assembly and inferring genomic location based on synteny with the related silk moth, B. mori. Using this method, we have been able to assess dosage compensation in an organism without having to first construct a complete reference genome, saving considerable time and resources. This study presents an effective strategy for significantly increasing the number of species for which transcriptome analysis in general, and dosage compensation in particular, can be determined. However, this may not be a satisfactory solution in all taxa if a representative genome has not yet been completed or synteny is low. Although this approach is subject to biases of orthology-based studies, it allows data points from multiple species to be gathered from homologous chromosomes to strengthen our understanding of dosage compensation in each clade and will help to resolve the status of dosage compensation in contested taxa. Future studies could also use RNA-seq barcoding information to track male and female individuals (Craig et al. 2008; Cronn et al. 2008) and then verify Z-linkage by comparing male and female heterozygosity levels.

Gene Expression and Dosage Effects

In addition to an absence of global dosage compensation in Plodia, the actual dosage difference between the sexes is not uniform across all types of contigs. Our comparison of male and female expression levels for Z-linked contigs showed that dosage differences between the sexes are more pronounced for highly expressed contigs. Detailed quartile-based analysis of expression level showed that dosage effects for contigs with the lowest expression levels were not significantly different between the sexes and that there was a trend for greater dosage effect between the sexes as expression level increased. This demonstrates the multiplying effect of Z copy number in highly expressed contigs when a global dosage compensation mechanism is absent. The largest dosage differences are observed in highly expressed contigs, which is intriguing, particularly considering that it is believed that highly expressed genes are generally more functionally important (having larger fitness effects), as shown by their elevated levels of purifying selection, conservation, and slower evolution (Pál et al. 2001; Drummond et al. 2005; Liao and Zhang 2006; Resch et al. 2007). The large dosage differences in highly expressed contigs may also be indicative of the importance of sex-specific fitness effects, discussed further later.

Unlike Plodia, birds (Arnold et al. 2008) and Schistosoma (Vicoso and Bachtrog 2011), some male heterogametic systems possess global compensation mechanisms. For example, Drosophila males hypertranscribe the X chromosome to equalize expression levels, achieving parity between the sexes, as well as between the single male X and the diploid autosomes (Bachtrog et al. 2010). This means that male-biased X-linked genes in Drosophila must be expressed over and above the already hypertranscribed levels associated with dosage compensation. At high expression levels, the ability to hypertranscribe from the single X chromosome may become saturated, explaining the deficit of male-biased highly expressed genes on the Drosophila X chromosome (Vicoso and Charlesworth 2009). Applied to our data in Plodia, the higher dosage difference for highly expressed contigs suggests that the ability of regulatory networks to buffer dosage effects of Z-linked genes decreases with expression level, possibly because the transcriptional process becomes similarly saturated.

Sexual Conflict over Dosage Compensation

The lack of complete dosage compensation in Plodia raises questions about sex chromosome evolution and the role of genes on the sex chromosomes. The single copy of the Z chromosome possessed by females equates to on average slightly more than half the expression level of males, as observed for the majority of female Z-linked contigs in Plodia. The differing expression levels presumably translate into differences in protein levels, and because Z-linked loci interact with the remainder of the genome through genetic networks, dosage differences between males and females have the potential to affect a large proportion of autosomal loci. This difference in expression of the Z chromosome in female heterogametic species and the X chromosome in male heterogametic systems creates conflict between the sexes over optimal transcription level (Mank et al. 2011). Selection to hypertranscribe X- or Z-linked genes in the heterogametic sex can lead to overexpression in the homogametic sex, which for dosage-sensitive genes could be as harmful as underexpression. Dosage compensation effectively decouples the correlation between male and female transcription rates, thereby resolving sexual conflict related to Z chromosome dosage differences.

Sexual conflict over transcription rates is believed to have brought about X chromosome inactivation in therian mammals. Therian X inactivation, previously thought to be a mechanism of dosage compensation, is in fact the result of conflict over optimal transcription of X-linked dosage-sensitive genes (Pessia et al. 2012). In this clade, although a number of dosage-sensitive genes have been compensated to the levels of autosomal genes in both sexes, the majority of X-linked genes remain expressed at lower levels than the autosomal genes in both sexes. However, the status of dosage compensation in the eutherians remains contentious, with several recent conflicting studies (Xiong et al. 2010; Deng et al. 2011; Kharchenko et al. 2011; Julien et al. 2012; Lin et al. 2012).

In addition to Plodia, there are a number of species, including birds (Itoh et al. 2007; Mank 2009), Schistosoma (Vicoso and Bachtrog 2011), and platypus (Deakin et al. 2008) that only compensate a minority of X or Z chromosome genes. Birds and Schistosoma are female heterogametic, whereas platypus males have five XY sex chromosome pairs. The fitness effects of incomplete dosage compensation in these organisms are unclear, as sexual conflict over optimal transcription rates of Z- or X-linked genes is unresolved. It may be that the majority of uncompensated genes are insensitive to dosage effects, although the correlation between dosage effects and expression level that we observed would suggest otherwise. Alternatively, there may yet be unknown mechanisms of minimizing dosage effects in these systems.

Concluding Remarks

Here, we employed de novo transcriptome assembly methods, in conjunction with genomic positional location inference from the related B. mori reference genome, to demonstrate that the Indian meal moth, P. interpunctella, lacks complete Z chromosome dosage compensation. This approach expands the number of data points that can be obtained for each clade, increasing the scope for studies of sex chromosome dosage compensation and other transcriptomic analyses.

Our analysis shows that P. interpunctella lacks complete Z chromosome dosage compensation, showing approximately half the expression in females compared with males for Z-linked contigs, the largest dose effect related to sex chromosomes yet observed. This work suggests the entire Lepidopteran clade may lack complete Z chromosome dosage compensation, with previous studies on B. mori proving inconclusive (Zha et al. 2009; Walters and Hardcastle 2011), though studies on a range of species will be required to confirm this assertion. This finding adds further weight to the observation that complete sex chromosome dosage compensation may be less common in female- compared with male-heterogametic species (Mank 2009; Naurin et al. 2011), although this is still based on relatively few data points. It is hoped that the status of dosage compensation in an expanded range of species with old and highly diverged sex chromosomes will soon be determined, allowing for the synthesis and understanding of the evolution of complete versus incomplete sex chromosome dosage compensation and the sexual conflict over expression that underlies it.


Supplementary Material

Supplementary material S1 is available at Genome Biology and Evolution online (http://www.gbe.oxfordjournals.org/).


Acknowledgments

The authors thank JC Perry, MA Pointer, AE Wright, F Zimmer, and S Montgomery for helpful comments on this manuscript. They acknowledge the use of the UCL Unity SMP Facility, and associated support services, in the completion of this work. This work was supported by a Royal Society Wolfson Award to N.W. and the European Research Council under the Framework 7 Agreement to J.E.M. [grant agreement 260233].


Literature Cited
Adams MD,et al. The genome sequence of Drosophila melanogasterScienceYear: 20002872185219510731132
Altschul SF,Gish W,Miller W,Myers EF,Lipman DJ. Basic local alignment search toolJ Mol Biol.Year: 19902154034102231712
Arnold AP,Itoh Y,Melamed E. A bird's-eye view of sex chromosome dosage compensationAnnu Rev Genom Hum Genet.Year: 20089109127
Arunkumar KP,Mita K,Nagaraju J. The silkworm Z chromosome is enriched in testis-specific genesGeneticsYear: 200918249350119332883
Ashburner M,et al. Gene ontology: tool for the unification of biologyNat Genet.Year: 200025252910802651
Bachtrog D,Toda NRT,Lockton S. Dosage compensation and demasculinization of X chromosomes in DrosophilaCurr Biol.Year: 2010201476148120705467
Bachtrog D,et al. Are all sex chromosomes created equal?Trends Genet.Year: 20112735035721962970
Birchler JA,Riddle NC,Auger DL,Veitia RA. Dosage balance in gene regulation: biological implicationsTrends Genet.Year: 20052121922615797617
Birchler JA,Yao H,Chudalayandi S. Biological consequences of dosage dependent gene regulatory systemsBiochim Biophys Acta.Year: 2007176942242817276527
Castagné R,et al. The choice of the filtering method in microarrays affects the inference regarding dosage compensation of the active X-chromosomePLoS OneYear: 20116e2395621912656
Chambers JM,Cleveland WS,Kleiner B,Tukey PA. Graphical methods for data analysisYear: 1983Belmont (CA)Wadsworth & Brooks/Cole
Craig DW,et al. Identification of genetic variants using bar-coded multiplexed sequencingNat Methods.Year: 2008588789318794863
Cronn R,et al. Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technologyNucleic Acids Res.Year: 200836e12218753151
Deakin JE,Hore TA,Koina E,Graves JAM. The status of dosage compensation in the multiple X chromosomes of platypusPLoS Genet.Year: 20084e100014018654631
Deng XX,et al. Evidence for compensatory upregulation of expressed X-linked genes in mammals, Caenorhabditis elegans and Drosophila melanogasterNat Genet.Year: 2011431179118522019781
Drummond DA,et al. Why highly expressed proteins evolve slowlyProc Natl Acad Sci U S A.Year: 2005102143381434316176987
Eden E,et al. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene listsBMC BioinformaticsYear: 2009104819192299
Eden E,Lipson D,Yogev S,Yakhini Z. Discovering motifs in ranked lists of DNA sequencesPLoS Comput Biol.Year: 20073e3917381235
Ercan S,et al. X chromosome repression by localization of the C. elegans dosage compensation machinery to sites of transcription initiationNat Genet.Year: 20073940340817293863
Estivill X,Armengol L. Copy number variants and common disorders: filling the gaps and exploring complexity in genome-wide association studiesPLoS Genet.Year: 200731787179917953491
Grabherr MG,et al. Full-length transcriptome assembly from RNA-Seq data without a reference genomeNat Biotechnol.Year: 20112964465221572440
He X,et al. He et al. replyNat Genet.Year: 20114311711172
Heliconius Genome ConsortiumButterfly genome reveals promiscuous exchange of mimicry adaptations among speciesNatureYear: 2012487949822722851
Ingleby FC,Lewis Z,Wedell N. Level of sperm competition promotes evolution of male ejaculate allocation patterns in a mothAnim Behav.Year: 2010803743
Itoh Y,et al. Dosage compensation is less effective in birds than in mammalsJ Biol.Year: 20076217352797
Julien P,et al. Mechanisms and evolutionary patterns of mammalian and avian dosage compensationPLoS Biol.Year: 201210e100132822615540
Kharchenko PV,Xi RB,Park PJ. Evidence for dosage compensation between the X chromosome and autosomes in mammalsNat Genet.Year: 2011431167116922120048
Li B,Dewey C. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genomeBMC BioinformaticsYear: 20111232321816040
Liao B-Y,Zhang J. Low rates of expression profile divergence in highly expressed genes and tissue-specific genes during mammalian evolutionMol Biol Evol.Year: 2006231119112816520335
Lin F,Xing K,Zhang J,He X. Expression reduction in mammalian X chromosome evolution refutes Ohno’s hypothesis of dosage compensationProc Natl Acad Sci U S A.Year: 2012109117521175722753487
Lin H,et al. Relative overexpression of X-linked genes in mouse embryonic stem cells is consistent with Ohno's hypothesisNat Genet.Year: 2011431169117022120049
Lucchesi JC. Dosage compensation in DrosophilaAnnu Rev Genet.Year: 197372252374205904
Makino T,McLysaght A. Ohnologs in the human genome are dosage balanced and frequently associated with diseaseProc Natl Acad Sci U S A.Year: 20101079270927420439718
Mank JE. The W, X, Y, and Z of sex chromosome dosage compensationTrends Genet.Year: 20092522623319359064
Mank JE,Ellegren H. All dosage compensation is local: gene by gene regulation of sex-biased expression on the chicken Z chromosomeHeredityYear: 200910231232018985062
Mank JE,Hosken DJ,Wedell N. Some inconvenient truths about sex chromosome dosage compensation and the potential role of sexual conflictEvolutionYear: 2011652133214421790564
Marec F,Traut W. Synaptonemal complexes in female and male meiotic prophase of Ephestia kuehniella (Lepidoptera)HeredityYear: 199371394404
McDonel P,Jans J,Peterson BK,Meyer BJ. Clustered DNA motifs mark X chromosomes for repression by a dosage compensation complexNatureYear: 200644461461817122774
Mita K,et al. The genome sequence of silkworm, Bombyx moriDNA Res.Year: 200411273515141943
Morgan TH. No crossing over the male Drosophila of genes in the second and third pair of chromosomesBiol Bull.Year: 191426195204
Mortazavi A,Williams BA,McCue K,Schaeffer L,Wold B. Mapping and quantifying mammalian transcriptomes by RNA-SeqNat Methods.Year: 2008562162818516045
Naurin S,Hansson B,Hasselquist D,Kim YH,Bensch S. The sex-biased brain: sexual dimorphism in gene expression in two species of songbirdsBMC GenomicsYear: 2011123721235773
Ohno S. Sex chromosomes and sex linked genesYear: 1967Berlin (Germany)Springer-Verlag
Pál C,Papp B,Hurst LD. Highly expressed genes in yeast evolve slowlyGeneticsYear: 200115892793111430355
Papp B,Pal C,Hurst LD. Dosage sensitivity and the evolution of gene families in yeastNatureYear: 200342419419712853957
Pessia E,et al. Mammalian X chromosome inactivation evolved as a dosage-compensation mechanism for dosage-sensitive genes on the X chromosomeProc Natl Acad Sci U S A.Year: 20121095346535122392987
Phipson B,Smyth GK. Permutation P-values should never be zero. Calculating exact P-values when permutations are randomly drawnStat Appl Genet Mol Biol.Year: 20109 Article39.
Prince EG,Kirkland D,Demuth JP. Hyperexpression of the X chromosome in both sexes results in extensive female bias on X-linked genes in the flour beetleGenome Biol Evol.Year: 2010233634620624738
Redon R,et al. Global variation in copy number in the human genomeNatureYear: 200644444445417122850
Resch AM,et al. Widespread positive selection in synonymous sites of mammalian genesMol Biol Evol.Year: 2007241821183117522087
Sahara k,Yoshido A,Traut W. Sex chromosome evolution in moths and butterfliesChromosome Res.Year: 201220839422187366
Schlattl A,Anders S,Waszak SM,Huber W,Korbel JO. Relating CNVs to transcriptome data at fine resolution: assessment of the effect of variant size, type, and overlap with functional regionsGenome Res.Year: 2011212004201321862627
Shimomura M,et al. KAIKObase: an integrated silkworm genome database and data mining toolBMC GenomicsYear: 20091048619843344
Vicoso B,Bachtrog D. Lack of global dosage compensation in Schistosoma mansoni, a female-heterogametic parasiteGenome Biol Evol.Year: 2011323023521317157
Vicoso B,Charlesworth B. The deficit of male-biased genes on the D. melanogaster X chromosome is expression-dependent: a consequence of dosage compensation?J Mol Evol.Year: 20096857658319407921
Walters JR,Hardcastle TJ. Getting a full dose? Reconsidering sex chromosome dosage compensation in the silkworm, Bombyx moriGenome Biol Evol.Year: 2011349150421508430
Wang B. Achiasmy or heterochiasmy: does meiotic recombination occur in female Lepidoptera?J Res LepidopteraYear: 2011444345
Wolf JBW,Bryk J. General lack of global dosage compensation in ZZ/ZW systems? Broadening the perspective with RNA-SeqBMC GenomicsYear: 2011129121284834
Wright AE,Mank JE. Battle of the sexes: conflict over dosage-sensitive genes and the origin of X chromosome inactivationProc Natl Acad Sci U S A.Year: 20121095144514522440748
Xiong YY,et al. RNA sequencing shows no dosage compensation of the active X-chromosomeNat Genet.Year: 2010421043104721102464
Yoshido A,Bando H,Yasukochi H,Sahara K. The Bombyx mori karyotype and the assignment of linkage groupsGeneticsYear: 200517067568515802516
Zha X,et al. Dosage analysis of Z chromosome genes using microarray in silkworm, Bombyx moriInsect Biochem Mol.Year: 200935315321
Zhan S,Merlin C,Boore JL,Reppert SM. The monarch butterfly genome yields insights into long-distance migrationCellYear: 20111471171118522118469
Zhou J,Lemos B,Dopman EB,Hartl DL. Copy-number variation: the balance between gene dosage and expression in Drosophila melanogasterGenome Biol Evol.Year: 201131014102421979154

Figures

[Figure ID: evs086-F1]
FIG. 1.— 

FPKM levels of autosomal (green) and Z-linked (orange) contigs in sex-averaged pools. Genomic location is based on inferred synteny with Bombyx mori. Plot boxes represent the median and interquartile range of FPKM expression level. Box notches are a visual aid that are approximate to the 95% confidence interval of the median and signify strong evidence that the medians differ for sample pairs if the notches do not overlap (Chambers et al. 1983). The plot tails extend to the most extreme data point that is no more than 1.5× the interquartile range.



[Figure ID: evs086-F2]
FIG. 2.— 

Sex-average log2 FPKM values for autosomal (green) and Z-linked (orange) contigs based on inferred synteny with Bombyx. In both classes, a linear regression line has been fitted.



[Figure ID: evs086-F3]
FIG. 3.— 

FPKM-normalized average expression levels from male and female Plodia Z-linked contigs, identified by inferred synteny with Bombyx mori. Contigs were divided into four quartiles based on male expression level. Mann–Whitney U tests of significance of average Z-linked expression identified no significant difference between sexes in quartile 1 (P = 2.78 × 10−1) but significant differences in quartiles 2 (P = 4.41 × 10−3), 3 (P = 4.22 × 10−9), and 4 (P = 6.77 × 10−11). Plot boxes represent the median and interquartile range of FPKM expression level. Box notches are a visual aid that are approximate to the 95% confidence interval of the median and signify strong evidence that the medians differ for sample pairs if the notches do not overlap (Chambers et al. 1983). The plot tails extend to the most extreme data point that is no more than the 1.5× interquartile range.



Tables
[TableWrap ID: evs086-T1] Table 1 

Average FPKM Levels for Male and Female Plodia Z-linked and Autosomal Contigs


Male Female
Z 56.107 32.540
A 58.786 60.983
Z:A 0.954 0.534
MWU 8.96 × 10−1 2.2 × 10−16

evs086-TF1NOTE.—Z:A expression ratios for each sex and the P value for a Mann–Whitney U (MWU) test of significant difference between average Z and autosomal expression are listed.



Article Categories:
  • Research Article

Keywords: dosage compensation, Lepidoptera, sex chromosomes, de novo transcriptome assembly, orthology.

Previous Document:  Population diversity of ORFan genes in Escherichia coli.
Next Document:  Wegener's granulomatosis in a middle-aged woman presenting with dyspnea, rash, hemoptysis and recurr...