Application of the LUminometric Methylation Assay to ecological species: tissue quality requirements and a survey of DNA methylation levels in animals

The LUminometric Methylation Assay (LUMA) measures global DNA methylation. LUMA depends on digestion of DNA with methyl‐sensitive and methyl‐insensitive restriction enzymes, followed by pyrosequencing. Until recently, LUMA has been principally used for biomedical research. Here, we use chickens as a model to investigate sample quality issues relating to LUMA and then apply the method to ecological species. First, we assessed the effect of tissue storage conditions on DNA methylation values. This is an important consideration for ecological species because samples are not always ideally preserved and LUMA is sensitive to poor DNA quality. We found that good quality LUMA data could be obtained from chicken liver and brain tissues stored at 21 °C for at least 2 and 12 h, respectively. Longer storage times introduced nonspecific peaks to pyrograms which were associated with reduced DNA methylation. Repeatedly, freezing and thawing the tissues did not affect LUMA data. Second, we measured DNA methylation in 12 species representing five animal classes: amphibians (African and Western clawed frog), reptiles (green anole lizard), fish (yellow perch, goldfish, lake trout), mammals (American mink, polar bear, short‐beaked common dolphin, Atlantic white‐sided dolphin) and birds (chicken, Japanese quail). We saw a pattern of high DNA methylation in fish (84–87%), and intermediate levels in mammals (68–72%) and birds (52–71%). This pattern corresponds well with previous measures of DNA methylation generated by HPLC. Our data represent the first CpG methylation values to be reported in several species and provide a basis for studying patterns of epigenetic inheritance in an ecological context.     

LUMA (LUminometric Methylation Assay)--a high throughput method to the analysis of genomic DNA methylation

Changes in genomic DNA methylation are recognized as important events in normal and pathological cellular processes, contributing both to normal development and differentiation as well as cancer and other diseases. Here, we report a novel method to estimate genome-wide DNA methylation, referred to as LUminometric Methylation Assay (LUMA). The method is based on combined DNA cleavage by methylation-sensitive restriction enzymes and polymerase extension assay by Pyrosequencing. The method is quantitative, highly reproducible and easy to scale up. Since no primary modification of genomic DNA, such as bisulfite treatment, is needed, the total assay time is only 6 h. In addition, the assay requires only 200-500 ng of genomic DNA and incorporates an internal control to eliminate the problem of varying amounts of starting DNA. The accuracy and linearity of LUMA were verified by in vitro methylated lambda DNA. In addition, DNA methylation levels were assessed by LUMA in DNA methyltransferase knock-out cell lines and after treatment with the DNA methyltransferase inhibitor (5-AzaCytidine). The LUMA assay may provide a useful method to analyze genome-wide DNA methylation for a variety of physiological and pathological conditions including etiologic, diagnostic and prognostic aspects of cancer.     

Methyl CpG Binding Domain Ultra‐Sequencing: a novel method for identifying inter‐individual and cell‐type‐specific variation in DNA methylation

Experience‐dependent changes in DNA methylation can exert profound effects on neuronal function and behaviour. A single learning event can induce a variety of DNA modifications within the neuronal genome, some of which may be common to all individuals experiencing the event, whereas others may occur in a subset of individuals. Variations in experience‐induced DNA methylation may subsequently confer increased vulnerability or resilience to the development of neuropsychiatric disorders. However, the detection of experience‐dependent changes in DNA methylation in the brain has been hindered by the interrogation of heterogeneous cell populations, regional differences in epigenetic states and the use of pooled tissue obtained from multiple individuals. Methyl CpG Binding Domain Ultra‐Sequencing (MBD Ultra‐Seq) overcomes current limitations on genome‐wide epigenetic profiling by incorporating fluorescence‐activated cell sorting and sample‐specific barcoding to examine cell‐type‐specific CpG methylation in discrete brain regions of individuals. We demonstrate the value of this method by characterizing differences in 5‐methylcytosine (5mC) in neurons and non‐neurons of the ventromedial prefrontal cortex of individual adult C57BL/6 mice, using as little as 50 ng of genomic DNA per sample. We find that the neuronal methylome is characterized by greater CpG methylation as well as the enrichment of 5mC within intergenic loci. In conclusion, MBD Ultra‐Seq is a robust method for detecting DNA methylation in neurons derived from discrete brain regions of individual animals. This protocol will facilitate the detection of experience‐dependent changes in DNA methylation in a variety of behavioural paradigms and help identify aberrant experience‐induced DNA methylation that may underlie risk and resiliency to neuropsychiatric disease.     

Correlations in global DNA methylation measures in peripheral blood mononuclear cells and granulocytes

Alterations in global DNA methylation levels have been associated with chronic diseases. Despite the increase in the number of studies measuring markers of global methylation, few have adequately examined within-individual differences by source of DNA and whether within-individual differences by source of DNA differ by age, race and other lifestyle factors. We examined correlations between peripheral mononuclear cell (PBMC) and granulocyte DNA methylation levels measured by the luminometric methylation assay (LUMA), and in LINE-1, Sat2, and Alu by MethyLight and pyrosequencing, in the same individual in 112 women participating in The New York City Multiethnic Breast Cancer Project. Levels of DNA methylation of Sat2 by MethyLight (r = 0.57; P < 0.01) and LINE-1 by pyrosequencing (r = 0.30; P < 0.01) were correlated between PBMC and granulocyte DNA of the same individuals, but LUMA and Alu levels were not. The magnitude of the correlations for Sat2 and LINE-1 varied when stratified by selected demographic and lifestyle factors, although the study sample size limited our comparisons across subgroups. These results lend further support to the importance of considering the source of DNA in epidemiologic studies of white blood cell DNA methylation. Results from studies that combine individuals with different available DNA sources need to be interpreted with caution.     

The methylation patterns of the IGF2 and IGF2R genes in bovine spermatozoa are not affected by flow-cytometric sex sorting

The objectives of this study were to investigate the effect of sexing by flow cytometry on the methylation patterns of the IGF2 and IGF2R genes. Frozen‐thawed, unsorted, and sex‐sorted sperm samples from four Nellore bulls were used. Each ejaculate was separated into three fractions: non‐sexed (NS), sexed for X‐sperm (SX), and sexed for Y‐sperm (SY). Sperm were isolated from the extender, cryoprotectant, and other cell types by centrifugation on a 40:70% Percoll gradient, and sperm pellets were used for genomic DNA isolation. DNA was used for analyses of the methylation patterns by bisulfite sequencing. Methylation status of the IGF2 and IGF2R genes were evaluated by sequencing 195 and 147 individual clones, respectively. No global differences in DNA methylation were found between NS, SX, and SY groups for the IGF2 (P = 0.09) or IGF2R genes (P = 0.38). Very specific methylation patterns were observed in the 25th and 26th CpG sites in the IGF2R gene. representing higher methylation in NS than in the SX and SY groups compared with the other CpG sites. Further, individual variation in methylation patterns was found among bulls. In conclusion, the sex‐sorting procedure by flow cytometry did not affect the overall DNA methylation patterns of the IGF2 and IGF2R genes, although individual variation in their methylation patterns among bulls was observed. Mol. Reprod. Dev. 79:77–84, 2012. © 2011 Wiley Periodicals, Inc.     

Sexually divergent DNA methylation patterns with hippocampal aging

DNA methylation is a central regulator of genome function, and altered methylation patterns are indicative of biological aging and mortality. Age‐related cellular, biochemical, and molecular changes in the hippocampus lead to cognitive impairments and greater vulnerability to neurodegenerative disease that varies between the sexes. The role of hippocampal epigenomic changes with aging in these processes is unknown as no genome‐wide analyses of age‐related methylation changes have considered the factor of sex in a controlled animal model. High‐depth, genome‐wide bisulfite sequencing of young (3 month) and old (24 month) male and female mouse hippocampus revealed that while total genomic methylation amounts did not change with aging, specific sites in CG and non‐CG (CH) contexts demonstrated age‐related increases or decreases in methylation that were predominantly sexually divergent. Differential methylation with age for both CG and CH sites was enriched in intergenic and intronic regions and under‐represented in promoters, CG islands, and specific enhancer regions in both sexes, suggesting that certain genomic elements are especially labile with aging, even if the exact genomic loci altered are predominantly sex‐specific. Lifelong sex differences in autosomal methylation at CG and CH sites were also observed. The lack of genome‐wide hypomethylation, sexually divergent aging response, and autosomal sex differences at CG sites was confirmed in human data. These data reveal sex as a previously unappreciated central factor of hippocampal epigenomic changes with aging. In total, these data demonstrate an intricate regulation of DNA methylation with aging by sex, cytosine context, genomic location, and methylation level.     

Flowering and dwarfism induced by DNA demethylation in Pharbitis nil

Flowering and dwarfism induced by 5‐azacytidine and zebularine, which both cause DNA demethylation, were studied in a short‐day (SD) plant Pharbitis nil (synonym Ipomoea nil), var. Violet whose photoinduced flowering state does not last for a long period of time. The DNA demethylating reagents induced flowering under non‐inductive long‐day (LD) conditions. The flower‐inducing effect of 5‐azacytidine did not last for a long period of time, and the plants reverted to vegetative growth. The progeny of the plants that were induced to flower by DNA demethylation did not flower under the non‐inductive photoperiodic conditions. These results suggest that the flowering‐related genes were activated by DNA demethylation and then remethylated again in the progeny. The DNA demethylation also induced dwarfism. The dwarfism did not last for a long period of time, was not heritable and was overcome by gibberellin A₃ but not by t‐zeatin or kinetin. The change in the genome‐wide methylation state was examined by methylation‐sensitive amplified fragment length polymorphism (MS‐AFLP) analysis. The analysis detected many more polymorphic fragments between the DNA samples isolated from the cotyledons treated with SD than from the cotyledons under LD conditions, indicating that the DNA methylation state was altered by photoperiodic conditions. Seven LD‐specific fragments were extracted from the gel of the MS‐AFLP and were sequenced. One of these fragments was highly homologous with the genes encoding ribosomal proteins.     

Methylome Analysis and Epigenetic Changes Associated with Menarcheal Age

Reproductive factors have been linked to both breast cancer and DNA methylation, suggesting methylation as an important mechanism by which reproductive factors impact on disease risk. However, few studies have investigated the link between reproductive factors and DNA methylation in humans. Genome-wide methylation in peripheral blood lymphocytes of 376 healthy women from the prospective EPIC study was investigated using LUminometric Methylation Assay (LUMA). Also, methylation of 458877 CpG sites was additionally investigated in an independent group of 332 participants of the EPIC-Italy sub-cohort, using the Infinium HumanMethylation 450 BeadChip. Multivariate logistic regression and linear models were used to investigate the association between reproductive risk factors and genome wide and CpG-specific DNA methylation, respectively. Menarcheal age was inversely associated with global DNA methylation as measured with LUMA. For each yearly increase in age at menarche, the risk of having genome wide methylation below median level was increased by 32% (OR:1.32, 95%CI:1.14–1.53). When age at menarche was treated as a categorical variable, there was an inverse dose-response relationship with LUMA methylation levels (OR_{12–14vs.≤11 yrs}:1.78, 95%CI:1.01–3.17 and OR_{≥15vs.≤11 yrs}:4.59, 95%CI:2.04–10.33; P for trend<0.0001). However, average levels of global methylation as measured by the Illumina technology were not significantly associated with menarcheal age. In locus by locus comparative analyses, only one CpG site had significantly different methylation depending on the menarcheal age category examined, but this finding was not replicated by pyrosequencing in an independent data set. This study suggests a link between age at menarche and genome wide DNA methylation, and the difference in results between the two arrays suggests that repetitive element methylation has a role in the association. Epigenetic changes may be modulated by menarcheal age, or the association may be a mirror of other important changes in early life that have a detectable effect on both methylation levels and menarcheal age.     

BREX is a novel phage resistance system widespread in microbial genomes

The perpetual arms race between bacteria and phage has resulted in the evolution of efficient resistance systems that protect bacteria from phage infection. Such systems, which include the CRISPR‐Cas and restriction‐modification systems, have proven to be invaluable in the biotechnology and dairy industries. Here, we report on a six‐gene cassette in Bacillus cereus which, when integrated into the Bacillus subtilis genome, confers resistance to a broad range of phages, including both virulent and temperate ones. This cassette includes a putative Lon‐like protease, an alkaline phosphatase domain protein, a putative RNA‐binding protein, a DNA methylase, an ATPase‐domain protein, and a protein of unknown function. We denote this novel defense system BREX (Bacteriophage Exclusion) and show that it allows phage adsorption but blocks phage DNA replication. Furthermore, our results suggest that methylation on non‐palindromic TAGGAG motifs in the bacterial genome guides self/non‐self discrimination and is essential for the defensive function of the BREX system. However, unlike restriction‐modification systems, phage DNA does not appear to be cleaved or degraded by BREX, suggesting a novel mechanism of defense. Pan genomic analysis revealed that BREX and BREX‐like systems, including the distantly related Pgl system described in Streptomyces coelicolor, are widely distributed in ~10% of all sequenced microbial genomes and can be divided into six coherent subtypes in which the gene composition and order is conserved. Finally, we detected a phage family that evades the BREX defense, implying that anti‐BREX mechanisms may have evolved in some phages as part of their arms race with bacteria.     

DNA methylation stability in fish spermatozoa upon external constraint: Impact of fish hormonal stimulation and sperm cryopreservation

Highly differentiated mature spermatozoa carry not only genetic but also epigenetic information that is to be transmitted to the embryo. DNA methylation is one epigenetic actor associated with sperm nucleus compaction, gene silencing, and prepatterning of embryonic gene expression. Therefore, the stability of this mark toward reproductive biotechnologies is a major issue in animal production. The present work explored the impact of hormonal induction of spermiation and sperm cryopreservation in two cyprinids, the goldfish (Carassius auratus) and the zebrafish (Danio rerio), using LUminometric Methylation Assay (LUMA). We showed that while goldfish hormonal treatment did increase sperm production, it did not alter global DNA methylation of spermatozoa. Different sperm samples repeatedly collected from the same males for 2 months also showed the same global DNA methylation level. Similarly, global DNA methylation was not affected after cryopreservation of goldfish spermatozoa with methanol, whereas less efficient cryoprotectants (dimethylsulfoxide and 1,2‐propanediol) decreased DNA methylation. In contrast, cryopreservation of zebrafish spermatozoa with methanol induced a slight, but significant, increase in global DNA methylation. In the less compact nuclei, that is, goldfish fin somatic cells, cryopreservation did not change global DNA methylation regardless of the choice of cryoprotectant. To conclude, global DNA methylation is a robust parameter with respect to biotechnologies such as hormonal induction of spermiation and sperm cryopreservation, but it can be altered when the best sperm manipulation conditions are not met.     

Correlations in global DNA methylation measures in peripheral blood mononuclear cells and granulocytes

Alterations in global DNA methylation levels have been associated with chronic diseases. Despite the increase in the number of studies measuring markers of global methylation, few have adequately examined within-individual differences by source of DNA and whether within-individual differences by source of DNA differ by age, race and other lifestyle factors. We examined correlations between peripheral mononuclear cell (PBMC) and granulocyte DNA methylation levels measured by the luminometric methylation assay (LUMA), and in LINE-1, Sat2, and Alu by MethyLight and pyrosequencing, in the same individual in 112 women participating in The New York City Multiethnic Breast Cancer Project. Levels of DNA methylation of Sat2 by MethyLight (r = 0.57; P < 0.01) and LINE-1 by pyrosequencing (r = 0.30; P < 0.01) were correlated between PBMC and granulocyte DNA of the same individuals, but LUMA and Alu levels were not. The magnitude of the correlations for Sat2 and LINE-1 varied when stratified by selected demographic and lifestyle factors, although the study sample size limited our comparisons across subgroups. These results lend further support to the importance of considering the source of DNA in epidemiologic studies of white blood cell DNA methylation. Results from studies that combine individuals with different available DNA sources need to be interpreted with caution.     

DNA methylation of retrotransposons,DNA transposons and genes in sugar beet (Beta vulgaris L.)

The methylation of cytosines shapes the epigenetic landscape of plant genomes, coordinates transgenerational epigenetic inheritance, represses the activity of transposable elements (TEs), affects gene expression and, hence, can influence the phenotype. Sugar beet (Beta vulgaris ssp. vulgaris), an important crop that accounts for 30% of worldwide sugar needs, has a relatively small genome size (758 Mbp) consisting of approximately 485 Mbp repetitive DNA (64%), in particular satellite DNA, retrotransposons and DNA transposons. Genome‐wide cytosine methylation in the sugar beet genome was studied in leaves and leaf‐derived callus with a focus on repetitive sequences, including retrotransposons and DNA transposons, the major groups of repetitive DNA sequences, and compared with gene methylation. Genes showed a specific methylation pattern for CG, CHG (H = A, C, and T) and CHH sites, whereas the TE pattern differed, depending on the TE class (class 1, retrotransposons and class 2, DNA transposons). Along genes and TEs, CG and CHG methylation was higher than that of adjacent genomic regions. In contrast to the relatively low CHH methylation in retrotransposons and genes, the level of CHH methylation in DNA transposons was strongly increased, pointing to a functional role of asymmetric methylation in DNA transposon silencing. Comparison of genome‐wide DNA methylation between sugar beet leaves and callus revealed a differential methylation upon tissue culture. Potential epialleles were hypomethylated (lower methylation) at CG and CHG sites in retrotransposons and genes and hypermethylated (higher methylation) at CHH sites in DNA transposons of callus when compared with leaves.     

Comprehensive definition of genome features in Spirodela polyrhiza by high‐depth physical mapping and short‐read DNA sequencing strategies

Spirodela polyrhiza is a fast‐growing aquatic monocot with highly reduced morphology, genome size and number of protein‐coding genes. Considering these biological features of Spirodela and its basal position in the monocot lineage, understanding its genome architecture could shed light on plant adaptation and genome evolution. Like many draft genomes, however, the 158‐Mb Spirodela genome sequence has not been resolved to chromosomes, and important genome characteristics have not been defined. Here we deployed rapid genome‐wide physical maps combined with high‐coverage short‐read sequencing to resolve the 20 chromosomes of Spirodela and to empirically delineate its genome features. Our data revealed a dramatic reduction in the number of the rDNA repeat units in Spirodela to fewer than 100, which is even fewer than that reported for yeast. Consistent with its unique phylogenetic position, small RNA sequencing revealed 29 Spirodela‐specific microRNA, with only two being shared with Elaeis guineensis (oil palm) and Musa balbisiana (banana). Combining DNA methylation data and small RNA sequencing enabled the accurate prediction of 20.5% long terminal repeats (LTRs) that doubled the previous estimate, and revealed a high Solo:Intact LTR ratio of 8.2. Interestingly, we found that Spirodela has the lowest global DNA methylation levels (9%) of any plant species tested. Taken together our results reveal a genome that has undergone reduction, likely through eliminating non‐essential protein coding genes, rDNA and LTRs. In addition to delineating the genome features of this unique plant, the methodologies described and large‐scale genome resources from this work will enable future evolutionary and functional studies of this basal monocot family.     

Genome‐wide DNA methylation analysis using next‐generation sequencing to reveal candidate genes responsible for boar taint in pigs

Boar taint (BT) is an offensive flavor observed in non‐castrated male pigs that reduces the carcass price. Surgical castration effectively avoids the taint but is associated with animal welfare concerns. The functional annotation of farm animal genomes for understanding the biology of complex traits can be used in the selection of breeding animals to achieve favorable phenotypic outcomes. The characterization of pig epigenomes/methylation changes between animals with high and low BT and genome‐wide epigenetic markers that can predict BT are lacking. Reduced representation bisulfite sequencing of DNA methylation patterns based on next‐generation sequencing is an efficient technology to identify candidate epigenetic biomarkers associated with BT. Three different BT levels were analyzed using reduced representation bisulfite sequencing data to calculate the methylation levels of cytosine and guanine dinucleotide (CpG) sites. The co‐analysis of differentially methylated CpG sites identified by this study and differentially expressed genes identified by a previous study found 32 significant co‐located genes. The joint analysis of GO terms and pathways revealed that methylation and gene expression of seven candidate genes were associated with BT; in particular, FASN plays a key role in fatty acid biosynthesis, and PEMT might be involved in estrogen regulation and the development of BT. This study is the first to report the genome‐wide DNA methylation profiles of BT in pigs using next‐generation sequencing and summarize candidate genes associated with epigenetic markers of BT, which could contribute to the understanding of the functional biology of BT traits and selective breeding of pigs against BT based on epigenetic biomarkers.