Assessing the co-variability of DNA methylation across peripheral cells and tissues: Implications for the interpretation of findings in epigenetic epidemiology

Most epigenome-wide association studies (EWAS) quantify DNA methylation (DNAm) in peripheral tissues such as whole blood to identify positions in the genome where variation is statistically associated with a trait or exposure. As whole blood comprises a mix of cell types, it is unclear whether trait-associated DNAm variation is specific to an individual cellular population. We collected three peripheral tissues (whole blood, buccal epithelial and nasal epithelial cells) from thirty individuals. Whole blood samples were subsequently processed using fluorescence-activated cell sorting (FACS) to purify five constituent cell-types (monocytes, granulocytes, CD4⁺ T cells, CD8⁺ T cells, and B cells). DNAm was profiled in all eight sample-types from each individual using the Illumina EPIC array. We identified significant differences in both the level and variability of DNAm between different sample types, and DNAm data-derived estimates of age and smoking were found to differ dramatically across sample types from the same individual. We found that for the majority of loci variation in DNAm in individual blood cell types was only weakly predictive of variance in DNAm measured in whole blood, although the proportion of variance explained was greater than that explained by either buccal or nasal epithelial samples. Covariation across sample types was much higher for DNAm sites influenced by genetic factors. Overall, we observe that DNAm variation in whole blood is additively influenced by a combination of the major blood cell types. For a subset of sites, however, variable DNAm detected in whole blood can be attributed to variation in a single blood cell type providing potential mechanistic insight about EWAS findings. Our results suggest that associations between whole blood DNAm and traits or exposures reflect differences in multiple cell types and our data will facilitate the interpretation of findings in epigenetic epidemiology.     

The population biology of transposable elements

A transposable element can be defined as a DNA sequence capable of moving to new sites in the genome. Such DNA sequences have been described in a wide range of organisms. The evolutionary processes affecting transposable elements can thus be divided into two categories: changes in sequence and changes in genomic location. As with other types of evolutionary change, the nature of the evolutionary process will be reflected in the extent and type of genetic variation existing in wild populations. Quantitative models of the evolution of transposable element sequences and positions will be outlined, and related to relevant data. The extent to which models designed to describe obvious transposable elements such as the mobile sequences of Drosophila are also applicable to interspersed repetitive DNAs from other species will be discussed.     

Cancer‐associated neochromosomes: a novel mechanism of oncogenesis

Malignant tumours are often characterised by significant rearrangement of the genome. This may be visible in the form of a deranged karyotype with both loss and gain of DNA sequences extending from chromosomal regions to whole chromosomes. In several tumour types, however, gross genomic derangements are minimal, and tumour cells contain one or more additional (supernumerary) chromosomes that may be unrecognisable in terms of a single origin. In this review we term such chromosomes cancer‐associated neochromosomes (CaNCs). In the absence of other identified genomic abnormalities, and because the CaNC is a common feature of the cancer type, it is hypothesised that the genetic alterations required for cell transformation are contained within its structure. In this review, we discuss the potential impact of modern genomic technologies on our understanding of the nature and causes of CaNC formation, which is central to several cancer types, exemplified here by well‐differentiated liposarcoma.     

Abundance of plastid DNA insertions in nuclear genomes of rice and Arabidopsis

Pairwise comparison of whole plastid and draft nuclear genomic sequences of Arabidopsis thaliana and Oryza sativa L. ssp. indica shows that rice nuclear genomic sequences contain homologs of plastid DNA covering about 94 kb (83%) of plastid genome and including one or more full-length intact (without mutations resulting in premature stop codons) homologues of 26 known protein-coding (KPC) plastid genes. By contrast, only about 20 kb (16%) of chloroplast DNA, including a single intact plastid-derived KPC gene, is presented in the nucleus of A. thaliana. Sixteen rice plastid genes have at least one nuclear copy without any mutation or with only synonymous substitutions. Nuclear copies for other ten plastid genes contain both synonymous and non-synonymous substitutions. Multiple ESTs for 25 out of 26 KPC genes were also found, as well as putative promoters for some of them. The study of substitutions pattern shows that some of nuclear homologues of plastid genes may be functional and/or are under the pressure of the positive natural selection. The similar comparative analysis performed on rice chromosome 1 revealed 27 contigs containing plastid-derived sequences, totalling about 84 kb and covering two thirds of chloroplast DNA, with the intact nuclear copies of 26 different KPC genes. One of these contigs, AP003280, includes almost 57 kb (45%) of chloroplast genome with the intact copies of 22 KPC genes. At the same time, we observed that relative locations of homologues in plastid DNA and the nuclear genome are significantly different.     

Genomic profiles for human peripheral blood T cells, B cells, natural killer cells, monocytes, and polymorphonuclear cells: comparisons to ischemic stroke, migraine, and Tourette syndrome

Blood genomic profiling has been applied to disorders of the blood and various organ systems including brain to elucidate disease mechanisms and identify surrogate disease markers. Since most studies have not examined specific cell types, we performed a preliminary genomic survey of major blood cell types from normal individuals using microarrays. CD4+ T cells, CD8+ T cells, CD19+ B cells, CD56+ natural killer cells, and CD14+ monocytes were negatively selected using the RosetteSep antibody cocktail, while polymorphonuclear leukocytes were separated with density gradient media. Genes differentially expressed by each cell type were identified. To demonstrate the potential use of such cell subtype-specific genomic expression data, a number of the major genes previously reported to be regulated in ischemic stroke, migraine, and Tourette syndrome are shown to be associated with distinct cell populations in blood. These specific gene expression, cell-type-related profiles will need to be confirmed in larger data sets and could be used to study these and many other neurological diseases.     

LINE-1 methylation status of endogenous DNA double-strand breaks

DNA methylation and the repair of DNA double-strand breaks (DSBs) are important processes for maintaining genomic integrity. Although DSBs can be produced by numerous agents, they also occur spontaneously as endogenous DSBs (EDSBs). In this study, we evaluated the methylation status of EDSBs to determine if there is a connection between DNA methylation and EDSBs. We utilized interspersed repetitive sequence polymerase chain reaction (PCR), ligation-mediated PCR and combined bisulfite restriction analysis to examine the extent of EDSBs and methylation at long interspersed nuclear element-1 (LINE-1) sequences nearby EDSBs. We tested normal white blood cells and several cell lines derived from epithelial cancers and leukemias. Significant levels of EDSBs were detectable in all cell types. EDSBs were also found in both replicating and non-replicating cells. We found that EDSBs contain higher levels of methylation than the cellular genome. This hypermethylation is replication independent and the methylation was present in the genome at the location prior to the DNA DSB. The differences in methylation levels between EDSBs and the rest of the genome suggests that EDSBs are differentially processed, by production, end-modification, or repair, depending on the DNA methylation status.     

Mitochondrial pseudogenes are pervasive and often insidious in the snapping shrimp genus Alpheus

Here we show that multiple DNA sequences, similar to the mitochondrial cytochrome oxidase I (COI) gene, occur within single individuals in at least 10 species of the snapping shrimp genus Alpheus. Cloning of amplified products revealed the presence of copies that differed in length and (more frequently) in base substitutions. Although multiple copies were amplified in individual shrimp from total genomic DNA (gDNA), only one sequence was amplified from cDNA. These results are best explained by the presence of nonfunctional duplications of a portion of the mtDNA, probably located in the nuclear genome, since transfer into the nuclear gene would render the COI gene nonfunctional due to differences in the nuclear and mitochondrial genetic codes. Analysis of codon variation suggests that there have been 21 independent transfer events in the 10 species examined. Within a single animal, differences between the sequences of these pseudogenes ranged from 0.2% to 20.6%, and those between the real mtDNA and pseudogene sequences ranged from 0.2% to 18.8% (uncorrected). The large number of integration events and the large range of divergences between pseudogenes and mtDNA sequences suggest that genetic material has been repeatedly transferred from the mtDNA to the nuclear genome of snapping shrimp. Unrecognized pseudogenes in phylogenetic or population studies may result in spurious results, although previous estimates of rates of molecular evolution based on Alpheus sister taxa separated by the Isthmus of Panama appear to remain valid. Especially worrisome for researchers are those pseudogenes that are not obviously recognizable as such. An effective solution may be to amplify transcribed copies of protein-coding mitochondrial genes from cDNA rather than using genomic DNA.     

Patterns of Genomic Integration of Nuclear Chloroplast DNA Fragments in Plant Species

The transfer of organelle DNA fragments to the nuclear genome is frequently observed in eukaryotes. These transfers are thought to play an important role in gene and genome evolution of eukaryotes. In plants, such transfers occur from plastid to nuclear [nuclear plastid DNAs (NUPTs)] and mitochondrial to nuclear (nuclear mitochondrial DNAs) genomes. The amount and genomic organization of organelle DNA fragments have been studied in model plant species, such as Arabidopsis thaliana and rice. At present, publicly available genomic data can be used to conduct such studies in non-model plants. In this study, we analysed the amount and genomic organization of NUPTs in 17 plant species for which genome sequences are available. The amount and distribution of NUPTs varied among the species. We also estimated the distribution of NUPTs according to the time of integration (relative age) by conducting sequence similarity analysis between NUPTs and the plastid genome. The age distributions suggested that the present genomic constitutions of NUPTs could be explained by the combination of the rapidly eliminated deleterious parts and few but constantly existing less deleterious parts.     

Genomic Treasure Troves: Complete Genome Sequencing of Herbarium and Insect Museum Specimens

Unlocking the vast genomic diversity stored in natural history collections would create unprecedented opportunities for genome-scale evolutionary, phylogenetic, domestication and population genomic studies. Many researchers have been discouraged from using historical specimens in molecular studies because of both generally limited success of DNA extraction and the challenges associated with PCR-amplifying highly degraded DNA. In today''s next-generation sequencing (NGS) world, opportunities and prospects for historical DNA have changed dramatically, as most NGS methods are actually designed for taking short fragmented DNA molecules as templates. Here we show that using a standard multiplex and paired-end Illumina sequencing approach, genome-scale sequence data can be generated reliably from dry-preserved plant, fungal and insect specimens collected up to 115 years ago, and with minimal destructive sampling. Using a reference-based assembly approach, we were able to produce the entire nuclear genome of a 43-year-old Arabidopsis thaliana (Brassicaceae) herbarium specimen with high and uniform sequence coverage. Nuclear genome sequences of three fungal specimens of 22–82 years of age (Agaricus bisporus, Laccaria bicolor, Pleurotus ostreatus) were generated with 81.4–97.9% exome coverage. Complete organellar genome sequences were assembled for all specimens. Using de novo assembly we retrieved between 16.2–71.0% of coding sequence regions, and hence remain somewhat cautious about prospects for de novo genome assembly from historical specimens. Non-target sequence contaminations were observed in 2 of our insect museum specimens. We anticipate that future museum genomics projects will perhaps not generate entire genome sequences in all cases (our specimens contained relatively small and low-complexity genomes), but at least generating vital comparative genomic data for testing (phylo)genetic, demographic and genetic hypotheses, that become increasingly more horizontal. Furthermore, NGS of historical DNA enables recovering crucial genetic information from old type specimens that to date have remained mostly unutilized and, thus, opens up a new frontier for taxonomic research as well.     

Chromatin remodeling and the maintenance of genome integrity

DNA damage of any type is threatening for a cell. If lesions are left unrepaired, genomic instability can arise, faithful transmission of genetic information is greatly compromised eventually leading the cell to undergo apoptosis or carcinogenesis. In order to access/detect and repair these damages, repair factors must circumvent the natural repressive barrier of chromatin. This review will present recent progress showing the intricate link between chromatin, its remodeling and the DNA repair process. Several studies demonstrated that one of the first events following specific types of DNA damage is the phosphorylation of histone H2A. This mark or the damage itself are responsible for the association of chromatin-modifying complexes near damaged DNA. These complexes are able to change the chromatin structure around the wounded DNA in order to allow the repair machinery to gain access and repair the lesion. Chromatin modifiers include ATP-dependent remodelers such as SWI/SNF and Rad54 as well as histone acetyltransferases (HATs) like SAGA/NuA4-related complexes and p300/CBP, which have been shown to facilitate DNA accessibility and repair in different pathways leading to the maintenance of genome integrity.     

Transformation of primary rat kidney cells by DNA fragments of weakly oncogenic adenoviruses.

Primary cultures of baby rat kidney (BRK) cells were transformed by intact DNA and DNA fragments of weakly oncogenic human adenovirus types 3 and 7. The smallest fragment found to contain transforming activity was the left-terminal 4% endo R.HindIII fragment (for both adenovirus type 3 and 7 DNAs). The efficiency of transformation of this fragment was low, and no permanent cell line could be established. Left-terminal fragments ranging from 84 to 4,5% of the viral genome could all transform BRK cells with the same efficiency as intact viral DNA. A number of adenovirus type 7 DNA fragment-transformed lines were established and were found to contain persistent viral DNA sequences and adenovirus subgroup B-specific T antigen. Consequently, the transforming functions of adenovirus types 3 and 7 are located at the extreme left-hand end of the genome, and the minimum size for a DNA fragment with transforming activity is 1.0 X 10(6) daltons. These results do not rule out the possibility that viral genes located outside the transforming region may also influence transformation.     

Cloned microsatellite repeats differ between 4-base restriction endonucleases.

Simple sequence repeat (SSR) loci are an important marker type for population genetic studies despite the limitation that development of novel loci requires construction and screening of genomic DNA libraries. The common practice of size fractioning genomic DNA before cloning could lead to differential representation of SSR loci within genomic libraries. In addition, linkage mapping studies have shown that small numbers of SSR markers are not randomly distributed within the genomes from which they are isolated. From attempts to clone five SSR repeat sequences in two wild plant species we show that the numbers and repeat type of potential SSR markers depend on the restriction endonuclease used to sample the genome when constructing DNA libraries. This observation is consistent with unequal sampling of the genome by different restriction enzymes. However, as a group the five SSR repeat sequences are not associated with a given restriction enzyme, suggesting they are not clumped within the genome. Use of multiple restriction enzymes to construct DNA libraries may help ensure that cloned SSR loci are drawn from diverse locations in the genome, helping to meet the assumption of randomly located marker loci required for population genetic inferences.     

Differential subtraction chain, a method for identifying differences in genomic DNA and mRNA.

Identifying the genetic differences between two organisms or cell types has been a major goal in modern biomedical research. Recently, we developed a novel methodology that can rapidly identify the differences between two populations of DNA. This method, termed 'differential subtraction chain' (DSC), is based on a novel 'negative amplification' strategy that converts (amplifiable) tester sequences to counterpart (unamplifiable) driver sequences. The result is a double exponential elimination of amplifiable sequences in the testers, while preserving the sequences in the testers that have no counterpart in the drivers. We applied this methodology to the genome of a glioblastoma cell line. A homozygous deletion was rapidly identified. We extended this technique to identifying the unique sequences in mRNA. Two CDC25 transgene fragments were quickly identified in a cdc25B transgenic mouse. We also applied this methodology to systems with profound differences in mRNA expression. In a 'prostate epithelia subtracting blood cells' DSC reaction, a sample of unique gene fragments which are absent in the prostate but present in the blood were identified. Lastly, we detected rare (1 virus/100 cells) Herpes simplex virus type 2 (HSV-2) sequences in a tissue culture, indicating good sensitivity of this methodology. Overall, DSC represents a fast, efficient and sensitive method for identifying differences in genomic DNA and mRNA and can be easily applied in a variety of biological systems.     

Simple sequence repeats for genetic studies of alpaca

Trace sequences from the 2X alpaca genome sequencing effort were examined to identify simple sequence repeats (microsatellites) for genetic studies. A total of 6,685 repeat-containing sequences were downloaded from GenBank, processed, and assembled into contigs representing an estimated 4,278 distinct sequences. This sequence set contained 2,290 sequences of length > 100 nucleotides that contained microsatellites of length > or = 14 dinucleotide or 10 trinucleotide repeats with purity equal to 100%. An additional 13 sequences contained a GC microsatellite of length > or = 12 repeats (purity = 100%) were also obtained. Primer pairs for amplification of 1,516 putative loci are presented. Amplification of genomic DNA from alpaca and llama by PCR was demonstrated for 14 primer sets including one from each of the microsatellite repeat types. Comparative chromosomal location for the alpaca markers was predicted in the bovine genome by BLAT searches against assembly 4.0 of the bovine whole genome sequence. A total of 634 markers (41.8%) returned BLAT hits with score > 100 and Identity > 85%, with the majority assignable to unique locations. We show that microsatellites are abundant and easily identified within the alpaca genome sequence. These markers will provide a valuable resource for further genetic studies of the alpaca and related species.     

Patterns of DNA methylation are indistinguishable in different individuals over a wide range of human DNA sequences.

Patterns of DNA methylation at 5'-CCGG-3' and 5'-GCGC-3' sequences were determined in about 570 kb, equivalent to about 0.02% of the human genome, by using HpaII and HhaI restriction endonucleases, respectively, and randomly selected cosmid clones of human DNA as hybridization probes. Many of these human DNA sequences were of the repetitive type. The DNAs from human lymphocytes, from a mixture of all blood cells or from several established human cell lines (HeLa, KB, 293, or DEV) were included in these analyses. In the segments of the human genome investigated, the patterns of DNA methylation were characterized by often completely or partly methylated 5'-CCGG-3' or by partly methylated 5'-GCGC-3' sequences. Even among individuals of different genetic origins (East-Asian or Caucasian), these patterns of DNA methylation proved indistinguishable by the method applied. The cytokine-dependent stimulation of human lymphocytes to replicate in culture did not affect the stability of these patterns. In the same DNA sequences from several human cell lines, much lower levels of DNA methylation were observed. In human cell lines some of the investigated sequences were unmethylated. The results presented lend credence to the notion that the human genome exhibits highly cell type-specific patterns of DNA methylation which are often indistinguishable among different individuals even of different genetic backgrounds.     

Efficient procedure for purification of obligate intracellular Wolbachia pipientis and representative amplification of its genome by multiple-displacement amplification

Bacteria belonging to the genus Wolbachia are obligatory microendocytobionts that infect a variety of arthropods and a majority of filarial nematode species, where they induce reproductive alterations or establish a mutualistic symbiosis. Although two whole genome sequences of Wolbachia pipientis, for strain wMel from Drosophila melanogaster and strain wBm from Brugia malayi, have been fully completed and six other genome sequencing projects are ongoing (http://www.genomesonline.org/index.cgi?want=Prokaryotic+Ongoin), genetic analyses of these bacteria are still scarce, mainly due to the inability to cultivate them outside of eukaryotic cells. Usually, a large amount of host tissue (a thousand individuals, or about 10 g) is required in order to purify Wolbachia and extract its DNA, which is often recovered in small amounts and contaminated by host cell DNA, thus hindering genomic studies. In this report, we describe an efficient and reliable procedure to representatively amplify the Wolbachia genome by multiple-displacement amplification from limited infected host tissue (0.2 g or 2 x 10(7) cells). We obtained sufficient amounts (8 to 10 microg) of DNA of suitable quality for genomic studies, and we demonstrated that the amplified DNA contained all of the Wolbachia loci targeted. In addition, our data indicated that the genome of strain wRi, an obligatory endosymbiont of Drosophila simulans, shares a similar overall architecture with its relative strain wMel.     

Three-dimensional positioning of genes in mouse cell nuclei

To understand the regulation of the genome, it is necessary to understand its three-dimensional organization in the nucleus. We investigated the positioning of eight gene loci on four different chromosomes, including the β-globin gene, in mouse embryonic stem cells and in in vitro differentiated macrophages by fluorescence in situ hybridization on structurally preserved nuclei, confocal microscopy, and 3D quantitative image analysis. We found that gene loci on the same chromosome can significantly differ from each other and from their chromosome territory in their average radial nuclear position. Radial distribution of a given gene locus can change significantly between cell types, excluding the possibility that positioning is determined solely by the DNA sequence. For the set of investigated gene loci, we found no relationship between radial distribution and local gene density, as it was described for human cell nuclei. We did find, however, correlation with other genomic properties such as GC content and certain repetitive elements such as long terminal repeats or long interspersed nuclear elements. Our results suggest that gene density itself is not a driving force in nuclear positioning. Instead, we propose that other genomic properties play a role in determining nuclear chromatin distribution. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.