Plant-RRBS, a bisulfite and next-generation sequencing-based methylome profiling method enriching for coverage of cytosine positions

Background: Cytosine methylation in plant genomes is important for the regulation of gene transcription and transposon activity. Genome-wide methylomes are studied upon mutation of the DNA methyltransferases, adaptation to environmental stresses or during development. However, from basic biology to breeding programs, there is a need to monitor multiple samples to determine transgenerational methylation inheritance or differential cytosine methylation. Methylome data obtained by sodium hydrogen sulfite (bisulfite)-conversion and next-generation sequencing (NGS) provide genome- wide information on cytosine methylation. However, a profiling method that detects cytosine methylation state dispersed over the genome would allow high-throughput analysis of multiple plant samples with distinct epigenetic signatures. We use specific restriction endonucleases to enrich for cytosine coverage in a bisulfite and NGS-based profiling method, which was compared to whole-genome bisulfite sequencing of the same plant material. Methods: We established an effective methylome profiling method in plants, termed plant-reduced representation bisulfite sequencing (plant-RRBS), using optimized double restriction endonuclease digestion, fragment end repair, adapter ligation, followed by bisulfite conversion, PCR amplification and NGS. We report a performant laboratory protocol and a straightforward bioinformatics data analysis pipeline for plant-RRBS, applicable for any reference-sequenced plant species. Results: As a proof of concept, methylome profiling was performed using an Oryza sativa ssp. indica pure breeding line and a derived epigenetically altered line (epiline). Plant-RRBS detects methylation levels at tens of millions of cytosine positions deduced from bisulfite conversion in multiple samples. To evaluate the method, the coverage of cytosine positions, the intra-line similarity and the differential cytosine methylation levels between the pure breeding line and the epiline were determined. Plant-RRBS reproducibly covers commonly up to one fourth of the cytosine positions in the rice genome when using MspI-DpnII within a group of five biological replicates of a line. The method predominantly detects cytosine methylation in putative promoter regions and not-annotated regions in rice. Conclusions: Plant-RRBS offers high-throughput and broad, genome- dispersed methylation detection by effective read number generation obtained from reproducibly covered genome fractions using optimized endonuclease combinations, facilitating comparative analyses of multi-sample studies for cytosine methylation and transgenerational stability in experimental material and plant breeding populations

Complete chloroplast genome sequence of Holoparasite Cistanche Deserticola (Orobanchaceae) reveals gene loss and horizontal gene transfer from Its host Haloxylon Ammodendron (Chenopodiaceae)

The central function of chloroplasts is to carry out photosynthesis, and its gene content and structure are highly conserved across land plants. Parasitic plants, which have reduced photosynthetic ability, suffer gene losses from the chloroplast (cp) genome accompanied by the relaxation of selective constraints. Compared with the rapid rise in the number of cp genome sequences of photosynthetic organisms, there are limited data sets from parasitic plants. The authors report the complete sequence of the cp genome of Cistanche deserticola, a holoparasitic desert species belonging to the family Orobanchaceae

The Mitochondrial Genome of the Legume Vigna radiata and the Analysis of Recombination across Short Mitochondrial Repeats

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38–297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes

BORIS/CTCFL-mediated transcriptional regulation of the hTERT telomerase gene in testicular and ovarian tumor cells

Telomerase activity, not detectable in somatic cells but frequently activated during carcinogenesis, confers immortality to tumors. Mechanisms governing expression of the catalytic subunit hTERT, the limiting factor for telomerase activity, still remain unclear. We previously proposed a model in which the binding of the transcription factor CTCF to the two first exons of hTERT results in transcriptional inhibition in normal cells. This inhibition is abrogated, however, by methylation of CTCF binding sites in 85% of tumors. Here, we showed that hTERT was unmethylated in testicular and ovarian tumors and in derivative cell lines. We demonstrated that CTCF and its paralogue, BORIS/CTCFL, were both present in the nucleus of the same cancer cells and bound to the first exon of hTERT in vivo. Moreover, exogenous BORIS expression in normal BORIS-negative cells was sufficient to activate hTERT transcription with an increasing number of cell passages. Thus, expression of BORIS was sufficient to allow hTERT transcription in normal cells and to counteract the inhibitory effect of CTCF in testicular and ovarian tumor cells. These results define an important contribution of BORIS to immortalization during tumorigenesis

Global DNA Hypomethylation in Peripheral Blood Leukocytes as a Biomarker for Cancer Risk: A Meta-Analysis

BACKGROUND: Good biomarkers for early detection of cancer lead to better prognosis. However, harvesting tumor tissue is invasive and cannot be routinely performed. Global DNA methylation of peripheral blood leukocyte DNA was evaluated as a biomarker for cancer risk. METHODS: We performed a meta-analysis to estimate overall cancer risk according to global DNA hypomethylation levels among studies with various cancer types and analytical methods used to measure DNA methylation. Studies were systemically searched via PubMed with no language limitation up to July 2011. Summary estimates were calculated using a fixed effects model. RESULTS: The subgroup analyses by experimental methods to determine DNA methylation level were performed due to heterogeneity within the selected studies (p<0.001, I(2): 80%). Heterogeneity was not found in the subgroup of %5-mC (p = 0.393, I(2): 0%) and LINE-1 used same target sequence (p = 0.097, I(2): 49%), whereas considerable variance remained in LINE-1 (p<0.001, I(2): 80%) and bladder cancer studies (p = 0.016, I(2): 76%). These results suggest that experimental methods used to quantify global DNA methylation levels are important factors in the association study between hypomethylation levels and cancer risk. Overall, cancer risks of the group with the lowest DNA methylation levels were significantly higher compared to the group with the highest methylation levels [OR (95% CI): 1.48 (1.28-1.70)]. CONCLUSIONS: Global DNA hypomethylation in peripheral blood leukocytes may be a suitable biomarker for cancer risk. However, the association between global DNA methylation and cancer risk may be different based on experimental methods, and region of DNA targeted for measuring global hypomethylation levels as well as the cancer type. Therefore, it is important to select a precise and accurate surrogate marker for global DNA methylation levels in the association studies between global DNA methylation levels in peripheral leukocyte and cancer risk

The Evolution of Epigenetic Regulators CTCF and BORIS/CTCFL in Amniotes

CTCF is an essential, ubiquitously expressed DNA-binding protein responsible for insulator function, nuclear architecture, and transcriptional control within vertebrates. The gene CTCF was proposed to have duplicated in early mammals, giving rise to a paralogue called “brother of regulator of imprinted sites” (BORIS or CTCFL) with DNA binding capabilities similar to CTCF, but testis-specific expression in humans and mice. CTCF and BORIS have opposite regulatory effects on human cancer-testis genes, the anti-apoptotic BAG1 gene, the insulin-like growth factor 2/H19 imprint control region (IGF2/H19 ICR), and show mutually exclusive expression in humans and mice, suggesting that they are antagonistic epigenetic regulators. We discovered orthologues of BORIS in at least two reptilian species and found traces of its sequence in the chicken genome, implying that the duplication giving rise to BORIS occurred much earlier than previously thought. We analysed the expression of CTCF and BORIS in a range of amniotes by conventional and quantitative PCR. BORIS, as well as CTCF, was found widely expressed in monotremes (platypus) and reptiles (bearded dragon), suggesting redundancy or cooperation between these genes in a common amniote ancestor. However, we discovered that BORIS expression was gonad-specific in marsupials (tammar wallaby) and eutherians (cattle), implying that a functional change occurred in BORIS during the early evolution of therian mammals. Since therians show imprinting of IGF2 but other vertebrate taxa do not, we speculate that CTCF and BORIS evolved specialised functions along with the evolution of imprinting at this and other loci, coinciding with the restriction of BORIS expression to the germline and potential antagonism with CTCF

Mosaic Origins of a Complex Chimeric Mitochondrial Gene in Silene vulgaris

Chimeric genes are significant sources of evolutionary innovation that are normally created when portions of two or more protein coding regions fuse to form a new open reading frame. In plant mitochondria astonishingly high numbers of different novel chimeric genes have been reported, where they are generated through processes of rearrangement and recombination. Nonetheless, because most studies do not find or report nucleotide variation within the same chimeric gene, evolution after the origination of these chimeric genes remains unstudied. Here we identify two alleles of a complex chimera in Silene vulgaris that are divergent in nucleotide sequence, genomic position relative to other mitochondrial genes, and expression patterns. Structural patterns suggest a history partially influenced by gene conversion between the chimeric gene and functional copies of subunit 1 of the mitochondrial ATP synthase gene (atp1). We identified small repeat structures within the chimeras that are likely recombination sites allowing generation of the chimera. These results establish the potential for chimeric gene divergence in different plant mitochondrial lineages within the same species. This result contrasts with the absence of diversity within mitochondrial chimeras found in crop species

A Cell Cycle Role for the Epigenetic Factor CTCF-L/BORIS

CTCF is a ubiquitous epigenetic regulator that has been proposed as a master keeper of chromatin organisation. CTCF-like, or BORIS, is thought to antagonise CTCF and has been found in normal testis, ovary and a large variety of tumour cells. The cellular function of BORIS remains intriguing although it might be involved in developmental reprogramming of gene expression patterns. We here unravel the expression of CTCF and BORIS proteins throughout human epidermis. While CTCF is widely distributed within the nucleus, BORIS is confined to the nucleolus and other euchromatin domains. Nascent RNA experiments in primary keratinocytes revealed that endogenous BORIS is present in active transcription sites. Interestingly, BORIS also localises to interphase centrosomes suggesting a role in the cell cycle. Blocking the cell cycle at S phase or mitosis, or causing DNA damage, produced a striking accumulation of BORIS. Consistently, ectopic expression of wild type or GFP- BORIS provoked a higher rate of S phase cells as well as genomic instability by mitosis failure. Furthermore, downregulation of endogenous BORIS by specific shRNAs inhibited both RNA transcription and cell cycle progression. The results altogether suggest a role for BORIS in coordinating S phase events with mitosis

Long interspersed nuclear element-1 hypomethylation in cancer: biology and clinical applications

Epigenetic changes in long interspersed nuclear element-1s (LINE-1s or L1s) occur early during the process of carcinogenesis. A lower methylation level (hypomethylation) of LINE-1 is common in most cancers, and the methylation level is further decreased in more advanced cancers. Consequently, several previous studies have suggested the use of LINE-1 hypomethylation levels in cancer screening, risk assessment, tumor staging, and prognostic prediction. Epigenomic changes are complex, and global hypomethylation influences LINE-1s in a generalized fashion. However, the methylation levels of some loci are dependent on their locations. The consequences of LINE-1 hypomethylation are genomic instability and alteration of gene expression. There are several mechanisms that promote both of these consequences in cis. Therefore, the methylation levels of different sets of LINE-1s may represent certain phenotypes. Furthermore, the methylation levels of specific sets of LINE-1s may indicate carcinogenesis-dependent hypomethylation. LINE-1 methylation pattern analysis can classify LINE-1s into one of three classes based on the number of methylated CpG dinucleotides. These classes include hypermethylation, partial methylation, and hypomethylation. The number of partial and hypermethylated loci, but not hypomethylated LINE-1s, is different among normal cell types. Consequently, the number of hypomethylated loci is a more promising marker than methylation level in the detection of cancer DNA. Further genome-wide studies to measure the methylation level of each LINE-1 locus may improve PCR-based methylation analysis to allow for a more specific and sensitive detection of cancer DNA or for an analysis of certain cancer phenotypes