Large-scale survey of cytosine methylation of retrotransposons and the impact of readout transcription from long terminal repeats on expression of adjacent rice genes

DNA Methylation and Alternative Splicing Safeguard Genome and Transcriptome After a Retrotransposition Burst in Arabidopsis thaliana

Transposable elements (TEs) are major drivers of plant genome plasticity, but the immediate molecular consequences of new TE insertions remain poorly understood. In this study, we generated a wild-type Arabidopsis thaliana population with novel insertions of ONSEN retrotransposon to investigate early epigenomic and transcriptomic changes using whole-genome and cDNA nanopore sequencing. We found that novel ONSEN insertions were distributed non-randomly, with a strong preference for genic regions, particularly in chromatin enriched for H2A.Z, H3K27me3, and H3K4me2. Most full-length ONSEN insertions within genes were rapidly recognized and spliced out as new introns (intronization), thereby mitigating potential deleterious effects on transcript isoforms. In some cases, ONSEN insertions provided alternative transcription start or termination sites, generating novel transcript isoforms. Genome-wide methylation analysis revealed that new ONSEN copies were efficiently and precisely targeted by DNA methylation. Independently on the location of the original ONSEN element, the euchromatic and heterochromatic insertions display distinct methylation signatures, reflecting the action of different epigenetic pathways. In conclusion, our results demonstrate that DNA methylation and alternative splicing are effective control mechanisms safeguarding the plant genome and transcriptome integrity after retrotransposition burst.

Exploring the Effect of High-Energy Heavy Ion Beam on Rice Genome: Transposon Activation

High-energy heavy ion beams are a new type of physical mutagen that can produce a wide range of phenotypic variations. In order to understand the mechanism of high-energy heavy ion beams, we resequenced the whole genome of individual plants with obvious phenotypic variations in rice. The sequence alignment results revealed a large number of SNPs and InDels, as well as genetic variations related to grain type and heading date. The distribution of SNP and InDel on chromosomes is random, but they often occur in the up/downstream regions and the intergenic region. Mutagenesis can cause changes in transposons such as Dasheng, mPing, Osr13 and RIRE2, affecting the stability of the genome. This study obtained the major gene mutation types, discovered differentially active transposons, screened out gene variants related to phenotype, and explored the mechanism of high-energy heavy ion beam radiation on rice genes.

Epigenetic variation: A major player in facilitating plant fitness under changing environmental conditions

Recent research in plant epigenetics has increased our understanding of how epigenetic variability can contribute to adaptive phenotypic plasticity in natural populations. Studies show that environmental changes induce epigenetic switches either independently or in complementation with the genetic variation. Although most of the induced epigenetic variability gets reset between generations and is short-lived, some variation becomes transgenerational and results in heritable phenotypic traits. The short-term epigenetic responses provide the first tier of transient plasticity required for local adaptations while transgenerational epigenetic changes contribute to stress memory and help the plants respond better to recurring or long-term stresses. These transgenerational epigenetic variations translate into an additional tier of diversity which results in stable epialleles. In recent years, studies have been conducted on epigenetic variation in natural populations related to various biological processes, ecological factors, communities, and habitats. With the advent of advanced NGS-based technologies, epigenetic studies targeting plants in diverse environments have increased manifold to enhance our understanding of epigenetic responses to environmental stimuli in facilitating plant fitness. Taking all points together in a frame, the present review is a compilation of present-day knowledge and understanding of the role of epigenetics and its fitness benefits in diverse ecological systems in natural populations.

DNA methylation and histone modifications induced by abiotic stressors in plants

BACKGROUND

A review of research shows that methylation in plants is more complex and sophisticated than in microorganisms and animals. Overall, studies on the effects of abiotic stress on epigenetic modifications in plants are still scarce and limited to few species. Epigenetic regulation of plant responses to environmental stresses has not been elucidated. This study summarizes key effects of abiotic stressors on DNA methylation and histone modifications in plants.

DISCUSSION

Plant DNA methylation and histone modifications in responses to abiotic stressors varied and depended on the type and level of stress, plant tissues, age, and species. A critical analysis of the literature available revealed that 44% of the epigenetic modifications induced by abiotic stressors in plants involved DNA hypomethylation, 40% DNA hypermethylation, and 16% histone modification. The epigenetic changes in plants might be underestimated since most authors used methods such as methylation-sensitive amplification polymorphism (MSAP), High performance liquid chromatography (HPLC), and immunolabeling that are less sensitive compared to bisulfite sequencing and single-base resolution methylome analyses. More over, mechanisms underlying epigenetic changes in plants have not yet been determined since most reports showed only the level or/and distribution of DNA methylation and histone modifications.

CONCLUSIONS

Various epigenetic mechanisms are involved in response to abiotic stressors, and several of them are still unknown. Integrated analysis of the changes in the genome by omic approaches should help to identify novel components underlying mechanisms involved in DNA methylation and histone modifications associated with plant response to environmental stressors.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE-based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty-five transposon-inserted genes were sequenced and homology-based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE-associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Identification and characterization of differentially expressed genes in the rice root following exogenous application of spermidine during salt stress

Salinity is a major limiting factor in crop production. Exogenous spermidine (spd) effectively ameliorates salt injury, though the underlying molecular mechanism is poorly understood. We have used a suppression subtractive hybridization method to construct a cDNA library that has identified up-regulated genes from rice root under the treatment of spd and salt. Total 175 high-quality ESTs of about 100-500 bp in length with an average size of 200 bp are isolated, clustered and assembled into a collection of 62 unigenes. Gene ontology analysis using the KEGG pathway annotation database has classified the unigenes into 5 main functional categories and 13 subcategories. The transcripts abundance has been validated using Real-Time PCR. We have observed seven different types of post-translational modifications in the DEPs. 44 transmembrane helixes are predicted in 6 DEPs. This above information can be used as first-hand data for dissecting the administrative role of spd during salinity.

Draft genome sequence of Solanum aethiopicum provides insights into disease resistance, drought tolerance, and the evolution of the genome

BACKGROUND

The African eggplant (Solanum aethiopicum) is a nutritious traditional vegetable used in many African countries, including Uganda and Nigeria. It is thought to have been domesticated in Africa from its wild relative, Solanum anguivi. S. aethiopicum has been routinely used as a source of disease resistance genes for several Solanaceae crops, including Solanum melongena. A lack of genomic resources has meant that breeding of S. aethiopicum has lagged behind other vegetable crops.

RESULTS

We assembled a 1.02-Gb draft genome of S. aethiopicum, which contained predominantly repetitive sequences (78.9%). We annotated 37,681 gene models, including 34,906 protein-coding genes. Expansion of disease resistance genes was observed via 2 rounds of amplification of long terminal repeat retrotransposons, which may have occurred ∼1.25 and 3.5 million years ago, respectively. By resequencing 65 S. aethiopicum and S. anguivi genotypes, 18,614,838 single-nucleotide polymorphisms were identified, of which 34,171 were located within disease resistance genes. Analysis of domestication and demographic history revealed active selection for genes involved in drought tolerance in both "Gilo" and "Shum" groups. A pan-genome of S. aethiopicum was assembled, containing 51,351 protein-coding genes; 7,069 of these genes were missing from the reference genome.

CONCLUSIONS

The genome sequence of S. aethiopicum enhances our understanding of its biotic and abiotic resistance. The single-nucleotide polymorphisms identified are immediately available for use by breeders. The information provided here will accelerate selection and breeding of the African eggplant, as well as other crops within the Solanaceae family.

Identification of transposons near predicted lncRNA and mRNA pools of Prunus mume using an integrative transposable element database constructed from Rosaceae plant genomes

This study focused on the construction of a database of transposable elements (TEs) from Rosaceae plants, the third most economically important plant family in temperate regions, and its transcriptomics applications. The evolutionary effects of TEs on gene regulation have been explored, and TE insertions can be the molecular bases of changes in gene structure and function. However, a specific Rosaceae plant TE database (RPTEdb) is lacking. The genomes of several Rosaceae plants have been sequenced, providing the opportunity to mine TE data at a whole-genome level. Therefore, we constructed the RPTEdb, a collective and comprehensive database of 19,596 annotated TEs in the genomes of Rosaceae plants using previously described identification and annotation methods and published genome sequences. The user-friendly web-based database provides access to research tools through hyperlinks, including Browse, TE tree, tools, JBrowse, and search sections, and through the inputting of sequences on the main webpage. Next, we performed one advanced application in which TEs near predicted long non-coding RNA (lncRNA) and mRNA domains within white and red petal-tissue transcriptomes of Prunus mume 'Fuban Tiaozhi' were identified, revealing 16 TEs that overlapped or were near 16 differentially expressed lncRNA domains, and 54 TEs that overlapped or were near 54 differentially expressed mRNA domains, and the TEs' possible functions were also discussed. We believe that the RPTEdb will contribute to the understanding of TE roles in the structural, functional and evolutionary dynamics of Rosaceae plant genomes.

Analysis of somaclonal variation in transgenic and regenerated plants of Arabidopsis thaliana using methylation related metAFLP and TMD markers

We provide evidence that nucleotide sequence and methylation status changes occur in the Arabidopsis genome during in vitro tissue culture at a frequency high enough to represent an important source of variation. Somaclonal variation is a general consequence of the tissue culture process that has to be analyzed specifically when regenerated plants are obtained in any plant species. Currently, there are few studies about the variability comprising sequence changes and methylation status at the DNA level, generated by the culture of A. thaliana cells and tissues. In this work, two types of highly reproducible molecular markers, modified methylation sensitive AFLP (metAFLP) and transposon methylation display (TMD) have been used for the first time in this species to analyze the nucleotide and cytosine methylation changes induced by transformation and tissue culture protocols. We found significantly higher average methylation values (7.5%) in regenerated and transgenic plants when compared to values obtained from seed derived plants (3.2%) and that the main component of the somaclonal variation present in Arabidopsis clonal plants is genetic rather than epigenetic. However, we have found that the Arabidopsis regenerated and transgenic plants had a higher number of non-fully methylated sites flanking transposable elements than the control plants, and therefore, their mobilization can be facilitated. These data provide further evidence that changes in nucleotide sequence and methylation status occur in the Arabidopsis genome during in vitro tissue culture frequently enough to be an important source of variation in this species.

Insertion of a solo LTR retrotransposon associates with spur mutations in 'Red Delicious' apple (Malus × domestica)

Insertion of a solo LTR, which possesses strong bidirectional, stem-specific promoter activities, is associated with the evolution of a dwarfing apple spur mutation. Spur mutations in apple scions revolutionized global apple production. Since long terminal repeat (LTR) retrotransposons are tightly related to natural mutations, inter-retrotransposon-amplified polymorphism technique and genome walking were used to find sequences in the apple genome based on these LTRs. In 'Red Delicious' spur mutants, a novel, 2190-bp insertion was identified as a spur-specific, solo LTR (sLTR) located at the 1038th nucleotide of another sLTR, which was 1536 bp in length. This insertion-within-an-insertion was localized within a preexisting Gypsy-50 retrotransposon at position 3,762,767 on chromosome 4. The analysis of transcriptional activity of the two sLTRs (the 2190- and 1536-bp inserts) indicated that the 2190-bp sLTR is a promoter, capable of bidirectional transcription. GUS expression in the 2190-bp-sense and 2190-bp-antisense transgenic lines was prominent in stems. In contrast, no promoter activity from either the sense or the antisense strand of the 1536-bp sLTR was detected. From ~150 kb of DNA on each side of the 2190 bp, sLTR insertion site, corresponding to 300 kb of the 'Golden Delicious' genome, 23 genes were predicted. Ten genes had predicted functions that could affect shoot development. This first report, of a sLTR insertion associated with the evolution of apple spur mutation, will facilitate apple breeding, cloning of spur-related genes, and discovery of mechanisms behind dwarf habit.

LTR-retrotransposons in plants: Engines of evolution

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.

Profiling Transposable Elements and Their Epigenetic Effects in Non-model Species

Taking transposable elements into consideration in surveys of genetic and epigenetic variation remains challenging in species lacking a high-quality reference genome. Here, molecular techniques reducing genome complexity and specifically targeting restructuring and methylation changes in TE genome fractions are described. In particular, methyl-sensitive transposon display (MSTD) uses isoschizomers and PCR amplifications to assess the methylation environment of TE insertions. MSTD offers reliable insights into genome-wide epigenetic changes associated with TEs, especially when used together with similar techniques tracking random sequences.

Laser Irradiation-Induced DNA Methylation Changes Are Heritable and Accompanied with Transpositional Activation of mPing in Rice

DNA methylation is an integral component of the epigenetic code in most higher eukaryotes. Exploring the extent to which DNA methylation can be altered under a specific condition and its heritability is important for elucidating the biological functions of this epigenetic modification. Here, we conducted MSAP analysis of rice plants with altered phenotypes subsequent to a low-dose Nd³⁺YAG laser irradiation. We found that all four methylation patterns at the 5'-CCGG sites that are analyzable by MSAP showed substantial changes in the immediately treated M0 plants. Interestingly, the frequencies of hypo- and hypermethylation were of similar extents, which largely offset each other and render the total methylation levels unchanged. Further analysis revealed that the altered methylation patterns were meiotically heritable to at least the M2 generation but accompanied with further changes in each generation. The methylation changes and their heritability of the metastable epigenetic state were verified by bisulfite sequencing of portion of the retrotranspon, Tos17, an established locus for assessing DNA methylation liability in rice. Real-time PCR assay indicated that the expression of various methylation-related chromatin genes was perturbed, and a Pearson correlation analysis showed that many of these genes, especially two AGOs (AGO4-1 and AGO4-2), were significantly correlated with the methylation pattern alterations. In addition, excisions of a MITE transposon, mPing, occurred rampantly in the laser irradiated plants and their progenies. Together, our results indicate that heritable DNA methylation changes can be readily induced by low-dose laser irradiation, and which can be accompanied by transpostional activation of transposable elements.

Ty1-copia elements reveal diverse insertion sites linked to polymorphisms among flax (Linum usitatissimum L.) accessions

BACKGROUND

Initial characterization of the flax genome showed that Ty1-copia retrotransposons are abundant, with several members being recently inserted, and in close association with genes. Recent insertions indicate a potential for ongoing transpositional activity that can create genomic diversity among accessions, cultivars or varieties. The polymorphisms generated constitute a good source of molecular markers that may be associated with phenotype if the insertions alter gene activity. Flax, where accessions are bred mainly for seed nutritional properties or for fibers, constitutes a good model for studying the relationship of transpositional activity with diversification and breeding. In this study, we estimated copy number and used a type of transposon display known as Sequence-Specific Amplification Polymorphisms (SSAPs), to characterize six families of Ty1-copia elements across 14 flax accessions. Polymorphic insertion sites were sequenced to find insertions that could potentially alter gene expression, and a preliminary test was performed with selected genes bearing transposable element (TE) insertions.

RESULTS

Quantification of six families of Ty1-copia elements indicated different abundances among TE families and between flax accessions, which suggested diverse transpositional histories. SSAPs showed a high level of polymorphism in most of the evaluated retrotransposon families, with a trend towards higher levels of polymorphism in low-copy number families. Ty1-copia insertion polymorphisms among cultivars allowed a general distinction between oil and fiber types, and between spring and winter types, demonstrating their utility in diversity studies. Characterization of polymorphic insertions revealed an overwhelming association with genes, with insertions disrupting exons, introns or within 1 kb of coding regions. A preliminary test on the potential transcriptional disruption by TEs of four selected genes evaluated in three different tissues, showed one case of significant impact of the insertion on gene expression.

CONCLUSIONS

We demonstrated that specific Ty1-copia families have been active since breeding commenced in flax. The retrotransposon-derived polymorphism can be used to separate flax types, and the close association of many insertions with genes defines a good source of potential mutations that could be associated with phenotypic changes, resulting in diversification processes.

Methods for accurate quantification of LTR-retrotransposon copy number using short-read sequence data: a case study in Sorghum

Transposable elements (TEs) are ubiquitous in eukaryotic genomes and their mobility impacts genome structure and function in myriad ways. Because of their abundance, activity, and repetitive nature, the characterization and analysis of TEs remain challenging, particularly from short-read sequencing projects. To overcome this difficulty, we have developed a method that estimates TE copy number from short-read sequences. To test the accuracy of our method, we first performed an in silico analysis of the reference Sorghum bicolor genome, using both reference-based and de novo approaches. The resulting TE copy number estimates were strikingly similar to the annotated numbers. We then tested our method on real short-read data by estimating TE copy numbers in several accessions of S. bicolor and its close relative S. propinquum. Both methods effectively identify and rank similar TE families from highest to lowest abundance. We found that de novo characterization was effective at capturing qualitative variation, but underestimated the abundance of some TE families, specifically families of more ancient origin. Also, interspecific reference-based mapping of S. propinquum reads to the S. bicolor database failed to fully describe TE content in S. propinquum, indicative of recent TE activity leading to changes in the respective repetitive landscapes over very short evolutionary timescales. We conclude that reference-based analyses are best suited for within-species comparisons, while de novo approaches are more reliable for evolutionarily distant comparisons.

The Role of DNA Methylation in Xylogenesis in Different Tissues of Poplar

In trees, xylem tissues play a key role in the formation of woody tissues, which have important uses for pulp and timber production; also DNA methylation plays an important part in gene regulation during xylogenesis in trees. In our study, methylation-sensitive amplified polymorphism (MSAP) analysis was used to analyze the role cytosine methylation plays in wood formation in the commercially important tree species Populus tomentosa. This analysis compared the methylation patterns between xylem tissues (developing xylem and mature xylem) and non-xylem tissues (cambium, shoot apex, young leaf, mature leaf, phloem, root, male catkin, and female catkin) and found 10,316 polymorphic methylation sites. MSAP identified 132 candidate genes with the same methylation patterns in xylem tissues, including seven wood-related genes. The expression of these genes differed significantly between xylem and non-xylem tissue types (P < 0.01). This indicated that the difference of expression of specific genes with unique methylation patterns, rather than relative methylation levels between the two tissue types plays a critical role in wood biosynthesis. However, 46.2% of candidate genes with the same methylation pattern in vascular tissues (cambium, phloem, and developing xylem) did not have distinct expression patterns in xylem and non-xylem tissue. Also, bisulfite sequencing and transcriptome sequencing of MYB, NAC and FASCICLIN-LIKE AGP 13 revealed that the location of cytosine methylation in the gene might affect the expression of different transcripts from the corresponding gene. The expression of different transcripts that produce distinct proteins from a single gene might play an important role in the regulation of xylogenesis.

Structure and extent of DNA methylation-based epigenetic variation in wild emmer wheat (T. turgidum ssp. dicoccoides) populations

BACKGROUND

The genetic structure and differentiation of wild emmer wheat suggests that genetic diversity is eco-geographically structured. However, very little is known about the structure and extent of the heritable epigenetic variation and its influence on local adaptation in natural populations.

RESULTS

The structure and extent of the heritable methylation-based epigenetic variation were assessed within and among natural populations of Triticum turgidum ssp. dicoccoides. We used methylation sensitive amplified polymorphism (MSAP) and transposon methylation display (TMD) techniques, to assess the methylation status of random genomic CCGG sites and CCGG sites flanking transposable elements (TEs), respectively. Both techniques were applied to the DNA of 50 emmer accessions which were collected from five different geographically isolated regions. In order to ensure the assessment of heritable epigenetic variation, all accessions were grown under common garden conditions for two generations. In all accessions, the difference in methylation levels of CCGG sites, including CCGG sites that flanked TEs, were not statistically significant and relatively high, ranging between 46 and 76 %. The pattern of methylation was significantly different among accessions, such that clear and statistically significant population-specific methylation patterns were observed.

CONCLUSION

In this study, we have observed population-unique heritable methylation patterns in emmer wheat accessions originating from five geographically isolated regions. Our data indicate that methylation-based epigenetic diversity might be eco-geographically structured and might be partly determined by climatic and edaphic factors.

Short-term hybridisation activates Tnt1 and Tto1 Copia retrotransposons in wild tuber-bearing Solanum species

Interspecific hybridisation in tuber-bearing species of Solanum is a common phenomenon and represents an important source of variability, crucial for adaptation and speciation of potato species. In this regard, the effects of interspecific hybridisation on retrotransposon families present in the genomes, and their consequent effects on generation of genetic variability in wild tuber-bearing Solanum species, are poorly characterised. The aim of this study was to analyse the activity of retrotransposons in inter- and intraspecific hybrids between S. kurtzianum and S. microdontum, obtained by controlled crosses, and the effects on morphological, genetic and epigenetic variability. For genetic and epigenetic analysis, S-SAP (sequence-specific amplification polymorphism) and TMD (transposon methylation display) techniques were used, respectively, with specific primers for Tnt1 and Tto1 retrotransposon families (Order LTR, Superfamily Copia). The results indicate that at morphological level, interspecific hybrid genotypes differ from their parental species, whereas derived intraspecific hybrids do not. In both cases, we observed significant reductions in pollen grain viability, and a negative correlation with Tnt1 mobility. Both retrotransposons, Tto1 and Tnt1, were mobilised in the genotypes analysed, with mobility ranging from 0 to 7.8%. Furthermore, at the epigenetic level, demethylation was detected in the vicinity of Tnt1 and Tto1 in the hybrids compared with the parental genotypes. These patterns were positively correlated with the activity of the retrotransposons. The results suggest a possible mechanism through which hybridisation events generate genetic variability in tuber-bearing species of Solanum through retrotranposon activation.

Differential regulation of genes by retrotransposons in rice promoters

Rice genome harbors genes and promoters with retrotransposon insertions. There is very little information about their function. The effect of retrotransposon insertions in four rice promoter regions on gene regulation, was investigated using promoter-reporter gene constructs with and without retrotransposons. Differences in expression levels of gus and egfp reporter genes in forward orientation and rfp in reverse orientation were evaluated in rice plants with transient expression employing quantitative RT-PCR analysis, histochemical GUS staining, and eGFP and RFP fluorescent microscopy. The presence of SINE in the promoter 1 (P1) resulted in higher expression levels of the reporter genes, whereas the presence of LINE in P2 or gypsy LTR retrotransposon in P3 reduced expression of the reporter genes. Furthermore, the SINE in P1 acts as an enhancer in contrast with the LINE in P2 and the gypsy LTR retrotransposon in P3 which act as silencers. CTAA and CGG motifs in these retrotransposons are the likely candidates for the downregulation compared to TCTT motif (SINE) which is a candidate for the upregulation of gene expression. The effect of retrotransposons on gene regulation correlated with the earlier investigation of conservation patterns of these four retrotransposon insertions in several rice accessions implying their evolutionary significance.

Genome-wide analysis of DNA methylation in photoperiod- and thermo-sensitive male sterile rice Peiai 64S

Background

Epigenetic modifications play important roles in the regulation of plant development. DNA methylation is an important epigenetic modification that dynamically regulates gene expression during developmental processes. However, little studies have been reported about the methylation profiles of photoperiod- and thermo-sensitive genic male sterile (PTGMS) rice during the fertility transition.

Results

In this study, using methylated DNA immunoprecipitation sequencing (MeDIP-seq), the global DNA methylation patterns were compared in the rice PTGMS line PA64S under two different environments (different temperatures and day lengths). The profiling of the DNA methylation under two different phenotypes (sterility and fertility) revealed that hypermethylation was observed in PA64S (sterility), and 1258 differentially methylated regions (DMRs) were found between PA64S (sterility) and PA64S (fertility). Twenty differentially methylated genes of them were further validated through bisulfite sequencing, and four of these genes were analyzed by qRT-PCR. Especially, a differentially methylated gene (LOC_Os08g38210), which encoded transcription factor BIM2, is a component of brassinosteroid signaling in rice. The hypermethylated BIM2 gene may suppress some downstream genes in brassinosteroid signaling pathway, and thus affect the male fertility in PA64S.

Conclusions

The results presented here indicated that hypermethylation was observed in PA64S (sterility). Gene Ontology (GO) analysis and KEGG analysis revealed that flavone and flavonol biosynthrsis, circadian rhythm, photosynthesis and oxidative phosphorylation pathways were involved in sterility-fertility transition of PA64S.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1317-7) contains supplementary material, which is available to authorized users.

Nascent transcription affected by RNA polymerase IV in Zea mays

All eukaryotes use three DNA-dependent RNA polymerases (RNAPs) to create cellular RNAs from DNA templates. Plants have additional RNAPs related to Pol II, but their evolutionary role(s) remain largely unknown. Zea mays (maize) RNA polymerase D1 (RPD1), the largest subunit of RNA polymerase IV (Pol IV), is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs), and transcriptional regulation of specific alleles. Here, we define the nascent transcriptomes of rpd1 mutant and wild-type (WT) seedlings using global run-on sequencing (GRO-seq) to identify the broader targets of RPD1-based regulation. Comparisons of WT and rpd1 mutant GRO-seq profiles indicate that Pol IV globally affects transcription at both transcriptional start sites and immediately downstream of polyadenylation addition sites. We found no evidence of divergent transcription from gene promoters as seen in mammalian GRO-seq profiles. Statistical comparisons identify genes and TEs whose transcription is affected by RPD1. Most examples of significant increases in genic antisense transcription appear to be initiated by 3'-proximal long terminal repeat retrotransposons. These results indicate that maize Pol IV specifies Pol II-based transcriptional regulation for specific regions of the maize genome including genes having developmental significance.

Different tobacco retrotransposons are specifically modulated by the elicitor cryptogein and reactive oxygen species

Interactions of plant retrotransposons with different steps of biotic and abiotic stress-associated signaling cascades are still poorly understood. We perform here a finely tuned comparison of four tobacco retrotransposons (Tnt1, Tnt2, Queenti, and Tto1) responses to the plant elicitor cryptogein. We demonstrate that basal transcript levels in cell suspensions and plant leaves as well as the activation during the steps of defense signaling events are specific to each retrotransposon. Using antisense NtrbohD lines, we show that NtrbohD-dependent reactive oxygen species (ROS) production might act as negative regulator of retrotransposon activation.

Roles, and establishment, maintenance and erasing of the epigenetic cytosine methylation marks in plants

Heritable information in plants consists of genomic information in DNA sequence and epigenetic information superimposed on DNA sequence. The latter is in the form of cytosine methylation at CG, CHG and CHH elements (where H = A, T orC) and a variety of histone modifications in nucleosomes. The epialleles arising from cytosine methylation marks on the nuclear genomic loci have better heritability than the epiallelic variation due to chromatin marks. Phenotypic variation is increased manifold by epiallele comprised methylomes. Plants (angiosperms) have highly conserved genetic mechanisms to establish, maintain or erase cytosine methylation from epialleles. The methylation marks in plants fluctuate according to the cell/tissue/organ in the vegetative and reproductive phases of plant life cycle. They also change according to environment. Epialleles arise by gain or loss of cytosine methylation marks on genes. The changes occur due to the imperfection of the processes that establish and maintain the marks and on account of spontaneous and stress imposed removal of marks. Cytosine methylation pattern acquired in response to abiotic or biotic stress is often inherited over one to several subsequent generations.Cytosine methylation marks affect physiological functions of plants via their effect(s) on gene expression levels. They also repress transposable elements that are abundantly present in plant genomes. The density of their distribution along chromosome lengths affects meiotic recombination rate, while their removal increases mutation rate. Transposon activation due to loss of methylation causes rearrangements such that new gene regulatory networks arise and genes for microRNAs may originate. Cytosine methylation dynamics contribute to evolutionary changes. This review presents and discusses the available evidence on origin, removal and roles of cytosine methylation and on related processes, such as RNA directed DNA methylation, imprinting, paramutation and transgenerational memory in plants.

Tissue culture-induced heritable genomic variation in rice, and their phenotypic implications

BACKGROUND

Somaclonal variation generally occurs in plants regenerated from tissue culture. However, fundamental issues regarding molecular characteristics, mutation rates and mutation spectra of plant somatic variation as well as their phenotypic relevance have been addressed only recently. Moreover, these studies have reported highly discrepant results in different plant species and even in the same plant genotype.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated heritable genomic variation induced by tissue culture in rice by whole genome re-sequencing of an extensively selfed somaclonal line (TC-reg-2008) and its wild type (WT) donor (cv. Hitomebore). We computed the overall mutation rate, single nucleotide polymorphisms (SNPs), small scale insertions/deletions (Indels) and mobilization of transposable elements (TEs). We assessed chromosomal distribution of the various types of genomic variations, tested correlations between SNPs and Indels, and examined concomitancy between TE activity and its cytosine methylation states. We also performed gene ontology (GO) analysis of genes containing nonsynonymous mutations and large-effect mutations, and assayed effects of the genomic variations on phenotypes under both normal growing condition and several abiotic stresses. We found that heritable somaclonal genomic variation occurred extensively in rice. The genomic variations distributed non-randomly across each of the 12 rice chromosomes, and affected a large number of functional genes. The phenotypic penetrance of the genomic variations was condition-dependent.

CONCLUSIONS/SIGNIFICANCE

Tissue culture is a potent means to generate heritable genetic variations in rice, which bear distinct difference at least in space (chromosomal distribution) from those occurred under natural settings. Our findings have provided new information regarding the mutation rate and spectrum as well as chromosomal distribution pattern of somaclonal variation in rice. Our data also suggest that rice possesses a strong capacity to canalize genetic variations under normal growing conditions to maintain phenotypic robustness, which however can be released by certain abiotic stresses to generate variable phenotypes.

Transposable elements: powerful contributors to angiosperm evolution and diversity

Transposable elements (TEs) are a dominant feature of most flowering plant genomes. Together with other accepted facilitators of evolution, accumulating data indicate that TEs can explain much about their rapid evolution and diversification. Genome size in angiosperms is highly correlated with TE content and the overwhelming bulk (>80%) of large genomes can be composed of TEs. Among retro-TEs, long terminal repeats (LTRs) are abundant, whereas DNA-TEs, which are often less abundant than retro-TEs, are more active. Much adaptive or evolutionary potential in angiosperms is due to the activity of TEs (active TE-Thrust), resulting in an extraordinary array of genetic changes, including gene modifications, duplications, altered expression patterns, and exaptation to create novel genes, with occasional gene disruption. TEs implicated in the earliest origins of the angiosperms include the exapted Mustang, Sleeper, and Fhy3/Far1 gene families. Passive TE-Thrust can create a high degree of adaptive or evolutionary potential by engendering ectopic recombination events resulting in deletions, duplications, and karyotypic changes. TE activity can also alter epigenetic patterning, including that governing endosperm development, thus promoting reproductive isolation. Continuing evolution of long-lived resprouter angiosperms, together with genetic variation in their multiple meristems, indicates that TEs can facilitate somatic evolution in addition to germ line evolution. Critical to their success, angiosperms have a high frequency of polyploidy and hybridization, with resultant increased TE activity and introgression, and beneficial gene duplication. Together with traditional explanations, the enhanced genomic plasticity facilitated by TE-Thrust, suggests a more complete and satisfactory explanation for Darwin's "abominable mystery": the spectacular success of the angiosperms.

The bright side of transposons in crop evolution

The past decades have revealed an unexpected yet prominent role of so-called 'junk DNA' in the regulation of gene expression, thereby challenging our view of the mechanisms underlying phenotypic evolution. In particular, several mechanisms through which transposable elements (TEs) participate in functional genome diversity have been depicted, bringing to light the 'TEs bright side'. However, the relative contribution of those mechanisms and, more generally, the importance of TE-based polymorphisms on past and present phenotypic variation in crops species remain poorly understood. Here, we review current knowledge on both issues, and discuss how analyses of massively parallel sequencing data combined with statistical methodologies and functional validations will help unravelling the impact of TEs on crop evolution in a near future.

Detecting epigenetic effects of transposable elements in plants

Transposable elements (TE) represent a major fraction of eukaryotic genomes and play many roles in plant epigenetics. In this chapter, we describe the use of Sequence-Specific Amplified Polymorphism (SSAP) as a reliable Transposon Display technique applicable for use in many plant species. We also discuss the interpretation of SSAP data and associated risks. This technique has potential to allow rapid screening of plant populations, especially in nonmodel or wild species.

Copy number variation of transposable elements in Triticum-Aegilops genus suggests evolutionary and revolutionary dynamics following allopolyploidization

Here, we report on copy number variation of transposable elements and on the genome-specific proliferation in wheat. In addition, we report on revolutionary and evolutionary dynamics of transposons. Wheat is a valuable model for understanding the involvement of transposable elements (TEs) in speciation as wheat species (Triticum-Aegilops group) have diverged from a common ancestor, have undergone two events of speciation through allopolyploidy, and contain a very high fraction of TEs. However, an unbiased genome-wide examination of TE variation among these species has not been conducted. Our research utilized quantitative real time PCR to assess the relative copy numbers of 16 TE families in various Triticum and Aegilops species. We found (1) high variation and genome-specificity of TEs in wheat species, suggesting they were active throughout the evolution of wheat, (2) neither Ae. searsii nor Ae. speltoides by themselves can be the only contributors of the B genome to wheat, and (3) nonadditive changes in TE quantities in polyploid wheat. This study indicates the apparent involvement of large TEs in creating genetic variation in revolutionary and evolutionary scales following allopolyploidization events, presumably assisting in the diploidization of homeologous chromosomes.

Genome-wide analysis of short interspersed nuclear elements SINES revealed high sequence conservation, gene association and retrotranspositional activity in wheat

Short interspersed nuclear elements (SINEs) are non-autonomous non-LTR retroelements that are present in most eukaryotic species. While SINEs have been intensively investigated in humans and other animal systems, they are poorly studied in plants, especially in wheat (Triticum aestivum). We used quantitative PCR of various wheat species to determine the copy number of a wheat SINE family, termed Au SINE, combined with computer-assisted analyses of the publicly available 454 pyrosequencing database of T. aestivum. In addition, we utilized site-specific PCR on 57 Au SINE insertions, transposon methylation display and transposon display on newly formed wheat polyploids to assess retrotranspositional activity, epigenetic status and genetic rearrangements in Au SINE, respectively. We retrieved 3706 different insertions of Au SINE from the 454 pyrosequencing database of T. aestivum, and found that most of the elements are inserted in A/T-rich regions, while approximately 38% of the insertions are associated with transcribed regions, including known wheat genes. We observed typical retrotransposition of Au SINE in the second generation of a newly formed wheat allohexaploid, and massive hypermethylation in CCGG sites surrounding Au SINE in the third generation. Finally, we observed huge differences in the copy numbers in diploid Triticum and Aegilops species, and a significant increase in the copy numbers in natural wheat polyploids, but no significant increase in the copy number of Au SINE in the first four generations for two of three newly formed allopolyploid species used in this study. Our data indicate that SINEs may play a prominent role in the genomic evolution of wheat through stress-induced activation.

Retrotransposon replication in plants

Retrotransposons comprise the bulk of large plant genomes, replicating via an RNA intermediate whereby the original, integrated element remains in place. Of the two main orders, the LTR retrotransposons considerably outnumber the LINEs. LINEs integrate into target sites simultaneously with the RNA transcript being copied into cDNA by target-primed reverse transcription. LTR retrotransposon replication is basically equivalent to the intracellular phase of retroviral life cycles. The envelope gene giving extracellular mobility to retroviruses is in fact widespread in plants and their retrotransposons. Evolutionary analyses of the retrotransposons and retroviruses suggest that both form an ancient monophyletic group. The particular adaptations of LTR retrotransposons to plant life cycles enabling their success remain to be clarified.

Arguments for standardizing transposable element annotation in plant genomes

Whole genome sequence assemblies have been generated for many plants. Annotation of transposable elements (TEs), which constitute the major proportion of genomes and play a significant role in epigenome alterations under stress, has not been given equal importance to that of genes. In this opinion article, we argue that the lack of focus dedicated to the fine-scale characterization of repeat fractions and the absence of consistent methods for their annotation impede our ability to critically understand the influence of TEs on the epigenome with implications in gene expression and non-Mendelian inheritance. Major structural changes occur over an evolutionary time scale. However, epigenetic regulation mediated by TEs can happen in a single generation, thus emphasizing the need for their standardized annotation.

Co-evolution of plant LTR-retrotransposons and their host genomes

Transposable elements (TEs), particularly, long terminal repeat retrotransposons (LTR-RTs), are the most abundant DNA components in all plant species that have been investigated, and are largely responsible for plant genome size variation. Although plant genomes have experienced periodic proliferation and/or recent burst of LTR-retrotransposons, the majority of LTR-RTs are inactivated by DNA methylation and small RNA-mediated silencing mechanisms, and/or were deleted/truncated by unequal homologous recombination and illegitimate recombination, as suppression mechanisms that counteract genome expansion caused by LTR-RT amplification. LTR-RT DNA is generally enriched in pericentromeric regions of the host genomes, which appears to be the outcomes of preferential insertions of LTR-RTs in these regions and low effectiveness of selection that purges LTR-RT DNA from these regions relative to chromosomal arms. Potential functions of various TEs in their host genomes remain blurry; nevertheless, LTR-RTs have been recognized to play important roles in maintaining chromatin structures and centromere functions and regulation of gene expressions in their host genomes.

Genomic change, retrotransposon mobilization and extensive cytosine methylation alteration in Brassica napus introgressions from two intertribal hybridizations

Hybridization and introgression represent important means for the transfer and/or de novo origination of traits and play an important role in facilitating speciation and plant breeding. Two sets of introgression lines in Brassica napus L. were previously established by its intertribal hybridizations with two wild species and long-term selection. In this study, the methods of amplified fragment length polymorphisms (AFLP), sequence-specific amplification polymorphism (SSAP) and methylation-sensitive amplified polymorphism (MSAP) were used to determine their genomic change, retrotransposon mobilization and cytosine methylation alteration in these lines. The genomic change revealed by the loss or gain of AFLP bands occurred for ∼10% of the total bands amplified in the two sets of introgressions, while no bands specific for wild species were detected. The new and absent SSAP bands appeared for 9 out of 11 retrotransposons analyzed, with low frequency of new bands and their total percentage of about 5% in both sets. MSAP analysis indicated that methylation changes were common in these lines (33.4-39.8%) and the hypermethylation was more frequent than hypomethylation. Our results suggested that certain extents of genetic and epigenetic alterations were induced by hybridization and alien DNA introgression. The cryptic mechanism of these changes and potential application of these lines in breeding were also discussed.

How important are transposons for plant evolution?

For decades, transposable elements have been known to produce a wide variety of changes in plant gene expression and function. This has led to the idea that transposable element activity has played a key part in adaptive plant evolution. This Review describes the kinds of changes that transposable elements can cause, discusses evidence that those changes have contributed to plant evolution and suggests future strategies for determining the extent to which these changes have in fact contributed to plant adaptation and evolution. Recent advances in genomics and phenomics for a range of plant species, particularly crops, have begun to allow the systematic assessment of these questions.

Genome-wide analysis of Stowaway-like MITEs in wheat reveals high sequence conservation, gene association, and genomic diversification

The diversity and evolution of wheat (Triticum-Aegilops group) genomes is determined, in part, by the activity of transposable elements that constitute a large fraction of the genome (up to 90%). In this study, we retrieved sequences from publicly available wheat databases, including a 454-pyrosequencing database, and analyzed 18,217 insertions of 18 Stowaway-like miniature inverted-repeat transposable element (MITE) families previously characterized in wheat that together account for approximately 1.3 Mb of sequence. All 18 families showed high conservation in length, sequence, and target site preference. Furthermore, approximately 55% of the elements were inserted in transcribed regions, into or near known wheat genes. Notably, we observed significant correlation between the mean length of the MITEs and their copy number. In addition, the genomic composition of nine MITE families was studied by real-time quantitative polymerase chain reaction analysis in 40 accessions of Triticum spp. and Aegilops spp., including diploids, tetraploids, and hexaploids. The quantitative polymerase chain reaction data showed massive and significant intraspecific and interspecific variation as well as genome-specific proliferation and nonadditive quantities in the polyploids. We also observed significant differences in the methylation status of the insertion sites among MITE families. Our data thus suggest a possible role for MITEs in generating genome diversification and in the establishment of nascent polyploid species in wheat.

Identification of differentially methylated regions during vernalization revealed a role for RNA methyltransferases in bolting

Sugar beet (Beta vulgaris altissima) is a biennial root crop with an absolute requirement for cold exposure to bolt and flower, a process called vernalization. Global DNA methylation variations have been reported during vernalization in several plants. However, few genes targeted by DNA methylation during vernalization have been described. The objectives of this study were to identify differentially methylated regions and to study their involvement in bolting induction and tolerance. Restriction landmark genome scanning was applied to DNA from shoot apical meristems of sugar beet genotypes, providing a direct quantitative epigenetic assessment of several CG methylated genes without prior knowledge of gene sequence. Several differentially methylated regions exhibiting variations of gene-body DNA methylation and expression during cold exposure and/or between genotypes were identified, including an AROGENATE DEHYDRATASE and two RNA METHYLCYTOSINE TRANSFERASE sequences. One RNA METHYLCYTOSINE TRANSFERASE sequence displayed gene-body hypermethylation and activation of expression, while the other was hypomethylated and inhibited by cold exposure. Global RNA methylation and phenolic compound levels changed during cold exposure in a genotype-dependent way. The use of methyl RNA immunoprecipitation of total RNA and reverse transcription-PCR analysis revealed mRNA methylation in a vernalized bolting-resistant genotype for the FLOWERING LOCUS 1 gene, a repressor of flowering. Finally, Arabidopsis mutants for RNA METHYLCYTOSINE TRANSFERASE and AROGENATE DEHYDRATASE were shown to exhibit, under different environmental conditions, early or late bolting phenotypes, respectively. Overall, the data identified functional targets of DNA methylation during vernalization in sugar beet, and it is proposed that RNA methylation and phenolic compounds play a role in the floral transition.

Mobilization of Stowaway-like MITEs in newly formed allohexaploid wheat species

Transposable elements (TEs) dominate the genetic capacity of most eukaryotes, especially plants, where they can account for up to 90 % of the genome, such as in wheat. The relationship between TEs and their hosts and the role of TEs in organismal biology are poorly understood. In this study, we have applied next generation sequencing, together with a transposon display technique in order to test whether a Stowaway-like MITE, termed Minos, transposes following allopolyploidization events in wheat. We have generated a 454-pyrosequencing database of Minos-specific amplicons (transposon display products) from a newly formed wheat allohexaploid and its parental lines and retrieved hundreds of novel MITE insertions in the allohexaploid. Clear mobilization of Minos was also seen by site-specific PCR analysis and sequence validation. In addition, using real-time qPCR analysis we observed an insignificant change in the relative quantity of Minos from the expected value of merging the two parental genomes, indicating that, despite its activation, no significant burst in Minos copy number can be seen in the newly formed allohexaploid. Interestingly, we found that CCGG sites surrounding Minos underwent massive hypermethylation following the allohexaploidization process. Our data suggest that MITEs have maintained their capacity for activity throughout the evolution of wheat and might be epigenetically deregulated in the first generations following allopolyploidization.

Epigenetic interplay between mouse endogenous retroviruses and host genes

Background

Transposable elements are often the targets of repressive epigenetic modifications such as DNA methylation that, in theory, have the potential to spread toward nearby genes and induce epigenetic silencing. To better understand the role of DNA methylation in the relationship between transposable elements and genes, we assessed the methylation state of mouse endogenous retroviruses (ERVs) located near genes.

Results

We found that ERVs of the ETn/MusD family show decreased DNA methylation when near transcription start sites in tissues where the nearby gene is expressed. ERVs belonging to the IAP family, however, are generally heavily methylated, regardless of the genomic environment and the tissue studied. Furthermore, we found full-length ETn and IAP copies that display differential DNA methylation between their two long terminal repeats (LTRs), suggesting that the environment surrounding gene promoters can prevent methylation of the nearby LTR. Spreading from methylated ERV copies to nearby genes was rarely observed, with the regions between the ERVs and genes apparently acting as a boundary, enriched in H3K4me3 and CTCF, which possibly protects the unmethylated gene promoter. Furthermore, the flanking regions of unmethylated ERV copies harbor H3K4me3, consistent with spreading of euchromatin from the host gene toward ERV insertions.

Conclusions

We have shown that spreading of DNA methylation from ERV copies toward active gene promoters is rare. We provide evidence that genes can be protected from ERV-induced heterochromatin spreading by either blocking the invasion of repressive marks or by spreading euchromatin toward the ERV copy.

Salt stress induced variation in DNA methylation pattern and its influence on gene expression in contrasting rice genotypes

Background

Salinity is a major environmental factor limiting productivity of crop plants including rice in which wide range of natural variability exists. Although recent evidences implicate epigenetic mechanisms for modulating the gene expression in plants under environmental stresses, epigenetic changes and their functional consequences under salinity stress in rice are underexplored. DNA methylation is one of the epigenetic mechanisms regulating gene expression in plant’s responses to environmental stresses. Better understanding of epigenetic regulation of plant growth and response to environmental stresses may create novel heritable variation for crop improvement.

Methodology/Principal Findings

Methylation sensitive amplification polymorphism (MSAP) technique was used to assess the effect of salt stress on extent and patterns of DNA methylation in four genotypes of rice differing in the degree of salinity tolerance. Overall, the amount of DNA methylation was more in shoot compared to root and the contribution of fully methylated loci was always more than hemi-methylated loci. Sequencing of ten randomly selected MSAP fragments indicated gene-body specific DNA methylation of retrotransposons, stress responsive genes, and chromatin modification genes, distributed on different rice chromosomes. Bisulphite sequencing and quantitative RT-PCR analysis of selected MSAP loci showed that cytosine methylation changes under salinity as well as gene expression varied with genotypes and tissue types irrespective of the level of salinity tolerance of rice genotypes.

Conclusions/Significance

The gene body methylation may have an important role in regulating gene expression in organ and genotype specific manner under salinity stress. Association between salt tolerance and methylation changes observed in some cases suggested that many methylation changes are not “directed”. The natural genetic variation for salt tolerance observed in rice germplasm may be independent of the extent and pattern of DNA methylation which may have been induced by abiotic stress followed by accumulation through the natural selection process.

Analysis of copy-number variation, insertional polymorphism, and methylation status of the tiniest class I (TRIM) and class II (MITE) transposable element families in various rice strains

Transposable elements (TEs) dominate the genetic capacity of most eukaryotes, especially plants, where they may compose up to 90% of the genome. Many studies, both in plants and animals reported that in fact non-autonomous elements that have lost their protein-coding sequences and became miniature elements were highly associated with genes, and showed a high level of transpositional activity such as mPing family in rice. In this study, we have investigated in detail the copy number, insertional polymorphism and the methylation status of the tiniest LTR retrotransposon family, termed TRIM, in nine rice strains, in comparison with mPing. While TRIM showed similar copy numbers (average of 79 insertions) in all the nine rice strains, the copy number of mPing varied dramatically (ranging from 6 to 203 insertions) in the same strains. Site-specific PCR analysis revealed that ~58% of the TRIM elements have identical insertion sites among the nine rice strains, while none of the mPing elements (100% polymorphism) have identical insertion sites in the same strains. Finally, over 65% of the TRIM insertion sites were cytosine methylated in all nine rice strains, while the level of the methylated mPing insertion sites ranged between 43 and 81.5%. The findings of this study indicate that unlike mPing, TRIM is most probably a fossil TE family in rice. In addition, the data shows that there might be a strong correlation between TE methylation and copy number.

Dynamic DNA cytosine methylation in the Populus trichocarpa genome: tissue-level variation and relationship to gene expression

Background

DNA cytosine methylation is an epigenetic modification that has been implicated in many biological processes. However, large-scale epigenomic studies have been applied to very few plant species, and variability in methylation among specialized tissues and its relationship to gene expression is poorly understood.

Results

We surveyed DNA methylation from seven distinct tissue types (vegetative bud, male inflorescence [catkin], female catkin, leaf, root, xylem, phloem) in the reference tree species black cottonwood (Populus trichocarpa). Using 5-methyl-cytosine DNA immunoprecipitation followed by Illumina sequencing (MeDIP-seq), we mapped a total of 129,360,151 36- or 32-mer reads to the P. trichocarpa reference genome. We validated MeDIP-seq results by bisulfite sequencing, and compared methylation and gene expression using published microarray data. Qualitative DNA methylation differences among tissues were obvious on a chromosome scale. Methylated genes had lower expression than unmethylated genes, but genes with methylation in transcribed regions ("gene body methylation") had even lower expression than genes with promoter methylation. Promoter methylation was more frequent than gene body methylation in all tissues except male catkins. Male catkins differed in demethylation of particular transposable element categories, in level of gene body methylation, and in expression range of genes with methylated transcribed regions. Tissue-specific gene expression patterns were correlated with both gene body and promoter methylation.

Conclusions

We found striking differences among tissues in methylation, which were apparent at the chromosomal scale and when genes and transposable elements were examined. In contrast to other studies in plants, gene body methylation had a more repressive effect on transcription than promoter methylation.

Bioinformatics and genomic analysis of transposable elements in eukaryotic genomes

A major portion of most eukaryotic genomes are transposable elements (TEs). During evolution, TEs have introduced profound changes to genome size, structure, and function. As integral parts of genomes, the dynamic presence of TEs will continue to be a major force in reshaping genomes. Early computational analyses of TEs in genome sequences focused on filtering out "junk" sequences to facilitate gene annotation. When the high abundance and diversity of TEs in eukaryotic genomes were recognized, these early efforts transformed into the systematic genome-wide categorization and classification of TEs. The availability of genomic sequence data reversed the classical genetic approaches to discovering new TE families and superfamilies. Curated TE databases and their accurate annotation of genome sequences in turn facilitated the studies on TEs in a number of frontiers including: (1) TE-mediated changes of genome size and structure, (2) the influence of TEs on genome and gene functions, (3) TE regulation by host, (4) the evolution of TEs and their population dynamics, and (5) genomic scale studies of TE activity. Bioinformatics and genomic approaches have become an integral part of large-scale studies on TEs to extract information with pure in silico analyses or to assist wet lab experimental studies. The current revolution in genome sequencing technology facilitates further progress in the existing frontiers of research and emergence of new initiatives. The rapid generation of large-sequence datasets at record low costs on a routine basis is challenging the computing industry on storage capacity and manipulation speed and the bioinformatics community for improvement in algorithms and their implementations.