Genomic, regulatory and epigenetic mechanisms underlying duplicated gene evolution in the natural allotetraploid Oryza minuta
© Sui et al.; licensee BioMed Central Ltd. 2014
Received: 25 August 2013
Accepted: 30 December 2013
Published: 6 January 2014
Polyploid species contribute to Oryza diversity. However, the mechanisms underlying gene and genome evolution in Oryza polyploids remain largely unknown. The allotetraploid Oryza minuta, which is estimated to have formed less than one million years ago, along with its putative diploid progenitors (O. punctata and O. officinalis), are quite suitable for the study of polyploid genome evolution using a comparative genomics approach.
Here, we performed a comparative study of a large genomic region surrounding the Shattering4 locus in O. minuta, as well as in O. punctata and O. officinalis. Duplicated genomes in O. minuta have maintained the diploid genome organization, except for several structural variations mediated by transposon movement. Tandem duplicated gene clusters are prevalent in the Sh4 region, and segmental duplication followed by random deletion is illustrated to explain the gene gain-and-loss process. Both copies of most duplicated genes still persist in O. minuta. Molecular evolution analysis suggested that these duplicated genes are equally evolved and mostly manipulated by purifying selection. However, cDNA-SSCP analysis revealed that the expression patterns were dramatically altered between duplicated genes: nine of 29 duplicated genes exhibited expression divergence in O. minuta. We further detected one gene silencing event that was attributed to gene structural variation, but most gene silencing could not be related to sequence changes. We identified one case in which DNA methylation differences within promoter regions that were associated with the insertion of one hAT element were probably responsible for gene silencing, suggesting a potential epigenetic gene silencing pathway triggered by TE movement.
Our study revealed both genetic and epigenetic mechanisms involved in duplicated gene silencing in the allotetraploid O. minuta.
Rice is one of the most important crops for human consumption, as it feeds more than half of the world’s population. To facilitate rice improvement, wild relative species in the genus Oryza have been employed as excellent genetic resources for rice breeding and genetic modification[1, 2].
The genus Oryza, comprising two cultivated and approximately 22 wild species, is classified into 10 genome types, including six diploids (AA, BB, CC, EE, FF and GG) and four allotetraploids (BBCC, CCDD, HHJJ and KKLL)[3, 4]. To better exploit the superior wild rice genetic resources, a robust analysis of phylogeny among Oryza species was performed several years ago. Subsequently, efforts have focused on deciphering the evolutionary relationship among diploid Oryza species[5–7].
However, within the genus Oryza, almost one-third of rice species are considered to be allotetraploids, representing a large part of species diversity present in this genus. Elucidating the evolutionary history of allotetraploids in Oryza is highly important for obtaining a complete understanding of the evolution of Oryza. Unfortunately, the potential progenitors of only a few species with the BBCC genome type have been identified[4, 8]. Of these species, O. minuta was selected as the representative species of BBCC in the Oryza Map Alignment Project (OMAP). Comparative genomics resources are available, including a high quality bacterial artificial chromosome (BAC) library, BAC end sequences and a BAC-based physical map. The genome donors of O. minuta were identified as diploid O. punctata (BB) and O. officinalis (CC), although some studies suggest that an extinct Asian BB genome carrier is the direct genome donor. Several studies have deduced the molecular timing of BBCC formation, indicating that the formation of the allotetraploid occurred ~0.3 to 0.6 million years ago (Mya)[8, 10, 11].
Polyploidy and the consequences of duplicated genomes have been extensively studied in some model species[12–17]. However, few studies have investigated microstructural variations using a comparative genomics approach[18–20]. The evolutionary fate of duplicated genes has also been well-studied, and expression analysis has often allowed gene silencing to be detected[21–25]. However, the genetic and epigenetic regulatory pathways of gene expression divergence are still largely unknown.
Shattering4 (Sh4), a major quantitative trait locus responsible for rice grain shattering, which encodes a transcription factor with an MYB3 DNA binding domain, plays an important role in the establishment of the abscission layer. An amino acid substitution in the Sh4 coding region affects the normal development of abscission between the grain and the pedicel and further reduces grain shattering. Human selection of this single substitution promoted the domestication of rice from wild species. To deepen our understanding of the evolution of O. minuta, we conducted a comparative genomic analysis of a genomic segment surrounding the Sh4 locus among O. minuta, O. punctata and O. officinalis.
Sequencing and annotation of the Sh4-orthologous regions
Sequenced BAC clones from the Sh4 region
BAC length (bp)
Total sequence length (bp)*
O. minuta (BB)
O. minuta (BB)
O. minuta (BB)
O. minuta (CC)
O. minuta (CC)
O. minuta (CC)
Summary of the ratio of LTRs and solo LTRs between the allotetraploid and diploids
Pair-wise genome comparison
Intact solo LTR
Intact LTR_RT: Solo LTR (All)
BB VS BBCC_BB
Only present in BBCC_BB
Only present in BB
BBCC_CC VS CC
Only present in BBCC_CC
Only present in CC
BBCC_BB VS BBCC_CC
Only present in BBCC_BB
Only present in BBCC_CC
Structural variation after allotetraploid formation
Duplicated gene evolution in O. minuta
To investigate duplicated gene evolution, we chose genes that were covered by BAC sequencing in all of the genomes, and therefore genes before 17 and after 53 were excluded from the following analysis (Additional file1: Table S2). The Sh4 regions contain 41 sets of orthologous genes (from Gene 17 to 53). Among these genes, only one deletion of Gene 29–2 was observed (represented with purple bars in Figure 1) in BBCC_CC, and both copies of the other duplicated genes were maintained. We also identified four genes that were putatively pseudogenized in the polyploid genome, which formed after polyploidy (Additional file1: Figure S5).
We estimated the molecular evolutionary rate for duplicated genes and examined whether these genes had different rates of evolution after polyploid formation (Additional file1: Table S5). The results from the relative rate tests suggest that most genes do not exhibit obviously different rates between the polyploid and diploids. To determine which type of selection was acting upon the duplicated genes, we calculated the ratio of non-synonymous (Ka) to synonymous (Ks) substitutions of protein coding sequences. Most duplicated genes were under purifying selection (Ka/Ks < 1), indicating that these genes are still strongly controlled after polyploidy (Additional file1: Table S6). Similar patterns of Ks distribution between BB-BC_B and CC-BC_C also indicate that the evolutionary rates of duplicated genes were not obviously different (Additional file1: Figure S6). We used duplicated genes to deduce the molecular timing of O. minuta divergence, and the results suggest that O. minuta was formed approximately 0.8–1.0 MYA (Additional file1: Table S6), which is older than a previous estimate obtained by examining the Moc1 region (~0.4 MYA).
Gene expression divergence of duplicated genes
To investigate whether genetic or epigenetic mechanisms regulate these genes, we first examined the genomic sequences of silenced copies of these nine genes to determine whether or how the gene structures were destroyed. The silencing in only one gene (Gene 43, whose CC subgenome copy was silenced) could be attributed to sequence variations (LTR insertion). In addition, we examined the expression patterns of two other genes (Gene 22 and Gene 26) that were identified as pseudogenes and found that both copies of the genes could be transcribed (coexpressed) (Additional file1: Table S7). In summary, of the nine genes that exhibited expression silencing, only one was under genetic regulation. This result implies that other mechanisms, such as epigenetic regulation, may be involved in gene silencing.
Although several studies on synthetic and natural polyploids have provided evidence for rapid loss and gain of genomic segments (including genes) and extensive genomic reshuffling[14, 16, 25, 32, 33], sequence comparisons in the Sh4 genomic region, combined with the results from analysis of Adh1 and Moc1[8, 11], suggest that the natural Oryza allotetraploid O. minuta is perhaps a relatively stable polyploid. Such stabilization is supported by the presence of conserved genes, intergenic regions and even shared TEs between polyploid and diploid genomes. However, confirming this notion would require additional investigations of other larger segments or even whole genomes.
The regulation of duplicated gene expression in polyploids has been well-studied in several model species, but few of these studies have correlated expression divergence with sequence variations[10, 23, 34–36]. In this study, we found that two pairs of duplicated genes annotated as pseudogenes (22 and 26) could be coexpressed in a cDNA-SSCP assay. Analysis of cDNA sequences has indicated that transcripts from pseudogenes are non-functional. Therefore, coexpression of homoelogous genes revealed by microarray analysis (or other methods) does not guarantee that both gene copies are treated equally by the genome, as mRNA sequence variations cannot be detected using these approaches. Notably, the biological significance of pseudogenes has recently been examined, especially pseudogenes that can be transcribed[37, 38]. These findings suggest that pseudogenes can probably evolve from being buried in huge genomes to becoming new, functional elements, implying that pseudogenization can lead to neofunctionalization.
Although genome-wide experimental data are required to prove that the epigenetic pathway is the universal mechanism for gene silencing, it is currently more difficult to test epigenetic markers than genetic sequences, especially when dealing with polyploids. Nonetheless, we still propose that DNA methylation-controlled gene silencing is a prevalent mechanism, based on several observations. First, TEs have been reported to regulate gene expression for numerous genes[39, 40]. For instance, a SINE element inserted upstream of the FWA coding sequence caused this gene to be epigenetically silenced in vegetative tissues of Arabidopsis. Second, numerous TEs are likely to be inserted adjacent to genes in the BBCC genome. TEs can affect gene expression, but the effects of TEs decrease with increasing distance between the TEs and genes. We investigated the distribution of several types of DNA-type TEs in the japonica genome. These TEs have the greatest potential to regulate nearby genes due to their preferential insertion near or within genes. Three types of TEs, including hAT, Stowaway and Tourist, are all present in thousands of copies (6,728, 49,810 and 40,092, respectively), and the mean distance of these TEs to genes is within 2 kb (data not shown). Therefore, we postulate that a considerable number of potential epigenetic triggers have also been buried within the BBCC genome. More importantly, according to previous reports and the results of the current study, TE insertion is not sufficient to initiate the silencing pathway; siRNA or DNA methylation is essential for initiating this program. For example, the same hAT elements were found in regulatory regions of the FLC gene in two Arabidopsis ecotypes, Landsberg erecta (Ler) and Columbia (Col). However, the roles of these two hAT elements in regulating FLC gene expression are quite different, as ~24-nt siRNA was found at higher levels in Ler than in Col; this siRNA can mediate DNA methylation and gene silencing in the Ler ecotype. Recently, the global effects of TEs on gene expression were investigated in Arabidopsis and its close relative. Genome-wide analysis indicated that TEs can affect the expression of nearby genes, especially when the TEs are epigenetically modified (methylated or siRNA targeted)[42, 44]. We also calculated the methylation rates of the above three types of TEs in the japonica genome and found that methylated copies account for over 80% of these TEs, providing further support for the potential role of silenced TEs in gene regulation. Here, we used the japonica genome to represent O. minuta. Future studies should focus on comprehensive analysis of the O. minuta genome to help elucidate the epigenetic regulatory pathway on a genome-wide scale.
By integrating comparative genomic tools, gene expression and epigenetic analysis, our study comprehensively demonstrates how duplicated genomes and genes evolve within the Sh4 region in O. minuta. We found that duplicated genes are under both genetic and epigenetic regulation, and DNA methylation is proposed as a potentially important regulatory mechanism for gene silencing.
Plant materials and BAC library
Seeds and seedlings of Oryza minuta (Accession No.101141), O. punctata (Accession No.105690) and O. officinalis (Accession No.100896) were obtained from the International Rice Resource Institute (IRRI, Philippines). High-density BAC library filters and BAC clones for three Oryza species were purchased from the Arizona Genomics Institute (USA).
BAC identification and sequencing
BAC clones covering the orthologous regions of the japonica Sh4 genome segment were identified by screening Oryza genomic BAC libraries following the method described by Lu et al.. Initial selections were conducted using two unique probes (designed with two japonica gene models, LOC_Os04g57350 and LOC_Os04g57600, which are located upstream and downstream of the Shattering 4 [LOC_Os04g57530] gene locus, respectively), to hybridize to high-density filters containing three Oryza genomic BAC libraries. Combined with physical map positions in Finger Printed Contigs (FPC), a total of 98 positive BAC clones were identified. All screened BAC clones were digested with Hind III, size-selected by electrophoresis and transferred onto nylon filters for Southern blot analysis. For diploid genomes, eight additional probes were used to identify the BAC clones, which maximized the orthologous region coverage and minimized the gaps between consecutive BAC clones. To distinguish between the subgenomes of O. minuta, the digested map of each BAC was compared with that of O. punctata and O. officinalis. The tetrapolyploid BAC clones were divided into two groups (each from one parental genome) based on their Hind III digestion patterns. Ten BAC clones of three Oryza species were sequenced with an ABI 3730 automated sequencer (Table 1). Specifically, Sanger reads were assembled from each BAC into contigs, and BAC sequences were then merged from the same genome by identifying overlapping sequences (Table 1). Orthologous genome regions in O. brachyantha (FF) were recently generated.
A comparative gene annotation approach was taken to identify gene models in the BB, CC and BBCC genomes. Before using the japonica Sh4 genome annotation as a reference gene model, transposon-related gene models and hypothetical genes without cDNA, ESTs and homologous proteins were excluded. Gene structures were confirmed using full-length japonica cDNA. To annotate non-japonica genomes, four wild rice genomes were initially repeat-masked and predicted using FGENESH (http://linux1.softberry.com). All predicted genes were aligned with japonica cDNA and proteins from rice, Sorghum, Brachypodium and Maize. Gene models without cDNA or protein supports were excluded. All gene models were manually refined based on the AA-BB-BBCC-CC multiple alignment framework. Conserved gene structures were modified based on japonica gene models. RT-PCR was used to detect gene structure if great variation existed between japonica genes and non-japonica genes. Multiple sequence comparisons were performed with CLUSTALW. Structural variation was detected using Artemis Comparison Tool (ACT).
To identify transposon elements, RepeatMasker was initially used to annotate the TEs and other repeat sequences, followed by manual analysis (http://www.repeatmasker.org/). LTR-type transposons were also predicted with LTR_STRUC, LTR_Finder and LTRharvest to complement the results from RepeatMasker. Intact structures and other TE signatures such as target site duplication (TSD), terminal inverted repeats (TIR), polypurine tracts (PPT), primer binding sites (PBS) and long terminal repeats (LTR) were manually identified using Dotter software and ACT. The insertion time of each LTR was estimated using the baseml program in PAML at a mutation rate of 1.3 × 10-8 per site per year.
Genomic structural variations, such as inversions, were also detected by ACT, followed by thorough analysis of their boundary sequences. Pack-MuLE elements were annotated manually, and to determine the captured genome sequence in the wild rice genome, the homologous sequences were searched against the japonica reference genome.
Molecular evolution analysis of duplicated genes
To determine the type of selection acting upon duplicated genes, Ka and Ks values were calculated for duplicated genes in O. minuta using the baseml program with the pairwise model in PAML version 4.6. Alignments for coding sequences of duplicated genes were conducted with CLUSTALW. Divergence times of duplicated genes were calculated with a synonymous substitution rate of 6.5 × 10-9 substitutions per site per year. Relative evolutionary rates of duplicated genes were estimated using the Tajima relative rate test implemented in MEGA v5.2.
Total RNA samples were isolated from six different tissues of three Oryza species, including mature leaves, mature roots, young leaves, young roots, flowers and mixed panicles. Reverse transcription was performed on mixed total RNA from all tissues, and cDNA products were used to amplify orthologous genes. SSCP analysis was then conducted with the Bio-RAD Dcode™ system following the standard protocol (Bio-RAD, USA). The SSCP results revealed the gene expression patterns of duplicated genes in O. minuta, including both coexpressed and gene silencing. All primer sequences are listed in Additional file2.
Methylation-specific PCR analysis
DNA was extracted from mature O. punctata, O. officinalis and O. minuta leaves, digested with restriction enzymes and treated with sodium bisulfite to convert the unmethylated cytosine residues to uracil. A set of primers was designed to amplify the genomic region of each gene from the end of the first exon to ~2 kb upstream of the gene; the size of each PCR product was approximately 150–300 bp. Primer sequences are listed in Additional file2. The products were recovered, cloned and sequenced for DNA methylation analysis.
Availability of supporting data
Sh4 genome sequences for wild rices can be downloaded from NCBI at http://www.ncbi.nlm.nih.gov/nuccore/HQ827834 (BB), http://www.ncbi.nlm.nih.gov/nuccore/HQ827835 (BC_B), http://www.ncbi.nlm.nih.gov/nuccore/HQ827836 (BC_C) and http://www.ncbi.nlm.nih.gov/nuccore/HQ827837 (CC). O. sativa ssp. japonica sequences can be downloaded from http://rice.plantbiology.msu.edu/. FF genome can be downloaded from http://www.gramene.org/Oryza_brachyantha/Info/Index.
This work was supported by the National Natural Science Foundation of China (grants # 31171231 and 30770143) to MSC.
- Vaughan DA, Morishima H, Kadowaki K: Diversity in the Oryza genus. Curr Opin Plant Biol. 2003, 6 (2): 139-146. 10.1016/S1369-5266(03)00009-8.PubMedView ArticleGoogle Scholar
- Wing RA, Ammiraju JS, Luo M, Kim H, Yu Y, Kudrna D, Goicoechea JL, Wang W, Nelson W, Rao K: The Oryza map alignment project: the golden path to unlocking the genetic potential of wild rice species. Plant Mol Biol. 2005, 59 (1): 53-62. 10.1007/s11103-004-6237-x.PubMedView ArticleGoogle Scholar
- Aggarwal R, Brar D, Nandi S, Huang N, Khush G: Phylogenetic relationships among Oryza species revealed by AFLP markers. Theor Appl Genet. 1999, 98 (8): 1320-1328. 10.1007/s001220051198.View ArticleGoogle 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. 1999, 96 (25): 14400-14405. 10.1073/pnas.96.25.14400.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhu Q, Ge S: Phylogenetic relationships among A-genome species of the genus Oryza revealed by intron sequences of four nuclear genes. New Phytol. 2005, 167 (1): 249-265. 10.1111/j.1469-8137.2005.01406.x.PubMedView ArticleGoogle Scholar
- Zhu Q, Zheng X, Luo J, Gaut BS, Ge S: Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice. Mol Biol Evol. 2007, 24 (3): 875-888.PubMedView ArticleGoogle Scholar
- Zou XH, Zhang FM, Zhang JG, Zang LL, Tang L, Wang J, Sang T, Ge S: Analysis of 142 genes resolves the rapid diversification of the rice genus. Genome Biol. 2008, 9 (3): R49-10.1186/gb-2008-9-3-r49.PubMed CentralPubMedView ArticleGoogle Scholar
- Ammiraju JSS, Fan C, Yu Y, Song X, Cranston KA, Pontaroli AC, Lu F, Sanyal A, Jiang N, Rambo T: Spatio-temporal patterns of genome evolution in allotetraploid species of the genus Oryza. Plant J. 2010, 63 (3): 430-442. 10.1111/j.1365-313X.2010.04251.x.PubMedView ArticleGoogle Scholar
- Kim H, Hurwitz B, Yu Y, Collura K, Gill N, SanMiguel P, Mullikin JC, Maher C, Nelson W, Wissotski M: Construction, alignment and analysis of twelve framework physical maps that represent the ten genome types of the genus Oryza. Genome Biol. 2008, 9 (2): R45-10.1186/gb-2008-9-2-r45.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang B, Ding Z, Liu W, Pan J, Li C, Ge S, Zhang D: Polyploid evolution in Oryza officinalis complex of the genus Oryza. BMC Evol Biol. 2009, 9 (1): 250-10.1186/1471-2148-9-250.PubMed CentralPubMedView ArticleGoogle Scholar
- Lu F, Ammiraju JS, Sanyal A, Zhang S, Song R, Chen J, Li G, Sui Y, Song X, Cheng Z: Comparative sequence analysis of MONOCULM1-orthologous regions in 14 Oryza genomes. Proc Natl Acad Sci. 2009, 106 (6): 2071-2076. 10.1073/pnas.0812798106.PubMed CentralPubMedView ArticleGoogle Scholar
- Feldman M, Liu B, Segal G, Abbo S, Levy AA, Vega JM: Rapid elimination of low-copy DNA sequences in polyploid wheat: a possible mechanism for differentiation of homoeologous chromosomes. Genetics. 1997, 147 (3): 1381-1387.PubMed CentralPubMedGoogle Scholar
- Soltis DE, Soltis PS: The dynamic nature of polyploid genomes. Proc Natl Acad Sci. 1995, 92 (18): 8089-8091. 10.1073/pnas.92.18.8089.PubMed CentralPubMedView ArticleGoogle Scholar
- Song K, Lu P, Tang K, Osborn TC: Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution. Proc Natl Acad Sci. 1995, 92 (17): 7719-7723. 10.1073/pnas.92.17.7719.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiong Z, Gaeta RT, Pires JC: Homoeologous shuffling and chromosome compensation maintain genome balance in resynthesized allopolyploid Brassica napus. Proc Natl Acad Sci. 2011, 108 (19): 7908-7913. 10.1073/pnas.1014138108.PubMed CentralPubMedView ArticleGoogle Scholar
- Chester M, Gallagher JP, Symonds VV, Da Silva AVC, Mavrodiev EV, Leitch AR, Soltis PS, Soltis DE: Extensive chromosomal variation in a recently formed natural allopolyploid species, Tragopogon miscellus (Asteraceae). Proc Natl Acad Sci. 2012, 109 (4): 1176-1181. 10.1073/pnas.1112041109.PubMed CentralPubMedView ArticleGoogle Scholar
- Gaeta RT, Pires JC, Iniguez-Luy F, Leon E, Osborn TC: Genomic changes in resynthesized Brassica napus and their effect on gene expression and phenotype. Plant Cell Online. 2007, 19 (11): 3403-3417. 10.1105/tpc.107.054346.View ArticleGoogle Scholar
- Comai L, Tyagi AP, Winter K, Holmes-Davis R, Reynolds SH, Stevens Y, Byers B: Phenotypic instability and rapid gene silencing in newly formed Arabidopsis allotetraploids. Plant Cell Online. 2000, 12 (9): 1551-1568.View ArticleGoogle Scholar
- Fiebig A, Kimport R, Preuss D: Comparisons of pollen coat genes across Brassicaceae species reveal rapid evolution by repeat expansion and diversification. Proc Natl Acad Sci. 2004, 101 (9): 3286-3291. 10.1073/pnas.0305448101.PubMed CentralPubMedView ArticleGoogle Scholar
- Grover CE, Kim HR, Wing RA, Paterson AH, Wendel JF: Microcolinearity and genome evolution in the AdhA region of diploid and polyploid cotton (Gossypium). Plant J. 2007, 50 (6): 995-1006. 10.1111/j.1365-313X.2007.03102.x.PubMedView ArticleGoogle 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. 2003, 100 (8): 4649-4654. 10.1073/pnas.0630618100.PubMed CentralPubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Pumphrey M, Bai J, Laudencia-Chingcuanco D, Anderson O, Gill BS: Nonadditive expression of homoeologous genes is established upon polyploidization in hexaploid wheat. Genetics. 2009, 181 (3): 1147-1157. 10.1534/genetics.108.096941.PubMed CentralPubMedView ArticleGoogle 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 CentralPubMedView ArticleGoogle Scholar
- Kashkush K, Feldman M, Levy AA: Gene loss, silencing and activation in a newly synthesized wheat allotetraploid. Genetics. 2002, 160 (4): 1651-1659.PubMed CentralPubMedGoogle Scholar
- Li C, Zhou A, Sang T: Rice domestication by reducing shattering. Science. 2006, 311 (5769): 1936-1939. 10.1126/science.1123604.PubMedView ArticleGoogle Scholar
- Henderson IR, Jacobsen SE: Epigenetic inheritance in plants. Nature. 2007, 447 (7143): 418-424. 10.1038/nature05917.PubMedView ArticleGoogle Scholar
- Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33: 245-254. 10.1038/ng1089.PubMedView ArticleGoogle Scholar
- Chan SW-L, Henderson IR, Jacobsen SE: Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat Rev Genet. 2005, 6 (5): 351-360. 10.1038/nrg1601.PubMedView ArticleGoogle Scholar
- Cokus SJ, Feng S, Zhang X, Chen Z, Merriman B, Haudenschild CD, Pradhan S, Nelson SF, Pellegrini M, Jacobsen SE: Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature. 2008, 452 (7184): 215-219. 10.1038/nature06745.PubMed CentralPubMedView ArticleGoogle Scholar
- Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010, 328 (5980): 916-919. 10.1126/science.1186366.PubMedView ArticleGoogle Scholar
- Liu B, Vega J, Feldman M: Rapid genomic changes in newly synthesized amphiploids of Triticum and Aegilops. II. Changes in low-copy coding DNA sequences. Genome. 1998, 41 (4): 535-542.PubMedView ArticleGoogle Scholar
- Liu B, Vega J, Segal G, Abbo S, Rodova M, Feldman M: Rapid genomic changes in newly synthesized amphiploids of Triticum and Aegilops. I. Changes in low-copy noncoding DNA sequences. Genome. 1998, 41 (2): 272-277.View ArticleGoogle Scholar
- Flagel L, Udall J, Nettleton D, Wendel J: Duplicate gene expression in allopolyploid Gossypium reveals two temporally distinct phases of expression evolution. BMC Biol. 2008, 6 (1): 16-10.1186/1741-7007-6-16.PubMed CentralPubMedView ArticleGoogle Scholar
- Hovav R, Udall JA, Chaudhary B, Rapp R, Flagel L, Wendel JF: Partitioned expression of duplicated genes during development and evolution of a single cell in a polyploid plant. Proc Natl Acad Sci. 2008, 105 (16): 6191-6195. 10.1073/pnas.0711569105.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim E-D, Chen ZJ: Unstable transcripts in arabidopsis allotetraploids are associated with nonadditive gene expression in response to abiotic and biotic stresses. PLoS One. 2011, 6 (8): e24251-10.1371/journal.pone.0024251.PubMed CentralPubMedView ArticleGoogle Scholar
- Pei B, Sisu C, Frankish A, Howald C, Habegger L, Mu XJ, Harte R, Balasubramanian S, Tanzer A, Diekhans M: The GENCODE pseudogene resource. Genome Biol. 2012, 13 (9): R51-10.1186/gb-2012-13-9-r51.PubMed CentralPubMedView ArticleGoogle Scholar
- Kalyana-Sundaram S, Kumar-Sinha C, Shankar S, Robinson DR, Wu Y-M, Cao X, Asangani IA, Kothari V, Prensner JR, Lonigro RJ: Expressed pseudogenes in the transcriptional landscape of human cancers. Cell. 2012, 149 (7): 1622-1634. 10.1016/j.cell.2012.04.041.PubMed CentralPubMedView ArticleGoogle Scholar
- Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8 (4): 272-285.PubMedView ArticleGoogle Scholar
- Lisch D: How important are transposons for plant evolution?. Nat Rev Genet. 2013, 14 (1): 49-61.PubMedView ArticleGoogle Scholar
- Kinoshita Y, Saze H, Kinoshita T, Miura A, Soppe WJ, Koornneef M, Kakutani T: Control of FWA gene silencing in Arabidopsis thaliana by SINE-related direct repeats. Plant J. 2006, 49 (1): 38-45. 10.1111/j.1365-313X.2006.02936.x.PubMedView ArticleGoogle Scholar
- Hollister JD, Gaut BS: Epigenetic silencing of transposable elements: a trade-off between reduced transposition and deleterious effects on neighboring gene expression. Genome Res. 2009, 19 (8): 1419-1428. 10.1101/gr.091678.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhai J, Liu J, Liu B, Li P, Meyers BC, Chen X, Cao X: Small RNA-directed epigenetic natural variation in Arabidopsis thaliana. PLoS Genet. 2008, 4 (4): e1000056-10.1371/journal.pgen.1000056.PubMed CentralPubMedView ArticleGoogle Scholar
- Hollister JD, Smith LM, Guo Y-L, Ott F, Weigel D, Gaut BS: Transposable elements and small RNAs contribute to gene expression divergence between Arabidopsis thaliana and Arabidopsis lyrata. Proc Natl Acad Sci. 2011, 108 (6): 2322-2327. 10.1073/pnas.1018222108.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen J, Huang Q, Gao D, Wang J, Lang Y, Liu T, Li B, Bai Z, Luis Goicoechea J, Liang C, et al: Whole-genome sequencing of Oryza brachyantha reveals mechanisms underlying Oryza genome evolution. Nat Commun. 2013, 4: 1595-PubMed CentralPubMedView ArticleGoogle Scholar
- Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis comparison tool. Bioinformatics. 2005, 21 (16): 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar
- McCarthy EM, McDonald JF: LTR_STRUC: a novel search and identification program for LTR retrotransposons. Bioinformatics. 2003, 19 (3): 362-367. 10.1093/bioinformatics/btf878.PubMedView ArticleGoogle Scholar
- Xu Z, Wang H: LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007, 35 (suppl 2): W265-W268.PubMed CentralPubMedView ArticleGoogle Scholar
- Ellinghaus D, Kurtz S, Willhoeft U: LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons. BMC Bioinforma. 2008, 9 (1): 18-10.1186/1471-2105-9-18.View ArticleGoogle Scholar
- Sonnhammer ELL, Durbin R: A dot-matrix program with dynamic threshold control suited for genomic DNA and protein sequence analysis. Gene. 1995, 167 (1): GC1-GC10.PubMedGoogle Scholar
- Ma J, Bennetzen JL: Rapid recent growth and divergence of rice nuclear genomes. Proc Natl Acad Sci U S A. 2004, 101 (34): 12404-12410. 10.1073/pnas.0403715101.PubMed CentralPubMedView 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
- Larkin M, Blackshields G, Brown N, Chenna R, McGettigan P, McWilliam H, Valentin F, Wallace I, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Gaut BS, Morton BR, McCaig BC, Clegg MT: Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL. Proc Natl Acad Sci. 1996, 93 (19): 10274-10279. 10.1073/pnas.93.19.10274.PubMed CentralPubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.