Differential regulation of a MYB transcription factor is correlated with transgenerational epigenetic inheritance of trichome density in Mimulus guttatus

Small but strong: herbivory by sap-feeding insect reduces plant progeny growth but enhances direct and indirect anti-herbivore defenses

The transmission of resistance traits to herbivores across subsequent generations is an important strategy employed by plants to enhance their fitness in environments with high herbivore pressure. However, our understanding of the impact of maternal herbivory on direct and indirect induced chemical defenses of progeny, as well as the associated costs, is currently limited to herbivory by leaf-chewing insects. In this study, we investigated the transgenerational effects of a sap-feeding insect, the green peach aphid Myzus persicae, on direct and indirect chemical defenses of bell pepper plants (Capsicum annuum), and whether the effects entail costs to plant growth. Aphid herbivory on parental plants led to a reduced number of seeds per fruit, which exhibited lower germination rates and produced smaller seedlings compared to those from non-infested parental plants. In contrast, the progeny of aphid-infested plants were less preferred as hosts by aphids and less suitable than the progeny of non-infested plants. This enhanced resistance in the progeny of aphid-infested plants coincided with elevated levels of both constitutive and herbivore-induced total phenolic compounds, compared to the progeny of non-infested plants. Furthermore, the progeny of aphid-infested plants emitted herbivore-induced plant volatiles (HIPVs) that were more attractive to the aphid parasitoid Aphidius platensis than those emitted by the progeny of non-infested plants. Our results indicate that herbivory by sap-feeding insect induces transgenerational resistance on progeny bell pepper plants, albeit at the expense of vegetative growth.

Characterisation of the R2R3 Myb subgroup 9 family of transcription factors in tomato

Tomato (Solanum lycopersicum) has many epidermal cell outgrowths including conical petal cells and multiple types of trichomes. These include the anther-specific trichome mesh which holds the anthers connate. The R2R3 Myb Subgroup 9 family of transcription factors is involved in development of epidermal cell outgrowths throughout the angiosperms. No previous study has examined all members of this transcription factor family in a single species. All 7 R2R3 Myb Subgroup 9 genes were isolated from tomato. They were ectopically expressed in tobacco to assess their ability to induce epidermal cell outgrowth. Endogenous expression patterns were examined by semi-quantitative RT-PCR at different stages of floral development relative to the development of anther trichomes. We report variation in the degree of epidermal cell outgrowth produced in transgenic tobacco by each ectopically expressed gene. Based on expression profile and ectopic activity, SlMIXTA-2 is likely involved in the production of leaf trichomes. SlMIXTA-2 is expressed most strongly in the leaves, and not expressed in the floral tissue. SlMYB17-2 is the best candidate for the regulation of the anther trichome mesh. SlMYB17-2 is expressed strongly in the floral tissue and produces a clear phenotype of epidermal cell outgrowths when ectopically expressed in tobacco. Analysis of the phenotypes of transgenic plants ectopically expressing all 7 genes has revealed the different extent to which members of the same transcription factor subfamily can induce cellular outgrowth.

The Potential of Natural Carotenoids-Containing Sericin of the Domestic Silkworm Bombyx mori

Sericin derived from the white cocoon of Bombyx mori has been attracting more attention for its utilization in food, cosmetics, and biomedicine. The potential health benefits of natural carotenoids for humans have also been well-established. Some rare strains of Bombyx mori (B. mori) produce yellow-red cocoons, which endow a potential of natural carotenoid-containing sericin. We hypothesized that natural carotenoid-containing sericin from yellow-red cocoons would exhibit better properties compared with white cocoon sericin. To investigate the physicochemical attributes of natural carotenoid-containing sericin, we bred two silkworm strains from one common ancestor, namely XS7 and XS8, which exhibited different cocoon colors as a result of the inconsistent distribution of lutein and β-carotene. Compared with white cocoon sericin, the interaction between carotenoids and sericin molecules in carotenoid-containing sericin resulted in a unique fluorescence emission at 530, 564 nm. The incorporation of carotenoids enhanced the antibacterial effect, anti-cancer ability, cytocompatibility, and antioxidant of sericin, suggesting potential wide-ranging applications of natural carotenoid-containing sericin as a biomass material. We also found differences in fluorescence characteristics, antimicrobial effects, anti-cancer ability, and antioxidants between XS7 and XS8 sericin. Our work for the first time suggested a better application potential of natural carotenoid-containing sericin as a biomass material than frequently used white cocoon sericin.

A eudicot MIXTA family ancestor likely functioned in both conical cells and trichomes

The MIXTA family of MYB transcription factors modulate the development of diverse epidermal features in land plants. This study investigates the evolutionary history and function of the MIXTA gene family in the early-diverging eudicot model lineage Thalictrum (Ranunculaceae), with R2R3 SBG9-A MYB transcription factors representative of the pre-core eudicot duplication and thus hereby referred to as "paleoMIXTA" (PMX). Cloning and phylogenetic analysis of Thalictrum paleoMIXTA (ThPMX) orthologs across 23 species reveal a genus-wide duplication coincident with a whole-genome duplication. Expression analysis by qPCR confirmed that the highest expression is found in carpels, while newly revealing high expression in leaves and nuanced differences between paralogs in representative polyploid species. The single-copy ortholog from the diploid species T. thalictroides (TthPMX, previously TtMYBML2), which has petaloid sepals with conical-papillate cells and trichomes on leaves, was functionally characterized by virus-induced gene silencing (VIGS), and its role in leaves was also assessed from heterologous overexpression in tobacco. Another ortholog from a species with conical-papillate cells on stamen filaments, TclPMX, was also targeted for silencing. Overexpression assays in tobacco provide further evidence that the paleoMIXTA lineage has the potential for leaf trichome function in a core eudicot. Transcriptome analysis by RNA-Seq on leaves of VIGS-treated plants suggests that TthPMX modulates leaf trichome development and morphogenesis through microtubule-associated mechanisms and that this may be a conserved pathway for eudicots. These experiments provide evidence for a combined role for paleoMIXTA orthologs in (leaf) trichomes and (floral) conical-papillate cells that, together with data from other systems, makes the functional reconstruction of a eudicot ancestor most likely as also having a combined function.

The relative impact of parental and current environment on plant transcriptomes depends on type of stress and genotype

Through developmental plasticity, an individual organism integrates influences from its immediate environment with those due to the environment of its parents. While both effects on phenotypes are well documented, their relative impact has been little studied in natural systems, especially at the level of gene expression. We examined this issue in four genotypes of the annual plant Persicaria maculosa by varying two key resources-light and soil moisture-in both generations. Transcriptomic analyses showed that the relative effects of parent and offspring environment on gene expression (i.e. the number of differentially expressed transcripts, DETs) varied both for the two types of resource stress and among genotypes. For light, immediate environment induced more DETs than parental environment for all genotypes, although the precise proportion of parental versus immediate DETs varied among genotypes. By contrast, the relative effect of soil moisture varied dramatically among genotypes, from 8-fold more DETs due to parental than immediate conditions to 10-fold fewer. These findings provide evidence at the transcriptomic level that the relative impacts of parental and immediate environment on the developing organism may depend on the environmental factor and vary strongly among genotypes, providing potential for the interplay of these developmental influences to evolve.

Herbivory induced methylation changes in the Lombardy poplar: A comparison of results obtained by epiGBS and WGBS

DNA cytosine methylation is an epigenetic mechanism involved in regulation of plant responses to biotic and abiotic stress and its ability to change can vary with the sequence context in which a cytosine appears (CpG, CHG, CHH, where H = Adenine, Thymine, Cytosine). Quantification of DNA methylation in model plant species is frequently addressed by Whole Genome Bisulfite Sequencing (WGBS), which requires a good-quality reference genome. Reduced Representation Bisulfite Sequencing (RRBS) is a cost-effective potential alternative for ecological research with limited genomic resources and large experimental designs. In this study, we provide for the first time a comprehensive comparison between the outputs of RRBS and WGBS to characterize DNA methylation changes in response to a given environmental factor. In particular, we used epiGBS (recently optimized RRBS) and WGBS to assess global and sequence-specific differential methylation after insect and artificial herbivory in clones of Populus nigra cv. 'italica'. We found that, after any of the two herbivory treatments, global methylation percentage increased in CHH, and the shift was detected as statistically significant only by epiGBS. As regards to loci-specific differential methylation induced by herbivory (cytosines in epiGBS and regions in WGBS), both techniques indicated the specificity of the response elicited by insect and artificial herbivory, together with higher frequency of hypo-methylation in CpG and hyper-methylation in CHH. Methylation changes were mainly found in gene bodies and intergenic regions when present at CpG and CHG and in transposable elements and intergenic regions at CHH context. Thus, epiGBS succeeded to characterize global, genome-wide methylation changes in response to herbivory in the Lombardy poplar. Our results support that epiGBS could be particularly useful in large experimental designs aimed to explore epigenetic changes of non-model plant species in response to multiple environmental factors.

Function and molecular mechanism of a poplar placenta limited MIXTA gene in regulating differentiation of plant epidermal cells

The placenta in fruits of most plants either desiccate and shrink as the fruits mature or develop further to form the fleshy tissues. In poplars, placental epidermal cells protrude collectively to produce catkin fibers. In this study, three carpel limited MIXTA genes, PdeMIXTA02, PdeMIXTA03, PdeMIXTA04, were find to specifically expressed in carpel immediately after pollination. Heterologous expression of the three genes in Arabidopsis demonstrated that PdeMIXTA04 significantly promoted trichomes density and could restore trichomes in the trichomeless mutant. By contrast, such functions were not observed with PdeMIXTA02, PdeMIXTA03. In situ hybridization revealed that PdeMIXTA04 was explicitly expressed in poplar placental epidermal cells. We also confirmed trichome-specific expression of the PdeMIXTA04 promoter. Multiple experimental proofs have confirmed the interaction between PdeMIXTA04, PdeMYC and PdeWD40, indicating PdeMIXTA04 functioned through the MYB-bHLH-WD40 ternary complex. Our work provided distinctive understanding of the molecular mechanism triggering differentiation of poplar catkins.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Genetic and epigenetic variation separately contribute to range expansion and local metalliferous habitat adaptation during invasions of Chenopodium ambrosioides into China

BACKGROUND AND AIMS

Invasive plants often colonize wide-ranging geographical areas with various local microenvironments. The specific roles of epigenetic and genetic variation during such expansion are still unclear. Chenopodium ambrosioides is a well-known invasive alien species in China that can thrive in metalliferous habitats. This study aims to comprehensively understand the effects of genetic and epigenetic variation on the successful invasion of C. ambrosioides.

METHODS

We sampled 367 individuals from 21 heavy metal-contaminated and uncontaminated sites with a wide geographical distribution in regions of China. We obtained environmental factors of these sampling sites, including 13 meteorological factors and the contents of four heavy metals in soils. Microsatellite markers were used to investigate the demographic history of C. ambrosioides populations in China. We also analysed the effect of epigenetic variation on metalliferous microhabitat adaptation using methylation-sensitive amplified polymorphism (MSAP) markers. A common garden experiment was conducted to compare heritable phenotypic variations among populations.

KEY RESULTS

Two distinct genetic clusters that diverged thousands of years ago were identified, suggesting that the eastern and south-western C. ambrosioides populations in China may have originated from independent introduction events without recombination. Genetic variation was shown to be a dominant determinant of phenotypic differentiation relative to epigenetic variation, and further affected the geographical distribution pattern of invasive C. ambrosioides. The global DNA unmethylation level was reduced in metalliferous habitats. Dozens of methylated loci were significantly associated with the heavy metal accumulation trait of C. ambrosioides and may contribute to coping with metalliferous microenvironments.

CONCLUSIONS

Our study of C. ambrosioides highlighted the dominant roles of genetic variation in large geographical range expansion and epigenetic variation in local metalliferous habitat adaptation.

Adverse effects of inbreeding on the transgenerational expression of herbivore-induced defense traits in Solanum carolinense

In addition to directly inducing physical and chemical defenses, herbivory experienced by plants in one generation can influence the expression of defensive traits in offspring. Plant defense phenotypes can be compromised by inbreeding, and there is some evidence that such adverse effects can extend to the transgenerational expression of induced resistance. We explored how the inbreeding status of maternal Solanum carolinense plants influenced the transgenerational effects of herbivory on the defensive traits and herbivore resistance of offspring. Manduca sexta caterpillars were used to damage inbred and outbred S. carolinense maternal plants and cross pollinations were performed to produced seeds from herbivore-damaged and undamaged, inbred and outbred maternal plants. Seeds were grown in the greenhouse to assess offspring defense-related traits (i.e., leaf trichomes, internode spines, volatile organic compounds) and resistance to herbivores. We found that feeding by M. sexta caterpillars on maternal plants had a positive influence on trichome and spine production in offspring and that caterpillar development on offspring of herbivore-damaged maternal plants was delayed relative to that on offspring of undamaged plants. Offspring of inbred maternal plants had reduced spine production, compared to those of outbred maternal plants, and caterpillars performed better on the offspring of inbred plants. Both herbivory and inbreeding in the maternal generation altered volatile emissions of offspring. In general, maternal plant inbreeding dampened transgenerational effects of herbivory on offspring defensive traits and herbivore resistance. Taken together, this study demonstrates that inducible defenses in S. carolinense can persist across generations and that inbreeding compromises transgenerational resistance in S. carolinense.

DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium hybrid

BACKGROUND

Flower pigment and shape are determined by the coordinated expression of a set of structural genes during flower development. R2R3-MYB transcription factors are known regulators of structural gene expression. The current study focused on two members of this large family of transcription factors that were predicted to have roles in pigment biosynthesis and organ shape development in orchids.

METHODS

Phylogenetic analysis was used to identify candidate Dendrobium catenatum R2R3-MYB (DcaMYB) sequences associated with pigment and cell shape development. Gene silencing of candidate DhMYBs in Dendrobium hybrid by direct application of dsRNA to developing flowers was followed by observation of gene expression level and flower phenotypes. Silencing of the structural gene chalcone synthase was used as a comparative control.

KEY RESULTS

Ten candidate flower-associated DcaMYBs were identified. Flowers treated with dsRNA of DhMYB22 and DhMYB60 sequences were less pigmented and had relatively low expression of anthocyanin biosynthetic genes (F3'H and DFR), lower total anthocyanin concentration and markedly lower levels of cyanidin-3-glucoside and cyanidin-3-rutinoside. Petals of DhMYB22-treated flowers and sepals of DhMYB60-treated flowers showed the greatest colour difference relative to the same organs in untreated flowers. DhMYB22-treated flowers had relatively narrow and constricted lips, while DhMYB60-treated flowers had narrow and constricted sepals. No significant difference in shape was observed for DhCHS-treated or untreated flowers.

CONCLUSIONS

Our results demonstrate that DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium. This is a first report of MYB regulation of floral organ shape in orchids.

Major niche transitions in Pooideae correlate with variation in photoperiodic flowering and evolution of CCT domain genes

The external cues that trigger timely flowering vary greatly across tropical and temperate plant taxa, the latter relying on predictable seasonal fluctuations in temperature and photoperiod. In the grass family (Poaceae) for example, species of the subfamily Pooideae have become specialists of the northern temperate hemisphere, generating the hypothesis that their progenitor evolved a flowering response to long days from a short-day or day-neutral ancestor. Sampling across the Pooideae, we found support for this hypothesis, and identified several secondary shifts to day-neutral flowering and one to short-day flowering in a tropical highland clade. To explain the proximate mechanisms for the secondary transition back to short-day-regulated flowering, we investigated the expression of CCT domain genes, some of which are known to repress flowering in cereal grasses under specific photoperiods. We found a shift in CONSTANS 1 and CONSTANS 9 expression that coincides with the derived short-day photoperiodism of our exemplar species Nassella pubiflora. This sets up the testable hypothesis that trans- or cis-regulatory elements of these CCT domain genes were the targets of selection for major niche shifts in Pooideae grasses.

Molecular Mechanisms of Plant Trichome Development

Plant trichomes, protrusions formed from specialized aboveground epidermal cells, provide protection against various biotic and abiotic stresses. Trichomes can be unicellular, bicellular or multicellular, with multiple branches or no branches at all. Unicellular trichomes are generally not secretory, whereas multicellular trichomes include both secretory and non-secretory hairs. The secretory trichomes release secondary metabolites such as artemisinin, which is valuable as an antimalarial agent. Cotton trichomes, also known as cotton fibers, are an important natural product for the textile industry. In recent years, much progress has been made in unraveling the molecular mechanisms of trichome formation in Arabidopsis thaliana, Gossypium hirsutum, Oryza sativa, Cucumis sativus, Solanum lycopersicum, Nicotiana tabacum, and Artemisia annua. Here, we review current knowledge of the molecular mechanisms underlying fate determination and initiation, elongation, and maturation of unicellular, bicellular and multicellular trichomes in several representative plants. We emphasize the regulatory roles of plant hormones, transcription factors, the cell cycle and epigenetic modifications in different stages of trichome development. Finally, we identify the obstacles and key points for future research on plant trichome development, and speculated the development relationship between the salt glands of halophytes and the trichomes of non-halophytes, which provides a reference for future studying the development of plant epidermal cells.

Transgenerational changes in pod maturation phenology and seed traits of Glycine soja infested by the bean bug Riptortus pedestris

Land plants have diverse defenses against herbivores. In some cases, plant response to insect herbivory may be chronological and even transgenerational. Feeding by various stink bugs, such as the bean bug Riptortus pedestris (Hemiptera: Alydidae), induce physiological changes in soybean, called as green stem syndrome, which are characterized by delayed senescence in stems, leaves, and pods. To investigate the plant response to the bean bug feeding in the infested generation and its offspring, we studied the effects of R. pedestris infestation on Glycine soja, the ancestral wild species of soybean. Field surveys revealed that the occurrence of the autumn R. pedestris generation coincided with G. soja pod maturation in both lowland and mountainous sites. Following infestation by R. pedestris, pod maturation was significantly delayed in G. soja. When G. soja seeds obtained from infested and non-infested plants were cultivated, the progeny of infested plants exhibited much earlier pod maturation and larger-sized seed production than that of control plants, indicating that R. pedestris feeding induced transgenerational changes. Because earlier seed maturity results in asynchrony with occurrence of R. pedestris, the transgenerational changes in plant phenology are considered to be an adaptive transgenerational and chronological defense for the plant against feeding by the stink bug.

Lifetime genealogical divergence within plants leads to epigenetic mosaicism in the shrub Lavandula latifolia (Lamiaceae)

Epigenetic mosaicism is a possible source of within-plant phenotypic heterogeneity, yet its frequency and developmental origin remain unexplored. This study examines whether extant epigenetic heterogeneity within Lavandula latifolia (Lamiaceae) shrubs reflects recent epigenetic modifications experienced independently by different plant parts or, alternatively, it is the cumulative outcome of a steady lifetime process. Leaf samples from different architectural modules (branch tips) were collected from three L. latifolia plants and characterized epigenetically by global DNA cytosine methylation and methylation state of methylation-sensitive amplified fragment-length polymorphism (MS-AFLP) markers. Epigenetic characteristics of modules were then assembled with information on the branching history of plants. Methods borrowed from phylogenetic research were used to assess genealogical signal of extant epigenetic variation and reconstruct within-plant genealogical trajectory of epigenetic traits. Plants were epigenetically heterogeneous, as shown by differences among modules in global DNA methylation and variation in the methylation states of 6 to 8% of MS-AFLP markers. All epigenetic features exhibited significant genealogical signal within plants. Events of epigenetic divergence occurred throughout the lifespan of individuals and were subsequently propagated by branch divisions. Internal epigenetic diversification of L. latifolia individuals took place steadily during their development, a process which eventually led to persistent epigenetic mosaicism.

An HD-ZIP-MYB complex regulates glandular secretory trichome initiation in Artemisia annua

Plant glandular secretory trichomes (GSTs) produce various specialized metabolites. Increasing GST density represents a strategy to enhance the yield of these chemicals; however, the gene regulatory network that controls GST initiation remains unclear. In a previous study of Artemisia annua L., we found that a HD-ZIP IV transcription factor, AaHD1, promotes GST initiation by directly regulating AaGSW2. Here, we identified two AaHD1-interacting transcription factors, namely AaMIXTA-like 2 (AaMYB16) and AaMYB5. Through the generation and characterization of transgenic plants, we found that AaMYB16 is a positive regulator of GST initiation, whereas AaMYB5 has the opposite effect. Notably, neither of them regulates GST formation independently. Rather, they act competitively, by interacting and modulating AaHD1 promoter binding activity. Additionally, the phytohormone jasmonic acid (JA) was shown to be associated with the AaHD1-AaMYB16/AaMYB5 regulatory network through transcriptional regulation via a JASMONATE-ZIM DOMAIN (JAZ) protein repressor. These results bring new insights into the mechanism of GST initiation through regulatory complexes, which appear to have similar functions in a range of vascular plant taxa.

Analysis and review of trichomes in plants

BACKGROUND

Trichomes play a key role in the development of plants and exist in a wide variety of species.

RESULTS

In this paper, it was reviewed that the structure and morphology characteristics of trichomes, alongside the biological functions and classical regulatory mechanisms of trichome development in plants. The environment factors, hormones, transcription factor, non-coding RNA, etc., play important roles in regulating the initialization, branching, growth, and development of trichomes. In addition, it was further investigated the atypical regulation mechanism in a non-model plant, found that regulating the growth and development of tea (Camellia sinensis) trichome is mainly affected by hormones and the novel regulation factors.

CONCLUSIONS

This review further displayed the complex and differential regulatory networks in trichome initiation and development, provided a reference for basic and applied research on trichomes in plants.

Changes in DNA Methylation in Arabidopsis thaliana Plants Exposed Over Multiple Generations to Gamma Radiation

Previous studies have found indications that exposure to ionising radiation (IR) results in DNA methylation changes in plants. However, this phenomenon is yet to be studied across multiple generations. Furthermore, the exact role of these changes in the IR-induced plant response is still far from understood. Here, we study the effect of gamma radiation on DNA methylation and its effect across generations in young Arabidopsis plants. A multigenerational set-up was used in which three generations (Parent, generation 1, and generation 2) of 7-day old Arabidopsis thaliana plants were exposed to either of the different radiation treatments (30, 60, 110, or 430 mGy/h) or to natural background radiation (control condition) for 14 days. The parental generation consisted of previously non-exposed plants, whereas generation 1 and generation 2 plants had already received a similar irradiation in the previous one or two generations, respectively. Directly after exposure the entire methylomes were analysed with UPLC-MS/MS to measure whole genome methylation levels. Whole genome bisulfite sequencing was used to identify differentially methylated regions (DMRs), including their methylation context in the three generations and this for three different radiation conditions (control, 30 mGy/h, and 110 mGy/h). Both intra- and intergenerational comparisons of the genes and transposable elements associated with the DMRs were made. Taking the methylation context into account, the highest number of changes were found for cytosines followed directly by guanine (CG methylation), whereas only limited changes in CHG methylation occurred and no changes in CHH methylation were observed. A clear increase in IR-induced DMRs was seen over the three generations that were exposed to the lowest dose rate, where generation 2 had a markedly higher number of DMRs than the previous two generations (Parent and generation 1). Counterintuitively, we did not see significant differences in the plants exposed to the highest dose rate. A large number of DMRs associated with transposable elements were found, the majority of them being hypermethylated, likely leading to more genetic stability. Next to that, a significant number of DMRs were associated with genes (either in their promoter-associated region or gene body). A functional analysis of these genes showed an enrichment for genes related to development as well as various stress responses, including DNA repair, RNA splicing, and (a)biotic stress responses. These observations indicate a role of DNA methylation in the regulation of these genes in response to IR exposure and shows a possible role for epigenetics in plant adaptation to IR over multiple generations.

Small DNA Methylation, Big Player in Plant Abiotic Stress Responses and Memory

DNA methylation is a conserved epigenetic mark that plays important roles in maintaining genome stability and regulating gene expression. As sessile organisms, plants have evolved sophisticated regulatory systems to endure or respond to diverse adverse abiotic environmental challenges, i.e., abiotic stresses, such as extreme temperatures (cold and heat), drought and salinity. Plant stress responses are often accompanied by changes in chromatin modifications at diverse responsive loci, such as 5-methylcytosine (5mC) and N 6-methyladenine (6mA) DNA methylation. Some abiotic stress responses are memorized for several hours or days through mitotic cell divisions and quickly reset to baseline levels after normal conditions are restored, which is referred to as somatic memory. In some cases, stress-induced chromatin marks are meiotically heritable and can impart the memory of stress exposure from parent plants to at least the next stress-free offspring generation through the mechanisms of transgenerational epigenetic inheritance, which may offer the descendants the potential to be adaptive for better fitness. In this review, we briefly summarize recent achievements regarding the establishment, maintenance and reset of DNA methylation, and highlight the diverse roles of DNA methylation in plant responses to abiotic stresses. Further, we discuss the potential role of DNA methylation in abiotic stress-induced somatic memory and transgenerational inheritance. Future research directions are proposed to develop stress-tolerant engineered crops to reduce the negative effects of abiotic stresses.

The MIXTA-LIKE transcription factor CsMYB6 regulates fruit spine and tubercule formation in cucumber

Cucumber fruit wart composed of tubercule and spine (trichome on fruit) is not only an important fruit quality trait in cucumber production, but also a well-studied model for plant cell-fate determination. The development of spine is closely related to the initiation and formation of tubercule. The spine differentiation regulator CsGL1 has been proved to be epistatic to the tubercule initiation factor CsTu, which is the only connection to be identified between spine and tubercule formations. Our previous studies found that the MIXTA-LIKE transcription factor CsMYB6 can suppress fruit spine initiation, which is independent of CsGL1. How the formation of spine and tubercule is regulated at the molecular level by CsMYB6 remains poorly understood. In this study, we characterized cucumber 35S:CsMYB6 transgenic plants, which displayed an obvious reduction in the number and size of fruit spines and tubecules. Molecular analyses showed that CsMYB6 directly interacted with the key spine formation factor CsTTG1 in regulating the formation of fruit spine, and CsTu in regulating the initiation of fruit tubercule, respectively. Based on these evidences, a novel regulatory network is proposed by which CsMYB6/CsTTG1 and CsMYB6/CsTu complexes play an important role in regulating epidermal development, including spine formation and tubercule initiation in cucumber.

Mimulus sRNAs Are Wound Responsive and Associated with Transgenerationally Plastic Genes but Rarely Both

Organisms alter development in response to environmental cues. Recent studies demonstrate that they can transmit this plasticity to progeny. While the phenotypic and transcriptomic evidence for this “transgenerational plasticity” has accumulated, genetic and developmental mechanisms remain unclear. Plant defenses, gene expression and DNA methylation are modified as an outcome of parental wounding in Mimulus guttatus. Here, we sequenced M. guttatus small RNAs (sRNA) to test their possible role in mediating transgenerational plasticity. We sequenced sRNA populations of leaf-wounded and control plants at 1 h and 72 h after damage and from progeny of wounded and control parents. This allowed us to test three components of an a priori model of sRNA mediated transgenerational plasticity—(1) A subset of sRNAs will be differentially expressed in response to wounding, (2) these will be associated with previously identified differentially expressed genes and differentially methylated regions and (3) changes in sRNA abundance in wounded plants will be predictive of sRNA abundance, DNA methylation, and/or gene expression shifts in the following generation. Supporting (1) and (2), we found significantly different sRNA abundances in wounded leaves; the majority were associated with tRNA fragments (tRFs) rather than small-interfering RNAs (siRNA). However, siRNAs responding to leaf wounding point to Jasmonic Acid mediated responses in this system. We found that different sRNA classes were associated with regions of the genome previously found to be differentially expressed or methylated in progeny of wounded plants. Evidence for (3) was mixed. We found that non-dicer sRNAs with increased abundance in response to wounding tended to be nearby genes with decreased expression in the next generation. Counter to expectations, we did not find that siRNA responses to wounding were associated with gene expression or methylation changes in the next generation and within plant and transgenerational sRNA plasticity were negatively correlated.

Understanding 'Non-genetic' Inheritance: Insights from Molecular-Evolutionary Crosstalk

Understanding the evolutionary and ecological roles of 'non-genetic' inheritance (NGI) is daunting due to the complexity and diversity of epigenetic mechanisms. We draw on insights from molecular and evolutionary biology perspectives to identify three general features of 'non-genetic' inheritance systems: (i) they are functionally interdependent with, rather than separate from, DNA sequence; (ii) precise mechanisms vary phylogenetically and operationally; and (iii) epigenetic elements are probabilistic, interactive regulatory factors and not deterministic 'epialleles' with defined genomic locations and effects. We discuss each of these features and offer recommendations for future empirical and theoretical research that implements a unifying inherited gene regulation (IGR) approach to studies of 'non-genetic' inheritance.

The transcription factor MML4_D12 regulates fiber development through interplay with the WD40-repeat protein WDR in cotton

In planta, a vital regulatory complex, MYB-basic helix-loop-helix (bHLH)-WD40 (MBW), is involved in trichome development and synthesis of anthocyanin and proanthocyanin in Arabidopsis. Usually, WD40 proteins provide a scaffold for protein-protein interaction between MYB and bHLH proteins. Members of subgroup 9 of the R2R3 MYB transcription factors, which includes MYBMIXTA-Like (MML) genes important for plant cell differentiation, are unable to interact with bHLH. In this study, we report that a cotton (Gossypium hirsutum) seed trichome or lint fiber-related GhMML factor, GhMML4_D12, interacts with a diverged WD40 protein (GhWDR) in a process similar to but different from that of the MBW ternary complex involved in Arabidopsis trichome development. Amino acids 250-267 of GhMML4_D12 and the first and third WD40 repeat domains of GhWDR determine their interaction. GhWDR could rescue Arabidopsis ttg1 to its wild type, confirming its orthologous function in trichome development. Our findings shed more light towards understanding the key role of the MML and WD40 families in plants and in the improvement of cotton fiber production.

A Scanning Electron Micrograph-based Resource for Identification of Loci Involved in Epidermal Development in Tomato: Elucidation of a New Function for the Mixta-like Transcription Factor in Leaves

The aerial epidermis of plants plays a major role in environmental interactions, yet the development of the cellular components of the aerial epidermis-trichomes, stomata, and pavement cells-is still not fully understood. We have performed a detailed screen of the leaf epidermis in two generations of the well-established Solanum lycopersicum cv M82 × Solanum pennellii ac. LA716 introgression line (IL) population using a combination of scanning electron microscopy (SEM) techniques. Quantification of trichome and stomatal densities in the ILs revealed four genomic regions with a consistently low trichome density. This study also found ILs with abnormal proportions of different trichome types and aberrant trichome morphologies. This work has led to the identification of new, unexplored genomic regions with roles in trichome formation in tomato. This study investigated one interval in IL2-6 in more detail and identified a new function for the transcription factor SlMixta-like in determining trichome patterning in leaves. This illustrates how these SEM images, publicly available to the research community, provide an important dataset for further studies on epidermal development in tomato and other species of the Solanaceae family.

Co-expression clustering across flower development identifies modules for diverse floral forms in Achimenes (Gesneriaceae)

BACKGROUND

Genetic pathways involved with flower color and shape are thought to play an important role in the development of flowers associated with different pollination syndromes, such as those associated with bee, butterfly, or hummingbird pollination. Because pollination syndromes are complex traits that are orchestrated by multiple genes and pathways, the gene regulatory networks have not been explored. Gene co-expression networks provide a systems level approach to identify important contributors to floral diversification.

METHODS

RNA-sequencing was used to assay gene expression across two stages of flower development (an early bud and an intermediate stage) in 10 species of Achimenes (Gesneriaceae). Two stage-specific co-expression networks were created from 9,503 orthologs and analyzed to identify module hubs and the network periphery. Module association with bee, butterfly, and hummingbird pollination syndromes was tested using phylogenetic mixed models. The relationship between network connectivity and evolutionary rates (d N/d S) was tested using linear models.

RESULTS

Networks contained 65 and 62 modules that were largely preserved between developmental stages and contained few stage-specific modules. Over a third of the modules in both networks were associated with flower color, shape, and pollination syndrome. Within these modules, several hub nodes were identified that related to the production of anthocyanin and carotenoid pigments and the development of flower shape. Evolutionary rates were decreased in highly connected genes and elevated in peripheral genes.

DISCUSSION

This study aids in the understanding of the genetic architecture and network properties underlying the development of floral form and provides valuable candidate modules and genes for future studies.

Transgenerational effects of parental light environment on progeny competitive performance and lifetime fitness

Plant and animal parents may respond to environmental conditions such as resource stress by altering traits of their offspring via heritable non-genetic effects. While such transgenerational plasticity can result in progeny phenotypes that are functionally pre-adapted to the inducing environment, it is unclear whether such parental effects measurably enhance the adult competitive success and lifetime reproductive output of progeny, and whether they may also adversely affect fitness if offspring encounter contrasting conditions. In glasshouse experiments with inbred genotypes of the annual plant Polygonum persicaria, we tested the effects of parental shade versus sun on (a) competitive performance of progeny in shade, and (b) lifetime reproductive fitness of progeny in three contrasting treatments. Shaded parents produced offspring with increased fitness in shade despite competition, as well as greater competitive impact on plant neighbours. Inherited effects of parental light conditions also significantly altered lifetime fitness: parental shade increased reproductive output for progeny in neighbour and understorey shade, but decreased fitness for progeny in sunny, dry conditions. Along with these substantial adaptive and maladaptive transgenerational effects, results show complex interactions between genotypes, parent environment and progeny conditions that underscore the role of environmental variability and change in shaping future adaptive potential. This article is part of the theme issue 'The role of plasticity in phenotypic adaptation to rapid environmental change'.

Fixation of genetic variation and optimization of gene expression: The speed of evolution in isolated lizard populations undergoing Reverse Island Syndrome

The ecological theory of island biogeography suggests that mainland populations should be more genetically divergent from those on large and distant islands rather than from those on small and close islets. Some island populations do not evolve in a linear way, but the process of divergence occurs more rapidly because they undergo a series of phenotypic changes, jointly known as the Island Syndrome. A special case is Reversed Island Syndrome (RIS), in which populations show drastic phenotypic changes both in body shape, skin colouration, age of sexual maturity, aggressiveness, and food intake rates. The populations showing the RIS were observed on islets nearby mainland and recently raised, and for this they are useful models to study the occurrence of rapid evolutionary change. We investigated the timing and mode of evolution of lizard populations adapted through selection on small islets. For our analyses, we used an ad hoc model system of three populations: wild-type lizards from the mainland and insular lizards from a big island (Capri, Italy), both Podarcis siculus siculus not affected by the syndrome, and a lizard population from islet (Scopolo) undergoing the RIS (called P. s. coerulea because of their melanism). The split time of the big (Capri) and small (Scopolo) islands was determined using geological events, like sea-level rises. To infer molecular evolution, we compared five complete mitochondrial genomes for each population to reconstruct the phylogeography and estimate the divergence time between island and mainland lizards. We found a lower mitochondrial mutation rate in Scopolo lizards despite the phenotypic changes achieved in approximately 8,000 years. Furthermore, transcriptome analyses showed significant differential gene expression between islet and mainland lizard populations, suggesting the key role of plasticity in these unpredictable environments.

Genetic and epigenetic regulation of phenotypic variation in invasive plants – linking research trends towards a unified framework

Phenotypic variation in the introduced range of an invasive species can be modified by genetic variation, environmental conditions and their interaction, as well as stochastic events like genetic drift. Recent studies found that epigenetic modifications may also contribute to phenotypic variation being independent of genetic changes. Despite gaining profound ecological insights from empirical studies, understanding the relative contributions of these molecular mechanisms behind phenotypic variation has received little attention for invasive plant species in particular.

This review therefore aimed at summarizing and synthesizing information on the genetic and epigenetic basis of phenotypic variation of alien invasive plants in the introduced range and their evolutionary consequences. Transgenerational inheritance of epigenetic modifications was highlighted focusing on its influence on microevolution of the invasive plant species. We presented a comprehensive account of epigenetic regulation of phenotypic variation and its role in plant invasion in the presence of reduced standing genetic variation, inbreeding depression and associated genomic events which have often been observed during introduction and range expansion of an invasive alien species. Finally, taking clues from the studies conducted so far, we proposed a unified framework of future experimental approaches to understand ecological and evolutionary aspects of phenotypic variation. This holistic approach, being aligned to the invasion process in particular (introduction-establishment-spread), was intended to understand the molecular mechanisms of phenotypic variation of an invasive species in its introduced range and to disentangle the effects of standing genetic variation and epigenetic regulation of phenotypic variation.

The role of plant epigenetics in biotic interactions

Contents Summary 731 I. Biotic interactions in the context of genetic, epigenetic and environmental diversity 731 II. Biotic interactions affect epigenetic configuration 732 III. Plant epigenetic configuration influences biotic interactions 733 IV. Epigenetic memory in the context of biotic interactions 734 V. Conclusions and future research 735 Acknowledgements 735 Author contributions 735 References 735 SUMMARY: Plants are hubs of a wide range of biotic interactions with mutualist and antagonist animals, microbes and neighboring plants. Because the quality and intensity of those relationships can change over time, a fast and reversible response to stress is required. Here, we review recent studies on the role of epigenetic factors such as DNA methylation and histone modifications in modulating plant biotic interactions, and discuss the state of knowledge regarding their potential role in memory and priming. Moreover, we provide an overview of strategies to investigate the contribution of epigenetics to environmentally induced phenotypic changes in an ecological context, highlighting possible transitions from whole-genome high-resolution analyses in plant model organisms to informative reduced representation analyses in genomically less accessible species.

Parental experience modifies the Mimulus methylome

BACKGROUND

Transgenerational plasticity occurs when the environmental experience of an organism modifies the growth and development of its progeny. Leaf damage in Mimulus guttatus exhibits transgenerational plasticity mediated through differential expression of hundreds of genes. The epigenetic mechanisms that facilitate this response have yet to be described.

RESULTS

We performed whole genome bisulfite sequencing in the progeny of genetically identical damaged and control plants and developed a pipeline to compare differences in the mean and variance of methylation between treatment groups. We find that parental damage increases the variability of CG and CHG methylation among progeny, but does not alter the overall mean methylation. Instead it has positive effects in some regions and negative in others. We find 3,396 CHH, 203 CG, and 54 CHG Differentially Methylated Regions (DMRs) ranging from tens to thousands of base pairs scattered across the genome. CHG and CHH DMRs tended to overlap with transposable elements. CG DMRs tended to overlap with gene coding regions, many of which were previously found to be differentially expressed.

CONCLUSIONS

Genome-wide increases in methylome variation suggest that parental conditions can increase epigenetic diversity in response to stress. Additionally, the potential association between CG DMRs and differentially expressed genes supports the hypothesis that differential methylation is a mechanistic component of transgenerational plasticity in M. guttatus.

A CsMYB6-CsTRY module regulates fruit trichome initiation in cucumber

Fruit epidermal features such as the number and size of trichomes or spines are important fruit quality traits in cucumber production. Little is known about the molecular mechanisms underlying fruit spine formation in cucumber. Here, we report functional characterization of the cucumber CsMYB6 gene, which encodes a MIXTA-like MYB transcription factor that plays an important role in regulating fruit trichome development. Spatial-temporal expression analyses revealed high-level expression of CsMYB6 in the epidermis of cucumber ovaries during fruit spine initiation, which was similar to the expression of CsTRY, a homolog of the Arabidopsis TRY gene that also plays a key role in trichome development. Overexpression of CsMYB6 and CsTRY in cucumber and Arabidopsis revealed that CsMYB6 and CsTRY act as negative regulators of trichome initiation in both species, and that CsMYB6 acted upstream of CsTRY in this process. CsMYB6 was found to bind to the three MYB binding sites inside the promoter region of CsTRY, and protein-protein interaction assays suggested that CsTRY also directly interacted with CsMYB6 protein. The results also revealed conserved and divergent roles of CsMYB6 and its Arabidopsis homolog AtMYB106 in trichome development. Collectively, our results reveal a novel mechanism in which the CsMYB6-CsTRY complex negatively regulates fruit trichome formation in cucumber.

Context-Dependent Developmental Effects of Parental Shade Versus Sun Are Mediated by DNA Methylation

Parental environment influences progeny development in numerous plant and animal systems. Such inherited environmental effects may alter offspring phenotypes in a consistent way, for instance when resource-deprived parents produce low quality offspring due to reduced maternal provisioning. However, because development of individual organisms is guided by both inherited and immediate environmental cues, parental conditions may have different effects depending on progeny environment. Such context-dependent transgenerational plasticity suggests a mechanism of environmental inheritance that can precisely interact with immediate response pathways, such as epigenetic modification. We show that parental light environment (shade versus sun) resulted in context-dependent effects on seedling development in a common annual plant, and that these effects were mediated by DNA methylation. We grew replicate parents of five highly inbred Polygonum persicaria genotypes in glasshouse shade versus sun and, in a fully factorial design, measured ecologically important traits of their isogenic seedling offspring in both environments. Compared to the offspring of sun-grown parents, the offspring of shade-grown parents produced leaves with greater mean and specific leaf area, and had higher total leaf area and biomass. These shade-adaptive effects of parental shade were pronounced and highly significant for seedlings growing in shade, but slight and generally non-significant for seedlings growing in sun. Based on both regression and covariate analysis, inherited effects of parental shade were not mediated by changes to seed provisioning. To test for a role of DNA methylation, we exposed replicate offspring of isogenic shaded and fully insolated parents to either the demethylating agent zebularine or to control conditions during germination, then raised them in simulated growth chamber shade. Partial demethylation of progeny DNA had no phenotypic effect on offspring of shaded parents, but caused offspring of sun-grown parents to develop as if their parents had been shaded, with larger leaves and greater total canopy area and biomass. These results contribute to the increasing body of evidence that DNA methylation can mediate transgenerational environmental effects, and show that such effects may contribute to nuanced developmental interactions between parental and immediate environments.

DNA Methylation and the Evolution of Developmental Complexity in Plants

All land plants so far examined use DNA methylation to silence transposons (TEs). DNA methylation therefore appears to have been co-opted in evolution from an original function in TE management to a developmental function (gene regulation) in both phenotypic plasticity and in normal development. The significance of DNA methylation to the evolution of developmental complexity in plants lies in its role in the management of developmental pathways. As such it is more important in fine tuning the presence, absence, and placement of organs rather than having a central role in the evolution of new organs. Nevertheless, its importance should not be underestimated as it contributes considerably to the range of phenotypic expression and complexity available to plants: the subject of the emerging field of epi-evodevo. Furthermore, changes in DNA methylation can function as a "soft" mutation that may be important in the early stages of major evolutionary novelty.

Functional Characterization of the Versatile MYB Gene Family Uncovered Their Important Roles in Plant Development and Responses to Drought and Waterlogging in Sesame

The MYB gene family constitutes one of the largest transcription factors (TFs) modulating various biological processes in plants. Although genome-wide analysis of this gene family has been carried out in some species, only three MYB members have been functionally characterized heretofore in sesame (Sesamum indicum L.). Here, we identified a relatively high number (287) of sesame MYB genes (SIMYBs) with an uncommon overrepresentation of the 1R-subfamily. A total of 95% of SIMYBs was mapped unevenly onto the 16 linkage groups of the sesame genome with 55 SIMYBs tandemly duplicated. In addition, molecular characterization, gene structure, and evolutionary relationships of SIMYBs were established. Based on the close relationship between sesame and Arabidopsis thaliana, we uncovered that the functions of SIMYBs are highly diverse. A total of 65% of SIMYBs were commonly detected in five tissues, suggesting that they represent key TFs modulating sesame growth and development. Moreover, we found that SIMYBs regulate sesame responses to drought and waterlogging, which highlights the potential of SIMYBs towards improving stress tolerance in sesame. This work presents a comprehensive picture of the MYB gene family in sesame and paves the way for further functional validation of the members of this versatile gene family.

Parallel loss of introns in the ABCB1 gene in angiosperms

BACKGROUND

The presence of non-coding introns is a characteristic feature of most eukaryotic genes. While the size of the introns, number of introns per gene and the number of intron-containing genes can vary greatly between sequenced eukaryotic genomes, the structure of a gene with reference to intron presence and positions is typically conserved in closely related species. Unexpectedly, the ABCB1 (ATP-Binding Cassette Subfamily B Member 1) gene which encodes a P-glycoprotein and underlies dwarfing traits in maize (br2), sorghum (dw3) and pearl millet (d2) displayed considerable variation in intron composition.

RESULTS

An analysis of the ABCB1 gene structure in 80 angiosperms revealed that the number of introns ranged from one to nine. All introns in ABCB1 underwent either a one-time loss (single loss in one lineage/species) or multiple independent losses (parallel loss in two or more lineages/species) with the majority of losses occurring within the grass family. In contrast, the structure of the closest homolog to ABCB1, ABCB19, remained constant in the majority of angiosperms analyzed. Using known phylogenetic relationships within the grasses, we determined the ancestral branch-points where the losses occurred. Intron 7, the longest intron, was lost in only a single species, Mimulus guttatus, following duplication of ABCB1. Semiquantitative PCR showed that the M. guttatus ABCB1 gene copy without intron 7 had significantly lower transcript levels than the gene copy with intron 7. We further demonstrated that intron 7 carried two motifs that were highly conserved across the monocot-dicot divide.

CONCLUSIONS

The ABCB1 gene structure is highly dynamic, while the structure of ABCB19 remained largely conserved through evolution. Precise removal of introns, preferential removal of smaller introns and presence of at least 2 bp of microhomology flanking most introns indicated that intron loss may have predominantly occurred through non-homologous end-joining (NHEJ) repair of double strand breaks. Lack of microhomology in the exon upstream of lost phase I introns was likely due to release of the selective constraint on the penultimate base (3rd base in codon) of the terminal codon by the splicing machinery. In addition to size, the presence of regulatory motifs will make introns recalcitrant to loss.

Developmental plasticity: re-conceiving the genotype

In recent decades, the phenotype of an organism (i.e. its traits and behaviour) has been studied as the outcome of a developmental 'programme' coded in its genotype. This deterministic view is implicit in the Modern Synthesis approach to adaptive evolution as a sorting process among genetic variants. Studies of developmental pathways have revealed that genotypes are in fact differently expressed depending on environmental conditions. Accordingly, the genotype can be understood as a repertoire of potential developmental outcomes or norm of reaction. Reconceiving the genotype as an environmental response repertoire rather than a fixed developmental programme leads to three critical evolutionary insights. First, plastic responses to specific conditions often comprise functionally appropriate trait adjustments, resulting in an individual-level, developmental mode of adaptive variation. Second, because genotypes are differently expressed depending on the environment, the genetic diversity available to natural selection is itself environmentally contingent. Finally, environmental influences on development can extend across multiple generations via cytoplasmic and epigenetic factors transmitted to progeny individuals, altering their responses to their own, immediate environmental conditions and, in some cases, leading to inherited but non-genetic adaptations. Together, these insights suggest a more nuanced understanding of the genotype and its evolutionary role, as well as a shift in research focus to investigating the complex developmental interactions among genotypes, environments and previous environments.

Comparative epigenetic and genetic spatial structure of the perennial herb Helleborus foetidus: Isolation by environment, isolation by distance, and functional trait divergence

PREMISE OF THE STUDY

Epigenetic variation can play a role in local adaptation; thus, there should be associations among epigenetic variation, environmental variation, and functional trait variation across populations. This study examines these relationships in the perennial herb Helleborus foetidus (Ranunculaceae).

METHODS

Plants from 10 subpopulations were characterized genetically (AFLP, SSR markers), epigenetically (MSAP markers), and phenotypically (20 functional traits). Habitats were characterized using six environmental variables. Isolation-by-distance (IBD) and isolation-by-environment (IBE) patterns of genetic and epigenetic divergence were assessed, as was the comparative explanatory value of geographical and environmental distance as predictors of epigenetic, genetic, and functional differentiation.

KEY RESULTS

Subpopulations were differentiated genetically, epigenetically, and phenotypically. Genetic differentiation was best explained by geographical distance, while epigenetic differentiation was best explained by environmental distance. Divergence in functional traits was correlated with environmental and epigenetic distances, but not with geographical and genetic distances.

CONCLUSIONS

Results are compatible with the hypothesis that epigenetic IBE and functional divergence reflected responses to environmental variation. Spatial analyses simultaneously considering epigenetic, genetic, phenotypic and environmental information provide a useful tool to evaluate the role of environmental features as drivers of natural epigenetic variation between populations.

Evolution of the miR5200-FLOWERING LOCUS T flowering time regulon in the temperate grass subfamily Pooideae

Flowering time is a carefully regulated trait controlled primarily through the action of the central genetic regulator, FLOWERING LOCUS T (FT). Recently it was demonstrated that a microRNA, miR5200, targets the end of the second exon of FT under short-day photoperiods in the grass subfamily Pooideae, thus preventing FT transcripts from reaching threshold levels under non-inductive conditions. Pooideae are an interesting group in that they rapidly diversified from the tropics into the northern temperate region during a major global cooling event spanning the Eocene-Oligocene transition. We hypothesize that miR5200 photoperiod-sensitive regulation of Pooideae flowering time networks assisted their transition into northern seasonal environments. Here, we test predictions derived from this hypothesis that miR5200, originally found in bread wheat and later identified in Brachypodium distachyon, (1) was present in the genome of the Pooideae common ancestor, (2) is transcriptionally regulated by photoperiod, and (3) is negatively correlated with FT transcript abundance, indicative of miR5200 regulating FT. Our results demonstrate that miR5200 did evolve at or around the base of Pooideae, but only acquired photoperiod-regulated transcription within the Brachypodium lineage. Based on expression profiles and previous data, we posit that the progenitor of miR5200 was co-regulated with FT by an unknown mechanism.

Repeated and diverse losses of corolla bilateral symmetry in the Lamiaceae

BACKGROUND AND AIMS

Independent evolution of derived complex characters provides a unique opportunity to assess whether and how similar genetic changes correlate with morphological convergence. Bilaterally symmetrical corollas have evolved multiple times independently from radially symmetrical ancestors and likely represent adaptations to attract specific pollinators. On the other hand, losses of bilateral corolla symmetry have occurred sporadically in various groups, due to either modification of bilaterally symmetrical corollas in late development or early establishment of radial symmetry.

METHODS

This study integrated phylogenetic, scanning electron microscopy (SEM)-based morphological, and gene expression approaches to assess the possible mechanisms underlying independent evolutionary losses of corolla bilateral symmetry.

KEY RESULTS

This work compared three species of Lamiaceae having radially symmetrical mature corollas with a representative sister taxon having bilaterally symmetrical corollas and found that each reaches radial symmetry in a different way. Higher core Lamiales share a common duplication in the CYCLOIDEA (CYC ) 2 gene lineage and show conserved and asymmetrical expression of CYC2 clade and RAD genes along the adaxial-abaxial floral axis in species having bilateral corolla symmetry. In Lycopus americanus , the development and expression pattern of La-CYC2A and La-CYC2B are similar to those of their bilaterally symmetrical relatives, whereas the loss of La-RAD expression correlates with a late switch to radial corolla symmetry. In Mentha longifolia , late radial symmetry may be explained by the loss of Ml-CYC2A , and by altered expression of two Ml-CYC2B and Ml-RAD genes . Finally, expanded expression of Cc-CYC2A and Cc-RAD strongly correlates with the early development of radially symmetrical corollas in Callicarpa cathayana .

CONCLUSIONS

Repeated losses of mature corolla bilateral symmetry in Lamiaceae are not uncommon, and may be achieved by distinct mechanisms and various changes to symmetry genes, including the loss of a CYC2 clade gene from the genome, and/or contraction, expansion or alteration of CYC2 clade and RAD -like gene expression.

A Whole-Transcriptome Approach to Evaluating Reference Genes for Quantitative Gene Expression Studies: A Case Study in Mimulus

While quantitative PCR (qPCR) is widely recognized as being among the most accurate methods for quantifying gene expression, it is highly dependent on the use of reliable, stably expressed reference genes. With the increased availability of high-throughput methods for measuring gene expression, whole-transcriptome approaches may be increasingly utilized for reference gene selection and validation. In this study, RNA-seq was used to identify a set of novel qPCR reference genes and evaluate a panel of traditional “housekeeping” reference genes in two species of the evolutionary model plant genus Mimulus. More broadly, the methods proposed in this study can be used to harness the power of transcriptomes to identify appropriate reference genes for qPCR in any study organism, including emerging and nonmodel systems. We find that RNA-seq accurately estimates gene expression means in comparison to qPCR, and that expression means are robust to moderate environmental and genetic variation. However, measures of expression variability were only in agreement with qPCR for samples obtained from a shared environment. This result, along with transcriptome-wide comparisons, suggests that environmental changes have greater impacts on expression variability than on expression means. We discuss how this issue can be addressed through experimental design, and suggest that the ever-expanding pool of published transcriptomes represents a rich and low-cost resource for developing better reference genes for qPCR.

Mechanisms to Mitigate the Trade-Off between Growth and Defense

Plants have evolved an array of defenses against pathogens. However, mounting a defense response frequently comes with the cost of a reduction in growth and reproduction, carrying critical implications for natural and agricultural populations. This review focuses on how costs are generated and whether and how they can be mitigated. Most well-characterized growth-defense trade-offs stem from antagonistic crosstalk among hormones rather than an identified metabolic expenditure. A primary way plants mitigate such costs is through restricted expression of resistance; this can be achieved through inducible expression of defense genes or by the concentration of defense to particular times or tissues. Defense pathways can be primed for more effective induction, and primed states can be transmitted to offspring. We examine the resistance (R) genes as a case study of how the toll of defense can be generated and ameliorated. The fine-scale regulation of R genes is critical to alleviate the burden of their expression, and the genomic organization of R genes into coregulatory modules reduces costs. Plants can also recruit protection from other species. Exciting new evidence indicates that a plant's genotype influences the microbiome composition, lending credence to the hypothesis that plants shape their microbiome to enhance defense.

Transgenerational effects alter plant defence and resistance in nature

Trichomes, or leaf hairs, are epidermal extensions that take a variety of forms and perform many functions in plants, including herbivore defence. In this study, I document genetically determined variation, within-generation plasticity, and a direct role of trichomes in herbivore defence for Mimulus guttatus. After establishing the relationship between trichomes and herbivory, I test for transgenerational effects of wounding on trichome density and herbivore resistance. Patterns of interannual variation in herbivore density and the high cost of plant defence makes plant-herbivore interactions a system in which transgenerational phenotypic plasticity (TPP) is apt to evolve. Here, I demonstrate that parental damage alters offspring trichome density and herbivore resistance in nature. Moreover, this response varies between populations. This is among the first studies to demonstrate that TPP contributes to variation in nature, and also suggests that selection can modify TPP in response to local conditions.

Genome-Wide Survey Reveals Transcriptional Differences Underlying the Contrasting Trichome Phenotypes of Two Sister Desert Poplars

Trichomes, which are widely used as an important diagnostic characteristic in plant species delimitation, play important roles in plant defense and adaptation to adverse environments. In this study, we used two sister poplar species, Populus pruinosa and Populus euphratica—which have, respectively, dense and sparse trichomes—to examine the genetic differences associated with these contrasting phenotypes. The results showed that 42 and 45 genes could be identified as candidate genes related to trichomes in P. pruinosa and P. euphratica, respectively; most of these genes possessed high degrees of diversification in their coding sequences, but they were similar in intron/exon structure in the two species. We also found that most of the candidate trichome genes were expressed at higher levels in P. pruinosa, which has dense trichomes, than in P. euphratica, where there are few trichomes. Based on analyses of transcriptional profiles, a total of 195 genes, including many transcription factors, were found to show distinct differences in expression. The results of gene function annotation suggested that the genes identified as having contrasting levels of expression level are mainly associated with trichome elongation, ATPase activity, and hormone transduction. Changes in the expression of these and other related genes with high sequence diversification may have contributed to the contrast in the pattern of trichome phenotypes between the two species.

DNA methylation mediates genetic variation for adaptive transgenerational plasticity

Environmental stresses experienced by individual parents can influence offspring phenotypes in ways that enhance survival under similar conditions. Although such adaptive transgenerational plasticity is well documented, its transmission mechanisms are generally unknown. One possible mechanism is environmentally induced DNA methylation changes. We tested this hypothesis in the annual plant Polygonum persicaria, a species known to express adaptive transgenerational plasticity in response to parental drought stress. Replicate plants of 12 genetic lines (sampled from natural populations) were grown in dry versus moist soil. Their offspring were exposed to the demethylating agent zebularine or to control conditions during germination and then grown in dry soil. Under control germination conditions, the offspring of drought-stressed parents grew longer root systems and attained greater biomass compared with offspring of well-watered parents of the same genetic lines. Demethylation removed these adaptive developmental effects of parental drought, but did not significantly alter phenotypic expression in offspring of well-watered parents. The effect of demethylation on the expression of the parental drought effect varied among genetic lines. Differential seed provisioning did not contribute to the effect of parental drought on offspring phenotypes. These results demonstrate that DNA methylation can mediate adaptive, genotype-specific effects of parental stress on offspring phenotypes.

Evidence for an Early Origin of Vernalization Responsiveness in Temperate Pooideae Grasses

The ability of plants to match their reproductive output with favorable environmental conditions has major consequences both for lifetime fitness and geographic patterns of diversity. In temperate ecosystems, some plant species have evolved the ability to use winter nonfreezing cold (vernalization) as a cue to ready them for spring flowering. However, it is unknown how important the evolution of vernalization responsiveness has been for the colonization and subsequent diversification of taxa within the northern and southern temperate zones. Grasses of subfamily Pooideae, including several important crops, such as wheat (Triticum aestivum), barley (Hordeum vulgare), and oats (Avena sativa), predominate in the northern temperate zone, and it is hypothesized that their radiation was facilitated by the early evolution of vernalization responsiveness. Predictions of this early origin hypothesis are that a response to vernalization is widespread within the subfamily and that the genetic basis of this trait is conserved. To test these predictions, we determined and reconstructed vernalization responsiveness across Pooideae and compared expression of wheat vernalization gene orthologs VERNALIZATION1 (VRN1) and VRN3 in phylogenetically representative taxa under cold and control conditions. Our results demonstrate that vernalization responsive Pooideae species are widespread, suggesting that this trait evolved early in the lineage and that at least part of the vernalization gene network is conserved throughout the subfamily. These results are consistent with the hypothesis that the evolution of vernalization responsiveness was important for the initial transition of Pooideae out of the tropics and into the temperate zone.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

Small interfering RNAs from bidirectional transcripts of GhMML3_A12 regulate cotton fiber development

Natural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development. Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development. The extremely low expression of GhMML3_A12 in N1 is associated with NAT production, driven by its 3' antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton.