The epigenetic footprint of poleward range‐expanding plants in apomictic dandelions

Epigenetic modifications, such as DNA methylation variation, can generate heritable phenotypic variation independent of the underlying genetic code. However, epigenetic variation in natural plant populations is poorly documented and little understood. Here, we test whether northward range expansion of obligate apomicts of the common dandelion (Taraxacum officinale) is associated with DNA methylation variation. We characterized and compared patterns of genetic and DNA methylation variation in greenhouse‐reared offspring of T. officinale that were collected along a latitudinal transect of northward range expansion in Europe. Genetic AFLP and epigenetic MS‐AFLP markers revealed high levels of local diversity and modest but significant heritable differentiation between sampling locations and between the southern, central and northern regions of the transect. Patterns of genetic and epigenetic variation were significantly correlated, reflecting the genetic control over epigenetic variation and/or the accumulation of lineage‐specific spontaneous epimutations, which may be selectively neutral. In addition, we identified a small component of DNA methylation differentiation along the transect that is independent of genetic variation. This epigenetic differentiation might reflect environment‐specific induction or, in case the DNA methylation variation affects relevant traits and fitness, selection of heritable DNA methylation variants. Such generated epigenetic variants might contribute to the adaptive capacity of individual asexual lineages under changing environments. Our results highlight the potential of heritable DNA methylation variation to contribute to population differentiation along ecological gradients. Further studies are needed using higher resolution methods to understand the functional significance of such natural occurring epigenetic differentiation.     

The population genomic basis of geographic differentiation in North American common ragweed (Ambrosia artemisiifolia L.)

Common ragweed (Ambrosia artemisiifolia L.) is an invasive, wind‐pollinated plant nearly ubiquitous in disturbed sites in its eastern North American native range and present across growing portions of Europe, Africa, Asia, and Australia. Phenotypic divergence between European and native‐range populations has been described as rapid evolution. However, a recent study demonstrated major human‐mediated shifts in ragweed genetic structure before introduction to Europe and suggested that native‐range genetic structure and local adaptation might fully explain accelerated growth and other invasive characteristics of introduced populations. Genomic differentiation that potentially influenced this structure has not yet been investigated, and it remains unclear whether substantial admixture during historical disturbance of the native range contributed to the development of invasiveness in introduced European ragweed populations. To investigate fine‐scale population genetic structure across the species' native range, we characterized diallelic SNP loci via a reduced‐representation genotyping‐by‐sequencing (GBS) approach. We corroborate phylogeographic domains previously discovered using traditional sequencing methods, while demonstrating increased power to resolve weak genetic structure in this highly admixed plant species. By identifying exome polymorphisms underlying genetic differentiation, we suggest that geographic differentiation of this important invasive species has occurred more often within pathways that regulate growth and response to defense and stress, which may be associated with survival in North America's diverse climatic regions.     

msap: a tool for the statistical analysis of methylation‐sensitive amplified polymorphism data

In this study msap, an R package which analyses methylation‐sensitive amplified polymorphism (MSAP or MS‐AFLP) data is presented. The program provides a deep analysis of epigenetic variation starting from a binary data matrix indicating the banding pattern between the isoesquizomeric endonucleases HpaII and MspI, with differential sensitivity to cytosine methylation. After comparing the restriction fragments, the program determines if each fragment is susceptible to methylation (representative of epigenetic variation) or if there is no evidence of methylation (representative of genetic variation). The package provides, in a user‐friendly command line interface, a pipeline of different analyses of the variation (genetic and epigenetic) among user‐defined groups of samples, as well as the classification of the methylation occurrences in those groups. Statistical testing provides support to the analyses. A comprehensive report of the analyses and several useful plots could help researchers to assess the epigenetic and genetic variation in their MSAP experiments. msap is downloadable from CRAN ( http://cran.r-project.org/ ) and its own webpage ( http://msap.r-forge.R-project.org/ ). The package is intended to be easy to use even for those people unfamiliar with the R command line environment. Advanced users may take advantage of the available source code to adapt msap to more complex analyses.     

Utilization of UPLC/Q‐TOF‐MS‐Based Metabolomics and AFLP‐Based Marker‐Assisted Selection to Facilitate/Assist Conventional Breeding of Polygala tenuifolia

As one of the most important traditional Chinese medicine, the quality of Polygala tenuifolia is difficult to control and a new method must be established to facilitate/assist the breeding of P. tenuifolia. In this study, UPLC/Q‐TOF‐MS‐based metabolomics analysis was performed to determine the chemical composition and screen metabolite biomarkers according to agronomic traits. A total of 29 compounds and 18 metabolite biomarkers were found. AFLP‐based marker‐assisted selection (MAS) was used to identify molecular marker bands and screen characteristic bands associated with specific agronomic traits. 184 bands and 76 characteristic AFLP bands were found. The correlation network between compounds and characteristic AFLP bands was built, so we may directly breed certain P. tenuifolia herbs with special agronomic traits (or characteristic AFLP bands), which exhibit specific pharmacological functions depending on the content of the active compounds. The proposed method of metabolomics coupled with MAS could facilitate/assist the breeding of P. tenuifolia.     

MSAP markers and global cytosine methylation in plants: a literature survey and comparative analysis for a wild‐growing species

Methylation of DNA cytosines affects whether transposons are silenced and genes are expressed, and is a major epigenetic mechanism whereby plants respond to environmental change. Analyses of methylation‐sensitive amplification polymorphism (MS‐AFLP or MSAP) have been often used to assess methyl‐cytosine changes in response to stress treatments and, more recently, in ecological studies of wild plant populations. MSAP technique does not require a sequenced reference genome and provides many anonymous loci randomly distributed over the genome for which the methylation status can be ascertained. Scoring of MSAP data, however, is not straightforward, and efforts are still required to standardize this step to make use of the potential to distinguish between methylation at different nucleotide contexts. Furthermore, it is not known how accurately MSAP infers genome‐wide cytosine methylation levels in plants. Here, we analyse the relationship between MSAP results and the percentage of global cytosine methylation in genomic DNA obtained by HPLC analysis. A screening of literature revealed that methylation of cytosines at cleavage sites assayed by MSAP was greater than genome‐wide estimates obtained by HPLC, and percentages of methylation at different nucleotide contexts varied within and across species. Concurrent HPLC and MSAP analyses of DNA from 200 individuals of the perennial herb Helleborus foetidus confirmed that methyl‐cytosine was more frequent in CCGG contexts than in the genome as a whole. In this species, global methylation was unrelated to methylation at the inner CG site. We suggest that global HPLC and context‐specific MSAP methylation estimates provide complementary information whose combination can improve our current understanding of methylation‐based epigenetic processes in nonmodel plants.     

Insulin‐like growth factors and their potential role in cardiac epigenetics

Cardiovascular disease (CVD) constitutes a major public health threat worldwide, accounting for 17.3 million deaths annually. Heart disease and stroke account for the majority of healthcare costs in the developed world. While much has been accomplished in understanding the pathophysiology, molecular biology and genetics underlying the diagnosis and treatment of CVD, we know less about the role of epigenetics and their molecular determinants. The impact of environmental changes and epigenetics in CVD is now emerging as critically important in understanding the origin of disease and the development of new therapeutic approaches to prevention and treatment. This review focuses on the emerging role of epigenetics mediated by insulin like‐growth factors‐I and ‐II in major CVDs such as heart failure, cardiac hypertrophy and diabetes.     

Fine‐scale population epigenetic structure in relation to gastrointestinal parasite load in red grouse (Lagopus lagopus scotica)

Epigenetic modification of cytosine methylation states can be elicited by environmental stresses and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their host, but there is little understanding in how parasites may influence host methylation states. Here, we estimate epigenetic diversity and differentiation among 21 populations of red grouse (Lagopus lagopus scotica) in north‐east Scotland and test for association of gastrointestinal parasite load (caecal nematode Trichostrongylus tenuis) with hepatic genome‐wide and locus‐specific methylation states. Following methylation‐sensitive AFLP (MSAP), 129 bands, representing 73 methylation‐susceptible and 56 nonmethylated epiloci, were scored across 234 individuals. The populations differed significantly in genome‐wide methylation levels and were also significantly epigenetically (F_SC = 0.0227; P < 0.001) and genetically (F_SC = 0.0058; P < 0.001) differentiated. Parasite load was not associated with either genome‐wide methylation levels or epigenetic differentiation. Instead, we found eight disproportionately differentiated epilocus‐specific methylation states (F_ST outliers) using bayescan software and significant positive and negative association of 35 methylation states with parasite load from bespoke generalized estimating equations (GEE), simple logistic regression (sam ) and Bayesian environmental analysis (bayenv 2). Following Sanger sequencing, genome mapping and geneontology (go ) annotation, some of these epiloci were linked to genes involved in regulation of cell cycle, signalling, metabolism, immune system and notably rRNA methylation, histone acetylation and small RNAs. These findings demonstrate an epigenetic signature of parasite load in populations of a wild bird and suggest intriguing physiological effects of parasite‐associated cytosine methylation.     

Honey bee colonies from different races show variation in defenses against the varroa mite in a ‘common garden’

Honey bee [Apis mellifera L. (Hymenoptera: Apidae)] genetic diversity may be the key to responding to novel health challenges faced by this important pollinator. In this study, we first compared colonies of four honey bee races, A. m. anatoliaca, A. m. carnica, A. m. caucasica, and A. m. syriaca from Turkey, with respect to honey storage, bee population size, and defenses against varroa. The mite Varroa destructor Anderson & Trueman (Acari: Varroidae) is an important pest of honey bee colonies. There are genetic correlates with two main defenses of bees against this parasite: hygienic behavior, or removing infested brood, and grooming, which involves shaking and swiping off mites and biting them. In the second part of this study, we examined the relationship of these two types of defenses, hygiene and grooming, and their correlation with infestation rates in 32 genetically diverse colonies in a ‘common garden’ apiary. Mite biting was found to be negatively correlated with mite infestation levels.     

Global patterns of population genetic differentiation in seed plants

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by F_ST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337‐species data set for seed plants from publications with data on F_ST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that F_ST values were higher for tropical, mixed‐mating, non‐woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to F_ST. Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.     

A DNA mini‐barcode for land plants

Small portions of the barcode region – mini‐barcodes – may be used in place of full‐length barcodes to overcome DNA degradation for samples with poor DNA preservation. 591,491,286 rbcL mini‐barcode primer combinations were electronically evaluated for PCR universality, and two novel highly universal sets of priming sites were identified. Novel and published rbcL mini‐barcode primers were evaluated for PCR amplification [determined with a validated electronic simulation (n = 2765) and empirically (n = 188)], Sanger sequence quality [determined empirically (n = 188)], and taxonomic discrimination [determined empirically (n = 30 472)]. PCR amplification for all mini‐barcodes, as estimated by validated electronic simulation, was successful for 90.2–99.8% of species. Overall Sanger sequence quality for mini‐barcodes was very low – the best mini‐barcode tested produced sequences of adequate quality (B₂₀ ≥ 0.5) for 74.5% of samples. The majority of mini‐barcodes provide correct identifications of families in excess of 70.1% of the time. Discriminatory power noticeably decreased at lower taxonomic levels. At the species level, the discriminatory power of the best mini‐barcode was less than 38.2%. For samples believed to contain DNA from only one species, an investigator should attempt to sequence, in decreasing order of utility and probability of success, mini‐barcodes F (rbcL1/rbcLB), D (F52/R193) and K (F517/R604). For samples believed to contain DNA from more than one species, an investigator should amplify and sequence mini‐barcode D (F52/R193).     

Genome‐scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high‐throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20–33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non‐model systems.     

Determination of DNA methylation associated with Acer rubrum (red maple) adaptation to metals: analysis of global DNA modifications and methylation‐sensitive amplified polymorphism

Red maple (Acer rubum), a common deciduous tree species in Northern Ontario, has shown resistance to soil metal contamination. Previous reports have indicated that this plant does not accumulate metals in its tissue. However, low level of nickel and copper corresponding to the bioavailable levels in contaminated soils in Northern Ontario causes severe physiological damages. No differentiation between metal‐contaminated and uncontaminated populations has been reported based on genetic analyses. The main objective of this study was to assess whether DNA methylation is involved in A. rubrum adaptation to soil metal contamination. Global cytosine and methylation‐sensitive amplified polymorphism (MSAP) analyses were carried out in A. rubrum populations from metal‐contaminated and uncontaminated sites. The global modified cytosine ratios in genomic DNA revealed a significant decrease in cytosine methylation in genotypes from a metal‐contaminated site compared to uncontaminated populations. Other genotypes from a different metal‐contaminated site within the same region appear to be recalcitrant to metal‐induced DNA alterations even ≥30 years of tree life exposure to nickel and copper . MSAP analysis showed a high level of polymorphisms in both uncontaminated (77%) and metal‐contaminated (72%) populations. Overall, 205 CCGG loci were identified in which 127 were methylated in either outer or inner cytosine. No differentiation among populations was established based on several genetic parameters tested. The variations for nonmethylated and methylated loci were compared by analysis of molecular variance (AMOVA). For methylated loci, molecular variance among and within populations was 1.5% and 13.2%, respectively. These values were low (0.6% for among populations and 5.8% for within populations) for unmethylated loci. Metal contamination is seen to affect methylation of cytosine residues in CCGG motifs in the A. rubrum populations that were analyzed.     

A novel approach to parasite population genetics: Experimental infection reveals geographic differentiation,recombination and host‐mediated population structure in Pasteuria ramosa,a bacterial parasite of Daphnia

The population structure of parasites is central to the ecology and evolution of host‐parasite systems. Here, we investigate the population genetics of Pasteuria ramosa, a bacterial parasite of Daphnia. We used natural P. ramosa spore banks from the sediments of two geographically well‐separated ponds to experimentally infect a panel of Daphnia magna host clones whose resistance phenotypes were previously known. In this way, we were able to assess the population structure of P. ramosa based on geography, host resistance phenotype and host genotype. Overall, genetic diversity of P. ramosa was high, and nearly all infected D. magna hosted more than one parasite haplotype. On the basis of the observation of recombinant haplotypes and relatively low levels of linkage disequilibrium, we conclude that P. ramosa engages in substantial recombination. Isolates were strongly differentiated by pond, indicating that gene flow is spatially restricted. Pasteuria ramosa isolates within one pond were segregated completely based on the resistance phenotype of the host—a result that, to our knowledge, has not been previously reported for a nonhuman parasite. To assess the comparability of experimental infections with natural P. ramosa isolates, we examined the population structure of naturally infected D. magna native to one of the two source ponds. We found that experimental and natural infections of the same host resistance phenotype from the same source pond were indistinguishable, indicating that experimental infections provide a means to representatively sample the diversity of P. ramosa while reducing the sampling bias often associated with studies of parasite epidemics. These results expand our knowledge of this model parasite, provide important context for the large existing body of research on this system and will guide the design of future studies of this host‐parasite system.     

Variation in DNA methylation transmissibility,genetic heterogeneity and fecundity‐related traits in natural populations of the perennial herb Helleborus foetidus

Inferences about the role of epigenetics in plant ecology and evolution are mostly based on studies of cultivated or model plants conducted in artificial environments. Insights from natural populations, however, are essential to evaluate the possible consequences of epigenetic processes in biologically realistic scenarios with genetically and phenotypically heterogeneous populations. Here, we explore associations across individuals between DNA methylation transmissibility (proportion of methylation‐sensitive loci whose methylation status persists unchanged after male gametogenesis), genetic characteristics (assessed with AFLP markers), seed size variability (within‐plant seed mass variance), and realized maternal fecundity (number of recently recruited seedlings), in three populations of the perennial herb Helleborus foetidus along a natural ecological gradient in southeastern Spain. Plants (sporophytes) differed in the fidelity with which DNA methylation was transmitted to descendant pollen (gametophytes). This variation in methylation transmissibility was associated with genetic differences. Four AFLP loci were significantly associated with transmissibility and accounted collectively for ~40% of its sample‐wide variance. Within‐plant variance in seed mass was inversely related to individual transmissibility. The number of seedlings recruited by individual plants was significantly associated with transmissibility. The sign of the relationship varied between populations, which points to environment‐specific, divergent phenotypic selection on epigenetic transmissibility. Results support the view that epigenetic transmissibility is itself a phenotypic trait whose evolution may be driven by natural selection, and suggest that in natural populations epigenetic and genetic variation are two intertwined, rather than independent, evolutionary factors.     

Two combinatorial patterns of telomere histone marks in plants with canonical and non‐canonical telomere repeats

Telomeres, nucleoprotein structures at the ends of linear eukaryotic chromosomes, are crucial for the maintenance of genome integrity. In most plants, telomeres consist of conserved tandem repeat units comprising the TTTAGGG motif. Recently, non‐canonical telomeres were described in several plants and plant taxons, including the carnivorous plant Genlisea hispidula (TTCAGG/TTTCAGG), the genus Cestrum (Solanaceae; TTTTTTAGGG), and plants from the Asparagales order with either a vertebrate‐type telomere repeat TTAGGG or Allium genus‐specific CTCGGTTATGGG repeat. We analyzed epigenetic modifications of telomeric histones in plants with canonical and non‐canonical telomeres, and further in telomeric chromatin captured from leaves of Nicotiana benthamiana transiently transformed by telomere CRISPR‐dCas9‐eGFP, and of Arabidopsis thaliana stably transformed with TALE_telo C‐3×GFP. Two combinatorial patterns of telomeric histone modifications were identified: (i) an Arabidopsis‐like pattern (A. thaliana, G. hispidula, Genlisea nigrocaulis, Allium cepa, Narcissus pseudonarcissus, Petunia hybrida, Solanum tuberosum, Solanum lycopersicum) with telomeric histones decorated predominantly by H3K9me2; (ii) a tobacco‐like pattern (Nicotiana tabacum, N. benthamiana, C. elegans) with a strong H3K27me3 signal. Our data suggest that epigenetic modifications of plant telomere‐associated histones are related neither to the sequence of the telomere motif nor to the lengths of the telomeres. Nor the phylogenetic position of the species plays the role; representatives of the Solanaceae family are included in both groups. As both patterns of histone marks are compatible with fully functional telomeres in respective plants, we conclude that the described specific differences in histone marks are not critical for telomere functions.     

Transgenerational plasticity as an important mechanism affecting response of clonal species to changing climate

In spite of the increasing number of studies on the importance of transgenerational plasticity for species response to novel environments, its effects on species ability to respond to climate change are still largely unexplored. We study the importance of transgenerational plasticity for response of a clonal species Festuca rubra. Individuals from four natural populations representing two levels of temperature and two levels of precipitation were cultivated in four growth chambers that simulate the temperature and precipitation of origin of the populations (maternal phase). Each population was represented in each growth chamber. After 6 months, single young ramets of these plants were reshuffled among the growth chambers and let to grow for additional 2 months (offspring phase). The results show that transgenerational effects (i.e., maternal phase conditions) significantly modify species response to novel climates, and the direction and intensity of the response depend on the climate of origin of the plants. For traits related to recourse acquisition, the conditions of maternal phase, either alone or in interaction mainly with climate of origin, had stronger effect than the conditions of cultivation. Overall, the maternal climate interacted more intensively with the climate of origin than with the offspring climate. The direction of the effect of the maternal climate was of different directions and intensities depending on plant origin and trait studied. The data demonstrated strong significant effects of conditions during maternal phase on species response to novel climates. These transgenerational affects were, however, not adaptive. Still, transgenerational plasticity may be an important driver of species response to novel conditions across clonal generations. These effects thus need to be carefully considered in future studies exploring species response to novel climates. This will also have strong effects on species performance under increasingly variable climates expected to occur with the climate change.