Global transcriptome and gene co-expression network analyses reveal regulatory and non-additive effects of drought and heat stress in grapevine

Despite frequent co-occurrence of drought and heat stress, the molecular mechanisms governing plant responses to these stresses in combination have not often been studied. This is particularly evident in non-model, perennial plants. We conducted large scale physiological and transcriptome analyses to identify genes and pathways associated with grapevine response to drought and/or heat stress during stress progression and recovery. We identified gene clusters with expression correlated to leaf temperature and water stress and five hub genes for the combined stress co-expression network. Several differentially expressed genes were common to the individual and combined stresses, but the majority were unique to the individual or combined stress treatments. These included heat-shock proteins, mitogen-activated kinases, sugar metabolizing enzymes, and transcription factors, while phenylpropanoid biosynthesis and histone modifying genes were unique to the combined stress treatment. Following physiological recovery, differentially expressed genes were found only in plants under heat stress, both alone and combined with drought. Taken collectively, our results suggest that the effect of the combined stress on physiology and gene expression is more severe than that of individual stresses, but not simply additive, and that epigenetic chromatin modifications may play an important role in grapevine responses to combined drought and heat stress.

Phylogenomics reveals extensive introgression and a case of mito-nuclear discordance in the killifish genus Kryptolebias

Introgression is a widespread evolutionary process leading to phylogenetic inconsistencies among distinct parts of the genomes, particularly between mitochondrial and nuclear-based phylogenetic reconstructions (e.g., mito-nuclear discordances). Here, we used mtDNA and genome-wide nuclear sites to provide the first phylogenomic-based hypothesis on the evolutionary relationships within the killifish genus Kryptolebias. In addition, we tested for evidence of past introgression in the genus given the multiple reports of undergoing hybridization between its members. Our mtDNA phylogeny generally agreed with the relationships previously proposed for the genus. However, our reconstruction based on nuclear DNA revealed an unknown lineage - Kryptolebias sp. 'ESP' - as the sister group of the self-fertilizing mangrove killifishes, K. marmoratus and K. hermaphroditus. All individuals sequenced of Kryptolebias sp. 'ESP' had the same mtDNA haplotype commonly observed in K. hermaphroditus, demonstrating a clear case of mito-nuclear discordance. Our analysis further confirmed extensive history of introgression between Kryptolebias sp. 'ESP' and K. hermaphroditus. Population genomics analyses indicate no current gene flow between the two lineages, despite their current sympatry and history of introgression. We also confirmed introgression between other species pairs in the genus that have been recently reported to form hybrid zones. Overall, our study provides a phylogenomic reconstruction covering most of the Kryptolebias species, reveals a new lineage hidden in a case of mito-nuclear discordance, and provides evidence of multiple events of ancestral introgression in the genus. These findings underscore the importance of investigating different genomic information in a phylogenetic framework, particularly in taxa where introgression is common as in the sexually diverse mangrove killifishes.

Tissue and regional expression patterns of dicistronic tRNA-mRNA transcripts in grapevine (Vitis vinifera) and their evolutionary co-appearance with vasculature in land plants

Transfer RNAs (tRNA) are crucial adaptor molecules between messenger RNA (mRNA) and amino acids. Recent evidence in plants suggests that dicistronic tRNA-like structures also act as mobile signals for mRNA transcripts to move between distant tissues. Co-transcription is not a common feature in the plant nuclear genome and, in the few cases where polycistronic transcripts have been found, they include non-coding RNA species, such as small nucleolar RNAs and microRNAs. It is not known, however, the extent to which dicistronic transcripts of tRNA and mRNAs are expressed in field-grown plants, or the factors contributing to their expression. We analysed tRNA-mRNA dicistronic transcripts in the major horticultural crop grapevine (Vitis vinifera) using a novel pipeline developed to identify dicistronic transcripts from high-throughput RNA-sequencing data. We identified dicistronic tRNA-mRNA in leaf and berry samples from 22 commercial vineyards. Of the 124 tRNA genes that were expressed in both tissues, 18 tRNA were expressed forming part of 19 dicistronic tRNA-mRNAs. The presence and abundance of dicistronic molecules was tissue and geographic sub-region specific. In leaves, the expression patterns of dicistronic tRNA-mRNAs significantly correlated with tRNA expression, suggesting that their transcriptional regulation might be linked. We also found evidence of syntenic genomic arrangements of tRNAs and protein-coding genes between grapevine and Arabidopsis thaliana, and widespread prevalence of dicistronic tRNA-mRNA transcripts among vascular land plants but no evidence of these transcripts in non-vascular lineages. This suggests that the appearance of plant vasculature and tRNA-mRNA occurred concurrently during the evolution of land plants.

Wine Terroir and the Soil Bacteria: An Amplicon Sequencing-Based Assessment of the Barossa Valley and Its Sub-Regions

A wines' terroir, represented as wine traits with regional distinctiveness, is a reflection of both the biophysical and human-driven conditions in which the grapes were grown and wine made. Soil is an important factor contributing to the uniqueness of a wine produced by vines grown in specific conditions. Here, we evaluated the impact of environmental variables on the soil bacteria of 22 Barossa Valley vineyard sites based on the 16S rRNA gene hypervariable region 4. In this study, we report that both dispersal isolation by geographic distance and environmental heterogeneity (soil plant-available P content, elevation, rainfall, temperature, spacing between row and spacing between vine) contribute to microbial community dissimilarity between vineyards. Vineyards located in cooler and wetter regions showed lower beta diversity and a higher ratio of dominant taxa. Differences in soil bacterial community composition were significantly associated with differences in fruit and wine composition. Our results suggest that environmental factors affecting wine terroir, may be mediated by changes in microbial structure, thus providing a basic understanding of how growing conditions affect interactions between plants and their soil bacteria.

Common garden experiment reveals altered nutritional values and DNA methylation profiles in micropropagated three elite Ghanaian sweet potato genotypes

Micronutrient deficiency is the cause of multiple diseases in developing countries. Staple crop biofortification is an efficient means to combat such deficiencies in the diets of local consumers. Biofortified lines of sweet potato (Ipomoea batata L. Lam) with enhanced beta-carotene content have been developed in Ghana to alleviate Vitamin A Deficiency. These genotypes are propagated using meristem micropropagation to ensure the generation of virus-free propagules. In vitro culture exposes micropropagated plants to conditions that can lead to the accumulation of somaclonal variation with the potential to generate unwanted aberrant phenotypes. However, the effect of micropropagation induced somaclonal variation on the production of key nutrients by field-grown plants has not been previously studied. Here we assessed the extent of in vitro culture induced somaclonal variation, at a phenotypic, compositional and genetic/epigenetic level, by comparing field-maintained and micropropagated lines of three elite Ghanaian sweet potato genotypes grown in a common garden. Although micropropagated plants presented no observable morphological abnormalities compared to field maintained lines, they presented significantly lower levels of iron, total protein, zinc, and glucose. Methylation Sensitive Amplification Polymorphism analysis showed a high level of in vitro culture induced molecular variation in micropropagated plants. Epigenetic, rather than genetic variation, accounts for most of the observed molecular variability. Taken collectively, our results highlight the importance of ensuring the clonal fidelity of the micropropagated biofortified lines in order to reduce potential losses in the nutritional value prior to their commercial release.

A Concise Review on Multi-Omics Data Integration for Terroir Analysis in Vitis vinifera

Vitis vinifera (grapevine) is one of the most important fruit crops, both for fresh consumption and wine and spirit production. The term terroir is frequently used in viticulture and the wine industry to relate wine sensory attributes to its geographic origin. Although, it can be cultivated in a wide range of environments, differences in growing conditions have a significant impact on fruit traits that ultimately affect wine quality. Understanding how fruit quality and yield are controlled at a molecular level in grapevine in response to environmental cues has been a major driver of research. Advances in the area of genomics, epigenomics, transcriptomics, proteomics and metabolomics, have significantly increased our knowledge on the abiotic regulation of yield and quality in many crop species, including V. vinifera. The integrated analysis of multiple 'omics' can give us the opportunity to better understand how plants modulate their response to different environments. However, 'omics' technologies provide a large amount of biological data and its interpretation is not always straightforward, especially when different 'omic' results are combined. Here we examine the current strategies used to integrate multi-omics, and how these have been used in V. vinifera. In addition, we also discuss the importance of including epigenomics data when integrating omics data as epigenetic mechanisms could play a major role as an intermediary between the environment and the genome.

Epigenetic-induced alterations in sex-ratios in response to climate change: An epigenetic trap?

We hypothesize that under the predicted scenario of climate change epigenetically mediated environmental sex determination could become an epigenetic trap. Epigenetically regulated environmental sex determination is a mechanism by which species can modulate their breeding strategies to accommodate environmental change. Growing evidence suggests that epigenetic mechanisms may play a key role in phenotypic plasticity and in the rapid adaptation of species to environmental change, through the capacity of organisms to maintain a non-genetic plastic memory of the environmental and ecological conditions experienced by their parents. However, inherited epigenetic variation could also be maladaptive, becoming an epigenetic trap. This is because environmental sex determination can alter sex ratios by increasing the survival of one of the sexes at the expense of negative fitness consequences for the other, which could lead not only to the collapse of natural populations, but also have an impact in farmed animal and plant species.

Meristem micropropagation of cassava (Manihot esculenta) evokes genome-wide changes in DNA methylation

There is great interest in the phenotypic, genetic and epigenetic changes associated with plant in vitro culture known as somaclonal variation. In vitro propagation systems that are based on the use of microcuttings or meristem cultures are considered analogous to clonal cuttings and so widely viewed to be largely free from such somaclonal effects. In this study, we surveyed for epigenetic changes during propagation by meristem culture and by field cuttings in five cassava (Manihot esculenta) cultivars. Principal Co-ordinate Analysis of profiles generated by methylation-sensitive amplified polymorphism revealed clear divergence between samples taken from field-grown cuttings and those recovered from meristem culture. There was also good separation between the tissues of field samples but this effect was less distinct among the meristem culture materials. Application of methylation-sensitive Genotype by sequencing identified 105 candidate epimarks that distinguish between field cutting and meristem culture samples. Cross referencing the sequences of these epimarks to the draft cassava genome revealed 102 sites associated with genes whose homologs have been implicated in a range of fundamental biological processes including cell differentiation, development, sugar metabolism, DNA methylation, stress response, photosynthesis, and transposon activation. We explore the relevance of these findings for the selection of micropropagation systems for use on this and other crops.

Epi-fingerprinting and epi-interventions for improved crop production and food quality

Increasing crop production at a time of rapid climate change represents the greatest challenge facing contemporary agricultural research. Our understanding of the genetic control of yield derives from controlled field experiments designed to minimize environmental variance. In spite of these efforts there is substantial residual variability among plants attributable to Genotype × Environment interactions. Recent advances in the field of epigenetics have revealed a plethora of gene control mechanisms that could account for much of this unassigned variation. These systems act as a regulatory interface between the perception of the environment and associated alterations in gene expression. Direct intervention of epigenetic control systems hold the enticing promise of creating new sources of variability that could enhance crop performance. Equally, understanding the relationship between various epigenetic states and responses of the crop to specific aspects of the growing environment (epigenetic fingerprinting) could allow for a more tailored approach to plant agronomy. In this review, we explore the many ways in which epigenetic interventions and epigenetic fingerprinting can be deployed for the improvement of crop production and quality.

Epigenetic rather than genetic factors may explain phenotypic divergence between coastal populations of diploid and tetraploid Limonium spp. (Plumbaginaceae) in Portugal

BACKGROUND

The genus Limonium Miller comprises annual and perennial halophytes that can produce sexual and/or asexual seeds (apomixis). Genetic and epigenetic (DNA methylation) variation patterns were investigated in populations of three phenotypically similar putative sexual diploid species (L. nydeggeri, L. ovalifolium, L. lanceolatum), one sexual tetraploid species (L. vulgare) and two apomict tetraploid species thought to be related (L. dodartii, L. multiflorum). The extent of morphological differentiation between these species was assessed using ten diagnostic morphometric characters.

RESULTS

A discriminant analysis using the morphometric variables reliably assigns individuals into their respective species groups. We found that only modest genetic and epigenetic differentiation was revealed between species by Methylation Sensitive Amplification Polymorphism (MSAP). However, whilst there was little separation possible between ploidy levels on the basis of genetic profiles, there was clear and pronounced interploidy discrimination on the basis of epigenetic profiles. Here we investigate the relative contribution of genetic and epigenetic factors in explaining the complex phenotypic variability seen in problematic taxonomic groups such as Limonium that operate both apomixis and sexual modes of reproduction.

CONCLUSIONS

Our results suggest that epigenetic variation might be one of the drivers of the phenotypic divergence between diploid and tetraploid taxa and discuss that intergenome silencing offers a plausible mechanistic explanation for the observed phenotypic divergence between these microspecies. These results also suggest that epigenetic profiling offer an additional tool to infer ploidy level in stored specimens and that stable epigenetic change may play an important role in apomict evolution and species recognition.

Pre-conditioning the epigenetic response to high vapor pressure deficit increases the drought tolerance of Arabidopsis thaliana

Epigenetic modification of the genome via cytosine methylation is a dynamic process that responds to changes in the growing environment. This modification can also be heritable. The combination of both properties means that there is the potential for the life experiences of the parental generation to modify the methylation profiles of their offspring and so potentially to “pre-condition” them to better accommodate abiotic conditions encountered by their parents. We recently identified high vapor pressure deficit (vpd)-induced DNA methylation at 2 gene loci in the stomatal development pathway and an associated reduction in leaf stomatal frequency.¹ Here, we test whether this epigenetic modification pre-conditioned parents and their offspring to the more severe water stress of periodic drought. We found that 3 generations of high vpd-grown plants were better able to withstand periodic drought stress over 2 generations. This resistance was not directly associated with de novo methylation of the target stomata genes, but was associated with the cmt3 mutant’s inability to maintain asymmetric sequence context methylation. If our finding applies widely, it could have significant implications for evolutionary biology and breeding for stressful environments.

Transgenerational, dynamic methylation of stomata genes in response to low relative humidity

Transgenerational inheritance of abiotic stress-induced epigenetic modifications in plants has potential adaptive significance and might condition the offspring to improve the response to the same stress, but this is at least partly dependent on the potency, penetrance and persistence of the transmitted epigenetic marks. We examined transgenerational inheritance of low Relative Humidity-induced DNA methylation for two gene loci in the stomatal developmental pathway in Arabidopsis thaliana and the abundance of associated short-interfering RNAs (siRNAs). Heritability of low humidity-induced methylation was more predictable and penetrative at one locus (SPEECHLESS, entropy ≤ 0.02; χ² < 0.001) than the other (FAMA, entropy ≤ 0.17; χ² ns). Methylation at SPEECHLESS correlated positively with the continued presence of local siRNAs (r² = 0.87; p = 0.013) which, however, could be disrupted globally in the progeny under repeated stress. Transgenerational methylation and a parental low humidity-induced stomatal phenotype were heritable, but this was reversed in the progeny under repeated treatment in a previously unsuspected manner.

Differential effect of three base modifications on DNA thermostability revealed by high resolution melting

High resolution melting (HRM) can detect and quantify the presence of 5-methylcytosine (5mC) in DNA samples, but the ability of HRM to diagnose other DNA modifications remains unexplored. The DNA bases N6-methyladenine and 5-hydroxymethylcytosine occur across almost all phyla. While their function remains controversial, their presence perturbs DNA structure. Such modifications could affect gene regulation, chromatin condensation and DNA packaging. Here, we reveal that DNA containing N6-methyladenine or 5-hydroxymethylcytosine exhibits reduced thermal stability compared to cytosine-methylated DNA. These thermostability changes are sufficiently divergent to allow detection and quantification by HRM analysis. Thus, we report that HRM distinguishes between sequence-identical DNA differing only in the modification type of one base. This approach is also able to distinguish between two DNA fragments carrying both N6-methyladenine and 5-methylcytosine but differing only in the distance separating the modified bases. This finding provides scope for the development of new methods to characterize DNA chemically and to allow for low cost screening of mutant populations of genes involved in base modification. More fundamentally, contrast between the thermostabilizing effects of 5mC on dsDNA compared with the destabilizing effects of N6-methyladenine (m6A) and 5-hydroxymethylcytosine (5hmC) raises the intriguing possibility of an antagonistic relationship between modification types with functional significance.

Low relative humidity triggers RNA-directed de novo DNA methylation and suppression of genes controlling stomatal development

Environmental cues influence the development of stomata on the leaf epidermis, and allow plants to exert plasticity in leaf stomatal abundance in response to the prevailing growing conditions. It is reported that Arabidopsis thaliana 'Landsberg erecta' plants grown under low relative humidity have a reduced stomatal index and that two genes in the stomatal development pathway, SPEECHLESS and FAMA, become de novo cytosine methylated and transcriptionally repressed. These environmentally-induced epigenetic responses were abolished in mutants lacking the capacity for de novo DNA methylation, for the maintenance of CG methylation, and in mutants for the production of short-interfering non-coding RNAs (siRNAs) in the RNA-directed DNA methylation pathway. Induction of methylation was quantitatively related to the induction of local siRNAs under low relative humidity. Our results indicate the involvement of both transcriptional and post-transcriptional gene suppression at these loci in response to environmental stress. Thus, in a physiologically important pathway, a targeted epigenetic response to a specific environmental stress is reported and several of its molecular, mechanistic components are described, providing a tractable platform for future epigenetics experiments. Our findings suggest epigenetic regulation of stomatal development that allows for anatomical and phenotypic plasticity, and may help to explain at least some of the plant's resilience to fluctuating relative humidity.

Direct detection and quantification of methylation in nucleic acid sequences using high-resolution melting analysis

High-resolution melting (HRM) analysis exploits the reduced thermal stability of DNA fragments that contain base mismatches to detect single nucleotide polymorphisms (SNPs). However, the capacity of HRM to reveal other features of DNA chemistry remains unexplored. DNA methylation plays a key role in regulating gene expression and is essential for normal development in many higher organisms. The presence of methylated bases perturbs the double-stranded DNA structure, although its effect on DNA thermal stability is largely unknown. Here, we reveal that methylated DNA has enhanced thermal stability and is sufficiently divergent from nonmethylated DNA to allow detection and quantification by HRM analysis. This approach reliably distinguishes between sequence-identical DNA differing only in the methylation of one base. The method also provides accurate discrimination between mixes of methylated and nonmethylated DNAs, allowing discrimination between DNA that is 1% and 0% methylated and also between 97.5% and 100% methylated. Thus, the method provides a new means of adjusting thermal optima for DNA hybridization and PCR-based techniques and to empirically measure the impact of DNA methylation marks on the thermostability of regulatory regions. In the longer term, it could enable the development of new techniques to quantify methylated DNA.

Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants

*Relatively little is known about the timing of genetic and epigenetic forms of somaclonal variation arising from callus growth. We surveyed for both types of change in cocoa (Theobroma cacao) plants regenerated from calli of various ages, and also between tissues from the source trees. *For genetic change, we used 15 single sequence repeat (SSR) markers from four source trees and from 233 regenerated plants. For epigenetic change, we used 386 methylation-sensitive amplified polymorphism (MSAP) markers on leaf and explant (staminode) DNA from two source trees and on leaf DNA from 114 regenerants. *Genetic variation within source trees was limited to one slippage mutation in one leaf. Regenerants were far more variable, with 35% exhibiting at least one mutation. Genetic variation initially accumulated with culture age but subsequently declined. MSAP (epigenetic) profiles diverged between leaf and staminode samples from source trees. Multivariate analysis revealed that leaves from regenerants occupied intermediate eigenspace between leaves and staminodes of source plants but became progressively more similar to source tree leaves with culture age. *Statistical analysis confirmed this rather counterintuitive finding that leaves of 'late regenerants' exhibited significantly less genetic and epigenetic divergence from source leaves than those exposed to short periods of callus growth.

Detection and quantification of in vitro-culture induced chimerism using simple sequence repeat (SSR) analysis in Theobroma cacao (L.)

Mutation rates are often elevated in plants regenerated from in vitro culture, giving rise to so-called 'somaclonal variation'. Detailed characterisation of mutation profiles that arise during culture should improve our understanding of processes influencing mutation and allow the selection of protocols yielding the fewest/least severe changes. Somatic mutations will usually produce genetic chimeras where unchanged alleles are retained by some cells. Such chimeras are difficult to detect but likely to form a significant proportion of any regenerant population. We present a simple protocol that enables the provisional diagnosis of both homogenous and chimeric mutants among large regenerant populations, together with a semi-quantitative means of estimating the proportion of mutant cells. The assay exploits consistent differential amplification of alternate simple sequence repeat alleles at heterozygous loci. Calibration of the relative amplification of alleles from two genotypes-and the synthetic chimeras created from them-revealed a strong linear relationship between 'peak heights' representing alternate alleles following capillary electrophoresis. The assay predicts chimeric composition to a reasonable level of confidence (+/-5%) so long as the infrequent allele exceeds 15% of the template. The system was applied to 233 regenerants of cocoa somatic embryogenesis and identified 72 (31%) putative chimeric mutants for slippage mutation or allele loss across two loci.