Genome-wide identification and expression profiling of DNA methyltransferase gene family in maize

DNA Methylation Dynamics in Response to Drought Stress in Crops

Drought is one of the most hazardous environmental factors due to its severe damage on plant growth, development and productivity. Plants have evolved complex regulatory networks and resistance strategies for adaptation to drought stress. As a conserved epigenetic regulation, DNA methylation dynamically alters gene expression and chromosome interactions in plants' response to abiotic stresses. The development of omics technologies on genomics, epigenomics and transcriptomics has led to a rapid increase in research on epigenetic variation in non-model crop species. In this review, we summarize the most recent findings on the roles of DNA methylation under drought stress in crops, including methylating and demethylating enzymes, the global methylation dynamics, the dual regulation of DNA methylation on gene expression, the RNA-dependent DNA methylation (RdDM) pathway, alternative splicing (AS) events and long non-coding RNAs (lnc RNAs). We also discuss drought-induced stress memory. These epigenomic findings provide valuable potential for developing strategies to improve crop drought tolerance.

Identification of the Gossypium hirsutum SDG Gene Family and Functional Study of GhSDG59 in Response to Drought Stress

SET-domain group histone methyltransferases (SDGs) are known to play crucial roles in plant responses to abiotic stress. However, their specific function in cotton's response to drought stress has not been well understood. This study conducted a comprehensive analysis of the SDG gene family in Gossypium hirsutum, identifying a total of 82 SDG genes. An evolutionary analysis revealed that the SDG gene family can be divided into eight subgroups. The expression analysis shows that some GhSDG genes are preferentially expressed in specific tissues, indicating their involvement in cotton growth and development. The transcription level of some GhSDG genes is induced by PEG, with GhSDG59 showing significant upregulation upon polyethylene glycol (PEG) treatment. Quantitative polymerase chain reaction (qPCR) analysis showed that the accumulation of transcripts of the GhSDG59 gene was significantly upregulated under drought stress. Further functional studies using virus-induced gene silencing (VIGS) revealed that silencing GhSDG59 reduced cotton tolerance to drought stress. Under drought conditions, the proline content, superoxide dismutase (SOD) and peroxidase (POD) enzyme activities in the GhSDG59-silenced plants were significantly lower than in the control plants, while the malondialdehyde (MDA) content was significantly higher. Transcriptome sequencing showed that silencing the GhSDG59 gene led to significant changes in the expression levels of 1156 genes. The KEGG enrichment analysis revealed that these differentially expressed genes (DEGs) were mainly enriched in the carbon metabolism and the starch and sucrose metabolism pathways. The functional annotation analysis identified known drought-responsive genes, such as ERF, CIPK, and WRKY, among these DEGs. This indicates that GhSDG59 is involved in the drought-stress response in cotton by affecting the expression of genes related to the carbon metabolism and the starch and sucrose metabolism pathways, as well as known drought-responsive genes. This analysis provides valuable information for the functional genomic study of SDGs and highlights potential beneficial genes for genetic improvement and breeding in cotton.

Genome-wide characterization of DNA methyltransferase family genes implies GhDMT6 improving tolerance of salt and drought on cotton

BACKGROUND

DNA methylation is an important epigenetic mode of genomic DNA modification and plays a vital role in maintaining epigenetic content and regulating gene expression. Cytosine-5 DNA methyltransferase (C5-MTase) are the key enzymes in the process of DNA methylation. However, there is no systematic analysis of the C5-MTase in cotton so far, and the function of DNMT2 genes has not been studied.

METHODS

In this study, the whole genome of cotton C5-MTase coding genes was identified and analyzed using a bioinformatics method based on information from the cotton genome, and the function of GhDMT6 was further validated by VIGS experiments and subcellular localization analysis.

RESULTS

33 C5-MTases were identified from three cotton genomes, and were divided into four subfamilies by systematic evolutionary analysis. After the protein domain alignment of C5-MTases in cotton, 6 highly conserved motifs were found in the C-terminus of 33 proteins involved in methylation modification, which indicated that C5-MTases had a basic catalytic methylation function. These proteins were divided into four classes based on the N-terminal difference, of which DNMT2 lacks the N-terminal regulatory domain. The expression of C5-MTases in different parts of cotton was different under different stress treatments, which indicated the functional diversity of cotton C5-MTase gene family. Among the C5-MTases, the GhDMT6 had a obvious up-regulated expression. After silencing GhDMT6 with VIGS, the phenotype of cotton seedlings under different stress treatments showed a significant difference. Compared with cotton seedlings that did not silence GhDMT6, cotton seedlings silencing GhDMT6 showed significant stress resistance.

CONCLUSION

The results show that C5-MTases plays an important role in cotton stress response, which is beneficial to further explore the function of DNMT2 subfamily genes.

Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses

DNA methylation is a highly conserved epigenetic modification involved in many biological processes, including growth and development, stress response, and secondary metabolism. DNA demethylase (DNA-deMTase) genes have been identified in some plant species; however, there are no reports on the identification and analysis of DNA-deMTase genes in Foxtail millet (Setaria italica L.). In this study, seven DNA-deMTases were identified in S. italica. These DNA-deMTase genes were divided into four subfamilies (DML5, DML4, DML3, and ROS1) by phylogenetic and gene structure analysis. Further analysis shows that the physical and chemical properties of these DNA-deMTases proteins are similar, contain the typical conserved domains of ENCO3c and are located in the nucleus. Furthermore, multiple cis-acting elements were observed in DNA-deMTases, including light responsiveness, phytohormone responsiveness, stress responsiveness, and elements related to plant growth and development. The DNA-deMTase genes are expressed in all tissues detected with certain tissue specificity. Then, we investigated the abundance of DNA-deMTase transcripts under abiotic stresses (cold, drought, salt, ABA, and MeJA). The results showed that different genes of DNA-deMTases were involved in the regulation of different abiotic stresses. In total, our findings will provide a basis for the roles of DNA-deMTase in response to abiotic stress.

Identification and virus-induced gene silencing (VIGS) analysis of methyltransferase affecting tomato (Solanum lycopersicum) fruit ripening

MAIN CONCLUSION

Six methyltransferase genes affecting tomato fruit ripening were identified through genome-wide screening, VIGS assay, and expression pattern analysis. The data provide the basis for understanding new mechanisms of methyltransferases. Fruit ripening is a critical stage for the formation of edible quality and seed maturation, which is finely modulated by kinds of factors, including genetic regulators, hormones, external signals, etc. Methyltransferases (MTases), important genetic regulators, play vital roles in plant development through epigenetic regulation, post-translational modification, or other mechanisms. However, the regulatory functions of numerous MTases except DNA methylation in fruit ripening remain limited so far. Here, six MTases, which act on different types of substrates, were identified to affect tomato fruit ripening. First, 35 MTase genes with relatively high expression at breaker (Br) stage of tomato fruit were screened from the tomato MTase gene database encompassing 421 genes totally. Thereafter, six MTase genes were identified as potential regulators of fruit ripening via virus-induced gene silencing (VIGS), including four genes with a positive regulatory role and two genes with a negative regulatory role, respectively. The expression of these six MTase genes exhibited diverse patterns during the fruit ripening process, and responded to various external ripening-related factors, including ethylene, 1-methylcyclopropene (1-MCP), temperature, and light exposure. These results help to further elaborate the biological mechanisms of MTase genes in tomato fruit ripening and enrich the understanding of the regulatory mechanisms of fruit ripening involving MTases, despite of DNA MTases.

Identification of watermelon H3K4 and H3K27 genes and their expression profiles during watermelon fruit development

BACKGROUND

The ClaH3K4s and ClaH3K27s gene families are subfamilies of the SET family, each with a highly conserved SET structure domain and a PHD structural domain. Both participate in histone protein methylation, which affects the chromosome structure and gene expression, and is essential for fruit growth and development.

METHODS AND RESULTS

In order to demonstrate the structure and expression characteristics of ClaH3K4s and ClaH3K27s in watermelon, members of the watermelon H3K4 and H3K27 gene families were identified, and their chromosomal localization, gene structure, and protein structural domains were analyzed. The phylogeny and covariance of the gene families with other species were subsequently determined, and the expression profiles were obtained by performing RNA-Seq and qRT-PCR. The watermelon genome had five H3K4 genes with 3207-8043 bp nucleotide sequence lengths and four H3K27 genes with a 1107-5499 bp nucleotide sequence. Synteny analysis revealed the close relationship between watermelon and cucumber, with the majority of members displaying a one-to-one covariance. Approximately half of the 'Hua-Jing 13 watermelon' ClaH3K4s and ClaH3K27s genes were expressed more in the late fruit development stages, while the changes were minimal for the remaining half. H3K4-2 expression was observed to be slightly greater on day 21 compared to other periods. Moreover, ClaH3K27-1 and ClaH3K27-2 were hardly expressed throughout the developing period, and ClaH3K27-4 exhibited the highest expression.

CONCLUSION

These results serve as a basis for further functional characterization of the H3K4 and H3K27 genes in the fruit development of watermelon.

Genome-wide association study of soybean (Glycine max [L.] Merr.) germplasm for dissecting the quantitative trait nucleotides and candidate genes underlying yield-related traits

Soybean (Glycine max [L.] Merr.) is one of the most significant crops in the world in terms of oil and protein. Owing to the rising demand for soybean products, there is an increasing need for improved varieties for more productive farming. However, complex correlation patterns among quantitative traits along with genetic interactions pose a challenge for soybean breeding. Association studies play an important role in the identification of accession with useful alleles by locating genomic sites associated with the phenotype in germplasm collections. In the present study, a genome-wide association study was carried out for seven agronomic and yield-related traits. A field experiment was conducted in 2015/2016 at two locations that include 155 diverse soybean germplasm. These germplasms were genotyped using SoySNP50K Illumina Infinium Bead-Chip. A total of 51 markers were identified for node number, plant height, pods per plant, seeds per plant, seed weight per plant, hundred-grain weight, and total yield using a multi-locus linear mixed model (MLMM) in FarmCPU. Among these significant SNPs, 18 were putative novel QTNs, while 33 co-localized with previously reported QTLs. A total of 2,356 genes were found in 250 kb upstream and downstream of significant SNPs, of which 17 genes were functional and the rest were hypothetical proteins. These 17 candidate genes were located in the region of 14 QTNs, of which ss715580365, ss715608427, ss715632502, and ss715620131 are novel QTNs for PH, PPP, SDPP, and TY respectively. Four candidate genes, Glyma.01g199200, Glyma.10g065700, Glyma.18g297900, and Glyma.14g009900, were identified in the vicinity of these novel QTNs, which encode lsd one like 1, Ergosterol biosynthesis ERG4/ERG24 family, HEAT repeat-containing protein, and RbcX2, respectively. Although further experimental validation of these candidate genes is required, several appear to be involved in growth and developmental processes related to the respective agronomic traits when compared with their homologs in Arabidopsis thaliana. This study supports the usefulness of association studies and provides valuable data for functional markers and investigating candidate genes within a diverse germplasm collection in future breeding programs.

Molecular cloning and functional characterization of the promoter of a novel Aspergillus flavus inducible gene (AhOMT1) from peanut

Peanut is an important oil and food legume crop grown in more than one hundred countries, but the yield and quality are often impaired by different pathogens and diseases, especially aflatoxins jeopardizing human health and causing global concerns. For better management of aflatoxin contamination, we report the cloning and characterization of a novel A. flavus inducible promoter of the O-methyltransferase gene (AhOMT1) from peanut. The AhOMT1 gene was identified as the highest inducible gene by A. flavus infection through genome-wide microarray analysis and verified by qRT-PCR analysis. AhOMT1 gene was studied in detail, and its promoter, fussed with the GUS gene, was introduced into Arabidopsis to generate homozygous transgenic lines. Expression of GUS gene was studied in transgenic plants under the infection of A. flavus. The analysis of AhOMT1 gene characterized by in silico assay, RNAseq, and qRT-PCR revealed minute expression in different organs and tissues with trace or no response to low temperature, drought, hormones, Ca2+, and bacterial stresses, but highly induced by A. flavus infection. It contains four exons encoding 297 aa predicted to transfer the methyl group of S-adenosyl-L-methionine (SAM). The promoter contains different cis-elements responsible for its expression characteristics. Functional characterization of AhOMT1P in transgenic Arabidopsis plants demonstrated highly inducible behavior only under A. flavus infection. The transgenic plants did not show GUS expression in any tissue(s) without inoculation of A. flavus spores. However, GUS activity increased significantly after inoculation of A. flavus and maintained a high level of expression after 48 hours of infection. These results provided a novel way for future management of peanut aflatoxins contamination through driving resistance genes in A. flavus inducible manner.

Genome-Wide Identification of Maize Protein Arginine Methyltransferase Genes and Functional Analysis of ZmPRMT1 Reveal Essential Roles in Arabidopsis Flowering Regulation and Abiotic Stress Tolerance

Histone methylation, as one of the important epigenetic regulatory mechanisms, plays a significant role in growth and developmental processes and stress responses of plants, via altering the methylation status or ratio of arginine and lysine residues of histone tails, which can affect the regulation of gene expression. Protein arginine methyltransferases (PRMTs) have been revealed to be responsible for histone methylation of specific arginine residues in plants, which is important for maintaining pleiotropic development and adaptation to abiotic stresses in plants. Here, for the first time, a total of eight PRMT genes in maize have been identified and characterized in this study, named as ZmPRMT1-8. According to comparative analyses of phylogenetic relationship and structural characteristics among PRMT gene family members from several representative species, all maize 8 PRMT proteins were categorized into three distinct subfamilies. Further, schematic structure and chromosome location analyses displayed evolutionarily conserved structure features and an unevenly distribution on maize chromosomes of ZmPRMT genes, respectively. The expression patterns of ZmPRMT genes in different tissues and under various abiotic stresses (heat, drought, and salt) were determined. The expression patterns of ZmPRMT genes indicated that they play a role in regulating growth and development and responses to abiotic stress. Eventually, to verify the biological roles of ZmPRMT genes, the transgenic Arabidopsis plants overexpressing ZmPRMT1 gene was constructed as a typical representative. The results demonstrated that overexpression of ZmPRMT1 can promote earlier flowering time and confer enhanced heat tolerance in transgenic Arabidopsis. Taken together, our results are the first to report the roles of ZmPRMT1 gene in regulating flowering time and resisting heat stress response in plants and will provide a vital theoretical basis for further unraveling the functional roles and epigenetic regulatory mechanism of ZmPRMT genes in maize growth, development and responses to abiotic stresses.

Uncovering the role of DNA methyltransferases in grapevine - Plasmopara viticola interaction: From genome-wide characterization to global methylation patterns

Epigenetic regulation has recently gained prominence in the field of plant-pathogen interactions, providing a deeper understanding of the molecular mechanisms associated with plant infection. In grapevine interaction with pathogens, epigenetic regulation still remains a black box. In this work, we characterized grapevine DNA methyltransferase gene family and identified nine DNA methyltransferases genes across eight grapevine chromosomes coding for 17 proteins. We also assessed the modulation of global cytosine methylation and gene expression levels of these genes with the aim of establishing a connection between DNA methylation and grapevine resistance towards downy mildew. Our results revealed that, in the incompatible interaction, an early hypomethylation, coupled with downregulation of DNMT and CMT genes occurs very early after pathogen inoculation. Additionally, the compatible interaction is characterized by a hypermethylation at 6hpi. A temporal delay is evident between the shifts in DNA methyltransferases gene expression in both compatible and incompatible interactions which in turn may be reflected in the global methylation percentage. Overall, we present the first evidence of an epigenetic regulation role in grapevine defense against P. viticola.

Exploration of the Potential Transcriptional Regulatory Mechanisms of DNA Methyltransferases and MBD Genes in Petunia Anther Development and Multi-Stress Responses

Cytosine-5 DNA methyltransferases (C5-MTases) and methyl-CpG-binding-domain (MBD) genes can be co-expressed. They directly control target gene expression by enhancing their DNA methylation levels in humans; however, the presence of this kind of cooperative relationship in plants has not been determined. A popular garden plant worldwide, petunia (Petunia hybrida) is also a model plant in molecular biology. In this study, 9 PhC5-MTase and 11 PhMBD proteins were identified in petunia, and they were categorized into four and six subgroups, respectively, on the basis of phylogenetic analyses. An expression correlation analysis was performed to explore the co-expression relationships between PhC5-MTases and PhMBDs using RNA-seq data, and 11 PhC5-MTase/PhMBD pairs preferentially expressed in anthers were identified as having the most significant correlations (Pearson’s correlation coefficients > 0.9). Remarkably, the stability levels of the PhC5-MTase and PhMBD pairs significantly decreased in different tissues and organs compared with that in anthers, and most of the selected PhC5-MTases and PhMBDs responded to the abiotic and hormonal stresses. However, highly correlated expression relationships between most pairs were not observed under different stress conditions, indicating that anther developmental processes are preferentially influenced by the co-expression of PhC5-MTases and PhMBDs. Interestingly, the nuclear localization genes PhDRM2 and PhMBD2 still had higher correlations under GA treatment conditions, implying that they play important roles in the GA-mediated development of petunia. Collectively, our study suggests a regulatory role for DNA methylation by C5-MTase and MBD genes in petunia anther maturation processes and multi-stress responses, and it provides a framework for the functional characterization of C5-MTases and MBDs in the future.

DNA methylation is involved in acclimation to iron-deficiency in rice (Oryza sativa)

Iron (Fe) is an essential micronutrient in plants, and Fe limitation significantly affects plant growth, yield and food quality. While many studies have reported the transcriptomic profile and pursue molecular mechanism in response to Fe limitation, little is known if epigenetic factors play a role in response to Fe-deficiency. In this study, whole-genome bisulfite sequencing analysis, high-throughput RNA-Seq of mRNA, small RNA and transposable element (TE) expression with root and shoot organs of rice seedlings under Fe-sufficient and Fe-deficient conditions were performed. The results showed that widespread hypermethylation, especially for the CHH context, occurred after Fe-deficiency. Integrative analysis of methylation and transcriptome revealed that the transcript abundance of Fe-deficiency-induced genes was negatively correlated with nearby TEs and positively with the 24-nucleotide siRNAs. The ability of methylation to affect the physiology and molecular response to Fe-deficiency was tested using an exogenous DNA methyltransferase inhibitor (5-azacytidine), and genetically using a mutant for domains rearranged methyltransferase 2 (DRM2), that lacks CHH methylation. Both approaches resulted in decreased growth and Fe content in rice plants. Thus, alterations in specific methylation patterns, directed by siRNAs, play an important role in acclimation of rice to Fe-deficient conditions. Furthermore, comparison with other reports suggests this may be a universal mechanism to acclimate to limited nutrient availability.

Genome-wide identification and transcriptional characterization of DNA methyltransferases conferring temperature-sensitive male sterility in wheat

Background

DNA methyltransferase (DMT) genes contribute to plant stress responses and development by de novo establishment and subsequent maintenance of DNA methylation during replication. The photoperiod and/or temperature-sensitive genic male sterile (P/TGMS) lines play an important role in hybrid seed production of wheat. However, only a few studies have reported on the effect of DMT genes on temperature-sensitive male sterility of wheat. Although DMT genes have been investigated in some plant species, the identification and analysis of DMT genes in wheat (Triticum aestivum L.) based on genome-wide levels have not been reported.

Results

In this study, a detailed overview of phylogeny of 52 wheat DMT (TaDMT) genes was presented. Homoeolog retention for TaDMT genes was significantly above the average retention rate for whole-wheat genes, indicating the functional importance of many DMT homoeologs. We found that the strikingly high number of TaDMT genes resulted mainly from the significant expansion of the TaDRM subfamily. Intriguingly, all 5 paralogs belonged to the wheat DRM subfamily, and we speculated that tandem duplications might play a crucial role in the TaDRM subfamily expansion. Through the transcriptional analysis of TaDMT genes in a TGMS line BS366 and its hybrids with the other six fertile lines under sterile and fertile conditions, we concluded that TaCMT-D2, TaMET1-B1, and TaDRM-U6 might be involved in male sterility in BS366. Furthermore, a correlation analysis showed that TaMET1-B1 might negatively regulate the expression of TaRAFTIN1A, an important gene for pollen development, so we speculated regarding an epigenetic regulatory mechanism underlying the male sterility of BS366 via the interaction between TaMET1-B1 and TaRAFTIN1A.

Conclusions

Our findings presented a detailed phylogenic overview of the DMT genes and could provide novel insights into the effects of DMT genes on TGMS wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07600-7.

Identification of virulence associated milRNAs and their bidirectional targets in Rhizoctonia solani and maize during infection

BACKGROUND

Anastomosis group 1 IA (AG1-IA) of Rhizoctonia solani is the major agent of banded leaf and sheath blight (BLSB) disease that causes severe yield loss in many worldwide crops. MicroRNAs (miRNAs) are ~ 22 nt non-coding RNAs that negatively regulate gene expression levels by mRNA degradation or translation inhibition. A better understanding of miRNA function during AG1-IA infection can expedite to elucidate the molecular mechanisms of fungi-host interactions.

RESULTS

In this study, we sequenced three small RNA libraries obtained from the mycelium of AG1-IA isolate, non-infected maize sheath and mixed maize sheath 3 days after inoculation. In total, 137 conserved and 34 novel microRNA-like small RNAs (milRNAs) were identified from the pathogen. Among these, one novel and 17 conserved milRNAs were identified as potential virulence-associated (VA) milRNAs. Subsequently, the prediction of target genes for these milRNAs was performed in both AG1-IA and maize, while functional annotation of these targets suggested a link to pathogenesis-related biological processes. Further, expression patterns of these virulence-associated milRNAs demonstrated that theyparticipate in the virulence of AG1-IA. Finally, regulation of one maize targeting gene, GRMZM2G412674 for Rhi-milRNA-9829-5p, was validated by dual-luciferase assay and identified to play a positive role in BLSB resistance in two maize mutants. These results suggest the global differentially expressed milRNAs of R. solani AG1-IA that participate in the regulation of target genes in both AG1-IA and maize to reinforce its pathogenicity.

CONCLUSIONS

Our data have provided a comprehensive overview of the VA-milRNAs of R. solani and identified that they are probably the virulence factors by directly interfered in host targeting genes. These results offer new insights on the molecular mechanisms of R.solani-maize interactions during the process of infection.

Genome-wide identification and analysis of DNA methyltransferase and demethylase gene families in Dendrobium officinale reveal their potential functions in polysaccharide accumulation

BACKGROUND

DNA methylation is a conserved and important epigenetic modification involved in the regulation of numerous biological processes, including plant development, secondary metabolism, and response to stresses. However, no information is available regarding the identification of cytosine-5 DNA methyltransferase (C5-MTase) and DNA demethylase (dMTase) genes in the orchid Dendrobium officinale.

RESULTS

In this study, we performed a genome-wide analysis of DoC5-MTase and DodMTase gene families in D. officinale. Integrated analysis of conserved motifs, gene structures and phylogenetic analysis showed that eight DoC5-MTases were divided into four subfamilies (DoCMT, DoDNMT, DoDRM, DoMET) while three DodMTases were divided into two subfamilies (DoDML3, DoROS1). Multiple cis-acting elements, especially stress-responsive and hormone-responsive ones, were found in the promoter region of DoC5-MTase and DodMTase genes. Furthermore, we investigated the expression profiles of DoC5-MTase and DodMTase in 10 different tissues, as well as their transcript abundance under abiotic stresses (cold and drought) and at the seedling stage, in protocorm-like bodies, shoots, and plantlets. Interestingly, most DoC5-MTases were downregulated whereas DodMTases were upregulated by cold stress. At the seedling stage, DoC5-MTase expression decreased as growth proceeded, but DodMTase expression increased.

CONCLUSIONS

These results provide a basis for elucidating the role of DoC5-MTase and DodMTase in secondary metabolite production and responses to abiotic stresses in D. officinale.

Genome-Wide Identification and Coexpression Network Analysis of DNA Methylation Pathway Genes and Their Differentiated Functions in Ginkgo biloba L

DNA methylation plays a vital role in diverse biological processes. DNA methyltransferases (DNMTs) genes and RNA-directed DNA methylation (RdDM)-related genes are key genes responsible for establishing and maintaining genome DNA methylation in plants. In the present study, we systematically identified nine GbDNMTs in Ginkgo biloba, including the three common families of GbMET1a/1b, GbCMT2, and GbDRMa/b/2a/2b/2c, and a fourth family—GbDNMT3—which is absent in most angiosperms. We also identified twenty RdDM-related genes, including four GbDCLs, six GbAGOs, and ten GbRDRs. Expression analysis of the genes showed the different patterns of individual genes, and 15 of 29 genes displayed expression change under five types of abiotic stress. Gene coexpression analysis and weighted gene co-expression network analysis (WGCNA) using 126 public transcriptomic datasets revealed that these genes were clustered into two groups. In group I, genes covered members from all six families which were preferentially expressed in the ovulate strobile and fruit. A gene ontology (GO) enrichment analysis of WGCNA modules indicated that group I genes were most correlated with the biological process of cell proliferation. Group II only consisted of RdDM-related genes, including GbDRMs, GbAGOs, and GbRDRs, but no GbDCLs, and these genes were specifically expressed in the cambium, suggesting that they may function in a dicer-like (DCL)-independent RdDM pathway in specific tissues. The gene module related to group II was most enriched in signal transduction, cell communication, and the response to the stimulus. These results demonstrate that gene family members could be conserved or diverged across species, and multi-member families in the same pathway may cluster into different modules to function differentially. The study provides insight into the DNA methylation genes and their possible functions in G. biloba, laying a foundation for the further study of DNA methylation in gymnosperms.

The Interplay between Toxic and Essential Metals for Their Uptake and Translocation Is Likely Governed by DNA Methylation and Histone Deacetylation in Maize

The persistent nature of lead (Pb) and cadmium (Cd) in the environment severely affects plant growth and yield. Conversely, plants acquire zinc (Zn) from the soil for their vital physiological and biochemical functions. However, the interplay and coordination between essential and toxic metals for their uptake and translocation and the putative underlying epigenetic mechanisms have not yet been investigated in maize. Here, we report that the presence of Zn facilitates the accumulation and transport of Pb and Cd in the aerial parts of the maize plants. Moreover, the Zn, Pb, and Cd interplay specifically interferes with the uptake and translocation of other divalent metals, such as calcium and magnesium. Zn, Pb, and Cd, individually and in combinations, differentially regulate the expression of DNA methyltransferases, thus alter the DNA methylation levels at the promoter of Zinc-regulated transporters, Iron-regulated transporter-like Protein (ZIP) genes to regulate their expression. Furthermore, the expression of histone deacetylases (HDACs) varies greatly in response to individual and combined metals, and HDACs expression showed a negative correlation with ZIP transporters. Our study highlights the implication of DNA methylation and histone acetylation in regulating the metal stress tolerance dynamics through Zn transporters and warns against the excessive use of Zn fertilizers in metal contaminated soils.

Bioinformatics and expression analysis of histone modification genes in grapevine predict their involvement in seed development, powdery mildew resistance, and hormonal signaling

BACKGROUND

Histone modification genes (HMs) play potential roles in plant growth and development via influencing gene expression and chromatin structure. However, limited information is available about HMs genes in grapes (Vitis vinifera L.).

RESULTS

Here, we described detailed genome-wide identification of HMs gene families in grapevine. We identified 117 HMs genes in grapevine and classified these genes into 11 subfamilies based on conserved domains and phylogenetic relationships with Arabidopsis. We described the genes in terms of their chromosomal locations and exon-intron distribution. Further, we investigated the evolutionary history, gene ontology (GO) analysis, and syntenic relationships between grapes and Arabidopsis. According to results 21% HMs genes are the result of duplication (tandem and segmental) events and all the duplicated genes have negative mode of selection. GO analysis predicted the presence of HMs proteins in cytoplasm, nucleus, and intracellular organelles. According to seed development expression profiling, many HMs grapevine genes were differentially expressed in seeded and seedless cultivars, suggesting their roles in seed development. Moreover, we checked the response of HMs genes against powdery mildew infection at different time points. Results have suggested the involvement of some genes in disease resistance regulation mechanism. Furthermore, the expression profiles of HMs genes were analyzed in response to different plant hormones (Abscisic acid, Jasmonic acid, Salicylic acid, and Ethylene) at different time points. All of the genes showed differential expression against one or more hormones.

CONCLUSION

VvHMs genes might have potential roles in grapevine including seed development, disease resistance, and hormonal signaling pathways. Our study provides first detailed genome-wide identification and expression profiling of HMs genes in grapevine.

Genome-wide identification and expression profiling of cytosine-5 DNA methyltransferases during drought and heat stress in wheat (Triticum aestivum)

DNA methylation is a potential epigenetic mechanism that regulates genome stability, development, and stress mitigation in plants. It is mediated by cytosine-5 DNA methyltransferases (C5-MTases). We identified 52 wheat C5-MTases; and based on domain structure and phylogenetics, these 52 C5-MTases were classified into four sub-families including MET, CMT, DRM and DNMT2; and were distributed on 18 chromosomes. Cis-acting regulatory elements analysis identified abiotic stress-responsive, phytohormone-responsive, development-related and light-related elements in the promoters of TaC5-MTases. We also examined the transcript abundance of TaC5-MTases in different tissues, developmental stages and under abiotic stresses. Notably, most of the TaC5-MTases (TaCMT2, TaCMT3b, TaCMT3c, TaMET1, TaDRM10, TaDNMT2) showed differential regulation of their transcript abundance during drought and heat stress. Overall, the above results provide significant insights into the expression and the probable functions of TaC5-MTases and will also expedite future research programs to explore the mechanisms of epigenetic regulation in wheat.

Systematic Analysis of the DNA Methylase and Demethylase Gene Families in Rapeseed (Brassica napus L.) and Their Expression Variations After Salt and Heat stresses

DNA methylation is a process through which methyl groups are added to the DNA molecule, thereby modifying the activity of a DNA segment without changing the sequence. Increasing evidence has shown that DNA methylation is involved in various aspects of plant growth and development via a number of key processes including genomic imprinting and repression of transposable elements. DNA methylase and demethylase are two crucial enzymes that play significant roles in dynamically maintaining genome DNA methylation status in plants. In this work, 22 DNA methylase genes and six DNA demethylase genes were identified in rapeseed (Brassica napus L.) genome. These DNA methylase and DNA demethylase genes can be classified into four (BnaCMTs, BnaMET1s, BnaDRMs and BnaDNMT2s) and three (BnaDMEs, BnaDML3s and BnaROS1s) subfamilies, respectively. Further analysis of gene structure and conserved domains showed that each sub-class is highly conserved between rapeseed and Arabidopsis. Expression analysis conducted by RNA-seq as well as qRT-PCR suggested that these DNA methylation/demethylation-related genes may be involved in the heat/salt stress responses in rapeseed. Taken together, our findings may provide valuable information for future functional characterization of these two types of epigenetic regulatory enzymes in polyploid species such as rapeseed, as well as for analyzing their evolutionary relationships within the plant kingdom.

Genome-wide identification and expression profiling of SET DOMAIN GROUP family in Dendrobium catenatum

BACKGROUND

Dendrobium catenatum, as a precious Chinese herbal medicine, is an epiphytic orchid plant, which grows on the trunks and cliffs and often faces up to diverse environmental stresses. SET DOMAIN GROUP (SDG) proteins act as histone lysine methyltransferases, which are involved in pleiotropic developmental events and stress responses through modifying chromatin structure and regulating gene transcription, but their roles in D. catenatum are unknown.

RESULTS

In this study, we identified 44 SDG proteins from D. catenatum genome. Subsequently, comprehensive analyses related to gene structure, protein domain organization, and phylogenetic relationship were performed to evaluate these D. catenatum SDG (DcSDG) proteins, along with the well-investigated homologs from the model plants Arabidopsis thaliana and Oryza sativa as well as the newly characterized 42 SDG proteins from a closely related orchid plant Phalaenopsis equestris. We showed DcSDG proteins can be grouped into eight distinct classes (I~VII and M), mostly consistent with the previous description. Based on the catalytic substrates of the reported SDG members mainly in Arabidopsis, Class I (E(z)-Like) is predicted to account for the deposition of H3K27me2/3, Class II (Ash-like) for H3K36me, Class III (Trx/ATX-like) for H3K4me2/3, Class M (ATXR3/7) for H3K4me, Class IV (Su (var)-like) for H3K27me1, Class V (Suv-like) for H3K9me, as well as class VI (S-ET) and class VII (RBCMT) for methylation of both histone and non-histone proteins. RNA-seq derived expression profiling showed that DcSDG proteins usually displayed wide but distinguished expressions in different tissues and organs. Finally, environmental stresses examination showed the expressions of DcASHR3, DcSUVR3, DcATXR4, DcATXR5b, and DcSDG49 are closely associated with drought-recovery treatment, the expression of DcSUVH5a, DcATXR5a and DcSUVR14a are significantly influenced by low temperature, and even 61% DcSDG genes are in response to heat shock.

CONCLUSIONS

This study systematically identifies and classifies SDG genes in orchid plant D. catenatum, indicates their functional divergence during the evolution, and discovers their broad roles in the developmental programs and stress responses. These results provide constructive clues for further functional investigation and epigenetic mechanism dissection of SET-containing proteins in orchids.

Genome-wide investigation and transcriptional analysis of cytosine-5 DNA methyltransferase and DNA demethylase gene families in tea plant (Camellia sinensis) under abiotic stress and withering processing

DNA methylation is a highly conserved epigenetic modification involved in many biological processes, including growth and development, stress response, and secondary metabolism. In the plant kingdom, cytosine-5 DNA methyltransferase (C5-MTase) and DNA demethylase (dMTase) genes have been identified in some plant species. However, to the best of our knowledge, no investigator has focused on the identification and analysis of C5-MTase and dMTase genes in tea plants (Camellia sinensis) based on genome-wide levels. In this study, eight CsC5-MTases and four dMTases were identified in tea plants. These CsC5-MTase genes were divided into four subfamilies, including CsMET, CsCMT, CsDRM and CsDNMT2. The CsdMTase genes can be classified into CsROS, CsDME and CsDML. Based on conserved domain analysis of these genes, the gene loss and duplication events occurred during the evolution of CsC5-MTase and CsdMTase. Furthermore, multiple cis-acting elements were observed in the CsC5-MTase and CsdMTase, including light responsiveness, phytohormone responsiveness, stress responsiveness, and plant growth and development-related elements. Then, we investigated the transcript abundance of CsC5-MTase and CsdMTase under abiotic stress (cold and drought) and withering processing (white tea and oolong tea). Notably, most CsC5-MTases, except for CsCMT1 and CsCMT2, were significantly downregulated under abiotic stress, while the transcript abundance of all four CsdMTase genes was significantly induced. Similarly, the same transcript abundance of CsC5-MTase and CsdMTase was found during withering processing of white tea and oolong tea, respectively. In total, our findings will provide a basis for the roles of CsC5-MTase and CsdMTase in response to abiotic stress and the potential functions of these two gene families in affecting tea flavor during tea withering processing.

Identification of DNA methyltransferases and demethylases in Solanum melongena L., and their transcription dynamics during fruit development and after salt and drought stresses

DNA methylation through the activity of cytosine-5-methyltransferases (C5-MTases) and DNA demethylases plays important roles in genome protection as well as in regulating gene expression during plant development and plant response to environmental stresses. In this study, we report on a genome-wide identification of six C5-MTases (SmelMET1, SmelCMT2, SmelCMT3a, SmelCMT3b, SmelDRM2, SmelDRM3) and five demethylases (SmelDemethylase_1, SmelDemethylase_2, SmelDemethylase_3, SmelDemethylase_4, SmelDemethylase_5) in eggplant. Gene structural characteristics, chromosomal localization and phylogenetic analyses are also described. The transcript profiling of both C5-MTases and demethylases was assessed at three stages of fruit development in three eggplant commercial F1 hybrids: i.e. 'Clara', 'Nite Lady' and 'Bella Roma', representative of the eggplant berry phenotypic variation. The trend of activation of C5-MTases and demethylase genes varied in function of the stage of fruit development and was genotype dependent. The transcription pattern of C5MTAses and demethylases was also assessed in leaves of the F1 hybrid 'Nite Lady' subjected to salt and drought stresses. A marked up-regulation and down-regulation of some C5-MTases and demethylases was detected, while others did not vary in their expression profile. Our results suggest a role for both C5-MTases and demethylases during fruit development, as well as in response to abiotic stresses in eggplant, and provide a starting framework for supporting future epigenetic studies in the species.

Natural variation of DNA methylation and gene expression may determine local adaptations of Scots pine populations

Long-lived conifers are vulnerable to climate change because classical evolutionary processes are slow in developing adaptive responses. Therefore, the capacity of a genotype to adopt different phenotypes is important. Gene expression is the primary mechanism that converts genome-encoded information into phenotypes, and DNA methylation is employed in the epigenetic regulation of gene expression. We investigated variations in global DNA methylation and gene expression between three Scots pine (Pinus sylvestris L.) populations located in northern and southern Finland using mature seeds. Gene expression levels were studied in six DNA methyltransferase (DNMT) genes, which were characterized in this study, and in 19 circadian clock genes regulating adaptive traits. In embryos, expression diversity was found for three DNMT genes, which maintain DNA methylation. The expression of two DNMT genes was strongly correlated with climate variables, which suggests a role for DNA methylation in local adaptation. For adaptation-related genes, expression levels showed between-population variation in 11 genes in megagametophytes and in eight genes in embryos, and many of these genes were linked to climate factors. Altogether, our results suggest that differential DNA methylation and gene expression contribute to local adaptation in Scots pine populations and may enhance the fitness of trees under rapidly changing climatic conditions.

Identification and characterization of the cytosine-5 DNA methyltransferase gene family in Salvia miltiorrhiza

Cytosine DNA methylation is highly conserved epigenetic modification involved in a wide range of biological processes in eukaryotes. It was established and maintained by cytosine-5 DNA methyltransferases (C5-MTases) in plants. Through genome-wide identification, eight putative SmC5-MTase genes were identified from the genome of Salvia miltiorrhiza, a well-known traditional Chinese medicine material and an emerging model medicinal plant. Based on conserved domains and phylogenetic analysis, eight SmC5-MTase genes were divided into four subfamilies, including MET, CMT, DRM and DNMT2. Genome-wide comparative analysis of the C5-MTase gene family in S. miltiorrhiza and Arabidopsis thaliana, including gene structure, sequence features, sequence alignment and conserved motifs, was carried out. The results showed conservation and divergence of the members of each subfamily in plants. The length of SmC5-MTase open reading frames ranges widely from 1,152 (SmDNMT2) to 5,034 bp (SmMET1). The intron number of SmC5-MTases varies between 7 (SmDRM1) and 20 (SmCMT1 and SmCMT2b). These features were similar to their counterparts from Arabidopsis. Sequence alignment and conserved motif analysis showed the existence of highly conserved and subfamily-specific motifs in the C5-MTases analyzed. Differential transcript abundance was detected for SmC5-MTases, implying genome-wide variance of DNA methylation in different organs and tissues. Transcriptome-wide analysis showed that the transcript levels of all SmC5-MTase genes was slightly changed under yeast extract and methyl jasmonate treatments. Six SmC5-MTases, including SmMET1, SmCMT1, SmCMT2a, SmCMT2b, SmCMT3 and SmDRM1, were salicylic acid-responsive, suggesting the involvement of SmC5-MTases in salicylic acid-dependent immunity. These results provide useful information for demonstrating the role of DNA methylation in bioactive compound biosynthesis and Dao-di herb formation in medicinal plants.

Temporal and spatial expression of genes involved in DNA methylation during reproductive development of sexual and apomictic Eragrostis curvula

Recent reports in model plant species have highlighted a role for DNA methylation pathways in the regulation of the somatic-to-reproductive transition in the ovule, suggesting that apomixis (asexual reproduction through seeds) likely relies on RdDM downregulation. Our aim was therefore to explore this hypothesis by characterizing genes involved in DNA methylation in the apomictic grass Eragrostis curvula. We explored floral transcriptomes to identify homologs of three candidate genes, for which mutations in Arabidopsis and maize mimic apomixis (AtAGO9/ZmAGO104, AtCMT3/ZmDMT102/ZmDMT105, and AtDDM1/ZmCHR106), and compared both their spatial and temporal expression patterns during reproduction in sexual and apomictic genotypes. Quantitative expression analyses revealed contrasting expression patterns for the three genes in apomictic vs sexual plants. In situ hybridization corroborated these results for two candidates, EcAGO104 and EcDMT102, and revealed an unexpected ectopic pattern for the AGO gene during germ line differentiation in apomicts. Although our data partially support previous results obtained in sexual plant models, they suggest that rather than an RdDM breakdown in the ovule, altered localization of AtAGO9/ZmAGO104 expression is required for achieving diplospory in E. curvula. The differences in the RdDM machinery acquired during plant evolution might have promoted the emergence of the numerous apomictic paths observed in plants.

Temporal and spatial expression of genes involved in DNA methylation during reproductive development of sexual and apomictic Eragrostis curvula

Recent reports in model plant species have highlighted a role for DNA methylation pathways in the regulation of the somatic-to-reproductive transition in the ovule, suggesting that apomixis (asexual reproduction through seeds) likely relies on RdDM downregulation. Our aim was therefore to explore this hypothesis by characterizing genes involved in DNA methylation in the apomictic grass Eragrostis curvula. We explored floral transcriptomes to identify homologs of three candidate genes, for which mutations in Arabidopsis and maize mimic apomixis (AtAGO9/ZmAGO104, AtCMT3/ZmDMT102/ZmDMT105, and AtDDM1/ZmCHR106), and compared both their spatial and temporal expression patterns during reproduction in sexual and apomictic genotypes. Quantitative expression analyses revealed contrasting expression patterns for the three genes in apomictic vs sexual plants. In situ hybridization corroborated these results for two candidates, EcAGO104 and EcDMT102, and revealed an unexpected ectopic pattern for the AGO gene during germ line differentiation in apomicts. Although our data partially support previous results obtained in sexual plant models, they suggest that rather than an RdDM breakdown in the ovule, altered localization of AtAGO9/ZmAGO104 expression is required for achieving diplospory in E. curvula. The differences in the RdDM machinery acquired during plant evolution might have promoted the emergence of the numerous apomictic paths observed in plants.

DNA Methylation Signatures of the Plant Chromomethyltransferases

DNA methylation in plants is traditionally partitioned into CG, CHG and CHH contexts (with H any nucleotide but G). By investigating DNA methylation patterns in trinucleotide contexts in four angiosperm species, we show that such a representation hides spatial and functional partitioning of different methylation pathways and is incomplete. CG methylation (mCG) is largely context-independent whereas, at CHG motifs, there is under-representation of mCCG in pericentric regions of A. thaliana and tomato and throughout the chromosomes of maize and rice. In A. thaliana the biased representation of mCCG in heterochromatin is related to specificities of H3K9 methyltransferase SUVH family members. At CHH motifs there is an over-representation of different variant forms of mCHH that, similarly to mCCG hypomethylation, is partitioned into the pericentric regions of the two dicots but dispersed in the monocot chromosomes. The over-represented mCHH motifs in A. thaliana associate with specific types of transposon including both class I and II elements. At mCHH the contextual bias is due to the involvement of various chromomethyltransferases whereas the context-independent CHH methylation in A. thaliana and tomato is mediated by the RNA-directed DNA methylation process that is most active in the gene-rich euchromatin. This analysis therefore reveals that the sequence context of the methylome of plant genomes is informative about the mechanisms associated with maintenance of methylation and the overlying chromatin structure.

Dense cytosine DNA methylation (mC) in eukaryotes is associated with closed chromatin and gene silencing. In plants it is well known that the sequence context of the mC (either mCG, mCHG or mCHH) provides a clue as to which of several mechanisms is involved but now, based on detailed analyses of the DNA methylome in wild type and mutants of four plant species, we reveal that there is additional information in the mC sequence context. Low mCCG and over-representation of mCAA and mCTA or mCAT in A. thaliana and tomato differentiates regions of the chromosomes near the centromere where methylation is dominated by chromomethyltransferases from the chromosome arms in which mCHH is context-independent and predominantly RNA-directed. Rice and maize have similar sequence context-dependent DNA methylation but the corresponding chromosome domains are not spatially separate as in the dicots. The discovery of the subcomponents of plant methylomes based on sequence context will allow greater resolution in past and future analyses of plant methylomes.

Variation of DNA methylation patterns associated with gene expression in rice (Oryza sativa) exposed to cadmium

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.