- Research article
- Open Access
CenH3 evolution in diploids and polyploids of three angiosperm genera
© Masonbrink et al.; licensee BioMed Central. 2014
- Received: 15 August 2014
- Accepted: 12 December 2014
- Published: 30 December 2014
Centromeric DNA sequences alone are neither necessary nor sufficient for centromere specification. The centromere specific histone, CenH3, evolves rapidly in many species, perhaps as a coevolutionary response to rapidly evolving centromeric DNA. To gain insight into CenH3 evolution, we characterized patterns of nucleotide and protein diversity among diploids and allopolyploids within three diverse angiosperm genera, Brassica, Oryza, and Gossypium (cotton), with a focus on evidence for diversifying selection in the various domains of the CenH3 gene. In addition, we compare expression profiles and alternative splicing patterns for CenH3 in representatives of each genus.
All three genera retain both duplicated CenH3 copies, while Brassica and Gossypium exhibit pronounced homoeologous expression level bias. Comparisons among genera reveal shared and unique aspects of CenH3 evolution, variable levels of diversifying selection in different CenH3 domains, and that alternative splicing contributes significantly to CenH3 diversity.
Since the N terminus is subject to diversifying selection but the DNA binding domains do not appear to be, rapidly evolving centromere sequences are unlikely to be the primary driver of CenH3 sequence diversification. At present, the functional explanation for the diversity generated by both conventional protein evolution in the N terminal domain, as well as alternative splicing, remains unexplained.
- Diversifying selection
- CenH3 evolution
- Alternative splicing
The centromere is a specific region of the eukaryotic chromosome that is the assembly point of the kinetochore, a group of proteins that act as a tether for microtubules during cell division. Although eukaryotic centromeres have highly conserved machinery for chromosome segregation, centromere sequences and binding proteins specific to centromeric chromatin are highly variable, even among closely related taxa ,. Specific retroelements and highly homogenized tandem repeats are common in the DNA of eukaryotic centromeres ,, though these too apparently diverge rapidly among closely related species ,,. Appropriate recruitment of these potentially co-evolving molecular components to the same site is enigmatic, as is their mechanism of spread and homogenization among chromosomes. Explanations for the apparent paradox between conservation of function but variability in sequence typically invoke an interplay between centromere function, centromere sequences, and epigenetic factors, such that DNA sequences per se become less functionally constrained ,,. Accordingly, it is generally accepted that centromeres are specified epigenetically .
The evolution of CenH3 is of particular interest due to its centrality in centromere specification and function . Unlike its highly conserved counterpart (histone H3), CenH3 has extensive sequence variability, particularly in two regions: (1) the non-canonical NH2 terminal tail, (2) the longer loop 1 region ,. DNA sequence diversity as well as alternative splicing can both play roles in generating diversity, which is thought to compensate for the fast-evolving centromeric DNA to ensure consistent centromeric function, although centromeric sequences are neither necessary nor sufficient for kinetochore assembly ,.
Given the twin observations of rapid sequence homogenization of centromeric repeats among chromosomes within species, yet rapid divergence among repeats between species, the question arises as to the fate of centromeric repeats, and CenH3 evolution, in allopolyploids. Polyploidy plays an integral role in the evolution of many organisms, particularly plants -. The consequences of polyploidy are often extensive -, and the resulting duplicated genes have myriad possible fates ,, including gene loss. Patterns of gene retention and loss following genome duplication have been extensively studied with respect to their broad classifications -. Because two presumably divergent suites of centromeric sequences become reunited into common nucleus at the time of allopolyploid formation, it is intriguing to investigate the subsequent evolutionary dynamics of the centromeric repeats as well as the now duplicated CenH3 sequences. In most modern diploids, CenH3 appears to have returned to single copy status following paleopolyploidy events , with a few notable exceptions . In contrast, recent allopolyploids often have multiple CenH3 gene copies -.
Here we evaluate the fate of CenH3 in allopolyploids from three divergent genera (Brassica, Oryza, and Gossypium) to address the question of whether the evolution of CenH3 is similar across a broad range of angiosperm taxa. The genera selected contain a diversity of allopolyploid species having either monophyletic or polyphyletic origins. The Brassica genus contains three diverse and widely cultivated diploid species (genomes designated A – C) and three allopolyploid species resulting from independent polyploidization events (BBCC, AACC, AABB) , while Oryza contains multiple diploid genome groups (designated A – G) and allopolyploids of diverse genomic origin (BBCC, CCDD, and HHJJ) . Gossypium includes 45 diploid species divided into genome groups designated A – G and K , as well as a single, monophyletic  polyploid clade (AD genome) containing 6 species. Previous work on CenH3 in Brassica and Oryza allopolyploids has focused on selection in specific regions of the gene and the relative expression of the retained homoeologs ,. Less is known about CenH3 and centromeric evolution in allopolyploid Gossypium. A prior study reported a centromeric gypsy-like retroelement (CRG) present in all centromeres of both the allotetraploid G. hirsutum (AD-genome) and the model progenitor D-genome diploid (G. raimondii), but absent from A-genome species . Neither the sequence nor the expression of duplicated CenH3 have been evaluated.
We characterize CenH3 sequence evolution on in three phylogenetically disparate angiosperm genera containing diploids and allopolyploids, and assess patterns of molecular evolution. We address whether allopolyploids retain duplicated copies following allopolyploid formation; the dynamics of sequence evolution of the duplicated, newly co-resident sequences; and the relative expression levels of homoeologous copies. In addition to reporting on CenH3 sequence evolution within and between genera, we describe novel patterns of alternative splicing in CenH3.
We cloned and sequenced the CenH3 genes from 7 diploid and 5 allotetraploid Gossypium species. The length of genomic sequences varied from 2565 to 2673 bp in the diploids and 2654 to 2673 bp in the polyploids (Additional file 1), although protein-coding length was consistent for all species, 492 bp. While length variation was not detected among cDNAs, we cannot account for length variation that occurs outside of our external primers in the first and last exons. The structure of CenH3, 7 exons and 6 introns, was conserved among all polyploid species, (diploid cDNA sequences were inferred from genomic sequences). As expected, CenH3 was largely conserved across the genus, with the majority of polymorphisms occurring in introns. Gossypium exiguum exhibited the greatest difference observed among the diploids (93 bp deletion in intron 6); the remainder of the polymorphisms were small (<10 bp), five of which were phylogenetically informative (Additional file 2).
Jukes Cantor corrected estimates of Ka/Ks within each genus, and between each genus and outgroup
Significant p-values associated with codons in Brassica and Oryza , as inferred from MEME analysis
CenH3evolution in allopolyploids
Given the single copy status observed in most diploid angiosperms sequenced to date, despite an evolutionary history, which encompasses multiple episodes of polyploidy, we wished to address whether or not CenH3 has been retained in duplicate following recent allopolyploidy events. In Gossypium, we found that each diploid representative had a single CenH3 gene and that all allopolyploid species had two, indicating retention of both parental copies, as reported for Brassica and Oryza ,,.
To address whether or not duplicated CenH3 sequences evolve independently of one another following allopolyploid formation, or if instead they are subjected to some form of sequence interaction or homogenization, we manually analyzed CenH3 sequences from all three genera. This lack of independence has been demonstrated for other homoeologous single copy genes in allopolyploids, most notably in Gossypium where the phenomenon was first described ,. CenH3 gene conversion was absent in all three genera, which is consistent with previous reports in Brassica and Oryza ,.
Genomic DNA sequences for CenH3 from Gossypium were used to construct a maximum likelihood tree, which concurs with the currently accepted phylogeny for the genus (Additional file 3). We resolved monophyletic clades for both the AT and DT homoeologs and their respective model diploid progenitors. Gene conversion was not detected in the homoeologs, further confirming that independent evolution of CenH3 homeologs occurred in this 1–2 MYD polyploid clade.
As mentioned above, there are various possible fates for genes duplicated via polyploidy. While both parental copies of CenH3 were retained in all studied allopolyploid species of Oryza, Brassica, and Gossypium, the transcriptional usage of each parental copy can vary from equivalent to complete silencing of one parental copy. To assess expression of CenH3 in Gossypium, we used three independent methods that allow us to assess the relative expression of homoeologs. The same tissue source was used for clone-counting and chromatogram measurements, while RNA-seq data from different sources was used to investigate other aspects of CenH3 expression. The RNA-seq data sources included (A2 vs D5) to determine the relative expression of CenH3 in the model diploid progenitors to the cotton polyploids, a synthetic hybrid between these two diploids (A2XD5F1), a synthetic polyploid (2_A2D1), and a domesticated and wild accession of AD1 (maxxa and yuc respectively).
Due to the sequence similarity between homoeologs, specific primers could not be designed for quantitative PCR, and thus we measured homoeolog expression bias with three separate methods, accompanied by different caveats. Due to the low number of SNPs between the homoeologs, a large number of RNA-seq reads could not be allocated to a particular homoeolog. We also cloned CenH3 cDNAs, which were counted to calculate relative expression of homoeologs, albeit with a smaller sample size (Additional file 4).
G. hirsutum (AD1) had the highest frequency of spliced transcripts at 44.4% (Additional file 5). Each species had a different proportion of each splicing variant; for example, the exon 137–142 deletion was present in 25% of G. hirsutum (AD1) clones, while it was present at 2.9-7.8% of other species. 17.6% of G. mustelinum (AD4) clones included the 45 bp intron insertion at position 99, which was only present 5.6-7.8% in the other species (Figure 5).
Despite its crucial role in centromere specification, many aspects of CenH3 evolution are poorly understood. Recent years have brought advances in our understanding of centromere epigenetics and evolution. For example, tandem repeats in many species have an evolutionary relationship with the CenH3 protein to reduce nucleosomal bending energy , multiple proteins interact with centromeric DNA to induce positive supercoiling of centromeric DNA , and CenPA (CenH3 in plants) provides the foundation for binding other kinetochore proteins .
CenH3sequence evolution in angiosperms
The present study extends our knowledge of the pace and process of CenH3 evolution by evaluating genomic and expression changes in three diverse angiosperm genera. Previous research in Brassica and Oryza showed that the CenH3 N terminus and CATD (loop 1 and α2 helix domains) sequences were under diversifying selection in lineage specific manners ,. To assess the generality of these findings in angiosperms, we reanalyzed the Oryza and Brassica sequences at a finer scale to identify specific regions of CenH3 that have actively diversified, and performed the first CenH3 sequence analysis for Gossypium. Across all three genera, diversification generally occurs in the N terminus, a result consistent with previous reports of rapid evolution in this domain, but contrary to reports of diversification in the CenPA Targeting Domain ,-.
Retention and expression of CenH3alleles in allopolyploids
In all three genera, both homoeologous copies are retained following genome doubling, demonstrating that restoration to single copy status, as widely observed among modern “diploid” plants, need not occur quickly following WGD events. Interestingly, expression of homoeologous CenH3 copies in Gossypium exhibited directional bias in all samples, although differences were not always statistically significant (Figure 4); this result contrasts with reports from Brassica and Oryza allopolyploids ,. In Brassica allotetraploids, a variety of CenH3 expression patterns were found for homoeologs, from a 2:1 ratio in an accession of B. juncea to complete B-genome CenH3 suppression in an accession of B. carinata . In allotetraploid O. minuta and O. alta, CenH3 expression is unbiased . The variation in expression profiles among allopolyploids in these three genera is notable and without an obvious explanation, although it seems likely that homoeolog expression levels reflect the unique genomic and evolutionary idiosyncrasies that characterize hybridization and genome doubling in each genus.
Alternative splicing of CenH3
In addition to non-synonymous evolution, protein diversity may be generated by alternative splicing. Considering the three genera collectively, it appears that alternative splicing frequently modifies the N terminus of CenH3 in Brassica and Gossypium, and this is the only domain modified in Oryza ,. Interestingly, the N terminus often cannot be aligned among closely related genera, yet it is necessary for centromeric deposition of CenH3 during meiosis in A. thaliana ,. The N terminus also interacts with kinetochore proteins in S. cerivisiae . Perhaps alternative splicing represents another means to generate the diversity of sequences necessary for CenH3 to target centromeres in meiosis, or a rapid defense that introduces new CenH3 proteins in response to increased centromere size. An alternative is that alternative splicing can also lead to a differing abundance of alternate transcripts between cells and tissues, which has implications in centromeric DNA if CenH3 and centromeric DNA are indeed coevolving.
By comparing the CenH3 sequences from three disparate angiosperm genera, we have gained insight into the rates and regions of evolution in this important protein. The most commonly mutated domain is the N terminus, which is also subject to alternative splicing, and contributes significantly to diversity in the N terminus of Brassica and Gossypium CenH3s. Alternative splicing is largely absent from the histone fold domain, even though loop 1 and α2 helix domains bind centromeric DNA, which therefore likely are subject to length constraints.
While CenH3 is considered a rapidly evolving protein, the N terminus is the only domain that is unalignable among closely related genera and is the most diverse domain. Roles have been attributed to the N terminus of CenH3, such as ubiquitin-mediated proteolysis  and it is indispensable for meiotic localization of CenH3, yet the domain is expendable in mitosis ,. In humans the N terminus is a hotspot for posttranslational modifications that interact with other centromeric proteins , and perhaps this is the case in plants as well. At present, the functional explanation for the diversity generated by both conventional protein evolution in the N terminal domain, as well as alternative splicing, remain unexplained.
The basis for the interest in CenH3 is its ability to bind to rapidly evolving centromeric repeats and yet still interact with conserved elements in the kinetochore. The dynamics of this relationship at the polyploid scale are increasingly complex due to the duplication of CenH3 and exposure to a new regime centromere repeats. The mono and polyphyletic origins of polyploids seem to influence the direction and level of expression bias between CenH3 homoeologs (Figure 4) and is most pronounced after hybridization, as is seen in the F1(A2D5) of Gossypium, but not following polyploidy 2(A2D1). Another compelling issue is obtaining diversity in these transcripts to respond to repeat evolution in the centromere, which is sparse in our tests for selection. Alternative splicing may be an another means to obtain this diversity, although neither selection nor alternative splicing modified the DNA binding domain of CenH3. Since the N terminus is the only domain modified by alternative splicing and subject to diversifying selection, but the DNA binding domains are not, rapidly evolving centromere sequences are unlikely to be the primary driver of CenH3 sequence diversification.
Leaves (2–4 cm) and leaf buds were collected from eight diploid cottons: G. arboreum cv. 101 [A2] (ISC427583), G. anomalum [B1] (ISC447893), G. robinsonii [C2] ISC451818), G. raimondii [D5] (ISC429440), G. stocksii [E1] (ISC447876), G. longicalyx [F1] (ISC418550), G. bickii [G1] (ISC414834), G. exiguum [K] (ISC416400), and from five allopolyploid cottons (G. hirsutum cv. TM1 [AD1] (ISC451819), G. barbadense cv. Pima S6 [AD2] (ISC451820), G. tomentosum 95 [AD3] (ISC451821), G. mustelinum local lab accession [AD4] (ISC429442), and G. darwinii PW45 [AD5] (ISC429431). For purposes of phylogenetic reconstruction, we included the outgroup species Gossypioides kirkii (ISC 418555) . All lines were grown in the Pohl Conservatory at Iowa State University and were used for both DNA and RNA extractions. DNA was extracted using the Qiagen DNeasy Plant Kit following the manufacturers recommended protocol. RNA was extracted using the Sigma Spectrum RNA Extraction Kit following the manufacturers recommended protocol with the following modifications: protocol A was followed at step 4, one wash each was performed for wash solutions I and II, and on-column DNA digestion was performed with the Sigma On-Column DNase I Digest Set.
CenH3 gene sequence data for G. raimondii was obtained from Phytozome ,, from which primers were designed (Additional file 6) for PCR amplification and sequencing from other species. PCR amplifications were performed using the manufacturers recommended reaction mixtures/cycling conditions and a melting temperature of 57°C. PCR amplicons from all diploid accessions (except G. exiguum) were cleaned via the Qiaquick PCR column cleanup (Qiagen) and sequenced with the amplification primers and a set of internal sequencing primers (Additional file 6).
PCR products from G. exiguum and all polyploid species were visualized on an Invitrogen E-gel to isolate bands, and cloned with the P-GEM-T Easy Vector ligation kit (Promega) and Top10 Competent Cells (Invitrogen) according to the recommended protocol. Clones were sequenced at the Iowa State DNA Sequencing Facility using both M13 primers and internal primers (Additional file 6).
Evaluation of selection
Jukes Cantor corrected estimates of Ka/Ks were measured using DNAsp. The MEME software package  accessed via the (http://www.datamonkey.org) server , was used to test for selection at the codon level. The automatic selection tool was used to determine the correct substitution models for Brassica (F81), Gossypium (F81), and Oryza (HKY85). The significance level cutoff was set at P = <0.05. We used Arabidopsis thaliana (24 million years divergence (MYD); Lysak et al. ), Gossypioides kirki (13.6 MYD; Cronn et al. ), and Brachypodium distachyon (46 MYD; Sanderson ) as the outgroups for Brassica (7.9 MYD; Jacquemin et al. ), Gossypium (5–10 MYD; Senchina et al. ), and Oryza (15 MYD; Sanderson ), respectively.
We used the SynMap tool of CoGe (http://genomeevolution.org/CoGe) , to identify blocks of syntenic orthologs to evaluate whole genome Ka/Ks for a representative species from each genus (Brassica rapa, Oryza sativa, Gossypium raimondi). The following parameters were used: BlastN, relative gene order, −D 50, −A 10, quota align merge –Dm 80, quota align with a ratio of coverage depth at 3:1 for (B. rapa: A. thaliana), 1:1 for (O. sativa: B. distachyon), and 6:1 for (G. raimondii: Theobroma cacao), overlap distance 40. Each species was compared to their previously described outgroups, except T. cacao (60 MYD)  was the outgroup for G. raimondii, since Gossypioides kirki lacks a sequenced genome.
Genomic DNA sequences were aligned using CLUSTALW  in BioEdit  and converted to NEXUS format using readal (http://trimal.cgenomics.org). The best fitting model of DNA sequence evolution was determined using the AIC and BIC, as calculated by jModelTest ,. Since both the GTR + Γ model and the HKY + Γ model were favored by AIC and BIC respectively, MEGA6 was used to build bootstrapped maximum likelihood trees with 100 replicates under both models ,. The log likelihood of the GTR + Γ tree was slightly higher and is reported here; however, both trees exhibited a similar topology.
Protein structure prediction
Secondary and tertiary protein structures were modeled using RaptorX -. This software compares alignments of the sample protein to other proteins with known structural information to determine a probable structure using statistics. A representative CenH3 sequence from each genus was modeled (G.raimondii (D5), Oryza australiensis (EE), and Brassica napus (AC8)) (Figure 1).
cDNA generation and sequencing of CenH3transcripts
Reverse transcription was performed using the Invitrogen SuperScript III First-Strand Synthesis System Kit using oligo dT primers, and CenH3 was amplified from the cDNA pool using primers that were designed from the 5’ and 3’ outermost exons (Additional file 6). PCR products were processed using the Qiaquick PCR Cleanup columns and sequenced with the amplification primers (Additional file 6).
Chromatogram-based expression estimates were calculated as described previously . At least three replicates were used for each tissue to permit standard error calculations and paired, two-tailed T tests were used to test for significance. Expression levels for the polyploid accessions were secondarily estimated with RNA-seq data and by cloning cDNA amplicons (as described above). The clones were randomly selected from each sample, sequenced, and then grouped by their subgenomic origin “AT” and “DT”. Since the samples should follow a binomial distribution, the null hypothesis for the rate of cloning each homeologous copy of CenH3 should be 0.5. To control for the FWER (Family-wise Error Rate) at α = 0.05, the Bonferroni correction was determine the significance.
Expression estimation via RNA-seq
To assess CenH3 gene expression, we analyzed previously generated leaf transcriptome data (SRA BioProject PRJNA171342)  for both model diploid parents, an F1 hybrid of G. arboreum (A2 genome) and G. raimondii (D5 genome), a colchicine doubled F1 hybrid of G. arboreum and G. thurberi (D1 genome), and two accessions of the allopolyploid (AD genome) G. hirsutum (G. hirsutum var yucatanense, a wild accession; G. hirsutum cv Maxxa, a domesticated accession). Raw reads were trimmed with sickle (https://github.com/najoshi/sickle), and mapped to the generated CenH3 sequences using GSNAP (batch 4, novel splicing on)  in conjunction with a CenH3-specific SNP index to efficiently map sequences from different species and subgenomes. The SNP indices were manually curated from Sanger sequencing of the cDNA’s and gene sequences. For sequences from the hybrid and polyploid, PolyCat  was used to partition A- and D-genome derived reads. The significance for homeolog bias was calculated using a paired, Student’s T-test with log2 transformation to ensure the normality in expression values.
Availability of supporting data
CenH3 gene sequence data have been submitted to GenBank. Accession numbers can be found in Additional file 1.
The authors thank Patrick Edger for comments and suggestions with Brassica, Guanjing Hu for help with the Gossypium experimental work, and Matt Hufford for assistance with the evolutionary analyses. Research was supported by the NSF Plant Genome Program and by Cotton Incorporated.
- Hirsch CD, Jiang J: Centromeres: sequences, structure, and biology plant genome diversity. Plant Genome Diversity. Edited by: Wendel JF, Greilhuber J, Dolezel J, Leitch IJ. 2012, Springer, Vienna, 59-70. 10.1007/978-3-7091-1130-7_4.View ArticleGoogle Scholar
- Henikoff S, Ahmad K, Malik HS: The centromere paradox: stable inheritance with rapidly evolving DNA. Science. 2001, 293 (5532): 1098-1102. 10.1126/science.1062939.View ArticlePubMedGoogle Scholar
- Willard HF, Waye JS: Chromosome-specific subsets of human alpha satellite DNA: Analysis of sequence divergence within and between chromosomal subsets and evidence for an ancestral pentameric repeat. J Mol Evol. 1987, 25 (3): 207-214. 10.1007/BF02100014.View ArticlePubMedGoogle Scholar
- Jiang J, Birchler JA, Parrott WA, Kelly Dawe R: A molecular view of plant centromeres. Trends Plant Sci. 2003, 8 (12): 570-575. 10.1016/j.tplants.2003.10.011.View ArticlePubMedGoogle Scholar
- Melters DP, Bradnam KR, Young HA, Telis N, May MR, Ruby JG, Sebra R, Peluso P, Eid J, Rank D: Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution. Genome Biol. 2013, 14 (1): R10-10.1186/gb-2013-14-1-r10.PubMed CentralView ArticlePubMedGoogle Scholar
- Meraldi P, McAinsh AD, Rheinbay E, Sorger PK: Phylogenetic and structural analysis of centromeric DNA and kinetochore proteins. Genome Biol. 2006, 7 (3): R23-10.1186/gb-2006-7-3-r23.PubMed CentralView ArticlePubMedGoogle Scholar
- Cleveland DW, Mao Y, Sullivan KF: Centromeres and kinetochores: from epigenetics to mitotic checkpoint signaling. Cell. 2003, 112 (4): 407-421. 10.1016/S0092-8674(03)00115-6.View ArticlePubMedGoogle Scholar
- Malik HS, Vermaak D, Henikoff S: Recurrent evolution of DNA-binding motifs in the Drosophila centromeric histone. Proc Natl Acad Sci U S A. 2002, 99 (3): 1449-1454. 10.1073/pnas.032664299.PubMed CentralView ArticlePubMedGoogle Scholar
- Talbert PB, Masuelli R, Tyagi AP, Comai L, Henikoff S: Centromeric localization and adaptive evolution of an Arabidopsis histone H3 variant. Plant Cell. 2002, 14 (5): 1053-1066. 10.1105/tpc.010425.PubMed CentralView ArticlePubMedGoogle Scholar
- Malik HS, Henikoff S: Adaptive evolution of CID, a centromere-specific histone in Drosophila . Genetics. 2001, 157 (3): 1293-1298.PubMed CentralPubMedGoogle Scholar
- Allshire RC, Karpen GH: Epigenetic regulation of centromeric chromatin: old dogs, new tricks?. Nat Rev Genet. 2008, 9 (12): 923-937. 10.1038/nrg2466.PubMed CentralView ArticlePubMedGoogle Scholar
- Song C, Liu S, Xiao J, He W, Zhou Y, Qin Q, Zhang C, Liu Y: Polyploid organisms. Sci China Life Sci. 2012, 55 (4): 301-311. 10.1007/s11427-012-4310-2.View ArticlePubMedGoogle Scholar
- Heslop-Harrison JS: Genome evolution: extinction, continuation or explosion?. Curr Opin Plant Biol. 2012, 15 (2): 115-121. 10.1016/j.pbi.2012.03.006.View ArticlePubMedGoogle Scholar
- Van de Peer Y, Maere S, Meyer A: The evolutionary significance of ancient genome duplications. Nat Rev Genet. 2009, 10 (10): 725-732. 10.1038/nrg2600.View ArticlePubMedGoogle Scholar
- Wendel JF: Genome evolution in polyploids. Plant Mol Biol. 2000, 42 (1): 225-249. 10.1023/A:1006392424384.View ArticlePubMedGoogle Scholar
- Soltis DE, Albert VA, Leebens-Mack J, Bell CD, Paterson AH, Zheng C, Sankoff D, de Pamphilis CW, Wall PK, Soltis PS: Polyploidy and angiosperm diversification. Am J Bot. 2009, 96 (1): 336-348. 10.3732/ajb.0800079.View ArticlePubMedGoogle Scholar
- Adams KL, Wendel JF: Polyploidy and genome evolution in plants. Curr Opin Plant Biol. 2005, 8 (2): 135-141. 10.1016/j.pbi.2005.01.001.View ArticlePubMedGoogle Scholar
- Adams KL, Cronn R, Percifield R, Wendel JF: Genes duplicated by polyploidy show unequal contributions to the transcriptome and organ-specific reciprocal silencing. Proc Natl Acad Sci U S A. 2003, 100 (8): 4649-4654. 10.1073/pnas.0630618100.PubMed CentralView ArticlePubMedGoogle Scholar
- Wendel JF, Flagel LE, Adams KL: Jeans, genes, and genomes: cotton as a model for studying polyploidy. Polyploidy and Genome Evolution. Edited by: Soltis PS, Soltis DE. 2012, Springer, New York, 181-207. 10.1007/978-3-642-31442-1_10.View ArticleGoogle Scholar
- Wendel JF, Doyle JJ: Polyploidy and evolution in plants. Curr Opin Plant Biol. 2005, 8: 135-141. 10.1016/j.pbi.2005.01.001.View ArticlePubMedGoogle Scholar
- Soltis DE, Soltis PS: Polyploidy: recurrent formation and genome evolution. Trends Eco Evol. 1999, 14 (9): 348-352. 10.1016/S0169-5347(99)01638-9.View ArticleGoogle Scholar
- Chen ZJ, Ni Z: Mechanisms of genomic rearrangements and gene expression changes in plant polyploids. Bioessays. 2006, 28 (3): 240-252. 10.1002/bies.20374.View ArticlePubMedGoogle Scholar
- Chen ZJ: Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu Rev Plant Biol. 2007, 58: 377-406. 10.1146/annurev.arplant.58.032806.103835.PubMed CentralView ArticlePubMedGoogle Scholar
- Adams KL, Wendel JF: Novel patterns of gene expression in polyploid plants. Trends Genet. 2005, 21 (10): 539-543. 10.1016/j.tig.2005.07.009.View ArticlePubMedGoogle Scholar
- Scannell DR, Byrne KP, Gordon JL, Wong S, Wolfe KH: Multiple rounds of speciation associated with reciprocal gene loss in polyploid yeasts. Nature. 2006, 440 (7082): 341-345. 10.1038/nature04562.View ArticlePubMedGoogle Scholar
- Town CD, Cheung F, Maiti R, Crabtree J, Haas BJ, Wortman JR, Hine EE, Althoff R, Arbogast TS, Tallon LJ: Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy. Plant Cell Online. 2006, 18 (6): 1348-1359. 10.1105/tpc.106.041665.View ArticleGoogle Scholar
- Ku H-M, Vision T, Liu J, Tanksley SD: Comparing sequenced segments of the tomato and Arabidopsis genomes: large-scale duplication followed by selective gene loss creates a network of synteny. Proc Natl Acad Sci U S A. 2000, 97 (16): 9121-9126. 10.1073/pnas.160271297.PubMed CentralView ArticlePubMedGoogle Scholar
- Doyle JJ, Flagel LE, Paterson AH, Rapp RA, Soltis DE, Soltis PS, Wendel JF: Evolutionary genetics of genome merger and doubling in plants. Ann Rev Genet. 2008, 42 (1): 443-461. 10.1146/annurev.genet.42.110807.091524.View ArticlePubMedGoogle Scholar
- Barker MS, Kane NC, Matvienko M, Kozik A, Michelmore RW, Knapp SJ, Rieseberg LH: Multiple paleopolyploidizations during the evolution of the Compositae reveal parallel patterns of duplicate gene retention after millions of years. Mol Biol Evol. 2008, 25 (11): 2445-2455. 10.1093/molbev/msn187.PubMed CentralView ArticlePubMedGoogle Scholar
- Hirsch CD, Wu Y, Yan H, Jiang J: Lineage-specific adaptive evolution of the centromeric protein CenH3 in diploid and allotetraploid Oryza species. Mol Biol Evol. 2009, 26 (12): 2877-2885. 10.1093/molbev/msp208.View ArticlePubMedGoogle Scholar
- Lermontova I, Schubert I: CenH3 for establishing and maintaining centromeres. In: Plant Centromere Biology. Oxford, UK: Wiley-Blackwell; 2013: 67–82.View ArticleGoogle Scholar
- Wang G, He Q, Liu F, Cheng Z, Talbert P, Jin W: Characterization of CenH3 proteins and centromere-associated DNA sequences in diploid and allotetraploid Brassica species. Chromosoma. 2011, 120 (4): 353-365. 10.1007/s00412-011-0315-z.View ArticlePubMedGoogle Scholar
- U. N: Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Japan J Bot. 1935, 7: 389-452.Google Scholar
- Ge S, Sang T, Lu B-R, Hong D-Y: Phylogeny of rice genomes with emphasis on origins of allotetraploid species. Proc Natl Acad Sci U S A. 1999, 96 (25): 14400-14405. 10.1073/pnas.96.25.14400.PubMed CentralView ArticlePubMedGoogle Scholar
- Grover C, Grupp K, Wanzek R, Wendel J: Assessing the monophyly of polyploid Gossypium species. Plant Syst Evol. 2012, 298 (6): 1177-1183. 10.1007/s00606-012-0615-7.View ArticleGoogle Scholar
- Luo S, Mach J, Abramson B, Ramirez R, Schurr R, Barone P, Copenhaver G, Folkerts O: The cotton centromere contains a Ty3-gypsy-like LTR retroelement. PLoS One. 2012, 7 (4): e35261-10.1371/journal.pone.0035261.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.View ArticlePubMedGoogle Scholar
- Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Pond SLK: Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012, 8 (7): e1002764-10.1371/journal.pgen.1002764.PubMed CentralView ArticlePubMedGoogle Scholar
- Hui L, Lu L, Heng Y, Qin R, Xing Y, Jin W: Expression of CENH3 alleles in synthesized allopolyploid Oryza species. J Genet Genomics. 2010, 37 (10): 703-711. 10.1016/S1673-8527(09)60088-6.View ArticleGoogle Scholar
- Salmon A, Flagel L, Ying B, Udall JA, Wendel JF: Homoeologous nonreciprocal recombination in polyploid cotton. New Phytol. 2010, 186 (1): 123-134. 10.1111/j.1469-8137.2009.03093.x.View ArticlePubMedGoogle Scholar
- Flagel L, Wendel J, Udall J: Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton. BMC Genomics. 2012, 13 (1): 302-10.1186/1471-2164-13-302.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang T, Talbert PB, Zhang W, Wu Y, Yang Z, Henikoff JG, Henikoff S, Jiang J: The CentO satellite confers translational and rotational phasing on CenH3 nucleosomes in rice centromeres. Proc Natl Acad Sci U S A. 2013, 110 (50): E4875-E4883. 10.1073/pnas.1319548110.PubMed CentralView ArticlePubMedGoogle Scholar
- Takeuchi K, Nishino T, Mayanagi K, Horikoshi N, Osakabe A, Tachiwana H, Hori T, Kurumizaka H, Fukagawa T: The centromeric nucleosome-like CENP–T–W–S–X complex induces positive supercoils into DNA. Nucleic Acids Res. 2014, 42 (3): 1644-1655. 10.1093/nar/gkt1124.PubMed CentralView ArticlePubMedGoogle Scholar
- Carroll CW, Milks KJ, Straight AF: Dual recognition of CENP-A nucleosomes is required for centromere assembly. J Cell Biol. 2010, 189 (7): 1143-1155. 10.1083/jcb.201001013.PubMed CentralView ArticlePubMedGoogle Scholar
- Vermaak D, Hayden HS, Henikoff S: Centromere targeting element within the histone fold domain of Cid. Mol Cell Biol. 2002, 22 (21): 7553-7561. 10.1128/MCB.22.21.7553-7561.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper JL, Henikoff S: Adaptive evolution of the histone fold domain in centromeric histones. Mol Biol Evol. 2004, 21 (9): 1712-1718. 10.1093/molbev/msh179.View ArticlePubMedGoogle Scholar
- Shelby RD, Vafa O, Sullivan KF: Assembly of CENP-A into centromeric chromatin requires a cooperative array of nucleosomal DNA contact sites. J Cell Biol. 1997, 136 (3): 501-513. 10.1083/jcb.136.3.501.PubMed CentralView ArticlePubMedGoogle Scholar
- Lermontova I, Schubert V, Fuchs J, Klatte S, Macas J, Schubert I: Loading of Arabidopsis centromeric histone CENH3 occurs mainly during G2 and requires the presence of the histone fold domain. Plant Cell. 2006, 18 (10): 2443-2451. 10.1105/tpc.106.043174.PubMed CentralView ArticlePubMedGoogle Scholar
- Black BE, Foltz DR, Chakravarthy S, Luger K, Woods VL, Cleveland DW: Structural determinants for generating centromeric chromatin. Nature. 2004, 430 (6999): 578-582. 10.1038/nature02766.View ArticlePubMedGoogle Scholar
- Lermontova I, Koroleva O, Rutten T, Fuchs J, Schubert V, Moraes I, Koszegi D, Schubert I: Knockdown of CenH3 in Arabidopsis reduces mitotic divisions and causes sterility by disturbed meiotic chromosome segregation. Plant J. 2011, 68 (1): 40-50. 10.1111/j.1365-313X.2011.04664.x.View ArticlePubMedGoogle Scholar
- Ravi M, Shibata F, Ramahi JS, Nagaki K, Chen C, Murata M, Chan SW: Meiosis-specific loading of the centromere-specific histone CENH3 in Arabidopsis thaliana . PLoS Genet. 2011, 7 (6): e1002121-10.1371/journal.pgen.1002121.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Baker RE, Keith KC, Harris K, Stoler S, Fitzgerald-Hayes M: The N terminus of the centromere H3-like protein Cse4p performs an essential function distinct from that of the histone fold domain. Mol Cell Biol. 2000, 20 (18): 7037-7048. 10.1128/MCB.20.18.7037-7048.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey AO, Panchenko T, Sathyan KM, Petkowski JJ, Pai P-J, Bai DL, Russell DH, Macara IG, Shabanowitz J, Hunt DF, Black BE, Foltz DR: Posttranslational modification of CENP-A influences the conformation of centromeric chromatin. Proc Natl Acad Sci U S A 2013, 110(29):11827–11832.PubMed CentralView ArticlePubMedGoogle Scholar
- Seelanan T, Schnabel A, Wendel JF: Congruence and consensus in the cotton tribe. Sys Bot. 1997, 22: 259-290. 10.2307/2419457.View ArticleGoogle Scholar
- Wang K, Wang Z, Li F, Ye W, Wang J, Song G, Yue Z, Cong L, Shang H, Zhu S: The draft genome of a diploid cotton Gossypium raimondii . Nat Genet. 2012, 44 (10): 1098-1103. 10.1038/ng.2371.View ArticlePubMedGoogle Scholar
- Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, Mitros T, Dirks W, Hellsten U, Putnam N: Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012, 40 (D1): D1178-D1186. 10.1093/nar/gkr944.PubMed CentralView ArticlePubMedGoogle Scholar
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2008, 36 (suppl 1): D25-D30.PubMed CentralPubMedGoogle Scholar
- Delport W, Poon AF, Frost SD, Pond SLK: Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics. 2010, 26 (19): 2455-2457. 10.1093/bioinformatics/btq429.PubMed CentralView ArticlePubMedGoogle Scholar
- Lysak MA, Koch MA, Pecinka A, Schubert I: Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005, 15 (4): 516-525. 10.1101/gr.3531105.PubMed CentralView ArticlePubMedGoogle Scholar
- Cronn RC, Small RL, Haselkorn T, Wendel JF: Rapid diversification of the cotton genus (Gossypium: Malvaceae) revealed by analysis of sixteen nuclear and chloroplast genes. Am J Botany. 2002, 89 (4): 707-725. 10.3732/ajb.89.4.707.View ArticleGoogle Scholar
- Sanderson M: A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997, 14 (12): 1218-10.1093/oxfordjournals.molbev.a025731.View ArticleGoogle Scholar
- Jacquemin J, Ammiraju JSS, Haberer G, Billheimer DD, Yu Y, Liu LC, Rivera LF, Mayer K, Chen M, Wing RA: Fifteen million years of evolution in the Oryza genus shows extensive gene family expansion. Mol Plant 2013, 7(4):642–656.View ArticlePubMedGoogle Scholar
- Senchina DS, Alvarez I, Cronn RC, Liu B, Rong J, Noyes RD, Paterson AH, Wing RA, Wilkins TA, Wendel JF: Rate variation among nuclear genes and the age of polyploidy in Gossypium . Mol Biol Evol. 2003, 20 (4): 633-643. 10.1093/molbev/msg065.View ArticlePubMedGoogle Scholar
- Lyons E, Freeling M: How to usefully compare homologous plant genes and chromosomes as DNA sequences. Plant J. 2008, 53 (4): 661-673. 10.1111/j.1365-313X.2007.03326.x.View ArticlePubMedGoogle Scholar
- Lyons E, Pedersen B, Kane J, Alam M, Ming R, Tang H, Wang X, Bowers J, Paterson A, Lisch D: Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with rosids. Plant Physiol. 2008, 148 (4): 1772-1781. 10.1104/pp.108.124867.PubMed CentralView ArticlePubMedGoogle Scholar
- Carvalho MR, Herrera FA, Jaramillo CA, Wing SL, Callejas R: Paleocene Malvaceae from northern South America and their biogeographical implications. Am J Bot. 2011, 98 (8): 1337-1355. 10.3732/ajb.1000539.View ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999, 95-98.Google Scholar
- Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012, 9 (8): 772-10.1038/nmeth.2109.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Kumar S, Nei M, Dudley J, Tamura K: MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9 (4): 299-306. 10.1093/bib/bbn017.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013, 30 (12): 2725-2729. 10.1093/molbev/mst197.PubMed CentralView ArticlePubMedGoogle Scholar
- Ma J, Wing RA, Bennetzen JL, Jackson SA: Plant centromere organization: a dynamic structure with conserved functions. Trends Genet. 2007, 23 (3): 134-139. 10.1016/j.tig.2007.01.004.View ArticlePubMedGoogle Scholar
- Källberg M, Wang H, Wang S, Peng J, Wang Z, Lu H, Xu J: Template-based protein structure modeling using the RaptorX web server. Nat Protocols. 2012, 7 (8): 1511-1522. 10.1038/nprot.2012.085.View ArticlePubMedGoogle Scholar
- Peng J, Xu J: A multiple-template approach to protein threading. Proteins. 2011, 79 (6): 1930-1939. 10.1002/prot.23016.PubMed CentralView ArticlePubMedGoogle Scholar
- Peng J, Xu J: Raptorx: exploiting structure information for protein alignment by statistical inference. Proteins. 2011, 79 (S10): 161-171. 10.1002/prot.23175.PubMed CentralView ArticlePubMedGoogle Scholar
- Adams KL, Percifield R, Wendel JF: Organ-specific silencing of duplicated genes in a newly synthesized cotton allotetraploidd. Genetics. 2004, 168 (4): 2217-2226. 10.1534/genetics.104.033522.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoo MJ, Szadkowski E, Wendel JF: Homoeolog expression bias and expression level dominance in allopolyploid cotton. Heredity. 2013, 110 (2): 171-180. 10.1038/hdy.2012.94.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu TD, Nacu S: Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010, 26 (7): 873-881. 10.1093/bioinformatics/btq057.PubMed CentralView ArticlePubMedGoogle Scholar
- Page JT, Gingle AR, Udall JA: PolyCat: a resource for genome categorization of sequencing reads from allopolyploid organisms. G3. 2013, 3 (3): 517-525. 10.1534/g3.112.005298.PubMed CentralView ArticlePubMedGoogle 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.