Methylation of microRNA genes regulates gene expression in bisexual flower development in andromonoecious poplar

Methylation of microRNA genes and its effect on secondary xylem development of stem in poplar

MicroRNAs (miRNAs) and DNA methylation are both vital regulators of gene expression. DNA methylation can affect the transcription of miRNAs, just like coding genes, through methylating the CpG islands in the gene regions of miRNAs. Although previous studies have shown that DNA methylation and miRNAs can each be involved in the process of wood formation, the relationship between the two has been relatively little studied in plant wood formation. Studies have shown that the second internode (IN2) (from top to bottom) of 3-month-old poplar trees can represent the primary stage of poplar stem development and IN8 can represent the secondary stage. There were also significant differences in DNA methylation patterns and miRNA expression patterns obtained from PS and SS. In this study, we first interactively analyzed methylation and miRNA sequencing data to identify 43 differentially expressed miRNAs regulated by differential methylation from the primary stage and secondary stage, which were found to be involved in multiple biological processes related to wood formation by enrichment analysis. In addition, six miRNA/target gene modules were finally identified as potentially involved in secondary xylem development of poplar stems through degradome sequencing and functional analysis. In conclusion, this study provides important reference information on the mechanism of interaction between different regulatory pathways of wood formation.

Genome-wide identification and comparative analysis of DNA methyltransferase and demethylase gene families in two ploidy Cyclocarya paliurus and their potential function in heterodichogamy

BACKGROUND

DNA methylation is one of the most abundant epigenetic modifications, which plays important roles in flower development, sex differentiation, and regulation of flowering time. Its pattern is affected by cytosine-5 DNA methyltransferase (C5-MTase) and DNA demethylase (dMTase). At present, there are no reports on C5-MTase and dMTase genes in heterodichogamous Cyclocarya paliurus.

RESULTS

In this study, 6 CpC5-MTase and 3 CpdMTase genes were identified in diploid (2n = 2 ×  = 32) C. paliurus, while 20 CpC5-MTase and 13 CpdMTase genes were identified in autotetraploid (2n = 4 ×  = 64). 80% of identified genes maintained relatively fixed positions on chromosomes during polyploidization. In addition, we found that some DRM subfamily members didn't contain the UBA domain. The transcript abundance of CpC5-MTase and CpdMTase in male and female flowers of two morphs (protandry and protogyny) from diploidy was analyzed. Results showed that all genes were significantly up-regulated at the stage of floral bud break (S2), but significantly down-regulated at the stage of flower maturation (S4). At S2, some CpC5-MTase genes showed higher expression levels in PG-M than in PG-F, whereas some CpdMTase genes showed higher expression levels in PA-M than in PA-F. In addition, these genes were significantly associated with gibberellin synthesis-related genes (e.g. DELLA and GID1), suggesting that DNA methylation may play a role in the asynchronous floral development process through gibberellin signal.

CONCLUSIONS

These results broaden our understanding of the CpC5-MTase and CpdMTase genes in diploid and autotetraploid C. paliurus, and provide a novel insight into regulatory mechanisms of DNA methylation in heterodichogamy.

The role of DNA methylation in the maintenance of phenotypic variation induced by grafting chimerism in Brassica

Grafting facilitates the interaction between heterologous cells with different genomes, resulting in abundant phenotypic variation, which provides opportunities for crop improvement. However, how grafting-induced variation occurs and is transmitted to progeny remains elusive. A graft chimera, especially a periclinal chimera, which has genetically distinct cell layers throughout the plant, is an excellent model to probe the molecular mechanisms of grafting-induced variation maintenance. Here we regenerated a plant from the T-cell layer of a periclinal chimera, TCC (where the apical meristem was artificially divided into three cell layers - from outside to inside, L1, L2, and L3; T = Tuber mustard, C = red Cabbage), named rTTT0 (r = regenerated). Compared with the control (rsTTT, s = self-grafted), rTTT0 had multiple phenotypic variations, especially leaf shape variation, which could be maintained in sexual progeny. Transcriptomes were analyzed and 58 phenotypic variation-associated genes were identified. Whole-genome bisulfite sequencing analyses revealed that the methylome of rTTT0 was changed, and the CG methylation level was significantly increased by 8.74%. In rTTT0, the coding gene bodies are hypermethylated in the CG context, while their promoter regions are hypomethylated in the non-CG context. DNA methylation changes in the leaf shape variation-associated coding genes, ARF10, IAA20, ROF1, and TPR2, were maintained for five generations of rTTT0. Interestingly, grafting chimerism also affected transcription of the microRNA gene (MIR), among which the DNA methylation levels of the promoters of three MIRs associated with leaf shape variation were changed in rTTT0, and the DNA methylation modification of MIR319 was maintained to the fifth generation of selfed progeny of rTTT0 (rTTT5). These findings demonstrate that DNA methylation of coding and non-coding genes plays an important role in heterologous cell interaction-induced variation formation and its transgenerational inheritance.

Methylation in DNA, histone, and RNA during flowering under stress condition: A review

Flowering is the most critical transition period in the whole lifecycle of plants, and it is a highly sensitive period to stress. New combinations of temperature, drought stress, carbon dioxide and other abiotic/biotic conditions resulting from contemporary climate change affect the flowering process. Plants have evolved several strategies to deal with environmental stresses, including epigenetic modifications. Numerous studies show that environmental stresses trigger methylation/demethylation during flowering to preserve/accelerate plant lifecycle. What's more, histone and DNA methylation can be induced to respond to stresses, resulting in changes of flowering gene expression and enhancing stress tolerance in plants. Furthermore, RNA methylation may influence stress-regulated flowering by regulating mRNA stability and antioxidant mechanism. Our review presents the involvement of methylation in stress-repressed and stress-induced flowering. The crosstalk between methylation and small RNAs, phytohormones and exogenous substances (such as salicylic acid, nitric oxide) during flowering under different stresses were discussed. The latest regulatory evidence of RNA methylation in stress-regulated flowering was collected for the first time. Meanwhile, the limited evidences of methylation in biotic stress-induced flowering were summarized. Thus, the review provides insights into understanding of methylation mechanism in stress-regulated flowering and makes use for the development of regulating plant flowering at epigenetic level in the future.

Transcriptome and methylome changes in two contrasting mungbean genotypes in response to drought stress

Background

Due to drought stress, the growth, distribution, and production of mungbean is severely restricted. Previous study combining physiological and transcriptomic data indicated different genotypes of mungbean exhibited variable responses when exposed to drought stress. Aside from the genetic variation, the modifications of environmentally induced epigenetics alterations on mungbean drought-stress responses were still elusive.

Results

In this study, firstly, we compared the drought tolerance capacity at seedling stage by detecting physiological parameters in two contrasting genotypes wild mungbean 61 and cultivar 70 in response to drought stress. We found that wild mungbean 61 showed lower level of MDA and higher levels of POD and CAT, suggesting wild mungbean 61 exhibited stronger drought resistance. Transcriptomic analysis indicated totally 2859 differentially expressed genes (DEGs) were detected when 70 compared with 61 (C70 vs C61), and the number increased to 3121 in the comparison of drought-treated 70 compared with drought-treated 61 (D70 vs D61). In addition, when drought-treated 61 and 70 were compared with their controls, the DEGs were 1117 and 185 respectively, with more down-regulated DEGs than up-regulated in D61 vs C61, which was opposite in D70 vs C70. Interestingly, corresponding to this, after drought stress, more hypermethylated differentially methylated regions (DMRs) in 61 were detected and more hypomethylated DMRs in 70 were detected. Further analysis suggested that the main variations between 61 and 70 existed in CHH methylation in promoter. Moreover, the preference of methylation status alterations in D61 vs C61 and D70 vs C70 also fell in CHH sequence context. Further analysis of the correlation between DMRs and DEGs indicated in both D61 vs C61 and D70 vs C70, the DMRs in gene body was significantly negatively correlated with DEGs.

Conclusions

The physiological parameters in this research suggested that wild mungbean 61 was more resistant to drought stress, with more hypermethylated DMRs and less hypomethylated DMRs after drought stress, corresponding to more down-regulated DEGs than up-regulated DEGs. Among the three DNA methylation contexts CG, CHG, and CHH, asymmetric CHH contexts were more dynamic and prone to be altered by drought stress and genotypic variations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08315-z.

Concerto on Chromatin: Interplays of Different Epigenetic Mechanisms in Plant Development and Environmental Adaptation

Epigenetic mechanisms such as DNA methylation, histone post-translational modifications, chromatin remodeling, and noncoding RNAs, play important roles in regulating plant gene expression, which is involved in various biological processes including plant development and stress responses. Increasing evidence reveals that these different epigenetic mechanisms are highly interconnected, thereby contributing to the complexity of transcriptional reprogramming in plant development processes and responses to environmental stresses. Here, we provide an overview of recent advances in understanding the epigenetic regulation of plant gene expression and highlight the crosstalk among different epigenetic mechanisms in making plant developmental and stress-responsive decisions. Structural, physical, transcriptional and metabolic bases for these epigenetic interplays are discussed.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

The Sw5a gene confers resistance to ToLCNDV and triggers an HR response after direct AC4 effector recognition

Several attempts have been made to identify antiviral genes against Tomato leaf curl New Delhi virus (ToLCNDV) and related viruses. This has led to the recognition of Ty genes (Ty1-Ty6), which have been successful in developing virus-resistant crops to some extent. Owing to the regular appearance of resistance-breaking strains of these viruses, it is important to identify genes related to resistance. In the present study, we identified a ToLCNDV resistance (R) gene, SlSw5a, in a ToLCNDV-resistant tomato cultivar, H-88-78-1, which lacks the known Ty genes. The expression of SlSw5a is controlled by the transcription factor SlMyb33, which in turn is regulated by microRNA159 (sly-miR159). Virus-induced gene silencing of either SlSw5a or SlMyb33 severely increases the disease symptoms and viral titer in leaves of resistant cultivar. Moreover, in SlMyb33-silenced plants, the relative messenger RNA level of SlSw5a was reduced, suggesting SlSw5a is downstream of the sly-miR159-SlMyb33 module. We also demonstrate that SlSw5a interacts physically with ToLCNDV-AC4 (viral suppressor of RNA silencing) to trigger a hypersensitive response (HR) and generate reactive oxygen species at infection sites to limit the spread of the virus. The "RTSK" motif in the AC4 C terminus is important for the interaction, and its mutation completely abolishes the interaction with Sw5a and HR elicitation. Overall, our research reports an R gene against ToLCNDV and establishes a connection between the upstream miR159-Myb33 module and its downstream target Sw5a to activate HR in the tomato, resulting in geminivirus resistance.

Transcriptomic and epigenomic remodeling occurs during vascular cambium periodicity in Populus tomentosa

Trees in temperate regions exhibit evident seasonal patterns, which play vital roles in their growth and development. The activity of cambial stem cells is the basis for regulating the quantity and quality of wood, which has received considerable attention. However, the underlying mechanisms of these processes have not been fully elucidated. Here we performed a comprehensive analysis of morphological observations, transcriptome profiles, the DNA methylome, and miRNAs of the cambium in Populus tomentosa during the transition from dormancy to activation. Anatomical analysis showed that the active cambial zone exhibited a significant increase in the width and number of cell layers compared with those of the dormant and reactivating cambium. Furthermore, we found that differentially expressed genes associated with vascular development were mainly involved in plant hormone signal transduction, cell division and expansion, and cell wall biosynthesis. In addition, we identified 235 known miRNAs and 125 novel miRNAs. Differentially expressed miRNAs and target genes showed stronger negative correlations than other miRNA/target pairs. Moreover, global methylation and transcription analysis revealed that CG gene body methylation was positively correlated with gene expression, whereas CHG exhibited the opposite trend in the downstream region. Most importantly, we observed that the number of CHH differentially methylated region (DMR) changes was the greatest during cambium periodicity. Intriguingly, the genes with hypomethylated CHH DMRs in the promoter were involved in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interactions during vascular cambium development. These findings improve our systems-level understanding of the epigenomic diversity that exists in the annual growth cycle of trees.

Differential DNA methylation and gene expression during development of reproductive and vegetative organs in Ilex species

Differential epigenetic (DNA cytosine methylation) and gene expression patterns were investigated in reproductive and vegetative organs from Ilex paraguariensis and I. dumosa, at distinct developmental stages. We aimed at contributing towards elucidating major molecular changes underlying the sexual differentiation processes which, in these dioecious species, are completely unknown. Simultaneously, as a first step towards the development of an early sexing system, we searched for promising molecular markers. This was assessed through Methylation Sensitive Amplified Polymorphism (MSAP) and Amplified Fragment Length Polymorphism on cDNA (cDNA-AFLP) techniques, applying discriminant multivariate analyses, and bioinformatic characterization of differential fragments. A significant positive correlation was found between epigenetic and indirect 'genetic' information for both species, indicating influence of the genetic background on the epigenetic variation. Higher epigenetic than genetic diversities were estimated. Our outcomes showed up to 1.86 times more representation of mCG subepiloci than mCCG in all organs sampled. Along the maturing stages of floral buds, the frequency of mCG evidenced an incremental trend, whereas mCCG and unmethylated conditions showed opposite tendencies. Reproductive and vegetative samples tended to cluster apart based on epigenetic patterns; at gene expression level, organs exhibited clear-cut distinctive patterns, nonetheless profiles of young leaves and floral primordia resemble. Epigenetic and expression data allowed discrimination of I. dumosa´s samples according to the gender of the donor; more elusive patterns were observed for I. paraguariensis. In total, 102 differentially methylated and expressed fragments were characterized bioinformatically. Forty-three were annotated in various functional categories; four candidate markers were validated through qPCR, finding statistical differences among organs but not among sexes. The methylation condition of epilocus C13m33 appears as indicative of gender in both species. Thirty-three organ-specific and 34 gender-specific methylated markers were discriminated and deserve further research, particularly those expressed in leaves. Our study contributes concrete candidate markers with potential for practical application.

Title: Hypermethylation of miRNA Genes During Nodule Development

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean (Glycine max) through the analysis of methylC-seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.

Evidences for a role of two Y-specific genes in sex determination in Populus deltoides

Almost all plants in the genus Populus are dioecious (i.e. trees are either male or female), but it is unknown whether dioecy evolved in a common ancestor or independently in different subgenera. Here, we sequence the small peritelomeric X- and Y-linked regions of P. deltoides chromosome XIX. Two genes are present only in the Y-linked region. One is a duplication of a non-Y-linked, female-specifically expressed response regulator, which produces siRNAs that block this gene's expression, repressing femaleness. The other is an LTR/Gypsy transposable element family member, which generates long non-coding RNAs. Overexpression of this gene in A. thaliana promotes androecium development. We also find both genes in the sex-determining region of P. simonii, a different poplar subgenus, which suggests that they are both stable components of poplar sex-determining systems. By contrast, only the duplicated response regulator gene is present in the sex-linked regions of P. davidiana and P. tremula. Therefore, findings in our study suggest dioecy may have evolved independently in different poplar subgenera.

Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Relationships between genome methylation, levels of non-coding RNAs, mRNAs and metabolites in ripening tomato fruit

Ripening of tomato fruit is a complex tightly orchestrated developmental process that involves multiple physiological and metabolic changes that render fruit attractive, palatable and nutritious. Ripening requires initiation, activation and coordination of key pathways at the transcriptional and post-transcriptional levels that lead to ethylene synthesis and downstream ripening events determining quality. We studied wild-type, Gr and r mutant fruits at the coding and non-coding transcriptomic, metabolomic and genome methylation levels. Numerous differentially expressed non-coding RNAs were identified and quantified and potential competing endogenous RNA regulation models were constructed. Multiple changes in gene methylation were linked to the ethylene pathway and ripening processes. A combined analysis of changes in genome methylation, long non-coding RNAs, circular RNAs, micro-RNAs and fruit metabolites revealed many differentially expressed genes (DEGs) with differentially methylated regions encoding transcription factors and key enzymes related to ethylene or carotenoid pathways potentially targeted by differentially expressed non-coding RNAs. These included ACO2 (targeted by MSTRG.59396.1 and miR396b), CTR1 (targeted by MSTRG.43594.1 and miR171b), ERF2 (targeted by MSTRG.183681.1), ERF5 (targeted by miR9470-3p), PSY1 (targeted by MSTRG.95226.7), ZISO (targeted by 12:66127788|66128276) and NCED (targeted by MSTRG.181568.2). Understanding the functioning of this intricate genetic regulatory network provides new insights into the underlying integration and relationships between the multiple events that collectively determine the ripe phenotype.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Analysis of Alternative Splicing and Alternative Polyadenylation in Populus alba var. pyramidalis by Single-Molecular Long-Read Sequencing

Poplars are worldwidely cultivated with ecologically and economically important value. Populus alba var. pyramidalis (= P. bolleana) is a main tree of the farmland shelter-belt system in the arid region of Northwest China due to its rapid growth, erect stems, and high biomass production. However, the full-length messenger RNA (mRNA) sequences and complete structure of P. alba var. pyramidalis remain unclear. In this study, using single-molecular real-time (SMRT) and next-generation high-throughput sequencing (NGS) platform, we sequenced transcripts from leaf, root, xylem, and phloem of P. alba var. pyramidalis, to obtain the full-length mRNA transcripts and annotate the complete structure. In total, 86,327 mapped full-length non-chimeric (FLNC) reads were identified, with 705 previously unannotated loci and 3,410 long noncoding RNAs (lncRNAs) and 174 fusion genes found. Alternative spicing (AS) events were detected in 7,536 genes, of which 4,652 genes had multiple AS events. A total of 10,213 alternative polyadenylation (APA) sites were identified, with two or more APA sites observed in 2,212 genes. Our transcriptome data provided the full-length sequences and gene isoforms of transcripts for P. alba var. pyramidalis, which will be helpful in improving our understanding for the genome annotation and gene structures of P. alba var. pyramidalis.

Linkage-linkage disequilibrium dissection of the epigenetic quantitative trait loci (epiQTLs) underlying growth and wood properties in Populus

Increasing evidence indicates that DNA methylation is heritable and serves as an essential marker contributing to phenotypic variation. Linkage-linkage disequilibrium mapping was used to decipher the epigenetic architecture underlying nine growth and wood property traits in a linkage population (550 F₁ progeny) and a natural population (435 unrelated individuals) of Populus using methylation-sensitive amplification polymorphism (MSAP)-based analysis. The interactions between genetic and epigenetic variants in the causative genes was further unveiled using expression quantitative trait methylation (eQTM) and nucleotide (eQTN) mapping strategies. A total of 163 epigenetic quantitative trait loci (epiQTLs; LOD ≥ 3.0), explaining 1.7-44.5% of phenotypic variations, were mapped to a high-resolution epigenetic map with 19 linkage groups, which was supported by the significant MSAP associations (P < 0.001) in the two populations. There were 23 causal genes involved in growth regulation and wood formation, whose markers were located in epiQTLs and associated with the same traits in both populations. Further eQTN and eQTM mapping showed that causal genetic and epigenetic variants within the 23 candidate genes may interact more in trans in gene expression and phenotype. The present study provides strategies for investigating epigenetic architecture and the interaction between genetic and epigenetic variants modulating complex traits in forest trees.

Integrated Methylome and Transcriptome Analysis between the CMS-D2 Line ZBA and Its Maintainer Line ZB in Upland Cotton

DNA methylation is an important epigenetic modification involved in multiple biological processes. Altered methylation patterns have been reported to be associated with male sterility in some plants, but their role in cotton cytoplasmic male sterility (CMS) remains unclear. Here, integrated methylome and transcriptome analyses were conducted between the CMS-D2 line ZBA and its near-isogenic maintainer line ZB in upland cotton. More methylated cytosine sites (mCs) and higher methylation levels (MLs) were found among the three sequence contexts in ZB compared to ZBA. A total of 4568 differentially methylated regions (DMRs) and 2096 differentially methylated genes (DMGs) were identified. Among the differentially expressed genes (DEGs) associated with DMRs (DMEGs), 396 genes were upregulated and 281 genes were downregulated. A bioinformatics analysis of these DMEGs showed that hyper-DEGs were significantly enriched in the “oxidative phosphorylation” pathway. Further qRT-PCR validation indicated that these hypermethylated genes (encoding the subunits of mitochondrial electron transport chain (ETC) complexes I and V) were all significantly upregulated in ZB. Our biochemical data revealed a higher extent of H₂O₂ production but a lower level of adenosine triphosphate (ATP) synthesis in CMS-D2 line ZBA. On the basis of the above results, we propose that disrupted DNA methylation in ZBA may disrupt the homeostasis of reactive oxygen species (ROS) production and ATP synthesis in mitochondria, triggering a burst of ROS that is transferred to the nucleus to initiate programmed cell death (PCD) prematurely, ultimately leading to microspore abortion. This study illustrates the important role of DNA methylation in cotton CMS.

The poplar pangenome provides insights into the evolutionary history of the genus

The genus Populus comprises a complex amalgam of ancient and modern species that has become a prime model for evolutionary and taxonomic studies. Here we sequenced the genomes of 10 species from five sections of the genus Populus, identified 71 million genomic variations, and observed new correlations between the single-nucleotide polymorphism-structural variation (SNP-SV) density and indel-SV density to complement the SNP-indel density correlation reported in mammals. Disease resistance genes (R genes) with heterozygous loss-of-function (LOF) were significantly enriched in the 10 species, which increased the diversity of poplar R genes during evolution. Heterozygous LOF mutations in the self-incompatibility genes were closely related to the self-fertilization of poplar, suggestive of genomic control of self-fertilization in dioecious plants. The phylogenetic genome-wide SNPs tree also showed possible ancient hybridization among species in sections Tacamahaca, Aigeiros, and Leucoides. The pangenome resource also provided information for poplar genetics and breeding.

Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica

The genus Populus is an effective model in tree genetics. This genus includes dioecious species and, recently, whole genome resequencing of P. trichocarpa and P. balsamifera enabled the identification of sex-linked regions and sex-associated single nucleotide polymorphisms (SNPs). These results created new opportunities to study sex determination in poplars. In the present work, we performed deep sequencing of genes encoding METHYLTRANSFERASE1 (MET1) and homolog of ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17), which are localized in a sex-linked region of Populus genome and contain a number of sex-associated SNPs. Amplicon libraries for 38 samples of P. × sibirica (19 males and 19 females) were sequenced on MiSeq Illumina (300 nt paired-end reads) and approximately 4000× coverage was obtained for each sample. In total, from 80 to 179 SNPs were detected in poplar individuals for MET1, and from 16 to 49 SNPs were detected for ARR17. We identified 17 sex-specific SNPs (11 in MET1 and 6 in ARR17) - they were present in all males but absent in all females. For identified sex-specific SNP sites, females were homozygous, while males were heterozygous. Moreover, colocation of sex-specific SNPs confirming the XY sex-determination system of poplars was revealed: in one allelic variant, males had the same nucleotides as females, while in the other, sex-specific SNPs were present. Based on the data obtained, we developed and successfully applied a high-resolution melting-based approach for sex identification in poplars. The developed molecular markers are useful for distinguishing between male and female poplars in scientific research and can also be applied to select male-only genotypes for use in city landscaping and production of paper, pulp, and biofuel.

The Effects of DNA Methylation Inhibition on Flower Development in the Dioecious Plant Salix Viminalis

DNA methylation, an important epigenetic modification, regulates the expression of genes and is therefore involved in the transitions between floral developmental stages in flowering plants. To explore whether DNA methylation plays different roles in the floral development of individual male and female dioecious plants, we injected 5-azacytidine (5-azaC), a DNA methylation inhibitor, into the trunks of female and male basket willow (Salix viminalis L.) trees before flower bud initiation. As expected, 5-azaC decreased the level of DNA methylation in the leaves of both male and female trees during floral development; however, it increased DNA methylation in the leaves of male trees at the flower transition stage. Furthermore, 5-azaC increased the number, length and diameter of flower buds in the female trees but decreased these parameters in the male trees. The 5-azaC treatment also decreased the contents of soluble sugars, starch and reducing sugars in the leaves of the female plants, while increasing them in the male plants at the flower transition stage; however, this situation was largely reversed at the flower development stage. In addition, 5-azaC treatment decreased the contents of auxin indoleacetic acid (IAA) in both male and female trees at the flower transition stage. These results indicate that hypomethylation in leaves at the flower transition stage promotes the initiation of flowering and subsequent floral growth in Salix viminalis, suggesting that DNA methylation plays a similar role in vegetative–reproductive transition and early floral development. Furthermore, methylation changes during the vegetative–reproductive transition and floral development were closely associated with the biosynthesis, metabolism and transportation of carbohydrates and IAA. These results provide insight into the epigenetic regulation of carbohydrate accumulation.

Sex and the flower - developmental aspects of sex chromosome evolution

Background

The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants.

Scope

This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed.

Conclusions

We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.

Integrative Analysis of Transcriptomic and Methylomic Data in Photoperiod-Dependent Regulation of Cucumber Sex Expression

The cucumber (Cucumis sativus) is characterized by its diversity and seasonal plasticity in sexual type. A long day length condition significantly decreased the cucumber female flower ratio by 17.7-52.9%, and the effect of photoperiod treatment is more significant under low temperature than under high temperature. Transcriptome analysis indicates that the photoperiod treatment preferentially significantly influenced flower development processes, particularly MADS-box genes in shoot apices. The long-day treatment resulted in predominantly transposable element (TE)- and gene-associated CHH-types of DNA methylation changes. Nevertheless, there was significant enrichment of CG- and CHG-types of DNA methylation changes nearing transcription start sites (TSSs)/transcription end sites (TESs) and gene bodies, respectively. Predominantly negative association between differentially methylated regions (DMRs) and differentially expressed genes (DEGs) were observed which implied epiregulation of DEGs. Two MADS-box genes that were significantly downregulated by long photoperiod showed significant hypermethylation in promoter regions that is essentially TE-rich. This study indicates MADS-box genes which are partially regulated by promoter methylation state may mediate photoperiod-dependent regulation of cucumber sex expression.

Intra-arterial transplantation of human bone marrow mesenchymal stem cells (hBMMSCs) improves behavioral deficits and alters gene expression in rodent stroke model

Stroke is a multi-factorial polygenic disease and is a major cause of death and adult disability. Administration of bone marrow stem cells protects ischemic rat brain by facilitating recovery of neurological functions. But the molecular mechanism of stem cells action and their effect on gene expression is not well explored. In this study, we have transplanted 1 × 10⁶ human bone marrow mesenchymal stem cells (hBMMSCs) in middle cerebral artery occluded (MCAo) adult male Wistar rats through intracarotid artery route at 24 h after surgery. Motor behavioral tests (rotarod and open field) were performed to assess the changes in motor functions at day 0 and day1, 4, 8 and 14. The expression of studied genes at mRNA and protein level was quantified by using Q-PCR and western blotting, respectively. Further, we have assessed the methylation pattern of promoter of these genes by using methylation-specific PCR. Data were analyzed statistically and correlated. A significant improvement in behavioral deficits was observed in stem cells treated group after 14th day post stroke. Significantly (p < 0.05) increased mRNA and protein levels of brain derived neurotrophic factor and ANP genes in hBMMSCs treated group along with decrease in methylation level at their promoter was observed. On the other hand, significantly decreased mRNA and protein level of TSP1 and WNK1 in hBMMSCs treated group was observed. In conclusion, hBMMSCs administration significantly improves the behavioral deficits by improving motor and locomotor coordination. The promoter of TSP1 and WNK1 genes was found to be hyper-methylated in hBMMSCs group resulting in their decreased expression while the promoter of ANP and brain derived neurotrophic factor was found to be hypo-methylated. This study might shed a light on how hBMMSCs affect the gene expression by modulating methylation status.

In Silico Analysis of Small RNAs Suggest Roles for Novel and Conserved miRNAs in the Formation of Epigenetic Memory in Somatic Embryos of Norway Spruce

Epigenetic memory in Norway spruce affects the timing of bud burst and bud set, vitally important adaptive traits for this long-lived forest species. Epigenetic memory is established in response to the temperature conditions during embryogenesis. Somatic embryogenesis at different epitype inducing (EpI) temperatures closely mimics the natural processes of epigenetic memory formation in seeds, giving rise to epigenetically different clonal plants in a reproducible and predictable manner, with respect to altered bud phenology. MicroRNAs (miRNAs) and other small non-coding RNAs (sRNAs) play an essential role in the regulation of plant gene expression and may affect this epigenetic mechanism. We used NGS sequencing and computational in silico methods to identify and profile conserved and novel miRNAs among small RNAs in embryogenic tissues of Norway spruce at three EpI temperatures (18, 23 and 28°C). We detected three predominant classes of sRNAs related to a length of 24 nt, followed by a 21–22 nt class and a third 31 nt class of sRNAs. More than 2100 different miRNAs within the prevailing length 21–22 nt were identified. Profiling these putative miRNAs allowed identification of 1053 highly expressed miRNAs, including 523 conserved and 530 novels. 654 of these miRNAs were found to be differentially expressed (DEM) depending on EpI temperature. For most DEMs, we defined their putative mRNA targets. The targets represented mostly by transcripts of multiple-repeats proteins, like TIR, NBS-LRR, PPR and TPR repeat, Clathrin/VPS proteins, Myb-like, AP2, etc. Notably, 124 DE miRNAs targeted 203 differentially expressed epigenetic regulators. Developing Norway spruce embryos possess a more complex sRNA structure than that reported for somatic tissues. A variety of the predicted miRNAs showed distinct EpI temperature dependent expression patterns. These putative EpI miRNAs target spruce genes with a wide range of functions, including genes known to be involved in epigenetic regulation, which in turn could provide a feedback process leading to the formation of epigenetic marks. We suggest that TIR, NBS and LRR domain containing proteins could fulfill more general functions for signal transduction from external environmental stimuli and conversion them into molecular response. Fine-tuning of the miRNA production likely participates in both developmental regulation and epigenetic memory formation in Norway spruce.

The association of changes in DNA methylation with temperature-dependent sex determination in cucumber

Cucumber (Cucumis sativus L.) is characterized by its diverse and flexible sexual types. Here, we evaluated the effect of low temperature (LT) exposure on cucumber femaleness under short-day conditions. Shoot apices were subjected to whole-genome bisulfate sequencing (WGBS), mRNA-seq, and sRNA-seq. The results showed that temperature had a substantial and global impact on transposable element (TE)-related small RNA-directed DNA methylation (RdDM) mechanisms, resulting in large amounts of CHH-type cytosine demethylation. In the cucumber genome, TEs are common in regions near genes that are also subject to DNA demethylation. TE-gene interactions showed very strong reactions to LT treatment, as nearly 80% of the differentially methylated regions (DMRs) were distributed in genic regions. Demethylation near genes led to the co-ordinated expression of genes and TEs. More importantly, genome-wide de novo methylation changes also resulted in small amounts of CG- and CHG-type DMRs. Methylation changes in CG-DMRs located <600 bp from the transcription start and end sites (TSSs/TESs) negatively correlated with transcription changes in differentially expressed genes (DEGs), probably indicating epiregulation. Ethylene is called the 'sex hormone' of cucumbers. We observed the up-regulation of ethylene biosynthesis-related CsACO3 and the down-regulation of an Arabidopsis RAP2.4-like ethylene-responsive (AP2/ERF) transcription factor, demonstrating the inferred epiregulation. Our study characterized the response of the apex methylome to LT and predicted the possible epiregulation of temperature-dependent sex determination (TSD) in cucumber.

Rosaceae Fruit Development, Ripening and Post-harvest: An Epigenetic Perspective

Rosaceae is a family with an extraordinary spectrum of fruit types, including fleshy peach, apple, and strawberry that provide unique contributions to a healthy diet for consumers, and represent an excellent model for studying fruit patterning and development. In recent years, many efforts have been made to unravel regulatory mechanism underlying the hormonal, transcriptomic, proteomic and metabolomic changes occurring during Rosaceae fruit development. More recently, several studies on fleshy (tomato) and dry (Arabidopsis) fruit model have contributed to a better understanding of epigenetic mechanisms underlying important heritable crop traits, such as ripening and stress response. In this context and summing up the results obtained so far, this review aims to collect the available information on epigenetic mechanisms that may provide an additional level in gene transcription regulation, thus influencing and driving the entire Rosaceae fruit developmental process. The whole body of information suggests that Rosaceae fruit could become also a model for studying the epigenetic basis of economically important phenotypes, allowing for their more efficient exploitation in plant breeding.

Site-Dependent Differences in DNA Methylation and Their Impact on Plant Establishment and Phosphorus Nutrition in Populus trichocarpa

The propagation via clonal stem cuttings is a frequent practice in tree plantations. Despite their clonal origin, the trees establish differently according to weather, temperature and nutrient availability, as well as the presence of various stresses. Here, clonal Populus trichocarpa (cv. Muhle Larson) cuttings from different sites were transferred into a common, fully nutrient supplied environment. Despite identical underlying genetics, stem cuttings derived from sites with lower phosphorus availability established worse, independent of phosphorus (P) level after transplantation. Differential growth of material from the sites was reflected in differences in the whole genome DNA methylome. Methylation differences were sequence context-dependent, but differentially methylated regions (DMRs) were apparently unrelated to P nutrition genes. Despite the undisputed negative general correlation of DNA promoter methylation with gene repression, only few of the top-ranked DMRs resulted in differential gene expression in roots or shoots. However, differential methylation was associated with site-dependent, different total amounts of microRNAs (miRNAs), with few miRNAs sequences directly targeted by differential methylation. Interestingly, in roots and shoots, the miRNA amount was dependent on the previous habitat and changed in roots in a habitat-dependent way under phosphate starvation conditions. Differentially methylated miRNAs, together with their target genes, showed P-dependent expression profiles, indicating miRNA expression differences as a P-related epigenetic modification in poplar. Together with differences in DNA methylation, such epigenetic mechanisms may explain habitat or seasonal memory in perennials and site-dependent growth performances.

Stable methylation of a non-coding RNA gene regulates gene expression in response to abiotic stress in Populus simonii

DNA methylation plays important roles in responses to environmental stimuli. However, in perennial plants, the roles of DNA methylation in stress-specific adaptions to different abiotic stresses remain unclear. Here, we present a systematic, comparative analysis of the methylome and gene expression in poplar under cold, osmotic, heat, and salt stress conditions from 3h to 24h. Comparison of the stress responses revealed different patterns of cytosine methylation in response to the four abiotic stresses. We isolated and sequenced 1376 stress-specific differentially methylated regions (SDMRs); annotation revealed that these SDMRs represent 1123 genes encoding proteins, 16 miRNA genes, and 17 long non-coding RNA (lncRNA) genes. The SDMR162 region, consisting of Psi-MIR396e and PsiLNCRNA00268512, is regulated by epigenetic pathways and we speculate that PsiLNCRNA00268512 regulates miR396e levels by acting as a target mimic. The ratios of methylated cytosine declined to ~35.1% after 1 month of recovery from abiotic stress and to ~15.3% after 6 months. Among methylated miRNA genes, only expression of the methylation-regulated gene MIRNA6445a showed long-term stability. Our data provide a strong basis for future work and improve our understanding of the effect of epigenetic regulation of non-coding RNA expression, which will enable in-depth functional analysis.

The Quest for Molecular Regulation Underlying Unisexual Flower Development

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled "Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber" published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.