Methylation of coding region alone inhibits gene expression in plant protoplasts

Genome-Wide Identification and Evolution Analysis of the CYP76 Subfamily in Rice (Oryza sativa)

The CYP76 subfamily, a member of the CYP superfamily, plays crucial roles in the biosynthesis of phytohormones in plants, involving biosynthesis of secondary metabolites, hormone signaling, and response to environmental stresses. Here, we conducted a genome-wide analysis of the CYP76 subfamily in seven AA genome species: Oryza sativa ssp. japonica, Oryza sativa ssp. indica, Oryza rufipogon, Oryza glaberrima, Oryza meridionalis, Oryza barthii, and Oryza glumaepatula. These were identified and classified into three groups, and it was found that Group 1 contained the largest number of members. Analysis of cis-acting elements revealed a large number of elements related to jasmonic acid and light response. The gene duplication analysis revealed that the CYP76 subfamily expanded mainly in SD/WGD and tandem forms and underwent strong purifying selection during evolution. Expression pattern analysis of OsCYP76s in various developmental stages revealed that the majority of OsCYP76s exhibit relatively restricted expression patterns in leaves and roots. We further analyzed the expression of CYP76s in O. sativa, japonica, and O. sativa, indica under cold, flooding, drought, and salt abiotic stresses by qRT-PCR. We found that OsCYP76-11 showed a huge increase in relative expression after drought and salt stresses. After flooding stress, OsiCYP76-4 showed a greater increase in expression compared to other genes. CYP76 in japonica and indica showed different response patterns to the same abiotic stresses, revealing functional divergence in the gene family during evolution; these may be the key genes responsible for the differences in tolerance to indica japonica. Our results provide valuable insights into the functional diversity and evolutionary history of the CYP76 subfamily and pave the way for the development of new strategies for improving stress tolerance and agronomic traits in rice.

Single-cell multi-omics sequencing and its applications in studying the nervous system

Single-cell sequencing has become one of the most powerful and popular techniques in dissecting molecular heterogeneity and modeling the cellular architecture of a biological system. During the past twenty years, the throughput of single-cell sequencing has increased from hundreds of cells to over tens of thousands of cells in parallel. Moreover, this technology has been developed from sequencing transcriptome to measure different omics such as DNA methylome, chromatin accessibility, and so on. Currently, multi-omics which can analyze different omics in the same cell is rapidly advancing. This work advances the study of many biosystems, including the nervous system. Here, we review current single-cell multi-omics sequencing techniques and describe how they improve our understanding of the nervous system. Finally, we discuss the open scientific questions in neural research that may be answered through further improvement of single-cell multi-omics sequencing technology.

Mitochondrial DHODH regulates hypoxia-inducible factor 1 expression in OTSCC

Oral tongue squamous cell carcinoma (OTSCC) was one of the most hypoxic tumors with unfavorable outcomes. Hypoxia-inducible factor-1 (HIF-1) signaling was associated with cancer proliferation, lymph node metastasis, angiogenesis and poor prognosis of OTSCC. Dihydroorotate dehydrogenase (DHODH) catalyzed the rate-limiting step in the de novo pyrimidine biosynthesis. The aim of the study was to explore the biological function of DHODH and investigate whether DHODH regulated HIF-1 signaling in OTSCC. Proliferation, migration and anoikis resistance were used to determine the function of DHODH. Western blot and luciferase activity assays were used to determine the regulatory role of DHODH on HIF-1. We found that increased DHODH expression was associated with advanced tumor stage and poorly differentiated tumor in head and neck cancer patients in The Cancer Genome Atlas (TCGA). DHODH enhanced the proliferation and aggressiveness of OTSCC. Moreover, DHODH prompted tumor growth and metastasis in vivo. DHODH promoted transcription, protein stability, and transactivation activity of HIF1A. DHODH-induced HIF1A upregulation in OTSCC can be reversed by reactive oxygen species (ROS) scavenger, indicating that DHODH enhanced HIF1A expression via ROS production. DHODH inhibitor suppressed DHODH-mediated ROS generation and HIF1A upregulation. Targeting DHODH using clinically available inhibitor, atovaquone, might provide a new strategy to treat OTSCC.

Genome-Wide DNA Methylation Profile in Jejunum Reveals the Potential Genes Associated With Paratuberculosis in Dairy Cattle

Paratuberculosis in cattle causes substantial economic losses to the dairy industry. Exploring functional genes and corresponding regulatory pathways related to resistance or susceptibility to paratuberculosis is essential to the breeding of disease resistance in cattle. Co-analysis of genome-wide DNA methylation and transcriptome profiles is a critically important approach to understand potential regulatory mechanism underlying the development of diseases. In this study, we characterized the profiles of DNA methylation of jejunum from nine Holstein cows in clinical, subclinical, and healthy groups using whole-genome bisulfite sequencing (WGBS). The average methylation level in functional regions was 29.95% in the promoter, 29.65% in the 5’ untranslated region (UTR), 68.24% in exons, 71.55% in introns, and 72.81% in the 3’ UTR. A total of 3,911, 4,336, and 4,094 differentially methylated genes (DMGs) were detected in clinical vs. subclinical, clinical vs. healthy, and subclinical vs. healthy comparative group, respectively. Gene ontology (GO) and analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that these DMGs were significantly enriched in specific biological processes related to immune response, such as Th1 and Th2 cell differentiation, wnt, TNF, MAPK, ECM-receptor interaction, cellular senescence, calcium, and chemokine signaling pathways (q value <0.05). The integration of information about DMGs, differentially expressed genes (DEGs), and biological functions suggested nine genes CALCRL, TNC, GATA4, CD44, TGM3, CXCL9, CXCL10, PPARG, and NFATC1 as promising candidates related to resistance/susceptibility to Mycobacterium avium subspecies paratuberculosis (MAP). This study reports on the high-resolution DNA methylation landscapes of the jejunum methylome across three conditions (clinical, subclinical, and healthy) in dairy cows. Our investigations integrated different sources of information about DMGs, DEGs, and pathways, enabling us to find nine functional genes that might have potential application in resisting paratuberculosis in dairy cattle.

Identification, methylation profiling, and expression analysis of stress-responsive cytochrome P450 genes in rice under abiotic and phytohormones stresses

The cytochrome P450 (CYP) is a large and complex eukaryotic gene superfamily with enzymatic activities involved in several physiological and regulatory processes. As an objective, an in-silico genome-wide DNA methylation (5mC) analysis was performed in rice (Oryza sativa cv. Zhonghua11), and the epigenetic role of CYPs in two abiotic stresses was observed. Being a stable representative mark, DNA-methylation alters the gene expression under stressful environmental conditions. Rice plants under salinity and drought stresses were analyzed through MeDIP-chip hybridization, and 14 unique genes of the CYP family were identified in the rice genome with varying degrees of methylation. The gene structure, promoter sequences, and phylogenetic analysis were performed. Furthermore, the responses of CYPs to various abiotic stresses, including salinity, drought, and cold were revealed. Similarly, the expression profile of potential CYPs was also investigated under various phytohormone stresses, which revealed the potential involvement of CYPs to hormone regulations. Overall, the current study provides evidence for CYP's stress regulation and fundamental for further characterization and understanding their epigenetic roles in gene expression regulation and environmental stress regulation in higher plants.

The genome of the Cauliflower mosaic virus, a plant pararetrovirus, is highly methylated in the nucleus

Cytosine methylation is an important defense against invasive DNAs. Here, cytosine methylation profiles of a plant pararetrovirus, Cauliflower mosaic virus (CaMV), were investigated. Nuclear CaMV DNA is highly methylated throughout the genome including at transcription regulatory regions, but the virion DNA is unmethylated. In vitro CG methylation of the viral 35S promoter reduces transcription from the downstream gene. Although nuclear CaMV DNA is highly methylated, its transcripts are accumulated in the nucleus. The data suggest that a small population of unmethylated viral genomes produced through reverse transcription are constantly delivered back to the nucleus. Small RNA profiles suggest that methylation of the CaMV DNA may be due to de novo methylation through 21-, 22-, and 24-nt small RNAs with adenines at their 5' terminus.

Heterologous Expression of SvMBD5 from Salix viminalis L. Promotes Flowering in Arabidopsis thaliana L

Methyl-CpG-binding domain (MBD) proteins have diverse molecular and biological functions in plants. Most studies of MBD proteins in plants have focused on the model plant Arabidopsis thaliana L. Here we cloned SvMBD5 from the willow Salix viminalis L. by reverse transcription-polymerase chain reaction (RT-PCR) and analyzed the structure of SvMBD5 and its evolutionary relationships with proteins in other species. The coding sequence of SvMBD5 is 645 bp long, encoding a 214 amino acid protein with a methyl-CpG-binding domain. SvMBD5 belongs to the same subfamily as AtMBD5 and AtMBD6 from Arabidopsis. Subcellular localization analysis showed that SvMBD5 is only expressed in the nucleus. We transformed Arabidopsis plants with a 35S::SvMBD5 expression construct to examine SvMBD5 function. The Arabidopsis SvMBD5-expressing line flowered earlier than the wild type. In the transgenic plants, the expression of FLOWERING LOCUS T and CONSTANS significantly increased, while the expression of FLOWERING LOCUS C greatly decreased. In addition, heterologously expressing SvMBD5 in Arabidopsis significantly inhibited the establishment and maintenance of methylation of CHROMOMETHYLASE 3 and METHYLTRANSFERASE 1, as well as their expression, and significantly increased the expression of the demethylation-related genes REPRESSOR OF SILENCING1 and DEMETER-LIKE PROTEIN3. Our findings suggest that SvMBD5 participates in the flowering process by regulating the methylation levels of flowering genes, laying the foundation for further studying the role of SvMBD5 in regulating DNA demethylation.

Spontaneous reactivation of a site-specifically placed transgene independent of copy number or DNA methylation

As millions of seeds are produced from a breeding line, the long-term stability of transgene expression is vital for commercial-scale production of seeds with transgenic traits. Transgenes can be silenced by epigenetic mechanisms, but reactivation of expression can occur as a result of treatment with chromatin modification inhibitors such as 5-azacytidine, from stress such as heat or UV-B, or in mutants that have acquired a defect in gene silencing. Previously, we targeted a gfp reporter gene into the tobacco (Nicotiana tabacum) genome by site-specific recombination but still found some silenced lines among independent integration events. One such line also had a second random copy and both copies showed DNA hypermethylation. To test whether removing the second copy would reactivate gfp expression, two T1 plants were backcrossed to the wild type. Whereas the silenced status was maintained in the progenies from one backcross, spontaneous partial reactivation of gfp expression was found among progenies from a second backcross. However, this reactivation did not correlate with loss of the second random copy or with a significant change in the pattern or amount of DNA hypermethylation. This finding supports the suggestion that gene reactivation does not necessarily involve loss of DNA homology or methylation.

Analysis of γ-Tocopherol methyl transferase3 promoter activity and study of methylation patterns of the promoter and its gene body

Soybeans are known for its good source of protein (40%), oil (20%) and also serve as a source of nutraceutical compounds including tocopherols (toc). To know the molecular basis of differential α-toc accumulation in two contrasting soybean genotypes: DS74 (low α-toc - 1.36 μg/g and total-toc -29.72 μg/g) and Bragg (high α-toc - 10.48 μg/g and total-toc 178.91 μg/g), the analysis of γ-TMT3 promoter activity and its methylation patterns were carried out. The sequencing results revealed nucleotide variation between Bragg:γ-TMT3-P and DS74:γ-TMT3-P, however none of the variations were found in core-promoter region or in cis-elements. The histochemical GUS assay revealed higher promoter activity of Bragg:γ-TMT3-P than that of DS74:γ-TMT3-P and correlated with significantly higher and lower (P < 0.05) expression of γ-TMT3 gene respectively. To know the molecular basis of differential accumulation of α-toc in these contrasting soybean genotypes, the DNA methylation pattern of γ-TMT3 gene body and its promoter was studied in both varieties. The results showed higher percentage (62.5%) of methylation in DS74:γ-TMT3-P than in Bragg:γ-TMT3-P (50%). Out of all the methylation sites in the promoter region, one of methylation site was found at CAAT box (-190 bp) of DS74:γ-TMT3-P. Further gene body methylation patterns revealed lowest % (40%) of CG methylation in DS74:γ-TMT3 gene as compared to Bragg:γ-TMT3 (64.2%). Thus our study revealed that, expression of γ-TMT3 gene was influenced by its promoter activity and methylation patterns in cis-elements of γ-TMT3 promoter and gene body. This study will help us to understand the possible role of methylation and promoter activity in determining the α-toc content in soybean seeds.

Identification of the expressome by machine learning on omics data

Accurate annotation of plant genomes remains complex due to the presence of many pseudogenes arising from whole-genome duplication-generated redundancy or the capture and movement of gene fragments by transposable elements. Machine learning on genome-wide epigenetic marks, informed by transcriptomic and proteomic training data, could be used to improve annotations through classification of all putative protein-coding genes as either constitutively silent or able to be expressed. Expressed genes were subclassified as able to express both mRNAs and proteins or only RNAs, and CG gene body methylation was associated only with the former subclass. More than 60,000 protein-coding genes have been annotated in the reference genome of maize inbred B73. About two-thirds of these genes are transcribed and are designated the filtered gene set (FGS). Classification of genes by our trained random forest algorithm was accurate and relied only on histone modifications or DNA methylation patterns within the gene body; promoter methylation was unimportant. Other inbred lines are known to transcribe significantly different sets of genes, indicating that the FGS is specific to B73. We accurately classified the sets of transcribed genes in additional inbred lines, arising from inbred-specific DNA methylation patterns. This approach highlights the potential of using chromatin information to improve annotations of functional genes.

5-Azacytidine promotes shoot regeneration during Agrobacterium-mediated soybean transformation

Agrobacterium-mediated soybean transformation has been greatly improved in recent years, however the transformation efficiency is still low and highly genotype-dependent when compared to other species. Here, we characterized seventeen soybean genotypes based on their genetic transformation efficiencies, i.e., high and low, during Agrobacterium-mediated transformation. To reveal the molecular basis of this transformation difference, we constructed a highly efficient transient transgene expression system using soybean cotyledon protoplasts and then assess the methylation levels of promoter and coding regions of an EYFP (enhanced yellow fluorescent protein) gene introduced into the protoplast cultures of various soybean genotypes using BSP (bisulfite sequencing PCR). Increased methylation was found to be associated with the considerably decreased transfection efficiency (as percentage of EYFP fluorescent protoplasts) in low-efficacy genotypes as compared with those in high-efficacy on three DAT (day after transfection). 5-Azacytidine (5-Azac), a demethylating reagent commonly applied in epigenetic researches, significantly improved the transient transfection efficiency and transgene expression level in low-efficiency genotypes. Furthermore, the shoot regeneration efficiency in low-efficiency genotypes was substantially increased by 5-Azac treatment in an Agrobacterium-mediated soybean transformation system. Taken together, we concluded that lower methylation level in transgene contributed to enhanced shoot regeneration in Agrobacterium-mediated soybean transformation.

Induction of Transgene Suppression in Plants via External Application of Synthetic dsRNA

Recent investigations show that exogenously applied small interfering RNAs (siRNA) and long double-stranded RNA (dsRNA) precursors can be taken up and translocated in plants to induce RNA interference (RNAi) in the plant or in its fungal pathogen. The question of whether genes in the plant genome can undergo suppression as a result of exogenous RNA application on plant surface is almost unexplored. This study analyzed whether it is possible to influence transcript levels of transgenes, as more prone sequences to silencing, in Arabidopsis genome by direct exogenous application of target long dsRNAs. The data revealed that in vitro synthesized dsRNAs designed to target the gene coding regions of enhanced green fluorescent protein (EGFP) or neomycin phosphotransferase II (NPTII) suppressed their transcript levels in Arabidopsis. The fact that, simple exogenous application of polynucleotides can affect mRNA levels of plant transgenes, opens new opportunities for the development of new scientific techniques and crop improvement strategies.

Physiological, Biochemical, Epigenetic and Molecular Analyses of Wheat (Triticum aestivum) Genotypes with Contrasting Salt Tolerance

Abiotic stress exerts significant impact on plant's growth, development, and productivity. Productivity of crop plants under salt stress is lagging behind because of our limited knowledge about physiological, biochemical, epigenetic, and molecular mechanisms of salt tolerance in plants. This study aimed to investigate physio-biochemical, molecular indices and defense responses of selected wheat cultivars to identify the most contrasting salt-responsive genotypes and the mechanisms associated with their differential responses. Physio-biochemical traits specifically membrane stability index, antioxidant potential, osmoprotectants and chlorophyll contents, measured at vegetative stage, were used for multivariate analysis to identify the most contrasting genotypes. Genetic and epigenetic analyses indicated the possible mechanisms associated with differential response of the wheat genotypes under salt stress. Better antioxidant potential, membrane stability, increased accumulation of osmolytes/phytophenolics, and higher K+/Na+ ratio under 200 mM NaCl stress identified Kharchia-65 to be the most salt-tolerant cultivar. By contrast, increased MDA level, reduced soluble sugar, proline, total chlorophyll, total phenolics contents, and lower antioxidant potential in HD-2329 marked it to be sensitive to the stress. Genetic and bioinformatics analyses of HKT1;4 of contrasting genotypes (Kharchia-65 and HD-2329) revealed deletions, transitions, and transversions resulting into altered structure, loss of conserved motifs (Ser-Gly-Gly-Gly and Gly-Arg) and function in salt-sensitive (HD-2329) genotype. Expression analysis of HKTs rationalized the observed responses. Epigenetic variations in cytosine methylation explained tissue- and genotype-specific differential expression of HKT2;1 and HKT2;3.

Salt-Induced Tissue-Specific Cytosine Methylation Downregulates Expression of HKT Genes in Contrasting Wheat (Triticum aestivum L.) Genotypes

Plants have evolved several strategies, including regulation of genes through epigenetic modifications, to cope with environmental stresses. DNA methylation is dynamically regulated through the methylation and demethylation of cytosine in response to environmental perturbations. High-affinity potassium transporters (HKTs) have accounted for the homeostasis of sodium and potassium ions in plants under salt stress. Wheat (Triticum aestivum L.) is sensitive to soil salinity, which impedes its growth and development, resulting in decreased productivity. The differential expression of HKTs has been reported to confer tolerance to salt stress in plants. In this study, we investigated variations in cytosine methylation and their effects on the expression of HKT genes in contrasting wheat genotypes under salt stress. We observed a genotype- and tissue-specific increase in cytosine methylation induced by NaCl stress that downregulated the expression of TaHKT2;1 and TaHKT2;3 in the shoot and root tissues of Kharchia-65, thereby contributing to its improved salt-tolerance ability. Although TaHKT1;4 was expressed only in roots and was downregulated under the stress in salt-tolerant genotypes, it was not regulated through variations in cytosine methylation. Thus, understanding epigenetic regulation and the function of HKTs would enable an improvement in salt tolerance and the development of salt-tolerant crops.

A Comparative Analysis of 5-Azacytidine- and Zebularine-Induced DNA Demethylation

The nonmethylable cytosine analogs, 5-azacytidine and zebularine, are widely used to inhibit DNA methyltransferase activity and reduce genomic DNA methylation. In this study, whole-genome bisulfite sequencing is used to construct maps of DNA methylation with single base pair resolution in Arabidopsis thaliana seedlings treated with each demethylating agent. We find that both inhibitor treatments result in nearly indistinguishable patterns of genome-wide DNA methylation and that 5-azacytidine had a slightly greater demethylating effect at higher concentrations across the genome. Transcriptome analyses revealed a substantial number of upregulated genes, with an overrepresentation of transposable element genes, in particular CACTA-like elements. This demonstrates that chemical demethylating agents have a disproportionately large effect on loci that are otherwise silenced by DNA methylation.

Comparative in silico analysis of SSRs in coding regions of high confidence predicted genes in Norway spruce (Picea abies) and Loblolly pine (Pinus taeda)

Background

Microsatellites or simple sequence repeats (SSRs) are DNA sequences consisting of 1–6 bp tandem repeat motifs present in the genome. SSRs are considered to be one of the most powerful tools in genetic studies. We carried out a comparative study of perfect SSR loci belonging to class I (≥20) and class II (≥12 and <20 bp) types located in coding regions of high confidence genes in Picea abies and Pinus taeda. SSRLocator was used to retrieve SSRs from the full length CDS of predicted genes in both species.

Results

Trimers were the most abundant motifs in class I followed by hexamers in Picea abies, while trimers and hexamers were equally abundant in Pinus taeda class I SSRs. Hexamers were most frequent within class II SSRs followed by trimers, in both species. Although the frequency of genes containing SSRs was slightly higher in Pinus taeda, SSR counts per Mbp for class I was similar in both species (P-value = 0.22); while for class II SSRs, it was significantly higher in Picea abies (P-value = 0.00009). AT-rich motifs were higher in abundance than the GC-rich motifs, within class II SSRs in both the species (P-values = 10⁻⁹ and 0). With reference to class I SSRs, AT-rich and GC-rich motifs were detected with equal frequency in Pinus taeda (P-value = 0.24); while in Picea abies, GC-rich motifs were detected with higher frequency than the AT-rich motifs (P-value = 0.0005).

Conclusions

Our study gives a comparative overview of the genome SSRs composition based on high confidence genes in the two recently sequenced and economically important conifers and, also provides information on functional molecular markers that can be applied in genetic studies in Pinus and Picea species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0304-y) contains supplementary material, which is available to authorized users.

Divergent DNA methylation patterns associated with gene expression in rice cultivars with contrasting drought and salinity stress response

DNA methylation is an epigenetic mechanism that play an important role in gene regulation in response to environmental conditions. The understanding of DNA methylation at the whole genome level can provide insights into the regulatory mechanisms underlying abiotic stress response/adaptation. We report DNA methylation patterns and their influence on transcription in three rice (Oryza sativa) cultivars (IR64, stress-sensitive; Nagina 22, drought-tolerant; Pokkali, salinity-tolerant) via an integrated analysis of whole genome bisulphite sequencing and RNA sequencing. We discovered extensive DNA methylation at single-base resolution in rice cultivars, identified the sequence context and extent of methylation at each site. Overall, methylation levels were significantly different in the three rice cultivars. Numerous differentially methylated regions (DMRs) among different cultivars were identified and many of which were associated with differential expression of genes important for abiotic stress response. Transposon-associated DMRs were found coupled to the transcript abundance of nearby protein-coding gene(s). Small RNA (smRNA) abundance was found to be positively correlated with hypermethylated regions. These results provide insights into interplay among DNA methylation, gene expression and smRNA abundance, and suggest a role in abiotic stress adaptation in rice.

AtROS1 overexpression provides evidence for epigenetic regulation of genes encoding enzymes of flavonoid biosynthesis and antioxidant pathways during salt stress in transgenic tobacco

In plants, epigenetic changes have been identified as regulators of developmental events during normal growth as well as environmental stress exposures. Flavonoid biosynthetic and antioxidant pathways play a significant role in plant defence during their exposure to environmental cues. The aim of this study was to unravel whether genes encoding enzymes of flavonoid biosynthetic and antioxidant pathways are under epigenetic regulation, particularly DNA methylation, during salt stress. For this, a repressor of silencing from Arabidopsis, AtROS1, was overexpressed in transgenic tobacco. Generated transgenics were evaluated to examine the influence of AtROS1 on methylation status of promoters as well as on coding regions of genes encoding enzymes of flavonoids biosynthesis and antioxidant pathways. Overexpression of AtROS1 increases the demethylation levels of both promoters as well as coding regions of genes encoding chalcone synthase, chalcone isomerase, flavanone 3-hydroxylase, flavonol synthase, dihydroflavonol 4-reductase, and anthocyanidin synthase of the flavonoid biosynthetic pathway, and glutathione S-transferase, ascorbate peroxidase, glutathione peroxidase, and glutathione reductase of the antioxidant pathway during control conditions. The level of demethylation was further increased at promoters as well as coding regions of these genes during salt-stress conditions. Transgenic tobacco overexpressing AtROS1 showed tolerance to salt stress that could have been due to the higher expression levels of the genes encoding enzymes of the flavonoid biosynthetic and antioxidant pathways. This is the first comprehensive study documenting the epigenetic regulation of flavonoid biosynthetic and antioxidant pathways during salt-stress exposure of plants.

Whole-genome bisulfite sequencing of multiple individuals reveals complementary roles of promoter and gene body methylation in transcriptional regulation

BACKGROUND

DNA methylation is an important type of epigenetic modification involved in gene regulation. Although strong DNA methylation at promoters is widely recognized to be associated with transcriptional repression, many aspects of DNA methylation remain not fully understood, including the quantitative relationships between DNA methylation and expression levels, and the individual roles of promoter and gene body methylation.

RESULTS

Here we present an integrated analysis of whole-genome bisulfite sequencing and RNA sequencing data from human samples and cell lines. We find that while promoter methylation inversely correlates with gene expression as generally observed, the repressive effect is clear only on genes with a very high DNA methylation level. By means of statistical modeling, we find that DNA methylation is indicative of the expression class of a gene in general, but gene body methylation is a better indicator than promoter methylation. These findings are general in that a model constructed from a sample or cell line could accurately fit the unseen data from another. We further find that promoter and gene body methylation have minimal redundancy, and either one is sufficient to signify low expression. Finally, we obtain increased modeling power by integrating histone modification data with the DNA methylation data, showing that neither type of information fully subsumes the other.

CONCLUSION

Our results suggest that DNA methylation outside promoters also plays critical roles in gene regulation. Future studies on gene regulatory mechanisms and disease-associated differential methylation should pay more attention to DNA methylation at gene bodies and other non-promoter regions.

Whole-genome bisulfite sequencing of multiple individuals reveals complementary roles of promoter and gene body methylation in transcriptional regulation

Background

DNA methylation is an important type of epigenetic modification involved in gene regulation. Although strong DNA methylation at promoters is widely recognized to be associated with transcriptional repression, many aspects of DNA methylation remain not fully understood, including the quantitative relationships between DNA methylation and expression levels, and the individual roles of promoter and gene body methylation.

Results

Here we present an integrated analysis of whole-genome bisulfite sequencing and RNA sequencing data from human samples and cell lines. We find that while promoter methylation inversely correlates with gene expression as generally observed, the repressive effect is clear only on genes with a very high DNA methylation level. By means of statistical modeling, we find that DNA methylation is indicative of the expression class of a gene in general, but gene body methylation is a better indicator than promoter methylation. These findings are general in that a model constructed from a sample or cell line could accurately fit the unseen data from another. We further find that promoter and gene body methylation have minimal redundancy, and either one is sufficient to signify low expression. Finally, we obtain increased modeling power by integrating histone modification data with the DNA methylation data, showing that neither type of information fully subsumes the other.

Conclusion

Our results suggest that DNA methylation outside promoters also plays critical roles in gene regulation. Future studies on gene regulatory mechanisms and disease-associated differential methylation should pay more attention to DNA methylation at gene bodies and other non-promoter regions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0408-0) contains supplementary material, which is available to authorized users.

Genome-wide mapping of DNA methylation in the human malaria parasite Plasmodium falciparum

Cytosine DNA methylation is an epigenetic mark in most eukaryotic cells that regulates numerous processes, including gene expression and stress responses. We performed a genome-wide analysis of DNA methylation in the human malaria parasite Plasmodium falciparum. We mapped the positions of methylated cytosines and identified a single functional DNA methyltransferase (Plasmodium falciparum DNA methyltransferase; PfDNMT) that may mediate these genomic modifications. These analyses revealed that the malaria genome is asymmetrically methylated and shares common features with undifferentiated plant and mammalian cells. Notably, core promoters are hypomethylated, and transcript levels correlate with intraexonic methylation. Additionally, there are sharp methylation transitions at nucleosome and exon-intron boundaries. These data suggest that DNA methylation could regulate virulence gene expression and transcription elongation. Furthermore, the broad range of action of DNA methylation and the uniqueness of PfDNMT suggest that the methylation pathway is a potential target for antimalarial strategies.

Electrochemical study on the effects of epigenetic cytosine methylation on anti-benzo[a]pyrene diol epoxide damage at TP53 oligomers

Anti-benzo[a]pyrene-r-7,t-8-dihydrodiol-t-9,10-epoxide (anti-BPDE) is a known carcinogen that damages DNA, and this damage is influenced by the DNA sequence and epigenetic factors. The influence of epigenetic cytosine methylation on the reaction with anti-BPDE at a known hotspot DNA damage site was studied electrochemically. Gold electrodes were modified with thiolated DNA oligomers spanning codons 270-276 of the TP53 gene. The oligomers exhibited 5-carbon cytosine methylation at the codon 273 location on the bound probe, the acquired complementary target, or both. Redox active diviologen compounds of the form C(12)H(25)V(2+)C(6)H(12)V(2+)C(12)H(25) (V(2+) = 4,4'-bipyridyl or viologen, C12-Viologen) were employed to detect anti-BPDE damage to DNA. DNA was exposed to racemic (±)- or enantiomerically pure (+)-anti-BPDE solutions followed by electrochemical interrogation in the presence of C12-Viologen. Background subtracted square wave voltammograms (SWV) showed the appearance of two peaks at approximately -0.38 V and -0.55 V vs Ag/AgCl upon anti-BPDE exposure. The acquired voltammetry is consistent with singly reduced C12-Viologen dimers bound at two different DNA environments, which arise from BPDE damage and are influenced by cytosine methylation and BPDE stereochemical considerations. UV spectroscopic and mass spectrometric methods employed to validate the electrochemical responses showed that (+)-anti-BPDE primarily adopts a minor groove bound orientation within the oligomers while selectively targeting the nontranscribed ssDNA sequence within the duplexes.

Polymorphism in exons CpG rich regions of the cyp17-II gene affecting its mRNA expression and reproductive endocrine levels in female Japanese flounder (Paralichthys olivaceus)

Cytochrome P450c17-II (cyp17-II) gene is an important factor affecting the growth, gonad differentiation and development, and other reproductive traits of fish. There are three CpG rich regions in the coding region of cyp17-II gene in Japanese flounder (Paralichthys olivaceus). The aim of this study was to understand whether mutations in exons of the cyp17-II gene occured at CpG sites, and mutations and methylation status of those CpG sites were involved in regulation of the expression level of cyp17-II gene and the reproductive endocrine of Japanese flounder. The results showed that three single nucleotide polymorphisms (SNPs) were identified. SNP1 [(c. G594A (p.Gly 188Arg)] located in exon 4 of L1 locus, and SNP2 (c.A939G) and SNP3 (c.C975T) of L2 locus located in CpG rich region of the exon 6 of cyp17-II gene. Furthermore, the A to G transition at 939bp position added a new methylation site to the cyp17-II coding region. According to multiple-comparison analysis, two loci (L1 and L2) were significantly associated with serum testosterone (T) level (P<0.05) and the expression of cyp17-II in ovary (P<0.01). Intriguingly, individuals with GG genotype of L2 locus containing eight CpG methylation sites had significantly lower serum testosterone level and cyp17-II mRNA expression than those with AA genotype containing seven CpG methylation sites. Moreover, the CpG site was highly methylated (≥77.8%) at 938 bp position of individuals with GG genotype of L2 locus. These implied that the mutation and methylation status of the coding region of cyp17-II could influence the gene expression and the reproductive endocrine levels in female Japanese flounder and L2 locus could be regarded as a candidate genetic or epigenetic marker for Japanese flounder breeding programs.

Demethylation of the coding region triggers the activation of the human testis-specific PDHA2 gene in somatic tissues

Human PDHA2 is a testis-specific gene that codes for the E(1)α subunit of Pyruvate Dehydrogenase Complex (PDC), a crucial enzyme system in cell energy metabolism. Since activation of the PDHA2 gene in somatic cells could be a new therapeutic approach for PDC deficiency, we aimed to identify the regulatory mechanisms underlying the human PDHA2 gene expression. Functional deletion studies revealed that the -122 to -6 promoter region is indispensable for basal expression of this TATA-less promoter, and suggested a role of an epigenetic program in the control of PDHA2 gene expression. Indeed, treatment of SH-SY5Y cells with the hypomethylating agent 5-Aza-2'-deoxycytidine (DAC) promoted the reactivation of the PDHA2 gene, by inducing the recruitment of the RNA polymerase II to the proximal promoter region and the consequent increase in PDHA2 mRNA levels. Bisulfite sequencing analysis revealed that DAC treatment induced a significant demethylation of the CpG island II (nucleotides +197 to +460) in PDHA2 coding region, while the promoter region remained highly methylated. Taken together with our previous results that show an in vivo correlation between PDHA2 expression and the demethylation of the CpG island II in testis germ cells, the present results show that internal methylation of the PDHA2 gene plays a part in its repression in somatic cells. In conclusion, our data support the novel finding that methylation of the PDHA2 coding region can inhibit gene transcription. This represents a key mechanism for absence of PDHA2 expression in somatic cells and a target for PDC therapy.

Site-specific methylation in gene coding region underlies transcriptional silencing of the Phytochrome A epiallele in Arabidopsis thaliana

DNA methylation in cytosine residues plays an important role in regulating gene expression. Densely methylated transgenes are often silenced. In contrast, several eukaryotic genomes express moderately methylated genes. These methylations are found in the CG context within the coding region (gene body). The role of gene body methylation in gene expression, however, is not clear. The Arabidopsis Phytochrome A epiallele, phyA', carries hypermethylation in several CG sites resident to the coding region. As a result, phyA' is transcriptionally silenced and confers strong mutant phenotype. Mutations in chromatin modification factors and RNAi genes failed to revert the mutant phenotype, suggesting the involvement of a distinct epigenetic mechanism associated with phyA' silencing. Using the forward genetics approach, a suppressor line, termed as suppressor of p hyA' silencing 1 (sps1), was isolated. Genetic and molecular analysis revealed that sps1 mutation reactivates the phyA' locus without altering its methylation density. However, hypomethylation at a specific CG site in exon 1 was consistently associated with the release of phyA' silencing. While gene underlying sps1 mutation is yet to be identified, microarray analysis suggested that its targets are the expressed genes or euchromatic loci in Arabidopsis genome. By identifying the association of phyA' silencing with the methylation of a specific CG site in exon 1, the present work shows that site-specific methylation confers greater effect on transcription than the methylation density within gene-body. Further, as the identified site (exon 1) is not critical for the promoter activity, transcription elongation rather than transcription initiation is likely to be affected by this site-specific CG methylation.

Development of useful recombinant promoter and its expression analysis in different plant cells using confocal laser scanning microscopy

BACKGROUND

Designing functionally efficient recombinant promoters having reduced sequence homology and enhanced promoter activity will be an important step toward successful stacking or pyramiding of genes in a plant cell for developing transgenic plants expressing desired traits(s). Also basic knowledge regarding plant cell specific expression of a transgene under control of a promoter is crucial to assess the promoter's efficacy.

METHODOLOGY/PRINCIPAL FINDINGS

We have constructed a set of 10 recombinant promoters incorporating different up-stream activation sequences (UAS) of Mirabilis mosaic virus sub-genomic transcript (MS8, -306 to +27) and TATA containing core domains of Figwort mosaic virus sub-genomic transcript promoter (FS3, -271 to +31). Efficacies of recombinant promoters coupled to GUS and GFP reporter genes were tested in tobacco protoplasts. Among these, a 369-bp long hybrid sub-genomic transcript promoter (MSgt-FSgt) showed the highest activity in both transient and transgenic systems. In a transient system, MSgt-FSgt was 10.31, 2.86 and 2.18 times more active compared to the CaMV35S, MS8 and FS3 promoters, respectively. In transgenic tobacco (Nicotiana tabaccum, var. Samsun NN) and Arabidopsis plants, the MSgt-FSgt hybrid promoter showed 14.22 and 7.16 times stronger activity compared to CaMV35S promoter respectively. The correlation between GUS activity and uidA-mRNA levels in transgenic tobacco plants were identified by qRT-PCR. Both CaMV35S and MSgt-FSgt promoters caused gene silencing but the degree of silencing are less in the case of the MSgt-FSgt promoter compared to CaMV35S. Quantification of GUS activity in individual plant cells driven by the MSgt-FSgt and the CaMV35S promoter were estimated using confocal laser scanning microscopy and compared.

CONCLUSION AND SIGNIFICANCE

We propose strong recombinant promoter MSgt-FSgt, developed in this study, could be very useful for high-level constitutive expression of transgenes in a wide variety of plant cells.

Nucleosome deposition and DNA methylation at coding region boundaries

BACKGROUND

Nucleosome deposition downstream of transcription initiation and DNA methylation in the gene body suggest that control of transcription elongation is a key aspect of epigenetic regulation.

RESULTS

Here we report a genome-wide observation of distinct peaks of nucleosomes and methylation at both ends of a protein coding unit. Elongating polymerases tend to pause near both coding ends immediately upstream of the epigenetic peaks, causing a significant reduction in elongation efficiency. Conserved features in underlying protein coding sequences seem to dictate their evolutionary conservation across multiple species. The nucleosomal and methylation marks are commonly associated with high sequence-encoded DNA-bending propensity but differentially with CpG density. As the gene grows longer, the epigenetic codes seem to be shifted from variable inner sequences toward boundary regions, rendering the peaks more prominent in higher organisms.

CONCLUSIONS

Recent studies suggest that epigenetic inhibition of transcription elongation facilitates the inclusion of constitutive exons during RNA splicing. The epigenetic marks we identified here seem to secure the first and last coding exons from exon skipping as they are indispensable for accurate translation.

Trans-generation inheritance of methylation patterns in a tobacco transgene following a post-transcriptional silencing event

We have studied the inheritance of the epigenetic state of tobacco transgenes whose expression was post-transcriptionally silenced by an invertedly repeated silencer locus. We show that, in hybrids, the coding region of the target neomycin phosphotransferase (nptII) gene was almost exclusively methylated at CG configurations, and dense non-CG methylation occurred in the 3' untranslated region. Homologous sequences in the silencer locus were heavily methylated at both CG and non-CG motifs. After segregation of the silencer locus, the CG methylation but not the non-CG methylation of the target genes was transmitted to the progeny. In the segregants, we observed an overall increase of CG methylation in the target genes, associated with a re-distribution from the 3' end of the coding region towards the middle. This pattern was inherited with some fluctuation for at least two additional generations in the absence of a detectable T-DNA-derived small RNA fraction. Thus CG methylation is not cleared during meiosis and may be inherited over generations without RNA signals being present. These epi-allelic variants re-expressed the reporter gene immediately after segregation of the trigger, showing that relatively dense CG methylation (approximately 60-80%) imprinted on most of the coding region (>500 bp) did not reduce expression compared with the parental non-methylated locus. We propose that the genic CG methylation seen in euchromatic regions of the genome may originate from ancient post-transcriptional gene silencing events as a result of adventitiously produced methylation-directing RNA molecules.

Global analysis of genetic, epigenetic and transcriptional polymorphisms in Arabidopsis thaliana using whole genome tiling arrays

Whole genome tiling arrays provide a high resolution platform for profiling of genetic, epigenetic, and gene expression polymorphisms. In this study we surveyed natural genomic variation in cytosine methylation among Arabidopsis thaliana wild accessions Columbia (Col) and Vancouver (Van) by comparing hybridization intensity difference between genomic DNA digested with either methylation-sensitive (HpaII) or -insensitive (MspI) restriction enzyme. Single Feature Polymorphisms (SFPs) were assayed on a full set of 1,683,620 unique features of Arabidopsis Tiling Array 1.0F (Affymetrix), while constitutive and polymorphic CG methylation were assayed on a subset of 54,519 features, which contain a 5′CCGG3′ restriction site. 138,552 SFPs (1% FDR) were identified across enzyme treatments, which preferentially accumulated in pericentromeric regions. Our study also demonstrates that at least 8% of all analyzed CCGG sites were constitutively methylated across the two strains, while about 10% of all analyzed CCGG sites were differentially methylated between the two strains. Within euchromatin arms, both constitutive and polymorphic CG methylation accumulated in central regions of genes but under-represented toward the 5′ and 3′ ends of the coding sequences. Nevertheless, polymorphic methylation occurred much more frequently in gene ends than constitutive methylation. Inheritance of methylation polymorphisms in reciprocal F1 hybrids was predominantly additive, with F1 plants generally showing levels of methylation intermediate between the parents. By comparing gene expression profiles, using matched tissue samples, we found that magnitude of methylation polymorphism immediately upstream or downstream of the gene was inversely correlated with the degree of expression variation for that gene. In contrast, methylation polymorphism within genic region showed weak positive correlation with expression variation. Our results demonstrated extensive genetic and epigenetic polymorphisms between Arabidopsis accessions and suggested a possible relationship between natural CG methylation variation and gene expression variation.

The functional expression of DNA sequence depends on the chromatin status. Epigenetic marks at specific loci could affect local chromatin accessibility, thus affect the gene activity of that loci. We applied an enzyme methylome approach to globally detect one type of epigenetic mark, cytosine methylation at CCGG restriction sites. Simultaneous transcriptional profiling allowed gene expression differences to be compared with DNA methylation differences, suggesting functional regulatory regions. Our method reveals natural variation in chromatin patterns which may underlie phenotypic variation.

High-resolution mapping of epigenetic modifications of the rice genome uncovers interplay between DNA methylation, histone methylation, and gene expression

We present high-resolution maps of DNA methylation and H3K4 di- and trimethylation of two entire chromosomes and two fully sequenced centromeres in rice (Oryza sativa) shoots and cultured cells. This analysis reveals combinatorial interactions between these epigenetic modifications and chromatin structure and gene expression. Cytologically densely stained heterochromatin had less H3K4me2 and H3K4me3 and more methylated DNA than the less densely stained euchromatin, whereas centromeres had a unique epigenetic composition. Most transposable elements had highly methylated DNA but no H3K4 methylation, whereas more than half of protein-coding genes had both methylated DNA and di- and/or trimethylated H3K4. Methylation of DNA but not H3K4 was correlated with suppressed transcription. By contrast, when both DNA and H3K4 were methylated, transcription was only slightly reduced. Transcriptional activity was positively correlated with the ratio of H3K4me3/H3K4me2: genes with predominantly H3K4me3 were actively transcribed, whereas genes with predominantly H3K4me2 were transcribed at moderate levels. More protein-coding genes contained all three modifications, and more transposons contained DNA methylation in shoots than cultured cells. Differential epigenetic modifications correlated to tissue-specific expression between shoots and cultured cells. Collectively, this study provides insights into the rice epigenomes and their effect on gene expression and plant development.

Transgene-induced silencing of Arabidopsis phytochrome A gene via exonic methylation

Transgene-induced promoter or enhancer methylation clearly retards gene activity. While exonic methylation of genes is frequently observed in the RNAi process, only sporadic evidence has demonstrated its definitive role in gene suppression. Here, we report the isolation of a transcriptionally suppressed epi-allele of the Arabidopsis thaliana phytochrome A gene (PHYA) termed phyA' that shows methylation only in symmetric CG sites resident in exonic regions. These exonic modifications confer a strong phyA mutant phenotype, characterized by elongated hypocotyls in seedlings grown under continuous far-red light. De-methylation of phyA' in the DNA methyl transferase I (met1) mutant background increased PHYA expression and restored the wild-type phenotype, confirming the pivotal role of exonic CG methylation in maintaining the altered epigenetic state. PHYA epimutation was apparently induced by a transgene locus; however, it is stably maintained following segregation. Chromatin immunoprecipitation assays revealed association with dimethyl histone H3 lysine 9 (H3K9me2), a heterochromatic marker, within the phyA' coding region. Therefore, transgene-induced exonic methylation can lead to chromatin alteration that affects gene expression, most likely through reduction in the transcription rate.

Mapping the genome landscape using tiling array technology

With the availability of complete genome sequences for a growing number of organisms, high-throughput methods for gene annotation and analysis of genome dynamics are needed. The application of whole-genome tiling microarrays for studies of global gene expression is providing a more unbiased view of the transcriptional activity within genomes. For example, this approach has led to the identification and isolation of many novel non-protein-coding RNAs (ncRNAs), which have been suggested to comprise a major component of the transcriptome that have novel functions involved in epigenetic regulation of the genome. Additionally, tiling arrays have been recently applied to the study of histone modifications and methylation of cytosine bases (DNA methylation). Surprisingly, recent studies combining the analysis of gene expression (transcriptome) and DNA methylation (methylome) using whole-genome tiling arrays revealed that DNA methylation regulates the expression levels of many ncRNAs. Further capture and integration of additional types of genome-wide data sets will help to illuminate additional hidden features of the dynamic genomic landscape that are regulated by both genetic and epigenetic pathways in plants.

DNA methylation dynamics in plant genomes

Cytosine bases are extensively methylated in the DNA of plant genomes. DNA methylation has been implicated in the silencing of transposable elements and genes, and loss of methylation can have severe consequences for the organism. The recent methylation profiling of the entire Arabidopsis genome has provided insight into the extent of DNA methylation and its functions in silencing and gene transcription. Patterns of DNA methylation are faithfully maintained across generations, but some changes in DNA methylation are observed in terminally differentiated tissues. Demethylation by a DNA glycosylase is required for the expression of imprinted genes in the endosperm and de novo methylation might play a role in the selective silencing of certain self-incompatibility alleles in the tapetum. Because DNA methylation patterns are faithfully inherited, changes in DNA methylation that arise somatically during the plant life cycle have the possibility of being propagated. Therefore, epimutations might be an important source of variation during plant evolution.

Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription

Cytosine methylation, a common form of DNA modification that antagonizes transcription, is found at transposons and repeats in vertebrates, plants and fungi. Here we have mapped DNA methylation in the entire Arabidopsis thaliana genome at high resolution. DNA methylation covers transposons and is present within a large fraction of A. thaliana genes. Methylation within genes is conspicuously biased away from gene ends, suggesting a dependence on RNA polymerase transit. Genic methylation is strongly influenced by transcription: moderately transcribed genes are most likely to be methylated, whereas genes at either extreme are least likely. In turn, transcription is influenced by methylation: short methylated genes are poorly expressed, and loss of methylation in the body of a gene leads to enhanced transcription. Our results indicate that genic transcription and DNA methylation are closely interwoven processes.

Novel relationships among DNA methylation, histone modifications and gene expression in Ascobolus

By studying Ascobolus strains methylated in various portions of the native met2 gene or of the hph transgene, we generalized our previous observation that methylation of the downstream portion of a gene promotes its stable silencing and triggers the production of truncated transcripts which rarely extend through the methylated region. In contrast, methylation of the promoter region does not promote efficient gene silencing. The chromatin state of met2 methylated strains was investigated after partial micrococcal nuclease (MNase) digestion. We show that MNase sensitive sites present along the unmethylated regions are no longer observed along the methylated ones. These chromatin changes are not resulting from the absence of transcription. They are associated, in both met2 and hph, with modifications of core histones corresponding, on the N terminus of histone H3, to an increase of dimethylation of lysine 9 and a decrease of dimethylation of lysine 4. Contrary to other organisms, these changes are independent of the transcriptional state of the genes, and furthermore, no decrease in acetylation of histone H4 is observed in silenced genes.

The cauliflower mosaic virus 35S promoter extends into the transcribed region

A 60-nucleotide region (S1) downstream of the transcription start site of the cauliflower mosaic virus 35S RNA can enhance gene expression. By using transient expression assays with plant protoplasts, this activity was shown to be at least partially due to the effect of transcriptional enhancers within this region. We identify sequence motifs with enhancer function, which are normally masked by the powerful upstream enhancers of the 35S promoter. A repeated CT-rich motif is involved both in enhancer function and in interaction with plant nuclear proteins. The S1 region can also enhance expression from heterologous promoters.

Hypermethylation of the wild-type ferrochelatase allele is closely associated with severe liver complication in a family with erythropoietic protoporphyria

Erythropoietic protoporphyria (EPP) is an inherited disorder of heme biosynthesis caused by cellular decreases in ferrochelatase (FECH) activity. Clinical expression of this disorder usually requires coinheritance of a mutant FECH allele and a normal FECH allele expressed at a low level. In this study, we investigated the methylation status of a normal, but poorly expressed, FECH gene in a single Japanese family with EPP. In this family, the proband died from liver failure, whereas the mother and sister exhibited overt EPP with mild liver dysfunction. A splicing mutation (IVS9+1g-->a) in the FECH gene, which produces a mutant FECH transcript lacking exon 9, was detected in the maternal allele of the proband and his sister. All subjects, including the father, who did not exhibit EPP, possessed the IVS3-48c/c genotype. This allele increases the proportion of aberrantly spliced mRNA, resulting in reduced FECH activity. Normal FECH transcripts were, however, detected in the mother and sister, but not in the proband. The CpG sites in the region from bases -78 to -31 were partially methylated in the proband and his father, but not in his mother or sister. Additionally, CpG methylation within this region reduced transcription of the FECH gene. These results suggest that whereas the combination of a maternal IVS9+1a allele and a paternal IVS3-48c allele results in overt EPP, CpG methylation of the FECH gene promoter, likely inherited from the father, increases the severity of EPP, leading to fatal liver failure, as seen in the proband.

Nonrandom intragenic variations in patterns of codon bias implicate a sequential interplay between transitional genetic drift and functional amino acid selection

Although most codon third bases appear to be functionless, the synonymous codons so defined exhibit a strikingly nonrandom distribution (codon bias) within human and other genes. To examine this phenomenon further, we generated a database of DNA sequences encoding human transmembrane cell-surface receptor proteins. Using this database we show here that the guanine and cytosine content of codon third bases (GC3) varies intragenically with the nature of the specified receptor domains (transmembrane > extracellular > intracellular domains; p < 0.001), the phenotype of the encoded amino acids (hydrophobic > hydrophilic > neutral amino acids; p < 0.001), and the receptor affiliation of the transmembrane (G-protein-coupled receptors > receptor tyrosine kinases; p < 0.001). Within gene regions specifying transmembrane domains, GC3 declines as domain functionality becomes redundant with increasing hydrophobicity (p < 0.001). Codons containing the second-base cytosine (XCZ, which encodes neutral amino acids) are selectively depleted of third-base adenine content (A3: XCA codons) when encoding transmembrane domain residues, consistent with positive selection for transitional mutation of XCG to XTG (which encodes hydrophobic amino acids) rather than to the synonymous XCA. Supporting this XCG --> XTG mechanism of codon bias, the G3:A3 ratio of codons specifying the transmembrane amino acid glycine (GGZ) is intermediate between that of its functional homolog alanine (GCZ) and that of hydrophobic valine (GTZ), even though the C3:T3 ratios are similar. Conversely, nearest-neighbor analysis of third bases 5' to codons specifying valine and leucine (CTZ) confirms a significant difference in C3:T3 but not G3:A3 ratios (i.e., C3/G1 --> T3/G1 > C3/A1; p < 0.001), consistent with the functionally advantageous retention of hydrophobic residues. These data raise the possibility that patterns of intragenic codon bias reflect a balance between negative and positive selection, suggesting in turn that analysis of codon third-base usage may help to predict the functional significance of encoded products.

A novel feature of microsatellites in plants: a distribution gradient along the direction of transcription

A computer-based analysis was conducted to assess the characteristics of microsatellites in transcribed regions of rice and Arabidopsis. In addition, two mammals were simultaneously analyzed for a comparative analysis. Our analyses confirmed a novel plant-specific feature in which there is a gradient in microsatellite density along the direction of transcription. It was also confirmed that pyrimidine-rich microsatellites are found intensively near the transcription start sites, specifically in the two plants, but not in the mammals. Our results suggest that microsatellites located at high frequency in the 5'-flanking regions of plant genes can potentially act as factors in regulating gene expression.

Reduced expression of GNA11 and silencing of MCT1 in human breast cancers

Alteration in the methylation status of a gene is often associated with its altered expression. Based on a genome scanning technique for differences in CpG methylations, methylation-sensitive representational difference analysis, DNA fragments hypermethylated in a human breast cancer were isolated. A DNA fragment was isolated from intron 1 of guanine-nucleotide-binding protein alpha-11 (GNA11). mRNA expression of GNA11 was shown to be decreased in 10 of 16 breast cancers by RT-PCR analysis, and the immunoreactivity of the GNA11 product, Galpha11 subunit of heterotrimeric G-protein, was observed to be reduced in 14 of the 16 cancers by immunohistochemistry. Methylation of a CpG island (CGI) in the 5' region of GNA11 or that of intron 1 did not show a clear correlation with its decreased expression. Another DNA fragment was isolated from a CGI in the 5' upstream region of monocarboxylate transporter 1 (MCT1), and was methylated in 4 of 20 breast cancers. The CGI was also methylated in a human breast cancer cell line, MDA-MB-231, and quantitative RT-PCR showed that its expression was almost lost in the cell line. By treatment of the cells with a demethylating agent, 5-aza-2'-deoxycytidine, the methylation was removed and the expression was restored. GNA11 is involved in signalling of gonadotropin-releasing hormone receptor, which negatively regulates cell growth. MCT1 is involved in cellular transportation of butyrate, which induces cellular differentiation. Downregulation of these two genes was suggested to be involved in human breast cancers.

Methylation of alpha-type embryonic globin gene alpha pi represses transcription in primary erythroid cells

The inverse relationship between expression and methylation of beta-type globin genes is well established. However, little is known about the relationship between expression and methylation of avian alpha-type globin genes. The embryonic alpha(pi)-globin promoter was unmethylated, and alpha(pi)-globin RNA was easily detected in 5-day chicken erythroid cells. A progressive methylation of the CpG dinucleotides in the alpha(pi) promoter associated with loss of expression of alpha(pi)-globin gene was seen during development in primary erythroid cells. A 315-bp alpha(pi)-globin promoter region was cloned in an expression construct (alpha(pi)pGL3E) containing a luciferase reporter gene and SV40 enhancer. The alpha(pi)pGL3E construct was transfected into primary erythroid cells derived from 5-day-old chicken embryos. Methylation of alpha(pi)pGL3E plasmid and alpha(pi)-globin promoter alone resulted in a 20-fold and 7-fold inhibition of expression, respectively. The fully methylated but not the unmethylated 315-bp alpha(pi)-globin gene promoter fragment formed a methyl cytosine-binding protein complex (MeCPC). Chromatin immunoprecipitation assays were combined with quantitative real-time polymerase chain reaction to assess histone acetylation associated with the alpha(pi)-globin gene promoter. Slight hyperacetylation of histone H3 but a marked hyperacetylation of histone H4 was seen in 5-day when compared with 14-day erythroid cells. These results demonstrate that methylation can silence transcription of an avian alpha-type embryonic globin gene in homologous primary erythroid cells, possibly by interacting with an MeCPC and histone deacetylase complex.

Inducible DNA demethylation mediated by the maize Suppressor-mutator transposon-encoded TnpA protein

Heritable epigenetic inactivation of the maize Suppressor-mutator (Spm) transposon is associated with promoter methylation, and its reversal is mediated by the transposon-encoded TnpA protein. We have developed an assay that permits demethylation of the Spm sequence to be controlled by inducing the expression of TnpA in plant cells. Using this assay, we show that demethylation is a rapid, active process. TnpA is a weak transcriptional activator, and deletions that abolish its transcriptional activity also eliminate its demethylation activity. We show that cell cycle and DNA synthesis inhibitors interfere with TnpA-mediated Spm demethylation. We further show that TnpA has a much lower affinity for fully methylated than for hemimethylated or unmethylated DNA fragments derived from Spm termini. Based on these observations, we suggest that TnpA binds to the postreplicative, hemimethylated Spm sequence and promotes demethylation either by creating an appropriate demethylation substrate or by itself participating in or recruiting a demethylase.

Reduced expression of the insulin-induced protein 1 and p41 Arp2/3 complex genes in human gastric cancers

Aberrantly methylated DNA fragments in a human gastric cancer were searched for by a genome-scanning method, methylation-sensitive-representational difference analysis (MS-RDA). Six DNA fragments flanked by CpG islands (CGIs) and hypermethylated in the cancer were isolated. Four of the 6 fragments possessed genes in their vicinities. Quantitative RT-PCR analysis of the 4 genes showed reduced expression of 2 genes in cancers: Insulin-induced protein 1 (INSIG1/CL-6) and p41 Arp2/3 complex (p41-Arc). As for INSIG1, a DNA fragment was derived from the edge of a CGI in the promoter region. The edge was methylated in 11 of 22 primary gastric cancers, whereas the center was not methylated in any cancer. INSIG1 expression was markedly reduced in 19 cancers, including the 11 cancers with the methylation. By 5-aza-2'-deoxycytidine treatment of 5 cell lines with the methylation of the edge, partial restoration of INSIG1 expression was detected only in 2 of them. These data indicated that, although the reduced INSIG1 expression in cancers was associated with the methylation at the edge of the CGI in the promoter region, the methylation was likely to be a secondary change. As for p41-Arc, a DNA fragment was derived from a CGI overlapping exon 8, and its methylation did not correlate with its expression. However, methylation of a CGI in the promoter region with a marked reduction of its expression was observed in 1 of the 22 primary cancers. INSIG1 and p41-Arc are known to be involved in cellular differentiation and morphology, respectively, and it was suggested that their reduced expressions might be involved in gastric cancer development or progression.

Methylation of promoter proximal-transcribed sequences of an embryonic globin gene inhibits transcription in primary erythroid cells and promotes formation of a cell type-specific methyl cytosine binding complex

The methylation pattern of a 248-base pair proximal transcribed region (rho248) of the avian embryonic rho-globin gene was found to correlate inversely with stage-specific expression in avian erythroid cells. In vitro methylation of the rho248 segment alone (in the absence of promoter methylation) resulted in a 5-fold inhibition of transcription in a transient transfection assay in primary erythroid cells in which the transfected gene is packaged into nucleosomal chromatin. This effect was observed if the rho248 segment was positioned adjacent to the promoter but not when it was located 2.7 kilobases downstream. Fully methylated but not unmethylated rho248 formed a novel cell type-specific methyl cytosine-binding protein complex (MeCPC) that contained methyl binding domain protein-2 (MBD-2) and histone deacetylase 1 proteins but differed from MeCP-1. The histone deacetylase inhibitor trichostatin A failed to relieve methylation-mediated repression of transcription from the rho-gene promoter, supporting the notion of the dominance of methylation over histone deacetylation in silencing through CpG-rich sequences at this locus. These data demonstrate that site-specific methylation of a vertebrate gene 5'-transcribed region alone at the exact CpGs that are methylated in vivo can suppress transcription in homologous primary cells and facilitate binding to a cell type-specific MeCPC.

Contribution of downstream promoter elements to transcriptional regulation of the rice tungro bacilliform virus promoter

Downstream sequences influence activity of the rice tungro bacilliform virus (RTBV) promoter in protoplasts derived from cultured rice cells. We previously identified a DNA element located between positions +50 and +90 relative to the transcription start site to which rice nuclear proteins bind. In this study, using DNA UV crosslinking assays, we show that two rice nuclear proteins bind specifically to this DNA element. We demonstrate that the DNA element enhances RTBV promoter activity in a copy number-dependent manner when transferred to a position upstream of the promoter. In addition, using electrophoretic mobility shift assays, we show that at least two novel nuclear proteins from rice cell suspension cultures bind to a subregion (from +50 to +59) of the DNA element and that a protein from rice root, but not shoot, nuclear extracts interacts with a perfect palindromic sequence motif located within the sequence +45 to +59. Furthermore, a position-dependent GAGA motif, present in three copies within downstream promoter sequences from +1 to +50, is involved in the regulation of RTBV promoter activity.

The late flowering phenotype of fwa mutants is caused by gain-of-function epigenetic alleles of a homeodomain gene

The transition to flowering in Arabidopsis thaliana is delayed in fwa mutant plants. FWA was identified by loss-of-function mutations in normally flowering revertants of the fwa mutant, and it encodes a homeodomain-containing transcription factor. The DNA sequence of wild-type and fwa mutant alleles was identical in the genomic region of FWA. Furthermore, the FWA gene is ectopically expressed in fwa mutants and silenced in mature wild-type plants. This silencing is associated with extensive methylation of two direct repeats in the 5' region of the gene. The late flowering phenotype, ectopic FWA expression, and hypomethylation of the repeats were also induced in the ddm1 hypomethylated background. Mechanisms for establishment and maintenance of the epigenetic mark on FWA are discussed.

Transgene silencing in monocots

Plant gene silencing was originally thought to be a quirk of transformation procedures, but is now recognized to be a facet of vitally important gene regulatory systems, present in all organisms. Monocot plants, especially the grasses, play a foremost role in the agricultural economy of all nations, and their biotechnological manipulation offers great potential for both developed and developing countries. Here, we review reported instances of transgene silencing in monocots and relate the processes of transcriptional and post-transcriptional gene silencing (TGS, PTGS) in perspective to the rapidly burgeoning knowledge of these phenomena in many organisms. Recent findings include the involvement of an RNA-dependent RNA polymerase and a nuclease in PTGS systems and the close relationship between methylation and chromatin structure in TGS events.

(TRANS)GENE SILENCING IN PLANTS: How Many Mechanisms?

Epigenetic silencing of transgenes and endogenous genes can occur at the transcriptional level (TGS) or at the posttranscriptional level (PTGS). Because they can be induced by transgenes and viruses, TGS and PTGS probably reflect alternative (although not exclusive) responses to two important stress factors that the plant's genome has to face: the stable integration of additional DNA into chromosomes and the extrachromosomal replication of a viral genome. TGS, which results from the impairment of transcription initiation through methylation and/or chromatin condensation, could derive from the mechanisms by which transposed copies of mobile elements and T-DNA insertions are tamed. PTGS, which results from the degradation of mRNA when aberrant sense, antisense, or double-stranded forms of RNA are produced, could derive from the process of recovery by which cells eliminate pathogens (RNA viruses) or their undesirable products (RNA encoded by DNA viruses). Mechanisms involving DNA-DNA, DNA-RNA, or RNA-RNA interactions are discussed to explain the various pathways for triggering (trans)gene silencing in plants.