International Journal of Genomics
Volume 2016 (2016), Article ID 4395153, 7 pages
http://dx.doi.org/10.1155/2016/4395153

Research Article

Differential Methylation of Genomic Regions Associated with Heteroblasty Detected by M&M Algorithm in the Nonmodel Species Eucalyptus globulus Labill.

Rodrigo Hasbún,¹ Carolina Iturra,¹ Soraya Bravo,² Boris Rebolledo-Jaramillo,³ and Luis Valledor⁴

¹Departamento de Silvicultura, Facultad de Ciencias Forestales, Universidad de Concepción, 4070386 Concepción, Chile
²Centro de Biotecnología Vegetal, Facultad de Ciencias Biológicas, Universidad Andrés Bello, 8370146 Santiago, Chile
³Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA 16802, USA
⁴Departamento de Biología de Organismos y Sistemas, Universidad de Oviedo, 33006 Oviedo, Spain

Received 24 December 2015; Revised 26 February 2016; Accepted 16 March 2016

Academic Editor: Cino Pertoldi

Copyright © 2016 Rodrigo Hasbún et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Epigenetic regulation plays important biological roles in plants, including timing of flowering and endosperm development. Little is known about the mechanisms controlling heterochrony (the change in the timing or rate of developmental events during ontogeny) in Eucalyptus globulus. DNA methylation has been proposed as a potential heterochrony regulatory mechanism in model species, but its role during the vegetative phase in E. globulus has not been explored. In order to investigate the molecular mechanisms governing heterochrony in E. globulus, we have developed a workflow aimed at generating high-resolution hypermethylome and hypomethylome maps that have been tested in two stages of vegetative growth phase: juvenile (6-month leaves) and adult (30-month leaves). We used the M&M algorithm, a computational approach that integrates MeDIP-seq and MRE-seq data, to identify differentially methylated regions (DMRs). Thousands of DMRs between juvenile and adult leaves of E. globulus were found. Although further investigations are required to define the loci associated with heterochrony/heteroblasty that are regulated by DNA methylation, these results suggest that locus-specific methylation could be major regulators of vegetative phase change. This information can support future conservation programs, for example, selecting the best methylomes for a determinate environment in a restoration project.

1. Introduction

Tree species are usually able to tolerate a wide range of environmental conditions. In fact, they are able to tolerate a wide range of environmental conditions and, in many cases, extreme seasonal changes [1]. Some of these organisms are able to manifest different phenotypes depending on the environment in which they grow. This phenomenon, called phenotypic plasticity, has been defined as a change in the phenotype expressed by a single genotype in different environments [2, 3]. Phenotypic plasticity evolves to maximize fitness when the environment is variable and increases with latitude [4–6].

Leaf heteroblasty is a significant and abrupt change in form and function that occurs over the maturation process (phase change from seedlings to reproductive individuals) of certain plants [7]. Characteristics commonly affected include leaf form, size, and arrangement. The earlier and later stages of leaf development are named juvenile and adult, respectively. In contrast to phenotypic plasticity, heteroblastic development does not depend on environmental cues. However, the timing or rate of heteroblastic changes, which can be referred to as a type of heterochrony—change of relative timing of events throughout development—can be modified by the environment [8]. Heterochrony has been implicated in plant evolution, because it can impact the physiology, phytochemistry, or resistance to pests and disease of certain plants [9, 10]; however it is a plastic response largely underexplored.

Phenotypic plasticity and heterochrony may interact to produce a pattern of variation in the leaf phenotype, even in organisms with little or no genetic diversity. Most endangered species have lower genetic diversity than related nonendangered species [11], and phenotypic plasticity and heterochrony together can increase the possibilities of adaption/persistence. Therefore, knowledge of the molecular mechanisms regulating both processes could open new alternatives to assist conservation programs.

Epigenetic regulation, in particular DNA methylation, plays an important biological role in plants, including timing of flowering and endosperm development [12]. Transgenerational inheritance of DNA methylation can mediate phenotypic plasticity via novel epialleles and phenotypes within populations/species [13]. In a review by Pascual et al. [3] it was shown that the coordination of genetic and epigenetic mechanisms mediated phenotypic variation in different plants. For example, in populations of Arabidopsis thaliana with experimental alteration of DNA methylation, the overall patterns of variability among the genotypes indicated that epigenetic changes could affect not only the short-term environmental responses of plants, but also the evolutionary potential of important traits and their plasticity [14]. Similarly, a recent study in invertebrates proposed that the absence of germline DNA methylation in genes involved in the response to fluctuating conditions facilitates phenotypic variation, which could contribute to increased adaptive potential [15]. In conifers it has been reported that environment influences a differential DNA methylation during embryogenesis, inducing differential priming of the embryos that causes differential capabilities to adapt to environment [16].

There is little evidence linking epigenetic regulation and heterochrony. Only a few studies focusing on epigenetic changes during leaf differentiation and development have been developed in Arabidopsis [17–19], rice [20], or pine [21] but none of these species are strongly heteroblastic. Environmental cues are perceived as input signals for the microRNA156/SQUAMOSA promoter-binding protein-like (SPL or SBP box) module and act as a quantitative developmental clock of phase transitions in plants [22, 23]. The evidence shows that sugars promote vegetative phase change through their effect on miR156 [24], but other endogenous factors could play additional important roles. The same authors proposed that heritable epigenetic modification of the miR156 precursor and/or additional chromatin structure alterations could regulate heteroblasty [24].

Eucalyptus globulus, a tree of Australian origin but cultivated worldwide, is strongly heteroblastic with clear differences between its juvenile and adult leaves [25]. Jordan et al. [26] found that genetic association of the timing of vegetative phase change with growth rate ranged from positive to negative at different sites. Early phase change may be favored in warm, wet environments to reduce damage produced by leaf fungi, but it may also be favored on exposed dry sites to increase form or plant structure by which it is protected from desiccation [27, 28]. Genome-wide DNA methylation maps of many model organisms have been reported, but in nonmodel organisms like Eucalyptus spp. the methylation patterns remain poorly studied. These types of maps can be applied to a wide range of biological problems, using the analysis of methylation differences between ecotypes or individuals within species [29–31].

The aim of this study is to set up the required methodology and assess the epigenetic changes (hypermethylation and hypomethylation) related to heteroblasty in E. globulus. Specifically, we want to detect differentially methylated regions (DMRs) between juvenile and adult leaves of E. globulus. DMRs are stretches of genomic DNA that have different DNA methylation patterns. Their natural variation could guide the conservation management of the species or the selection of individuals with potentially adequate methylomes—set of modifications of nucleic acid methylation in the genome of an organism—for a discrete environmental condition. We expect that the results and experience from this work could be used for the discovery of key regulators of heterochrony in future studies, which could be used to assist conservation programs of threatened species.

2. Materials and Methods

To investigate the molecular mechanisms governing heterochrony in Eucalyptus globulus, we applied the M&M algorithm [32] to identify DMRs related to heteroblasty of vegetative growth. We used the previously sequenced X46 clone (Mininco SA JGI Project ID: 401875). Ten ramets of the X46 clone were produced by cuttings, and a genetic trial was established in the commune of Renaico, province of Malleco, region of La Araucanía, Chile (latitude −37.67, longitude −72.59). Juvenile leaf material from five plants was harvested after 6 months at nodes 8 to 10, mixed, and stored at −70°C. Plants were grown until the vegetative phase change was evident (after more than 2 years of growth, average 45 nodes), and adult leaf material was collected. Genomic DNA was extracted using the DNeasy Plant Mini Kit (QIAGEN Inc.).

The methylation profiles were determined by DNA sequencing of enriched genomic libraries: (i) hypermethylome (the methylated part of the genome) using immunoprecipitation of methylated DNA (MeDIP-seq) and (ii) hypomethylome (the nonmethylated part of the genome) based on restriction enzymes sensitive to methylation (MRE-seq).

2.1. MeDIP-seq and MRE-seq Library Generation and Sequencing

MeDIP-seq is a large-scale purification technique used to enrich libraries for methylated DNA sequences. It consists of isolating methylated DNA fragments via an antibody raised against 5-methylcytosine (5mC). MRE-seq utilizes a combination of methyl sensitive enzymes to enrich libraries for unmethylated DNA sequences.

Libraries were generated as previously described in Li et al. [32], with minor modifications. For MeDIP-seq, 3500 ng of isolated DNA was sonicated using 26 pulses of 10 s ON/20 s OFF (Sonic Dismembrator model 100, Fisher Scientific) to a fragment size of 100–500 bp, end processed, and ligated to paired-end adapters using NEXTflex PCR free DNA Sequencing Kit (Bioo Scientific). After size selection of 166–566 bp using Agencourt AMPure XP (Beckman Coulter), DNA was heat denatured and then immunoprecipitated using Methylated DNA IP Kit (Zymo Research), using a mouse monoclonal anti-5-methylcytosine antibody according to manufacturer’s instructions. DNA was then purified with Agencourt AMPure XP (Beckman Coulter) and eluted in 25 mL resuspension buffer (10 mM Tris-HCl, pH 8.5). DNA was amplified by 12 cycles of PCR with the standard Illumina index primers. For MRE-seq, three digestion reactions (HpaII, Acil, and Hin6I; Fermentas) were performed in parallel, each with 600 ng of DNA. Ten units of enzyme (except Acil, which uses five units) were initially incubated with DNA for 3 h, and then additional five units of enzyme were added to the digestion for a total of 6 h of digestion time. Digested DNA from the different reactions was combined and purified using ChiP DNA Clean & Concentrator™ (Zymo Research). The purified DNA was end processed and ligated to single-end adapters using NEXTflex PCR free DNA Sequencing Kit (Bioo Scientific). After size selection (166–566 bp) with Agencourt AMPure XP (Beckman Coulter), the DNA was amplified by PCR for six cycles.

MeDIP and MRE libraries were sent to the DNA Sequencing Facility of the Biotechnology Center at the University of Wisconsin. Samples were sequenced on an Illumina HiSeq machine, yielding a total of 204 million MeDIP-seq reads and 236 million MRE-seq reads. The reads were mapped to the latest Eucalyptus grandis genome assembly (v2.0) [33], using BWA-MEM Li [34] with default settings.

2.2. Use of M&M Algorithm to Detect DMRs

We used an algorithm named M&M [32] that integrates data from both MeDIP-seq and MRE-seq to detect DMRs. M&M is available as an R package called “methylMnM.” Briefly, M&M integrates MeDIP-seq and MRE-seq by dynamically scaling, normalizing, and combining the datasets and provides exact value and -value for different sample comparison. The coverage of MeDIP and MRE sequencing data and genomic CpG information were calculated for each 2000 bp genomic bin. Before applying the M&M method, we generated two input files: (1) CpG sites of each window and (2) MRE-CpG sites of each window. To generate file (1) we used the script fasta2bed.py (Computational Genomics Analysis Tools CGAT 0.2.3), and to calculate (2) we used a script created specifically for this project (discussed at https://www.biostars.org/p/86480/). DMRs between developmental stages were identified using the M&M algorithm with default parameters in the R environment (version v.2.12.1). Briefly, the coverage of MeDIP and MRE sequencing data and genomic CpG information were calculated for each 2000 bp genomic bin. Scaffolds were excluded from the analysis. DMRs with a -value of were selected for analysis. Several statistical functions are implemented in the methylMnM package: “MnM.test()” to calculate the probability that the methylation levels of the two samples within each bin were different, “MnM.value()” to estimate -values based on all the values, and “MnM.select-DMR()” to select significant DMRs based on a cutoff of -value < . The output files contained genomic locations of statistically significant DMRs and their MeDIP-seq and MRE-seq values (in RPKM), as well as values and -values. The absolute values of genomic regions are negative log₁₀-transformed -values. If the value is negative, it represents hypermethylation in the vegetative juvenile sample and hypomethylation in the vegetative adult; if the value is positive, it represents hypomethylation in the vegetative juvenile sample and hypermethylation in the vegetative adult. Mapping results and detected DMRs were visually inspected with the GenomeView software [35] using default parameters.

3. Results and Discussion

3.1. Sequencing Results

For each developmental stage, we constructed one sequencing library using two complementary technologies: MeDIP-seq and MRE-seq. These libraries were sequenced to generate 446 million reads in total (Table 1), of which 387 million were mapped to the E. grandis genome covering at least 50%. We note that MeDIP-seq data has less mapping efficiency. Li et al. [32] explained this because relatively more MeDIP-seq reads are derived from repetitive regions of the genome, which are often heavily methylated. Some reads from repetitive regions cannot be mapped uniquely. In line with expectations MRE-seq scores were inversely correlated with MeDIP-seq scores.

Table 1: Summary of mapping statistic from MeDIP-seq and MRE-seq libraries of juvenile and adult leaves.

Total fragments (Tf)
Mapped fragments
% of Tf
Uniquely mapped fragments
% of Tf
Nonredundant uniquely mapped fragments
Genome coverage (%)

MeDIP-seq-juvenile

85.4

49.1

52.3

MeDIP-seq-adult

82.0

35.4

38.3

MRE-seq-juvenile

89.5

69.4

57.0

MRE-seq-adult

88.8

70.6

49.2

Total

86.7

58.3

The two sequencing signals covered nonoverlapping regions as expected, but some adjacent genomics regions showed overlapping signals (Figure 1).

Figure 1: Epigenome view (GenomeView software) of a 138.7 kb region of chromosome 07 (54,020,001–54,100,001) of Eucalyptus grandis reference genome. MeDIP-seq libraries and MRE-seq libraries covered largely nonoverlapping and short overlapping regions.

3.2. Identification of DMRs

The main aim of this work was the determination of the local DNA methylation changes between two developmental phases (i.e., vegetative juvenile and vegetative adult). We only considered CpG methylation and applied the M&M algorithm to our data to identify DMRs. This allowed us to find a total of 1090 putative DMRs (-value < ) between the two developmental stages. More than 70% of DMRs showed increasing DNA methylation levels from vegetative growth phases juvenile to adult. These preliminary results suggest that locus-specific methylation patterns could be an important feature of vegetative heteroblasty control. This increment in the number of methylated regions agrees with the reduction in the number of genes and proteins that is observed in mature leaves (when they reach their full physiological competence) compared to leaves in proliferative stage [36, 37] which have been also related to an increase of global DNA methylation in conifers [21, 38] and Arabidopsis [39].

In future work, candidate DMRs related to genes involved in plastic responses will be validated using more genotypes of E. globulus with heterochronic responses in the expression of heteroblasty. Using quantitative techniques for quantifying methylation (e.g., bisulfite/sequencing of specific genomic locus) we hope to find DMRs that could explain or predict heterochronic responses. Knowledge of genomic loci that regulate heterochrony in E. globulus will improve our understanding of molecular mechanism of this process and can support future conservation programs, for example, selecting the best methylomes for a specific environment in a restoration project. This complete workflow could be easily applied to other nonmodel species with the only requirement of a reference genome. Given that the cost of DNA sequencing is falling [40] and bioinformatics facilities for assembly of massive sequences are increasing quickly [41], generating a reference genome for almost any species is becoming more affordable. The latter also is valid for threatened species.

The introduction of a workflow allowing the study of the epigenetic regulation at genome level will provide new insight that goes further than the mere description of the genes that are regulated during plant development or environmental adaption. For instance, having the possibility of tracking epigenetic changes in the form of metastable epialleles—alleles that are variably expressed in genetically identical individuals due to epigenetic modifications established during early development and are thought to be particularly vulnerable to environmental influences—will allow us to add a new source of variation, which is particularly important for threatened species, which usually exhibit a very low genetic diversity [42, 43]. Although these studies can be done routinely in model species [44, 45], it is particularly difficult to adapt wet and in silico methodologies for nonmodel like trees, since tissues are rich in polyphenols and other contaminant molecules requiring optimizations (e.g., [46–48]) and at the same time bioinformatics pipelines (gene prediction, characterization, and annotation algorithms) still require improvements [49]. In this study we describe a workflow aiming at achieving results comparable to those obtained in Arabidopsis in terms of quality of libraries and quantity of reads. This procedure will allow future studies employing a higher number of samples and experimental situations that permit addressing important issues related to conservation and management of threatened species such as restoration and translocation. Furthermore, the integration of this regulatory layer together with other omic levels has been proved to be useful for explaining adaptive divergences (see Meijón et al. [50] for an example of the power of high-throughput approaches to explore natural variation in tree species).

4. Conclusion

These preliminary results suggest that locus-specific methylation process could be an important feature of vegetative heteroblasty control in E. globulus. The workflow set up in this project opens a promising future for discovering DNA methylation patterns among different tissue types, cell types, and individuals that will help us to explain phenotypic plasticity and adaption capabilities through the basis of a divergent epigenetic regulation. High-throughput epigenomic technology and analytical tools used in this study could be applied to population-based studies of nonmodel plants but principal challenges are experimental design, data analysis, and interpretation of results. The implementation of epigenetic fingerprinting as a support tool in restoration and conservation projects of threatened species requires the discovery of loci involved in adaptive variation.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

The authors thank MININCO Forest Company for providing plant material and land for the field trial. María Francisca Beltrán, Paulina Rivas, Macarena Arellano, Giannina Espinoza, Joseline Valdes, and Daniela Fernandez are gratefully acknowledged for providing technical assistance. The authors thank the DNA Sequencing Facility of the Biotechnology Center at the University of Wisconsin for providing sequencing and support services. Funding for this project was provided by CONICYT-Chile (Project FONDECYT INICIACION no. 11110505).

References

T. W. Crowther, H. B. Glick, K. R. Covey et al., “Mapping tree density at a global scale,” Nature, vol. 525, no. 7568, pp. 201–205, 2015.

L. Gratani, “Plant phenotypic plasticity in response to environmental factors,” Advances in Botany, vol. 2014, Article ID 208747, 17 pages, 2014.

J. Pascual, M. J. Cañal, B. Correia et al., “Can epigenetics help forest plants to adapt to climate change?” in Epigenetics in Plants of Agronomic Importance: Fundamentals and Applications, R. Alvarez-Venegas, C. De la Peña, and J. A. Casas-Mollano, Eds., pp. 125–146, Springer, Cham, Switzerland, 2014.

A. A. Agrawal, “Phenotypic plasticity in the interactions and evolution of species,” Science, vol. 294, no. 5541, pp. 321–326, 2001.

S. Bahrndorff, V. Loeschcke, C. Pertoldi, C. Beier, and M. Holmstrup, “The rapid cold hardening response of Collembola is influenced by thermal variability of the habitat,” Functional Ecology, vol. 23, no. 2, pp. 340–347, 2009.

J. Overgaard, T. N. Kristensen, K. A. Mitchell, and A. A. Hoffmann, “Thermal tolerance in widespread and tropical Drosophila species: does phenotypic plasticity increase with latitude?” The American Naturalist, vol. 178, no. 1, pp. S80–S96, 2011.

G. Zotz, K. Wilhelm, and A. Becker, “Heteroblasty—a review,” The Botanical Review, vol. 77, no. 2, pp. 109–151, 2011.

S. J. Gould, Ontogeny and Phylogeny, Harvard University Press, Cambridge, Mass, USA, 1st edition, 1977.

P. Li and M. O. Johnston, “Heterochrony in plant evolutionary studies through the twentieth century,” Botanical Review, vol. 66, no. 1, pp. 57–88, 2000.

C. J. Hudson, J. S. Freeman, R. C. Jones et al., “Genetic control of heterochrony in Eucalyptus globulus,” G3: Genes, Genomes, Genetics, vol. 4, no. 7, pp. 1235–1245, 2014.

R. Frankham, J. D. Ballou, and D. A. Briscoe, Introduction to Conservation Genetics, Cambridge University Press, Cambridge, UK, 1st edition, 2002.

E. J. Finnegan, W. J. Peacock, and E. S. Dennis, “DNA methylation, a key regulator of plant development and other processes,” Current Opinion in Genetics and Development, vol. 10, no. 2, pp. 217–223, 2000.

Y.-Y. Zhang, M. Fischer, V. Colot, and O. Bossdorf, “Epigenetic variation creates potential for evolution of plant phenotypic plasticity,” New Phytologist, vol. 197, no. 1, pp. 314–322, 2013.

O. Bossdorf, D. Arcuri, C. L. Richards, and M. Pigliucci, “Experimental alteration of DNA methylation affects the phenotypic plasticity of ecologically relevant traits in Arabidopsis thaliana,” Evolutionary Ecology, vol. 24, no. 3, pp. 541–553, 2010.

S. B. Roberts and M. R. Gavery, “Is there a relationship between DNA methylation and phenotypic plasticity in invertebrates?” Frontiers in Physiology, vol. 2, pp. 116–120, 2012.

I. A. Yakovlev, Y. Lee, B. Rotter et al., “Temperature-dependent differential transcriptomes during formation of an epigenetic memory in Norway spruce embryogenesis,” Tree Genetics and Genomes, vol. 10, no. 2, pp. 355–366, 2014.

M. Lafos, P. Kroll, M. L. Hohenstatt, F. L. Thorpe, O. Clarenz, and D. Schubert, “Dynamic regulation of H3K27 trimethylation during arabidopsis differentiation,” PLoS Genetics, vol. 7, no. 4, Article ID e1002040, 2011.

M. Lodha, C. F. Marco, and M. C. P. Timmermans, “The ASYMMETRIC LEAVES complex maintains repression of KNOX homeobox genes via direct recruitment of Polycomb-repressive complex2,” Genes and Development, vol. 27, no. 6, pp. 596–601, 2013.

J. A. Brusslan, A. M. Rus Alvarez-Canterbury, N. U. Nair, J. C. Rice, M. J. Hitchler, and M. Pellegrini, “Genome-wide evaluation of histone methylation changes associated with leaf senescence in Arabidopsis,” PLoS ONE, vol. 7, no. 3, Article ID e33151, 2012.

P. Sui, J. Jin, S. Ye et al., “H3K36 methylation is critical for brassinosteroid-regulated plant growth and development in rice,” Plant Journal, vol. 70, no. 2, pp. 340–347, 2012.

L. Valledor, J. Pascual, M. Meijón, M. Escandón, and M. J. Cañal, “Conserved epigenetic mechanisms could play a key role in regulation of photosynthesis and development-related genes during needle development of Pinus radiata,” PLoS ONE, vol. 10, no. 5, Article ID e0126405, 2015.

H. Wang and H. Wang, “The miR156/SPL module, a regulatory hub and versatile toolbox, gears up crops for enhanced agronomic traits,” Molecular Plant, vol. 8, no. 5, pp. 677–688, 2015.

A. Pulido and P. Laufs, “Co-ordination of developmental processes by small RNAs during leaf development,” Journal of Experimental Botany, vol. 61, no. 5, pp. 1277–1291, 2010.

L. Yang, M. Xu, Y. Koo, J. He, and R. Scott Poethig, “Sugar promotes vegetative phase change in Arabidopsis thaliana by repressing the expression of MIR156A and MIR156C,” eLife, vol. 2013, no. 2, Article ID e00260, 2013.

C. J. Hudson, J. S. Freeman, R. C. Jones et al., “Genetic control of heterochrony in Eucalyptus globulus,” G3: Genes∣Genomes∣Genetic, vol. 4, no. 7, pp. 1235–1245, 2014.

G. J. Jordan, B. M. Potts, P. Chalmers, and R. J. E. Wiltshire, “Quantitative genetic evidence that the timing of vegetative phase change in Eucalyptus globulus ssp. globulus is an adaptive trait,” Australian Journal of Botany, vol. 48, no. 5, pp. 561–567, 2000.

G. W. Dutkowski and B. M. Potts, “Geographic patterns of genetic variation in Eucalyptus globulus ssp. globulus and a revised racial classification,” Australian Journal of Botany, vol. 47, no. 2, pp. 237–263, 1999.

G. J. Jordan, B. M. Potts, and R. J. E. Wiltshire, “Strong, independent, quantitative genetic control of the timing of vegetative phase change and first flowering in Eucalyptus globulus ssp. globulus,” Heredity, vol. 83, no. 2, pp. 179–187, 1999.

N. A. Eckardt, “Deep sequencing maps the maize epigenomic landscape,” Plant Cell, vol. 21, no. 4, pp. 1024–1026, 2009.

F. Roudier, F. K. Teixeira, and V. Colot, “Chromatin indexing in Arabidopsis: an epigenomic tale of tails and more,” Trends in Genetics, vol. 25, no. 11, pp. 511–517, 2009.

X. Li, J. Zhu, F. Hu et al., “Single-base resolution maps of cultivated and wild rice methylomes and regulatory roles of DNA methylation in plant gene expression,” BMC Genomics, vol. 13, article 300, 2012.

D. Li, B. Zhang, X. Xing, and T. Wang, “Combining MeDIP-seq and MRE-seq to investigate genome-wide CpG methylation,” Methods C, vol. 72, pp. 29–40, 2015.

J. Bartholomé, E. Mandrou, A. Mabiala et al., “High-resolution genetic maps of Eucalyptus improve Eucalyptus grandis genome assembly,” New Phytologist, vol. 206, no. 4, pp. 1283–1296, 2015.

H. Li and R. Durbin, “Fast and accurate short read alignment with Burrows-Wheeler transform,” Bioinformatics, vol. 25, no. 14, pp. 1754–1760, 2009.

T. Abeel, T. Van Parys, Y. Saeys, J. Galagan, and Y. Van de Peer, “GenomeView: a next-generation genome browser,” Nucleic Acids Research, vol. 40, no. 2, p. e12, 2012.

L. Valledor, J. V. Jorrín, J. L. Rodríguez et al., “Combined proteomic and transcriptomic analysis identifies differentially expressed pathways associated to Pinus radiata needle maturation,” Journal of Proteome Research, vol. 9, no. 8, pp. 3954–3979, 2010.

M. Andriankaja, S. Dhondt, S. DeBodt et al., “Exit from proliferation during leaf development in Arabidopsis thaliana: a not-so-gradual process,” Developmental Cell, vol. 22, no. 1, pp. 64–78, 2012.

L. Valledor, M. Meijón, R. Hasbún, M. Jesús Cañal, and R. Rodríguez, “Variations in DNA methylation, acetylated histone H4, and methylated histone H3 during Pinus radiata needle maturation in relation to the loss of in vitro organogenic capability,” Journal of Plant Physiology, vol. 167, no. 5, pp. 351–357, 2010.

E. J. Finnegan, W. J. Peacock, and E. S. Dennis, “Reduced DNA methylation in Arabidopsis thaliana results in abnormal plant development,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 16, pp. 8449–8454, 1996.

K. A. Wetterstrand, “DNA Sequencing Costs: Data from the NHGRI Genome Sequencing Program (GSP),” 2016, http://www.genome.gov/sequencingcosts/

R. Ekblom and J. B. W. Wolf, “A field guide to whole-genome sequencing, assembly and annotation,” Evolutionary Applications, vol. 7, no. 9, pp. 1026–1042, 2014.

N. J. Ouborg, P. Vergeer, and C. Mix, “The rough edges of the conservation genetics paradigm for plants,” Journal of Ecology, vol. 94, no. 6, pp. 1233–1248, 2006.

R. Frankham, “Challenges and opportunities of genetic approaches to biological conservation,” Biological Conservation, vol. 143, no. 9, pp. 1919–1927, 2010.

M. J. Ziller, H. Gu, F. Müller et al., “Charting a dynamic DNA methylation landscape of the human genome,” Nature, vol. 500, no. 7463, pp. 477–481, 2013.

S. R. Eichten, R. Briskine, J. Song et al., “Epigenetic and genetic influences on DNA methylation variation in maize populations,” Plant Cell, vol. 25, no. 8, pp. 2783–2797, 2013.

R. Hasbún, L. Valledor, J. L. Rodríguez et al., “HPCE quantification of 5-methyl-2′deoxycytidine in genomic DNA: methodological optimization for chestnut and other woody species,” Plant Physiology and Biochemistry, vol. 46, no. 8-9, pp. 815–822, 2008.

L. Valledor, R. Hasbún, R. Rodríguez, and M. J. Cañal, “RNA-Free DNA extraction protocol from pinus tissues for molecular biology or HPCE/HPLC analyses,” Journal of Plant Biochemistry and Biotechnology, vol. 18, no. 2, pp. 229–232, 2009.

R. Hasbún, C. Iturra, P. Moraga et al., “An efficient and reproducible protocol for production of AFLP markers in tree genomes using fluorescent capillary detection,” Tree Genetics and Genomes, vol. 8, no. 4, pp. 925–931, 2012.

P. Radivojac, W. T. Clark, T. R. Oron et al., “A large-scale evaluation of computational protein function prediction,” Nature Methods, vol. 10, no. 3, pp. 221–227, 2013.

M. Meijón, I. Feito, M. Oravec et al., “Exploring natural variation of Pinus pinaster Aiton using metabolomics: is it possible to identify the region of origin of a pine from its metabolites?” Molecular Ecology, vol. 25, no. 4, pp. 959–976, 2016.