Environmental and genealogical effects on DNA methylation in a widespread apomictic dandelion lineage.

DNA methylation in plant genomes occurs in different sequences and genomic contexts that have very different properties. DNA methylation that occurs in CG (mCG) sequence context shows transgenerational stability and high epimutation rate, and can thus provide genealogical information at short time scales. However, due to meta-stability and because mCG variants may arise due to other factors than epimutation, such as environmental stress exposure, it is not clear how well mCG captures genealogical information at micro-evolutionary time scales. Here, we analysed DNA methylation variation between accessions from a geographically widespread, apomictic common dandelion (Taraxacum officinale) lineage when grown experimentally under different light conditions. Using a reduced-representation bisulphite sequencing approach, we show that the light treatment induced differentially methylated cytosines (DMCs) in all sequence contexts, with a bias towards transposable elements. Accession differences were associated mainly with DMCs in CG context. Hierarchical clustering of samples based on total mCG profiles revealed a perfect clustering of samples by accession identity, irrespective of light conditions. Using microsatellite information as a benchmark of genetic divergence within the clonal lineage, we show that genetic divergence between accessions correlates strongly with overall mCG profiles. However, our results suggest that environmental effects that do occur in CG context may produce a heritable signal that partly dilutes the genealogical signal. Our study shows that methylation information in plants can be used to reconstruct micro-evolutionary genealogy, providing a useful tool in systems that lack genetic variation such as clonal and vegetatively propagated plants.

Inheritance of DNA methylation differences in the mangrove Rhizophora mangle.

The capacity to respond to environmental challenges ultimately relies on phenotypic variation which manifests from complex interactions of genetic and nongenetic mechanisms through development. While we know something about genetic variation and structure of many species of conservation importance, we know very little about the nongenetic contributions to variation. Rhizophora mangle is a foundation species that occurs in coastal estuarine habitats throughout the neotropics where it provides critical ecosystem functions and is potentially threatened by anthropogenic environmental changes. Several studies have documented landscape-level patterns of genetic variation in this species, but we know virtually nothing about the inheritance of nongenetic variation. To assess one type of nongenetic variation, we examined the patterns of DNA sequence and DNA methylation in maternal plants and offspring from natural populations of R. mangle from the Gulf Coast of Florida. We used a reduced representation bisulfite sequencing approach (epi-genotyping by sequencing; epiGBS) to address the following questions: (a) What are the levels of genetic and epigenetic diversity in natural populations of R. mangle? (b) How are genetic and epigenetic variation structured within and among populations? (c) How faithfully is epigenetic variation inherited? We found low genetic diversity but high epigenetic diversity from natural populations of maternal plants in the field. In addition, a large portion (up to ~25%) of epigenetic differences among offspring grown in common garden was explained by maternal family. Therefore, epigenetic variation could be an important source of response to challenging environments in the genetically depauperate populations of this foundation species.

Genotyping by sequencing for estimating relative abundances of diatom taxa in mock communities.

BACKGROUND: Diatoms are present in all waters and are highly sensitive to pollution gradients. Therefore, they are ideal bioindicators for water quality assessment. Current indices used in these applications are based on identifying diatom species and counting their abundances using traditional light microscopy. Several molecular techniques have been developed to help automate different steps of this process, but obtaining reliable estimates of diatom community composition and species abundance remains challenging. RESULTS: Here, we evaluated a recently developed quantification method based on Genotyping by Sequencing (GBS) for the first time in diatoms to estimate the relative abundances within a species complex. For this purpose, a reference database comprised of thousands of genomic DNA clusters was generated from cultures of Nitzschia palea. The sequencing reads from calibration and mock samples were mapped against this database for parallel quantification. We sequenced 25 mock diatom communities containing up to five taxa per sample in different abundances. Taxon abundances in these communities were also quantified by a diatom expert using manual counting of cells on light microscopic slides. The relative abundances of strains across mock samples were over- or under-estimated by the manual counting method, and a majority of mock samples had stronger correlations using GBS. Moreover, one previously recognized putative hybrid had the largest number of false positive detections demonstrating the limitation of the manual counting method when morphologically similar and/or phylogenetically close taxa are analyzed. CONCLUSIONS: Our results suggest that GBS is a reliable method to estimate the relative abundances of the N. palea taxa analyzed in this study and outperformed traditional light microscopy in terms of accuracy. GBS provides increased taxonomic resolution compared to currently available quantitative molecular approaches, and it is more scalable in the number of species that can be analyzed in a single run. Hence, this is a significant step forward in developing automated, high-throughput molecular methods specifically designed for the quantification of [diatom] communities for freshwater quality assessments.

msGBS: A new high-throughput approach to quantify the relative species abundance in root samples of multispecies plant communities.

Plant interactions are as important belowground as aboveground. Belowground plant interactions are however inherently difficult to quantify, as roots of different species are difficult to disentangle. Although for a couple of decades molecular techniques have been successfully applied to quantify root abundance, root identification and quantification in multispecies plant communities remains particularly challenging. Here we present a novel methodology, multispecies genotyping by sequencing (msGBS), as a next step to tackle this challenge. First, a multispecies meta-reference database containing thousands of gDNA clusters per species is created from GBS derived High Throughput Sequencing (HTS) reads. Second, GBS derived HTS reads from multispecies root samples are mapped to this meta-reference which, after a filter procedure to increase the taxonomic resolution, allows the parallel quantification of multiple species. The msGBS signal of 111 mock-mixture root samples, with up to 8 plant species per sample, was used to calculate the within-species abundance. Optional subsequent calibration yielded the across-species abundance. The within- and across-species abundances highly correlated (R2 range 0.72-0.94 and 0.85-0.98, respectively) to the biomass-based species abundance. Compared to a qPCR based method which was previously used to analyse the same set of samples, msGBS provided similar results. Additional data on 11 congener species groups within 105 natural field root samples showed high taxonomic resolution of the method. msGBS is highly scalable in terms of sensitivity and species numbers within samples, which is a major advantage compared to the qPCR method and advances our tools to reveal hidden belowground interactions.

The cytosine methylation landscape of spring barley revealed by a new reduced representation bisulfite sequencing pipeline, WellMeth.

Patterns and level of cytosine methylation vary widely among plant species and are associated with genome size as well as the proportion of transposons and other repetitive elements in the genome. We explored epigenetic patterns and diversity in a representative proportion of the spring barley (Hordeum vulgare L.) genome across several commercial and historical cultivars. This study adapted a genotyping-by-sequencing (GBS) approach for the detection of methylated cytosines in genomic DNA. To analyze the data, we developed WellMeth, a complete pipeline for analysis of reduced representation bisulfite sequencing. WellMeth enabled quantification of context-specific DNA methylation at the single-base resolution as well as identification of differentially methylated sites (DMCs) and regions (DMRs). On average, DNA methylation levels were significantly higher than what is commonly observed in many plants species, reaching over 10-fold higher levels than those in Arabidopsis thaliana (L.) Heynh. in the CHH methylation. Preferential methylation was observed within and at the edges of long-terminal repeats (LTR) retrotransposons Gypsy and Copia. From a pairwise comparison of cultivars, numerous DMRs could be identified of which more than 5,000 were conserved within the analyzed set of barley cultivars. The subset of regions overlapping with genes showed enrichment in gene ontology (GO) categories associated with chromatin and cellular structure and organization. A significant correlation between genetic and epigenetic distances suggests that a considerable portion of methylated regions is under strict genetic control in barley. The data presented herein represents the first step in efforts toward a better understanding of genome-level structural and functional aspects of methylation in barley.

Early root overproduction not triggered by nutrients decisive for competitive success belowground.

BACKGROUND: Theory predicts that plant species win competition for a shared resource by more quickly preempting the resource in hotspots and by depleting resource levels to lower concentrations than its competitors. Competition in natural grasslands largely occurs belowground, but information regarding root interactions is limited, as molecular methods quantifying species abundance belowground have only recently become available. PRINCIPAL FINDINGS: In monoculture, the grass Festuca rubra had higher root densities and a faster rate of soil nitrate depletion than Plantago lanceolata, projecting the first as a better competitor for nutrients. However, Festuca lost in competition with Plantago. Plantago not only replaced the lower root mass of its competitor, but strongly overproduced roots: with only half of the plants in mixture than in monoculture, Plantago root densities in mixture were similar or higher than those in its monocultures. These responses occurred equally in a nutrient-rich and nutrient-poor soil layer, and commenced immediately at the start of the experiment when root densities were still low and soil nutrient concentrations high. CONCLUSIONS/SIGNIFICANCE: Our results suggest that species may achieve competitive superiority for nutrients by root growth stimulation prior to nutrient depletion, induced by the presence of a competitor species, rather than by a better ability to compete for nutrients per se. The root overproduction by which interspecific neighbors are suppressed independent of nutrient acquisition is consistent with predictions from game theory. Our results emphasize that root competition may be driven by other mechanisms than is currently assumed. The long-term consequences of these mechanisms for community dynamics are discussed.

The roles of ethylene, auxin, abscisic acid, and gibberellin in the hyponastic growth of submerged Rumex palustris petioles.

Rumex palustris responds to complete submergence with upward movement of the younger petioles. This so-called hyponastic response, in combination with stimulated petiole elongation, brings the leaf blade above the water surface and restores contact with the atmosphere. We made a detailed study of this differential growth process, encompassing the complete range of the known signal transduction pathway: from the cellular localization of differential growth, to the hormonal regulation, and the possible involvement of a cell wall loosening protein (expansin) as a downstream target. We show that hyponastic growth is caused by differential cell elongation across the petiole base, with cells on the abaxial (lower) surface elongating faster than cells on the adaxial (upper) surface. Pharmacological studies and endogenous hormone measurements revealed that ethylene, auxin, abscisic acid (ABA), and gibberellin regulate different and sometimes overlapping stages of hyponastic growth. Initiation of hyponastic growth and (maintenance of) the maximum petiole angle are regulated by ethylene, ABA, and auxin, whereas the speed of the response is influenced by ethylene, ABA, and gibberellin. We found that a submergence-induced differential redistribution of endogenous indole-3-acetic acid in the petiole base could play a role in maintenance of the response, but not in the onset of hyponastic growth. Since submergence does not induce a differential expression of expansins across the petiole base, it is unlikely that this cell wall loosening protein is the downstream target for the hormones that regulate the differential cell elongation leading to submergence-induced hyponastic growth in R. palustris.