Temporal genomic evolution of bird sex chromosomes
© Wang et al.; licensee BioMed Central Ltd. 2014
Received: 26 October 2014
Accepted: 20 November 2014
Published: 12 December 2014
Sex chromosomes exhibit many unusual patterns in sequence and gene expression relative to autosomes. Birds have evolved a female heterogametic sex system (male ZZ, female ZW), through stepwise suppression of recombination between chrZ and chrW. To address the broad patterns and complex driving forces of Z chromosome evolution, we analyze here 45 newly available bird genomes and four species’ transcriptomes, over their course of recombination loss between the sex chromosomes.
We show Z chromosomes in general have a significantly higher substitution rate in introns and synonymous protein-coding sites than autosomes, driven by the male-to-female mutation bias (‘male-driven evolution’ effect). Our genome-wide estimate reveals that the degree of such a bias ranges from 1.6 to 3.8 among different species. G + C content of third codon positions exhibits the same trend of gradual changes with that of introns, between chrZ and autosomes or regions with increasing ages of becoming Z-linked, therefore codon usage bias in birds is probably driven by the mutational bias. On the other hand, Z chromosomes also evolve significantly faster at nonsynonymous sites relative to autosomes (‘fast-Z’ evolution). And species with a lower level of intronic heterozygosities tend to evolve even faster on the Z chromosome. Further analysis of fast-evolving genes’ enriched functional categories and sex-biased expression patterns support that, fast-Z evolution in birds is mainly driven by genetic drift. Finally, we show in species except for chicken, gene expression becomes more male-biased within Z-linked regions that have became hemizygous in females for a longer time, suggesting a lack of global dosage compensation in birds, and the reported regional dosage compensation in chicken has only evolved very recently.
In conclusion, we uncover that the sequence and expression patterns of Z chromosome genes covary with their ages of becoming Z-linked. In contrast to the mammalian X chromosomes, such patterns are mainly driven by mutational bias and genetic drift in birds, due to the opposite sex-biased inheritance of Z vs. X.
KeywordsAvian genome Evolutionary strata Male-driven evolution Fast-Z evolution Dosage compensation
Sex chromosomes of birds and mammals originated independently from different pairs of ancestral autosomes ,, and have experienced stepwise suppressions of recombination in parallel ,. Since recombination was lost at different time points, each affected region shows a distinctive level of sequence divergence between sex chromosomes from adjacent regions, manifested as a pattern of ‘evolutionary strata’ -. As a result of recombination loss, mammalian Y and avian W chromosomes are specific to one sex and usually have lost most functional genes, rendering their counterpart (X and Z chromosomes) hemizygous and asymmetrically transmitted between sexes. These features subject the sex chromosomes to special mutational and selective regimes apart from autosomes, and may select for evolution of dosage compensation on X/Z in response to the unbalanced gene dosages between sex chromosomes and autosomes. Indeed, X chromosomes have evolved numerous unusual properties in their sequence substitution and gene expression relative to autosomes ,. Previous genome-wide studies in mammals consistently found that the X chromosome has a lower level of interspecific sequence divergence at neutral sites relative to autosomes ,, due to the effect of ‘male-driven’ evolution. This refers to the fact that males usually undergo many more rounds of germ line cell divisions, thus accumulate more mutations per time unit than females in most animals . The X chromosome spends less time in males, therefore it is expected to have a lower mutation rate than autosomes and Y chromosome. On the other hand, it also tends to experience more effective adaptive evolution and show a faster evolution rate, since its hemizygous state in male would more readily expose the male recessive beneficial mutations to positive selection (‘fast-X’ evolution) ,. Indeed, an accelerated rate of protein sequence evolution has been observed on human and chimpanzee X chromosomes ,, and sperm proteins show an even more pronounced pattern of ‘fast-X’ evolution when being X-linked, suggesting they are under strong positive selection ,. Finally, the gene content of X chromosomes of both Drosophila and mammals has become ‘demasculinized’ (deficient for testis genes) or ‘feminized’ (enriched for ovary genes) ,. Sexual antagonism seems to be a direct cause, given that the X chromosome spends more time in females and is expected to accumulate mutations favorable to females at the expense of males . However, later studies challenged the necessity of invoking such an explanation, and showed meiotic sex chromosome inactivation (MSCI) or dosage compensation is sufficient to explain the demasculinized X-linked gene content in mammal  or Drosophila ,.
These complex forces acting on the X chromosome are also expected to shape the evolution of Z chromosome, but with very different outcomes. Opposite to the X chromosome, the Z chromosome is more often transmitted in males than in females. And recent studies found that the chicken Z does not seem to undergo MSCI and only has evolved regional instead of chromosome-wide dosage compensation like mammalian X chromosomes do -. These features are likely to produce a distinctive sequence and expression pattern of Z chromosome comparing to that of X chromosome and autosomes. Although ZW sex systems are widely distributed in birds, reptiles, fishes and butterflies, most of our current knowledge on evolution of Z chromosome comes from studies of chicken, which is the first species having genome sequenced among the ZW species. Previous analyses of some intronic fragments show chicken and other birds have a significantly higher level of neutral sequence divergence on Z chromosome than autosomes, indicating a higher Z-linked mutation rate caused by the male-driven evolution effect. This leads to an estimate of male-to-female mutation rate ratio (α) ranging from 1.7 to 6.5 in birds ,. Chicken Z-linked genes have also been found to evolve significantly faster in their protein coding sequences than the autosomal genes (‘fast-Z’ effect) . In contrast to the pattern of ‘fast-X’ evolution, Z-linked chicken female-biased or non-biased genes do not show a more pronounced pattern of fast evolution than male-biased genes, suggesting a non-adaptive underlying driving force . Theoretical simulations suggest demographic and mating system factors would impact the X/Z chromosomes more than autosomes in their effective population size (N e) ,,. The relative N e of Z-linked genes vs. autosomal genes can be further decreased from the neutral expectation of 3/4 due to for example, widespread sexual selection acting on male birds. Such an effect of genetic drift can also produce a pattern of ‘fast-Z’ evolution, through fixation of excessive slightly deleterious mutations on the Z chromosome. This alternative explanation remains to be tested with more empirical data, with multiple species of varying relative N e values of Z chromosome, so that to establish its generality. In addition, a general comparison between Z chromosomes vs. autosomes would not only offer a powerful independent test on all the hypotheses proposed for X chromosome evolution, but also provide important insights into the distinctive evolution mode of Z chromosome under different inheritance and gene regulation programs .
We have recently reconstructed the evolutionary history of bird Z/W sex chromosomes, based on analyses of a total of 17 species’ genomes . To our surprise, many species besides the known cases of Palaeognathae birds (ostrich and tinamou) , degenerate very slowly in their W chromosomes. This is in great contrast to the mammalian Ys or chicken W with few functional genes left ,, whose boundaries of evolutionary strata have never been determined due to the scarcity of Y/W linked genes and sequences. We are allowed to infer the ancestral karyotype of bird sex chromosomes with the newly available ostrich genome, and precisely demarcate the evolutionary strata along the Z chromosome of chicken as well as other bird species . Within each stratum, genes have been truly Z-linked for a different time span. This provides a unique system for us to investigate the temporal evolution of sex chromosome for the first time. In this work, we greatly expand our analyses to cover the high-quality genomes of a total of 45 bird species , and transcriptomes of 4 bird species that have become recently available ,-. Comparing to the previous analyses in chicken ,,, the studied species span the entire avian phylogeny, representing the most comprehensive analyses of bird sex chromosomes up to date. With these species’ recently resolved phylogeny based on the genomic dataset , we aim to address the following questions: first, how are the sequence and gene expression patterns different between the Z chromosome and autosomes and what are the main causes of such differences? Second, what is the temporal dynamics of such differences accompanied by the stepwise recombination loss between avian sex chromosomes? Finally, there have been continuous efforts of quantifying the life history traits of birds ever since the times of Darwin, which show tremendous variations between species, and frequently between sexes of the same species. The abundant genome resources available for so many bird species provide us great opportunities to associate their trait variations to the patterns of genome evolution, so that to illuminate the possible interaction between the two.
Results and discussion
Bird species show different degrees of male-driven evolution
Many factors affecting numbers of germ line cell divisions or DNA damage may underlie the variation of spontaneous mutation rate and α among species. Their associations have been extensively studied in plants  and mammals ,, but not in birds at a genome-wide level. We further examine a diverse array of life-history traits of the studied 45 species, and their association with lineage-specific synonymous substitution rates (dS) and α. We find that the median value of lineage specific synonymous substitution rates of autosomal genes, as an approximation of species’ mutation rate, exhibits a significant (P < 0.05, Pearson’s correlation test and Chi-square test) negative correlation with species’ life span, generation time, body mass and basal metabolic rate (Figure 2B-E), after controlling for the phylogenetic dependence. This pattern is consistent with the previous findings in mammals ,. Interestingly, we haven’t found a significant correlation between these traits vs. α, as a recent study did analyzing 32 mammalian genomes  (Additional file 1: Figure S2). This could result from the high variation of both spermatogenesis and oogenesis processes among different bird species (e.g., seasonal vs. continual breeders) , which may respectively influence the male and female mutation rates toward different directions.
Synonymous sites evolve faster in the older strata of Z chromosome
Sex chromosomes and autosomes differ not only in their spontaneous mutation rate, but also in their rate of evolution within coding regions, which directly reflect the different intensities of selection and genetic drift. We further study 306 Z-linked and 5280 macro-chromosome orthologous gene pairs’ coding sequences, which together comprise about 40% of the entire avian gene repertoire. We perform the comparison between the these two chromosome sets (Additional file 3: Table S2), as well as between chrZ vs. chr5 and chr4, to further control for the variation of chromosomal size and gene density within macrochromosomes . If the effect of codon usage bias (CUB) is weak, synonymous sites are expected to evolve nearly neutrally, thus show a pattern of higher substitution rate on the Z-linked loci than autosomal loci similarly as at intronic sites. We found a significant difference of synonymous substitutions rates (dS) between the Z-linked and autosomal loci in only 11 out of 45 species (P < 0.05, Wilcoxon test, Additional file 3: Table S2), and the pattern of significance depends on the macrochromosome chosen (Figure 1B), implying some degree of natural selection or mutational bias is impacting these species’ synonymous sites. We also find loci that have become Z-linked for a longer time (e.g., genes of S0) tend to have a higher dS value (Figure 1B) and this trend is consistent across all the studied species (Additional file 1: Figure S3). Further inspection of different measurement of CUB, including Frequency of Optimal Codons (FOP) (Figure 1C) and codon bias index (CBI) consistently confirms Z-linked genes have a significantly higher degree of CUB (P < 0.05, Wilcoxon test) compared to the genes on macro-chromosomes, in almost all the species studied (Additional file 1: Figure S4, Additional file 3: Table S2), and genes that have recently become Z-linked tend to have a higher level of CUB (Figure 1C). Ostrich is the only exception, because 2/3rd of its Z chromosome is still recombining with the W chromosome as a pseudoautosomal region and evolves like an autosome ,.
This pattern is in accordance with a higher CUB found on the X chromosomes of C. elegans, Drosophila and human ,, and may reflect a more efficient purifying selection on the hemizygous X- or Z-linked loci, where recessive deleterious mutations are more readily exposed to natural selection . However, genetic drift and mutational bias probably have more important contribution. After recombination was suppressed on the Z chromosome in female, its effective population size becomes smaller than autosomes, and it maybe further reduced due to the variation in male mating success . Genetic drift thus would fix excessive slightly deleterious mutations in synonymous sites, due to Hill-Robertson interference among linked loci. This effect explains the trend of dS among evolutionary strata with different ages of becoming Z-linked, of which ‘older’ Z-linked genes tend to be fixed with more slightly deleterious mutations. In addition, we also find the G + C content at the third codon position (GC3) shows a similar trend with that of intron (GCi), between Z/autosomes and different strata (Figure 1D-E) across species (Additional file 1: Figure S5-6). The positive correlation between the G + C content and synonymous substitution rate is consistent with the pattern that has been previously reported in repetitive sequences of chicken . These results together argue against selection is shaping codon usage bias and evolution of synonymous sites in birds, and suggest mutational bias probably plays a dominant role .
Nonsynonymous sites evolve faster in the younger strata of Z chromosome
Analogous to the ‘fast-X’ evolution patterns reported in Drosophila and human, fast-Z evolution may result from a more efficient fixation of recessive mutations which are beneficial to females . This predicts an enrichment of functional categories that are related to the female reproduction process or sex-biased gene expression among genes that are fast-evolving on the Z chromosome. To test this, we compare the Gene Ontology (GO) terms in which candidate genes undergoing adaptive evolution on the Z chromosome and autosomes are enriched. However, we haven’t found any GO terms directly related to sexual reproduction processes (Additional file 4: Table S3) for the Z-linked genes. Instead, macro- or micro-chromosomes are enriched for fast-evolving genes with GO terms like ‘female gonad development’ (GO: 8585), ‘ovulation cycle process’ (GO: 22602) and ‘regulation of reproductive process’ (GO: 2000241) etc. Consistent with the previous analysis , sex-biased expression pattern also doesn’t seem to affect the evolutionary rate (measured by dN/dS ratio) of Z-linked genes. We examine the brain and liver sex-biased gene expression level in both ostrich and chicken (also see below), and study their correlations with evolutionary rates with Z and autosome linked genes in separate (Additional file 1: Figure S9). In contrast to the pattern found on X chromosomes of Drosophila  and mammals , female-biased or non-biased genes do not necessarily evolve significantly faster than male-biased genes on the Z chromosome.
These results are consistent with the previous findings on Z-linked genes of chicken , and snakes , but not silkmoths . They suggest the alternative explanation for the fast-Z evolution, i.e., fixation of slightly deleterious mutations by genetic drift -. The effective population size of chrZ (N ez) is three quarters of that of autosomes (N eA), under the hypothesis of random mating and equal numbers of males and females within population. However, this hypothesis is usually invalid due to the variation in male mating success caused by sexual selection or social mating system . This will result in a much smaller ratio of N ez/N eA. As shown by previous theoretical simulation, this ratio has great impact on the degree of fast-X or fast-Z evolution . Indeed, when we compare the neutral polymorphism level approximated by intron SNP densities on the Z chromosome and autosomes, we found the Z chromosomes show a polymorphism level that is significantly (P < 0.001, Wilcoxon test) lower than three quarters of that of autosomes (Figure 3C). And the reduction in neutral polymorphism level of Z shows a significant negative correlation (P = 0.0037, Pearson’s correlation test) with its increase of lineage-specific evolutionary rate relative to autosomes (Figure 3D), providing direct evidence that the fast-Z evolution of birds is mainly driven by genetic drift. Under this scenario, majority of excessive mutations fixed on the Z chromosome are slightly deleterious. So a larger dN or dN/dS value in younger strata (Figure 3B-C) is probably not caused by recent adaptive evolution. Instead, it may reflect older strata are undergoing stronger purifying selection against slightly deleterious mutations, since they have become hemizygous in female at an earlier time point and exposed recessive deleterious mutations to natural selection for a longer time span. We have recently shown W chromosomes of many bird species degenerate slowly after recombination was suppressed . Therefore, in younger strata like S2, recessive deleterious mutations may still be sheltered on the Z chromosomes by the residual W-linked genes, which probably accounts for a higher dN value in the younger strata.
Finally, to test the correlation between male mating success vs. the fast-Z evolution in birds, we study several life history traits that indicate or impact the strength of male sexual selection including residual testis mass, social mating systems, color dichromatism and tail sexual dimorphism . None of these traits show significant correlation with the ratio of evolutionary rate between Z chromosome and autosomes (Additional file 1: Figure S10). This could due to a lack of enough sample species. For example, most studied species are characterized as monogamous species; therefore it is difficult to evaluate the impact of different social mating systems. Also, life history traits can rapidly turn over between species while our measurement of fast-Z evolution is a result of long-term evolution.
Gene expression become more male-biased over time on the Z chromosome
Z chromosomes are transmitted in a male-biased manner, and may be subjected to selection for dosage compensation, after they suppressed recombination with W and became hemizygous in female. The former predicts more male-biased gene expression on the Z chromosome, which is counteracted by the latter. To examine relative importance of the two, we finally investigate male vs. female expression ratio of Z-linked genes across four species (ostrich, emu, chicken and zebra finch) representing the major groups of avian species over their course of recombination suppression.
Due to the lack of global dosage compensation, Z chromosomes show a significantly lower expression level relative to other chromosomes (P < 0.05, Wilcoxon test, Additional file 1: Figure S12) in female tissues and are generally male-biasedly expressed. However, we have not found a significant difference of overall testis gene expression level between the Z chromosome and autosomes in chicken (P = 0.699, Wilcoxon test, Additional file 1: Figure S12), suggesting the male-biased expression of Z chromosome is simply a result of lack of dosage compensation, rather than being driven by sexual antagonistic selection due to its male-biased transmission. This is consistent with the enriched GO results shown above, i.e., there’s no excessive Z-linked male-related genes driven by adaptive evolution. However, this pattern needs to be validated in more bird species without any form of dosage compensation to confirm its generality.
A complex suite of forces is driving the evolution of sex chromosomes apart from their autosomal ancestors or other autosomes within the same genome ,. One great advantage of studying Z chromosome is it parses the factors of maleness and hemizygosity combined on the X chromosome. We find in this work, most sequence and expression patterns of Z chromosome that we observed across different bird species can be explained by its male-biased transmission, including its higher mutation rate, but not a higher adaptive evolution rate. We find the effective population size of Z chromosome is significantly reduced than that of autosomes under the neutral expectation (Figure 3D), and such a reduction probably results in fixation of excessive slightly deleterious mutations by genetic drift on the Z chromosome, producing a pattern of ‘fast-Z’ evolution across different bird species. By contrast, variation in male mating success that reduces the relative effective population size of Z would instead increase that of X, thus the efficiency of natural selection. This may explain why it is easier to detect adaptive evolution as a cause of ‘fast-X’ evolution ,.
We also show the local dosage compensation in chicken  may have evolved very recently, and global dosage compensation probably never evolved on the Z chromosome at the common ancestor of birds. Although it is still unclear why female heterogametic systems generally lack global dosage compensation , we find regions that lost recombination more recently in birds don’t immediately reduce their gene expression in female to half of the male (Figure 4). This can be attributed to both the slow degeneration of W-linked genes  and transcriptional buffering  of Z-linked genes, which together might reduce the selective pressure for dosage compensation.
Estimating the degree of male-driven evolution
The sequencing, assembly and annotation of the genomes are described in . We used the orthologous intronic alignments generated by SATé  for all sequenced bird species, after removing ambiguous alignment regions . We calculated the pairwise substitution distance between certain species vs. chicken by the baseml program in the PAML package . We used the median value of Z-linked and macrochromosome linked introns’ substitution distance, and the approach of  to estimate the α value for each species, using the equation α = (3Z/A-2)/(4-3Z/A) (‘Z/A’ stands for the divergence ratio of Z-linked vs. autosome linked loci). We further estimated the confidence intervals of α by the nonparametric bootstrapping method described in , with the bootstrapping process repeated 1000 times.
Life history traits for all the studied bird species used for the regression analyses were collected from the published work. Basal metabolic rate data is from ; residual testis size, body mass, clutch size information is from ,; life span and generation time data is originally collected from Encyclopedia of Life website (www.eol.org) then manually curated by Frank Rheindt at National University of Singapore. We used BayesTraits (http://www.evolution.rdg.ac.uk/BayesTraits.html), and the recently resolved phylogenetic tree of the studied 45 species  to account for the phylogenetic dependence in our regression analyses . In brief, we first tested whether a specific trait’s evolutionary change is fully dependent on the phylogenetic structure (λ = 1 implemented in BayesTraits) by maximum likelihood ratio test. Then we used Continuous module of BayesTraits to test the significance of the correlation if certain trait shows the dependence. Otherwise, we performed generalized linear model analysis to test the significance of correlation between multiple trait data vs. α.
Demarcate the evolutionary strata
We divided all Z-linked genes into groups with different ages of becoming Z-linked based on the information of avian evolutionary strata that we inferred recently . In brief, we first inferred the sex of the sequenced bird samples by comparing the sequencing read depth level between Z chromosome vs. autosomes. Those show similar level of depths are recognized as male and confirmed with sample record, and those show reduction in depth on the Z chromosome are recognized as female. We assembled fragments of W-linked sequences for all the species with female sequenced by SOAPdenovo , and then aligned and ordered these fragments along the Z chromosomes using chicken or ostrich chromosome sequence as reference. Since suppressed recombination between the avian Z and W chromosomes occurred in a stepwise manner ,, each instance of recombination suppression would render the affected region show a distinctive pattern of Z/W pairwise sequence divergence from the nearby regions. We scanned the Z/W sequence divergence with a stepwise window, and determined the boundaries of strata when significant differences of divergence level between adjacent windows were detected. The age of each stratum was further inferred by constructing gene trees using Z and W linked gametolog sequences of multiple species, if available. A gene tree showing the clustering of Z/W gametologs of the same species suggests the residing stratum emerged after the speciation event; a clustering pattern of Z gametologs of different species separated from W gametologs suggests the stratum emerged at the ancestor of the species . The oldest stratum where the earliest recombination suppression between Z/W occurred with the highest Z/W divergence level is termed strata zero (S0), and then S1, S2 etc. They also form a gradient of the time span for which each has become truly Z-linked. For species with completely degraded W chromosomes (e.g., chicken), we inferred their strata information through the orthologous relationship with their related species that has the information.
We hypothesized the preferred codons are generally conserved among different bird species with few exceptions, as reported among Drosophila  and plant species . We used CodonW (http://codonw.sourceforge.net/), and the frequency of chicken optimal codon indices  to study the codon usage bias of all the other avian species’ genes.
To calculate gene evolutionary rate, we used multiple sequence alignments of the orthologous genes’ coding regions as input derived from  for codeml package of PAML. The alignments have been removed for ambiguously aligned regions by GBLOCK  and manual inspection. We calculated lineage-specific dN, dS and dN/dS ratios under the ‘free-ratio’ (model = 1) model, which assumes a different evolutionary rate along each lineage. We then filtered genes with dN or dS value higher 2 or lower than 0.001 according to a previous study that characterized the saturation of dS in birds , before comparing the entire sets of values between Z chromosome and autosomes. To identify the fast-evolving genes, we compared the two alternative models with one assuming all lineages evolve with the same rate (‘one-ratio’ , model = 0), and the other assuming one rate along the focus lineage and one different rate on all the other (background) lineages (‘two-ratio’ , model = 2). Genes evolve significantly faster than the background lineages according to maximum likelihood test are characterized as candidate genes undergoing positive selection. These genes of different species were pooled and finally subjected to GO enrichment analyses based on annotation of their chicken orthologs, by Cytoscape v.2.8 , using gene set of the entire chicken genome as background, with hypergeometric test and a significance cutoff of P < 0.01.
We aligned the genomic reads against each species genome by bowtie2  and then screened the alignment results by a mapping quality of 20. We used SOAPsnp  and snpEff  to identify SNP sites within intronic sequences of bird genomes, of which male sample have been sequenced. SNP sites with mapped read depths lower than 4 or higher than 100 were filtered; we also removed those clustered SNP sites whose distance from each other is lower than 5bp.
Gene expression analysis
We downloaded the RNA-seq reads of both sexes of ostrich, emu, chicken and zebra finch from NCBI SRA database (http://www.ncbi.nlm.nih.gov/sra), which were collected by previous studies ,,,. We aligned the reads to each species’ genome by tophat  after quality-trimming the reads, and calculated the normalized gene expression level FPKM (Fragments Per Kilobase of transcript per Million mapped reads) by cufflinks  for each gene. We annotated and ordered the genes of emu using ostrich chromosomal sequence and genes derived from  as reference. Then we scanned the male vs. female expression ratio along the Z chromosome of each species by a five gene non-overlapping window.
We used SOAPdenovo-trans  to assemble pooled reads from female ostrich and emu into transcriptomes. Then we searched chicken MHM locus’ homolog in ostrich genome or ostrich/emu transcriptome by BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
We thank Beatriz Vicoso, Doris Bachtrog and Melissa Sayres at University of California, Berkeley for inspiring discussions during the project, and Peter Dunn at University of Wisconsin, Milwaukee for sharing some of the avian life history trait data. This work is supported by the National Key Basic Research Program of China (973 Program: 2013CB945200) to G.Z. All the genomic data that have been used in this study have been deposited on http://phybirds.genomics.org.cn/download.jsp.
- Bull JJ: Evolution of Sex Determining Mechanisms. 1983, The Benjamin/Cummings Publishing Company, Inc., San FranciscoGoogle Scholar
- Bellott DW, Skaletsky H, Pyntikova T, Mardis ER, Graves T, Kremitzki C, Brown LG, Rozen S, Warren WC, Wilson RK, Page DC: Convergent evolution of chicken Z and human X chromosomes by expansion and gene acquisition. Nature. 2010, 466 (7306): 612-616. 10.1038/nature09172.PubMedPubMed CentralView ArticleGoogle Scholar
- Cortez D, Marin R, Toledo-Flores D, Froidevaux L, Liechti A, Waters PD, Grutzner F, Kaessmann H: Origins and functional evolution of Y chromosomes across mammals. Nature. 2014, 508 (7497): 488-493. 10.1038/nature13151.PubMedView ArticleGoogle Scholar
- Zhou Q, Zhang J, Bachtrog D, An N, Huang Q, Jarvis ED, Gilbert MTP, Zhang G: Sex chromosomes demonstrate complex evolutionary trajectories across bird taxa. Science 2014, in press.,Google Scholar
- Lahn BT, Page DC: Four evolutionary strata on the human X chromosome. Science. 1999, 286 (5441): 964-967. 10.1126/science.286.5441.964.PubMedView ArticleGoogle Scholar
- Nam K, Ellegren H: The chicken (Gallus gallus) Z chromosome contains at least three nonlinear evolutionary strata. Genetics. 2008, 180 (2): 1131-1136. 10.1534/genetics.108.090324.PubMedPubMed CentralView ArticleGoogle Scholar
- Roesti M, Moser D, Berner D: Recombination in the threespine stickleback genome–patterns and consequences. Mol Ecol. 2013, 22 (11): 3014-3027. 10.1111/mec.12322.PubMedView ArticleGoogle Scholar
- Bergero R, Qiu S, Forrest A, Borthwick H, Charlesworth D: Expansion of the pseudo-autosomal region and ongoing recombination suppression in the Silene latifolia sex chromosomes. Genetics. 2013, 194 (3): 673-686. 10.1534/genetics.113.150755.PubMedPubMed CentralView ArticleGoogle Scholar
- Natri HM, Shikano T, Merila J: Progressive recombination suppression and differentiation in recently evolved neo-sex chromosomes. Mol Biol Evol. 2013, 30 (5): 1131-1144. 10.1093/molbev/mst035.PubMedPubMed CentralView ArticleGoogle Scholar
- Vicoso B, Charlesworth B: Evolution on the X chromosome: unusual patterns and processes. Nat Rev Genet. 2006, 7 (8): 645-653. 10.1038/nrg1914.PubMedView ArticleGoogle Scholar
- Ellegren H: Sex-chromosome evolution: recent progress and the influence of male and female heterogamety. Nat Rev Genet. 2011, 12 (3): 157-166. 10.1038/nrg2948.PubMedView ArticleGoogle Scholar
- Xu K, Oh S, Park T, Presgraves DC, Yi SV: Lineage-specific variation in slow- and fast-X evolution in primates. Evolution. 2012, 66 (6): 1751-1761. 10.1111/j.1558-5646.2011.01556.x.PubMedView ArticleGoogle Scholar
- Lu J, Wu CI: Weak selection revealed by the whole-genome comparison of the X chromosome and autosomes of human and chimpanzee. Proc Natl Acad Sci U S A. 2005, 102 (11): 4063-4067. 10.1073/pnas.0500436102.PubMedPubMed CentralView ArticleGoogle Scholar
- Li WH, Yi S, Makova K: Male-driven evolution. Curr Opin Genet Dev. 2002, 12 (6): 650-656. 10.1016/S0959-437X(02)00354-4.PubMedView ArticleGoogle Scholar
- Charlesworth B, Coyne J, Barton N: The relative rates of evolution of sex chromosomes and autosomes. Am Nat. 1987, 130: 113-146. 10.1086/284701.View ArticleGoogle Scholar
- Meisel RP, Connallon T: The faster-X effect: integrating theory and data. Trends Genet. 2013, 29 (9): 537-544. 10.1016/j.tig.2013.05.009.PubMedPubMed CentralView ArticleGoogle Scholar
- Nielsen R, Bustamante C, Clark AG, Glanowski S, Sackton TB, Hubisz MJ, Fledel-Alon A, Tanenbaum DM, Civello D, White TJ, JS J, Adams MD, Cargill M: A scan for positively selected genes in the genomes of humans and chimpanzees. PLoS Biol 2005, 3(6):e170.,Google Scholar
- Hvilsom C, Qian Y, Bataillon T, Li Y, Mailund T, Salle B, Carlsen F, Li R, Zheng H, Jiang T, Jiang H, Jin X, Munch K, Hobolth A, Siegismund HR, Wang J, Schierup MH: Extensive X-linked adaptive evolution in central chimpanzees. Proc Natl Acad Sci U S A. 2012, 109 (6): 2054-2059. 10.1073/pnas.1106877109.PubMedPubMed CentralView ArticleGoogle Scholar
- Torgerson DG, Singh RS: Enhanced adaptive evolution of sperm-expressed genes on the mammalian X chromosome. Heredity. 2006, 96 (1): 39-44.PubMedGoogle Scholar
- Torgerson DG, Singh RS: Sex-linked mammalian sperm proteins evolve faster than autosomal ones. Mol Biol Evol. 2003, 20 (10): 1705-1709. 10.1093/molbev/msg193.PubMedView ArticleGoogle Scholar
- Sturgill D, Zhang Y, Parisi M, Oliver B: Demasculinization of X chromosomes in the Drosophila genus. Nature. 2007, 450 (7167): 238-241. 10.1038/nature06330.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang YE, Vibranovski MD, Landback P, Marais GA, Long M: Chromosomal redistribution of male-biased genes in mammalian evolution with two bursts of gene gain on the X chromosome. PLoS Biol 2010, 8(10):e1000494.,Google Scholar
- Wu CI, Xu EY: Sexual antagonism and X inactivation–the SAXI hypothesis. Trends Genet. 2003, 19 (5): 243-247. 10.1016/S0168-9525(03)00058-1.PubMedView ArticleGoogle Scholar
- Potrzebowski L, Vinckenbosch N, Marques AC, Chalmel F, Jegou B, Kaessmann H: Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes. PLoS Biol 2008, 6(4):e80.,Google Scholar
- 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. 2009, 68 (5): 576-583. 10.1007/s00239-009-9235-4.PubMedView ArticleGoogle Scholar
- Bachtrog D, Toda NR, Lockton S: Dosage compensation and demasculinization of X chromosomes in Drosophila. Curr Biol. 2010, 20 (16): 1476-1481. 10.1016/j.cub.2010.06.076.PubMedPubMed CentralView ArticleGoogle Scholar
- Guioli S, Lovell-Badge R, Turner JM: Error-prone ZW pairing and no evidence for meiotic sex chromosome inactivation in the chicken germ line. PLoS Genet 2012, 8(3):e1002560.,Google Scholar
- Itoh Y, Replogle K, Kim YH, Wade J, Clayton DF, Arnold AP: Sex bias and dosage compensation in the zebra finch versus chicken genomes: general and specialized patterns among birds. Genome Res. 2010, 20 (4): 512-518. 10.1101/gr.102343.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Schoenmakers S, Wassenaar E, Hoogerbrugge JW, Laven JS, Grootegoed JA, Baarends WM: Female meiotic sex chromosome inactivation in chicken. PLoS Genet 2009, 5(5):e1000466.,Google Scholar
- Ellegren H, Fridolfsson AK: Male-driven evolution of DNA sequences in birds. Nat Genet. 1997, 17 (2): 182-184. 10.1038/ng1097-182.PubMedView ArticleGoogle Scholar
- Axelsson E, Smith NG, Sundstrom H, Berlin S, Ellegren H: Male-biased mutation rate and divergence in autosomal, Z-linked and W-linked introns of chicken and Turkey. Mol Biol Evol. 2004, 21 (8): 1538-1547. 10.1093/molbev/msh157.PubMedView ArticleGoogle Scholar
- Mank JE, Axelsson E, Ellegren H: Fast-X on the Z: rapid evolution of sex-linked genes in birds. Genome Res. 2007, 17 (5): 618-624. 10.1101/gr.6031907.PubMedPubMed CentralView ArticleGoogle Scholar
- Mank JE, Nam K, Ellegren H: Faster-Z evolution is predominantly due to genetic drift. Mol Biol Evol. 2010, 27 (3): 661-670. 10.1093/molbev/msp282.PubMedView ArticleGoogle Scholar
- Mank JE, Vicoso B, Berlin S, Charlesworth B: Effective population size and the faster-X effect: empirical results and their interpretation. Evolution. 2010, 64 (3): 663-674. 10.1111/j.1558-5646.2009.00853.x.PubMedView ArticleGoogle Scholar
- Vicoso B, Charlesworth B: Effective population size and the faster-X effect: an extended model. Evolution. 2009, 63 (9): 2413-2426. 10.1111/j.1558-5646.2009.00719.x.PubMedView ArticleGoogle Scholar
- Tsuda Y, Nishida-Umehara C, Ishijima J, Yamada K, Matsuda Y: Comparison of the Z and W sex chromosomal architectures in elegant crested tinamou (Eudromia elegans) and ostrich (Struthio camelus) and the process of sex chromosome differentiation in palaeognathous birds. Chromosoma. 2007, 116 (2): 159-173. 10.1007/s00412-006-0088-y.PubMedView ArticleGoogle Scholar
- Ayers KL, Davidson NM, Demiyah D, Roeszler KN, Grutzner F, Sinclair AH, Oshlack A, Smith CA: RNA sequencing reveals sexually dimorphic gene expression before gonadal differentiation in chicken and allows comprehensive annotation of the W-chromosome. Genome Biol 2013, 14(3):R26.,Google Scholar
- Zhang G, Li C, Li Q, Wang Z, Jarvis E, Gilbert T: Comparative genomics across modern bird species reveal insights into avian genome evolution and adaptation. Science 2014, in press.,Google Scholar
- Julien P, Brawand D, Soumillon M, Necsulea A, Liechti A, Schutz F, Daish T, Grutzner F, Kaessmann H: Mechanisms and evolutionary patterns of mammalian and avian dosage compensation. PLoS Biol 2012, 10(5):e1001328.,Google Scholar
- Adolfsson S, Ellegren H: Lack of dosage compensation accompanies the arrested stage of sex chromosome evolution in ostriches. Mol Biol Evol. 2013, 30 (4): 806-810. 10.1093/molbev/mst009.PubMedPubMed CentralView ArticleGoogle Scholar
- Vicoso B, Kaiser VB, Bachtrog D: Sex-biased gene expression at homomorphic sex chromosomes in emus and its implication for sex chromosome evolution. Proc Natl Acad Sci U S A. 2013, 110 (16): 6453-6458. 10.1073/pnas.1217027110.PubMedPubMed CentralView ArticleGoogle Scholar
- Wright AE, Moghadam HK, Mank JE: Trade-off between selection for dosage compensation and masculinization on the avian Z chromosome. Genetics. 2012, 192 (4): 1433-1445. 10.1534/genetics.112.145102.PubMedPubMed CentralView ArticleGoogle Scholar
- Jarvis E, Warnow T, Mirarab S, Houde P, Cracraft J, Boussau B, Howard J, Zhang G, Gilbert T: Whole genome analyses resolve the early branches to the tree of life of modern birds. Science 2014, in press.,Google Scholar
- Nanda I, Schlegelmilch K, Haaf T, Schartl M, Schmid M: Synteny conservation of the Z chromosome in 14 avian species (11 families) supports a role for Z dosage in avian sex determination. Cytogenet Genome Res. 2008, 122 (2): 150-156. 10.1159/000163092.PubMedView ArticleGoogle Scholar
- Shetty S, Griffin DK, Graves JA: Comparative painting reveals strong chromosome homology over 80 million years of bird evolution. Chromosome Res. 1999, 7 (4): 289-295. 10.1023/A:1009278914829.PubMedView ArticleGoogle Scholar
- Derjusheva S, Kurganova A, Habermann F, Gaginskaya E: High chromosome conservation detected by comparative chromosome painting in chicken, pigeon and passerine birds. Chromosome Res. 2004, 12 (7): 715-723. 10.1023/B:CHRO.0000045779.50641.00.PubMedView ArticleGoogle Scholar
- Nanda I, Benisch P, Fetting D, Haaf T, Schmid M: Synteny conservation of chicken macrochromosomes 1-10 in different avian lineages revealed by cross-species chromosome painting. Cytogenet Genome Res. 2011, 132 (3): 165-181. 10.1159/000322358.PubMedView ArticleGoogle Scholar
- Nanda I, Karl E, Volobouev V, Griffin DK, Schartl M, Schmid M: Extensive gross genomic rearrangements between chicken and Old World vultures (Falconiformes : Accipitridae). Cytogenet Genome Res. 2006, 112 (3–4): 286-295. 10.1159/000089883.PubMedView ArticleGoogle Scholar
- Nishida C, Ishijima J, Kosaka A, Tanabe H, Habermann FA, Griffin DK, Matsuda Y: Characterization of chromosome structures of Falconinae (Falconidae, Falconiformes, Aves) by chromosome painting and delineation of chromosome rearrangements during their differentiation. Chromosome Res. 2008, 16 (1): 171-181. 10.1007/s10577-007-1210-6.PubMedView ArticleGoogle Scholar
- Burt DW: Origin and evolution of avian microchromosomes. Cytogenet Genome Res. 2002, 96 (1–4): 97-112. 10.1159/000063018.PubMedView ArticleGoogle Scholar
- Miyata T, Hayashida H, Kuma K, Mitsuyasu K, Yasunaga T: Male-driven molecular evolution: a model and nucleotide sequence analysis. Cold Spring Harb Symp Quant Biol. 1987, 52: 863-867. 10.1101/SQB.1987.052.01.094.PubMedView ArticleGoogle Scholar
- Carmichael AN, Fridolfsson AK, Halverson J, Ellegren H: Male-biased mutation rates revealed from Z and W chromosome-linked ATP synthase alpha-subunit (ATP5A1) sequences in birds. J Mol Evol. 2000, 50 (5): 443-447.PubMedGoogle Scholar
- Bartosch-Harlid A, Berlin S, Smith NG, Moller AP, Ellegren H: Life history and the male mutation bias. Evolution. 2003, 57 (10): 2398-2406. 10.1554/03-036.PubMedView ArticleGoogle Scholar
- Smith SA, Donoghue MJ: Rates of molecular evolution are linked to life history in flowering plants. Science. 2008, 322 (5898): 86-89. 10.1126/science.1163197.PubMedView ArticleGoogle Scholar
- Wilson Sayres MA, Venditti C, Pagel M, Makova KD: Do variations in substitution rates and male mutation bias correlate with life-history traits? A study of 32 mammalian genomes. Evolution. 2011, 65 (10): 2800-2815. 10.1111/j.1558-5646.2011.01337.x.PubMedView ArticleGoogle Scholar
- Nikolaev SI, Montoya-Burgos JI, Popadin K, Parand L, Margulies EH, Antonarakis SE, Program N: Life-history traits drive the evolutionary rates of mammalian coding and noncoding genomic elements. Proc Natl Acad Sci U S A. 2007, 104 (51): 20443-20448. 10.1073/pnas.0705658104.PubMedPubMed CentralView ArticleGoogle Scholar
- Jamieson BG: Reproductive Biology and Phylogeny of Birds. 2007, CRC Press, LondonGoogle Scholar
- International Chicken Genome Sequencing C: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432 (7018): 695-716. 10.1038/nature03154.View ArticleGoogle Scholar
- Zhou Q, Bachtrog D: Sex-specific adaptation drives early sex chromosome evolution in Drosophila. Science. 2012, 337 (6092): 341-345. 10.1126/science.1225385.PubMedPubMed CentralView ArticleGoogle Scholar
- Singh ND, Davis JC, Petrov DA: X-linked genes evolve higher codon bias in Drosophila and Caenorhabditis. Genetics. 2005, 171 (1): 145-155. 10.1534/genetics.105.043497.PubMedPubMed CentralView ArticleGoogle Scholar
- Charlesworth B: The effect of life-history and mode of inheritance on neutral genetic variability. Genet Res. 2001, 77 (2): 153-166. 10.1017/S0016672301004979.PubMedView ArticleGoogle Scholar
- Webster MT, Axelsson E, Ellegren H: Strong regional biases in nucleotide substitution in the chicken genome. Mol Biol Evol. 2006, 23 (6): 1203-1216. 10.1093/molbev/msk008.PubMedView ArticleGoogle Scholar
- Rao Y, Wu G, Wang Z, Chai X, Nie Q, Zhang X: Mutation bias is the driving force of codon usage in the Gallus gallus genome. DNA Res. 2011, 18 (6): 499-512. 10.1093/dnares/dsr035.PubMedPubMed CentralView ArticleGoogle Scholar
- Grath S, Parsch J: Rate of amino acid substitution is influenced by the degree and conservation of male-biased transcription over 50 myr of Drosophila evolution. Genome Biol Evol. 2012, 4 (3): 346-359. 10.1093/gbe/evs012.PubMedPubMed CentralView ArticleGoogle Scholar
- Vicoso B, Emerson JJ, Zektser Y, Mahajan S, Bachtrog D: Comparative sex chromosome genomics in snakes: differentiation, evolutionary strata, and lack of global dosage compensation. PLoS Biol 2013, 11(8):e1001643.,Google Scholar
- Sackton TB, Corbett-Detig RB, Nagaraju J, Vaishna L, Arunkumar KP, Hartl DL: Positve selection drives faster-Z evolution in silkmoths. Evolution. 2014, 68 (8): 2331-2342.PubMedPubMed CentralGoogle Scholar
- Dunn PO, Whittingham LA, Pitcher TE: Mating systems, sperm competition, and the evolution of sexual dimorphism in birds. Evolution. 2001, 55 (1): 161-175. 10.1111/j.0014-3820.2001.tb01281.x.PubMedView ArticleGoogle Scholar
- Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, Sinclair AH: The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature. 2009, 461 (7261): 267-271. 10.1038/nature08298.PubMedView ArticleGoogle Scholar
- Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, Searle S, White S, Vilella AJ, Fairley S, Heger A, Kong L, Ponting CP, Jarvis ED, Mello CV, Minx P, Lovell P, Velho TA, Ferris M, Balakrishnan CN, Sinha S, Blatti C, London SE, Li Y, Lin YC, George J, Sweedler J, Southey B, Gunaratne P, Watson M, et al: The genome of a songbird. Nature. 2010, 464 (7289): 757-762. 10.1038/nature08819.PubMedPubMed CentralView ArticleGoogle Scholar
- Stenberg P, Larsson J: Buffering and the evolution of chromosome-wide gene regulation. Chromosoma. 2011, 120 (3): 213-225. 10.1007/s00412-011-0319-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu K, Raghavan S, Nelesen S, Linder CR, Warnow T: Rapid and accurate large-scale coestimation of sequence alignments and phylogenetic trees. Science. 2009, 324 (5934): 1561-1564. 10.1126/science.1171243.PubMedView ArticleGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Miyata T, Hayashida H, Kuma K, Mitsuyasu K, Yasunaga T: Male-driven molecular evolution: a model and nucleotide sequence analysis. Cold Spring Harbor Symposia on Quantitative Biology. 1987, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, 863-867.Google Scholar
- McNab BK: Ecological factors affect the level and scaling of avian BMR. Comp Biochem Physiol A Mol Integr Physiol. 2009, 152 (1): 22-45. 10.1016/j.cbpa.2008.08.021.PubMedView ArticleGoogle Scholar
- Calhim S, Birkhead TR: Testes size in birds: quality versus quantity—assumptions, errors, and estimates. Behav Ecol. 2007, 18 (1): 271-275. 10.1093/beheco/arl076.View ArticleGoogle Scholar
- Pagel M: Inferring the historical patterns of biological evolution. Nature. 1999, 401 (6756): 877-884. 10.1038/44766.PubMedView ArticleGoogle Scholar
- Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y, Tang J, Wu G, Zhang H, Shi Y, Liu Y, Yu C, Wang B, Lu Y, Han C, Cheung DW, Yiu SM, Peng S, Xiaoqian Z, Liu G, Liao X, Li Y, Yang H, Wang J, Lam TW, Wang J: SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Giga Sci 2012, 1(1):18.,Google Scholar
- Charlesworth D, Charlesworth B, Marais G: Steps in the evolution of heteromorphic sex chromosomes. Heredity. 2005, 95 (2): 118-128. 10.1038/sj.hdy.6800697.PubMedView ArticleGoogle Scholar
- Vicario S, Moriyama EN, Powell JR: Codon usage in twelve species of Drosophila. BMC Evol Biol 2007, 7:226.,Google Scholar
- Qiu S, Bergero R, Zeng K, Charlesworth D: Patterns of codon usage bias in Silene latifolia. Mol Biol Evol. 2011, 28 (1): 771-780. 10.1093/molbev/msq251.PubMedView ArticleGoogle Scholar
- Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C, Storz JF, Antunes A, Greenwold MJ, Meredith RW, Odeen A, Cui J, Zhou Q, Xu L, Pan H, Wang Z, Jin L, Zhang P, Hu H, Yang W, Hu J, Xiao J, Yang Z, Liu Y, Xie Q, Yu H, Lian J, Wen P, Zhang F, Li H, et al: Comparative genomics reveals insights into avian genome evolution and adaptation. Science 2014, in press.,Google Scholar
- Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56 (4): 564-577. 10.1080/10635150701472164.PubMedView ArticleGoogle Scholar
- Axelsson E, Hultin-Rosenberg L, Brandstrom M, Zwahlen M, Clayton DF, Ellegren H: Natural selection in avian protein-coding genes expressed in brain. Mol Ecol. 2008, 17 (12): 3008-3017. 10.1111/j.1365-294X.2008.03795.x.PubMedView ArticleGoogle Scholar
- Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011, 27 (3): 431-432. 10.1093/bioinformatics/btq675.PubMedPubMed CentralView ArticleGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9 (4): 357-359. 10.1038/nmeth.1923.PubMedPubMed CentralView ArticleGoogle Scholar
- Li R, Li Y, Fang X, Yang H, Wang J, Kristiansen K, Wang J: SNP detection for massively parallel whole-genome resequencing. Genome Res. 2009, 19 (6): 1124-1132. 10.1101/gr.088013.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Cingolani P, Platts A, le Wang L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012, 6 (2): 80-92. 10.4161/fly.19695.PubMedPubMed CentralView ArticleGoogle Scholar
- Brawand D, Soumillon M, Necsulea A, Julien P, Csardi G, Harrigan P, Weier M, Liechti A, Aximu-Petri A, Kircher M, Albert FW, Zeller U, Khaitovich P, Grutzner F, Bergmann S, Nielsen R, Paabo S, Kaessmann H: The evolution of gene expression levels in mammalian organs. Nature. 2011, 478 (7369): 343-348. 10.1038/nature10532.PubMedView ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.PubMedPubMed CentralView ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-515. 10.1038/nbt.1621.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie Y, Wu G, Tang J, Luo R, Patterson J, Liu S, Huang W, He G, Gu S, Li S, Zhou X, Lam TW, Li Y, Xu X, Wong GK, Wang J: SOAPdenovo-Trans: de novo transcriptome assembly with short RNA-Seq reads. Bioinformatics. 2014, 30 (12): 1660-1666. 10.1093/bioinformatics/btu077.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.