Direct sequencing of haplotypes from diploid individuals through a modified emulsion PCR‐based single‐molecule sequencing approach

While standard DNA‐sequencing approaches readily yield genotypic sequence data, haplotype information is often of greater utility for population genetic analyses. However, obtaining individual haplotype sequences can be costly and time‐consuming and sometimes requires statistical reconstruction approaches that are subject to bias and error. Advancements have recently been made in determining individual chromosomal sequences in large‐scale genomic studies, yet few options exist for obtaining this information from large numbers of highly polymorphic individuals in a cost‐effective manner. As a solution, we developed a simple PCR‐based method for obtaining sequence information from individual DNA strands using standard laboratory equipment. The method employs a water‐in‐oil emulsion to separate the PCR mixture into thousands of individual microreactors. PCR within these small vesicles results in amplification from only a single starting DNA template molecule and thus a single haplotype. We improved upon previous approaches by including SYBR Green I and a melted agarose solution in the PCR, allowing easy identification and separation of individually amplified DNA molecules. We demonstrate the use of this method on a highly polymorphic estuarine population of the copepod Eurytemora affinis for which current molecular and computational methods for haplotype determination have been inadequate.     

Evolutionary maintenance of genomic diversity within arbuscular mycorrhizal fungi

Most organisms are built from a single genome. In striking contrast, arbuscular mycorrhizal fungi appear to maintain genomic variation within an individual fungal network. Arbuscular mycorrhizal fungi dwell in the soil, form mutualistic networks with plants, and bear multiple, potentially genetically diverse nuclei within a network. We explore, from a theoretical perspective, why such genetic diversity might be maintained within individuals. We consider selection acting within and between individual fungal networks. We show that genetic diversity could provide a benefit at the level of the individual, by improving growth in variable environments, and that this can stabilize genetic diversity even in the presence of nuclear conflict. Arbuscular mycorrhizal fungi complicate our understanding of organismality, but our findings offer a way of understanding such biological anomalies.     

Improved sensitivity in comparative genomic hybridization analysis of DNA heteroploid cell mixtures after pre-enrichment of subpopulations by fluorescence activated cell sorting.

Cytogenetic analysis of solid tumors with comparative genomic hybridization (CGH) is hampered by the dilution of DNA from individual tumor subpopulations with DNA from other cells. We investigated to what extent this dilution effect can be alleviated using fluorescence activated cell sorting (flow sorting) of experimental DNA heteroploid cell mixtures prior to CGH. From mixtures of normal lymphocytes with triploid K-562 cells the individual components were sorted according to stemline DNA content and processed by CGH in comparison with pure K-562 samples and the original mixtures. Compared with 30 autosome copy number imbalances found in pure K-562 samples, a mixture with 32% K-562 cells showed 16 imbalances, and none were detected in mixtures with 13% or 5% K-562 cells. In contrast, 29, 22 and 23 imbalances were detected in K-562 nuclei sorted from the 32%, 13% and 5% mixtures, respectively. This indicate that CGH analysis of flow sorted DNA aneuploid subpopulations enables a specific cytogenetic analysis of the individual subclones in a DNA heteroploid cell population.     

Allelic Skewing of DNA Methylation Is Widespread across the Genome

DNA methylation is assumed to be complementary on both alleles across the genome, although there are exceptions, notably in regions subject to genomic imprinting. We present a genome-wide survey of the degree of allelic skewing of DNA methylation with the aim of identifying previously unreported differentially methylated regions (DMRs) associated primarily with genomic imprinting or DNA sequence variation acting in cis. We used SNP microarrays to quantitatively assess allele-specific DNA methylation (ASM) in amplicons covering 7.6% of the human genome following cleavage with a cocktail of methylation-sensitive restriction enzymes (MSREs). Selected findings were verified using bisulfite-mapping and gene-expression analyses, subsequently tested in a second tissue from the same individuals, and replicated in DNA obtained from 30 parent-child trios. Our approach detected clear examples of ASM in the vicinity of known imprinted loci, highlighting the validity of the method. In total, 2,704 (1.5%) of our 183,605 informative and stringently filtered SNPs demonstrate an average relative allele score (RAS) change ≥0.10 following MSRE digestion. In agreement with previous reports, the majority of ASM (∼90%) appears to be cis in nature, and several examples of tissue-specific ASM were identified. Our data show that ASM is a widespread phenomenon, with >35,000 such sites potentially occurring across the genome, and that a spectrum of ASM is likely, with heterogeneity between individuals and across tissues. These findings impact our understanding about the origin of individual phenotypic differences and have implications for genetic studies of complex disease.     

Formation of homogeneous carbohydrate-lectin cross-linked precipitates from mixtures of D-galactose/N-acetyl-D-galactosamine-specific lectins and multiantennary galactosyl carbohydrates.

Quantitative precipitation studies have shown that the Man/Glc-specific lectin concanavalin A (ConA) forms homogeneous (homopolymeric) cross-linked precipitates with individual asparagine-linked oligomannose and bisected hybrid-type glycopeptides in the presence of binary mixtures of the carbohydrates [Bhattacharyya, L., Khan, M. I. & Brewer, C. F. (1988) Biochemistry 27, 8762-8767]. The results indicate that the ConA-glycopeptide precipitates are highly organized cross-linked lattices that are unique for each carbohydrate. Using similar techniques, the present study shows that the Gal-specific lectins from Erythrina indica and Ricinus communis (agglutinin I) form homogeneous cross-linked complexes with individual carbohydrates in binary mixtures of triantennary and tetraantennary complex-type oligosaccharides with terminal Gal residues. Conversely, binary mixtures of Gal/GalNAc-specific lectins from E. indica, Erythrina cristagalli, Erythrina flabelliformis, R. communis, soybean (Glycine max), and Wistaria floribunda (tetramer) in the presence of a naturally occurring or synthetic branched-chain oligosaccharide with terminal GalNAc or Gal residues provide evidence for the formation of separate cross-linked lattices between each lectin and the carbohydrate. The present results therefore demonstrate the formation of homogeneous lectin-carbohydrate cross-linked lattices in (a) a mixture of branched-chain complex-type oligosaccharides in the presence of a specific Gal/GalNAc-binding lectin, and (b) a mixture of lectins with similar physicochemical and carbohydrate binding properties in the presence of an oligosaccharide. These findings show that lectin-carbohydrate cross-linking interactions provide a high degree of specificity which may be relevant to their biological functions as receptors.     

Haplotype-Phased Synthetic Long Reads from Short-Read Sequencing

Next-generation DNA sequencing has revolutionized the study of biology. However, the short read lengths of the dominant instruments complicate assembly of complex genomes and haplotype phasing of mixtures of similar sequences. Here we demonstrate a method to reconstruct the sequences of individual nucleic acid molecules up to 11.6 kilobases in length from short (150-bp) reads. We show that our method can construct 99.97%-accurate synthetic reads from bacterial, plant, and animal genomic samples, full-length mRNA sequences from human cancer cell lines, and individual HIV env gene variants from a mixture. The preparation of multiple samples can be multiplexed into a single tube, further reducing effort and cost relative to competing approaches. Our approach generates sequencing libraries in three days from less than one microgram of DNA in a single-tube format without custom equipment or specialized expertise.     

Identification of point mutations in mixtures by capillary electrophoresis hybridization.

We have developed a rapid method for unambiguous identification and mutant fraction determination of individual mutants in mixtures of DNA sequence variants each differing by one or a few nucleotides. This method has applications to such diverse areas as interpretation of mutational spectra, screening of populations for polymorphisms and identification of species in environmental mixtures. In our approach, a mixture of unknown sequences labeled with a fluorescent dye is combined with a set of predetermined sequences (standards) representing the variants to be assayed. Labeling the standards with another dye allows the two sets of variants to be measured independently. Using constant denaturing capillary electrophoresis, the sequence variants are separated as individual peaks on the basis of differential melting equilibria. The unknown sequence variants are initially identified based on co-migration with particular standards. This preliminary identification is verified by hybridization of the unknown variants with the co-migrating standards within the capillary. We demonstrate the use of capillary electrophoresis hybridization to dissect complex mutational spectra of human cells in culture.     

Analysis of any point mutation in DNA. The amplification refractory mutation system (ARMS).

We have improved the "polymerase chain reaction" (PCR) to permit rapid analysis of any known mutation in genomic DNA. We demonstrate a system, ARMS (Amplification Refractory Mutation System), that allows genotyping solely by inspection of reaction mixtures after agarose gel electrophoresis. The system is simple, reliable and non-isotopic. It will clearly distinguish heterozygotes at a locus from homozygotes for either allele. The system requires neither restriction enzyme digestion, allele-specific oligonucleotides as conventionally applied, nor the sequence analysis of PCR products. The basis of the invention is that unexpectedly, oligonucleotides with a mismatched 3'-residue will not function as primers in the PCR under appropriate conditions. We have analysed DNA from patients with alpha 1-antitrypsin (AAT) deficiency, from carriers of the disease and from normal individuals. Our findings are in complete agreement with allele assignments derived by direct sequencing of PCR products.     

DNA-containing extracellular 50-nm particles in the ileal Peyer's patch of sheep

Follicles of the ileal Peyer's patch are sites of B cell proliferation and of diversification of the primary immunoglobulin repertoire in ruminants. We demonstrate here that 50-nm carbonic anhydrase-reactive particles released in the intercellular space in the follicle-associated epithelium of the ileal Peyer's patch of lambs contain DNA protected with a detergent-resistant membrane. We named these particles DiCAPs (DNA in carbonic anhydrase particles). DiCAPs can be purified from a suspension collected from ileal Peyer's patch follicles by sedimentation in a sucrose gradient. The DiCAP membrane is resistant to several ionic and non-ionic detergents alone, but can be disrupted by a combination of Triton X-100 and proteinase K. Differential nuclease treatment of purified DiCAPs indicates that they contain DNA. Digestion of DiCAP DNA with six-base pair restriction enzymes produces smears, suggesting that individual DiCAPs contain unique sequences. Nonetheless, the size of DiCAP DNA is smaller (approximately 16 kb) than that of lamb genomic DNA. Polymerase chain reaction and sequence analysis of DiCAP DNA reveals the presence of light and heavy chain variable genes as well as housekeeping genes. The data demonstrate the presence of DNA in these extracellular particles, and suggest a role of DiCAPs in transfer of DNA between cells within the ileal Peyer's patch. This raises the possibility of a novel form of communication between cells mediated by nucleic acids.     

Improving specificity of DNA hybridization-based methods

Methods based on DNA reassociation in solution with the subsequent PCR amplification of certain hybrid molecules, such as coincidence cloning and subtractive hybridization, all suffer from a common imperfection: cross-hybridization between various types of paralogous repetitive DNA fragments. Although the situation can be slightly improved by the addition of repeat-specific competitor DNA into the hybridization mixture, the cross-hybridization outcome is a significant number of background chimeric clones in resulting DNA libraries. In order to overcome this challenge, we developed a technique called mispaired DNA rejection (MDR), which utilizes a treatment of resulting reassociated DNA with mismatch-specific nucleases. We examined the MDR efficiency using cross-hybridization of complex, whole genomic mixtures derived from human and chimpanzee genomes, digested with frequent-cutter restriction enzyme. We show here that both single-stranded DNA-specific and mismatched double-stranded DNA-specific nucleases can be used for MDR separately or in combination, reducing the background level from 60 to 4% or lower. The technique presented here is of universal usefulness and can be applied to both cDNA and genomic DNA subtractions of very complex DNA mixtures. MDR is also useful for the genome-wide recovery of highly conserved DNA sequences, as we demonstrate by comparing human and pygmy marmoset genomes.     

Origin and molecular organization of supernumerary chromosomes of Prochilodus lineatus (Characiformes, Prochilodontidae) obtained by DNA probes

In Prochilodus lineatus B-chromosomes are visualized as reduced size extra elements identified as microchromosomes and are variable in morphology and number. We describe the specific total probe (B-chromosome probe) in P. lineatus obtained by chromosome microdissection and a whole genomic probe (genomic probe) from an individual without B-chromosome. The specific B-chromosome was scraped and processed to obtain DNA with amplification by DOP-PCR, and so did the genomic probe DNA. Fluorescence in situ hybridization using the B-chromosome probe labeled with dUTP-Tetramethyl-rhodamine and the genomic probe labeled with digoxigenin-FITC permitted to establish that in this species supernumerary chromosomes with varying number and morphology had different structure of chromatin when compared to that of the regular chromosomes or A complement, since only these extra elements were labeled in the metaphases. The present findings suggest that modifications in the chromatin structure of B-chromosomes to differentiate them from the A chromosomes could occur along their dispersion in the individuals of the population.     

A statistical approach for detecting genomic aberrations in heterogeneous tumor samples from single nucleotide polymorphism genotyping data

We describe a statistical method for the characterization of genomic aberrations in single nucleotide polymorphism microarray data acquired from cancer genomes. Our approach allows us to model the joint effect of polyploidy, normal DNA contamination and intra-tumour heterogeneity within a single unified Bayesian framework. We demonstrate the efficacy of our method on numerous datasets including laboratory generated mixtures of normal-cancer cell lines and real primary tumours.     

Zinc finger protein ZFP57 requires its co-factor to recruit DNA methyltransferases and maintains DNA methylation imprint in embryonic stem cells via its transcriptional repression domain

Previously, we discovered that ZFP57 is a maternal-zygotic effect gene, and it maintains DNA methylation genomic imprint at multiple imprinted regions in mouse embryos. Despite these findings, it remains elusive how DNA methyltransferases are targeted to the imprinting control regions to initiate and maintain DNA methylation imprint. To gain insights into these essential processes in genomic imprinting, we examined how ZFP57 maintains genomic DNA methylation imprint in mouse embryonic stem (ES) cells. Here we demonstrate that the loss of ZFP57 in mouse ES cells led to a complete loss of genomic DNA methylation imprint at multiple imprinted regions, similar to its role in mouse embryos. However, reintroduction of ZFP57 into Zfp57-null ES cells did not result in reacquisition of DNA methylation imprint, suggesting that the memory for genomic imprinting had been lost or altered in Zfp57-null ES cells in culture. Interestingly, ZFP57 and DNA methyltransferases could form complexes in the presence of KAP1/TRIM28/TIF1β when co-expressed in COS cells. We also found that the wild-type exogenous ZFP57 but not the mutant ZFP57 lacking the KRAB box that interacts with its co-factor KAP1/TRIM28/TIF1β could substitute for the endogenous ZFP57 in maintaining the DNA methylation imprint in ES cells. These results suggest that ZFP57 may recruit DNA methyltransferases to its target regions to maintain DNA methylation imprint, and this interaction is likely facilitated by KAP1/TRIM28/TIF1β.     

Bayesian hierarchically weighted finite mixture models for samples of distributions

Finite mixtures of Gaussian distributions are known to provide an accurate approximation to any unknown density. Motivated by DNA repair studies in which data are collected for samples of cells from different individuals, we propose a class of hierarchically weighted finite mixture models. The modeling framework incorporates a collection of k Gaussian basis distributions, with the individual-specific response densities expressed as mixtures of these bases. To allow heterogeneity among individuals and predictor effects, we model the mixture weights, while treating the basis distributions as unknown but common to all distributions. This results in a flexible hierarchical model for samples of distributions. We consider analysis of variance-type structures and a parsimonious latent factor representation, which leads to simplified inferences on non-Gaussian covariance structures. Methods for posterior computation are developed, and the model is used to select genetic predictors of baseline DNA damage, susceptibility to induced damage, and rate of repair.     

A new method of pebrine inspection of silkworm egg using multiprimer PCR

Using a mixture of several PCR primers, we evaluated whether multiprimer PCR is practically useful for the early and simultaneous detection of several kinds of microsporidia that cause silkworm pebrine. When genomic DNA extracted from silkworm eggs infected with Nosema bombycis was used as the DNA template, the specific DNA sequences were amplified by multiprimer PCR. In addition, similar results were obtained even when genomic DNA extracted from silkworms infected with N. bombycis was used as the DNA template. These findings suggest that multiprimer PCR using several primers designed for this study is practically useful for pebrine inspection of silkworm eggs.     

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.     

Genomic diversity,population structure,and migration following rapid range expansion in the Balsam Poplar,Populus balsamifera

Rapid range expansions can cause pervasive changes in the genetic diversity and structure of populations. The postglacial history of the Balsam Poplar, Populus balsamifera, involved the colonization of most of northern North America, an area largely covered by continental ice sheets during the last glacial maximum. To characterize how this expansion shaped genomic diversity within and among populations, we developed 412 SNP markers that we assayed for a range‐wide sample of 474 individuals sampled from 34 populations. We complemented the SNP data set with DNA sequence data from 11 nuclear loci from 94 individuals, and used coalescent analyses to estimate historical population size, demographic growth, and patterns of migration. Bayesian clustering identified three geographically separated demes found in the Northern, Central, and Eastern portions of the species’ range. These demes varied significantly in nucleotide diversity, the abundance of private polymorphisms, and population substructure. Most measures supported the Central deme as descended from the primary refuge of diversity. Both SNPs and sequence data suggested recent population growth, and coalescent analyses of historical migration suggested a massive expansion from the Centre to the North and East. Collectively, these data demonstrate the strong influence that range expansions exert on genomic diversity, both within local populations and across the range. Our results suggest that an in‐depth knowledge of nucleotide diversity following expansion requires sampling within multiple populations, and highlight the utility of combining insights from different data types in population genomic studies.     

Global mapping of transposon location

Transposable genetic elements are ubiquitous, yet their presence or absence at any given position within a genome can vary between individual cells, tissues, or strains. Transposable elements have profound impacts on host genomes by altering gene expression, assisting in genomic rearrangements, causing insertional mutations, and serving as sources of phenotypic variation. Characterizing a genome's full complement of transposons requires whole genome sequencing, precluding simple studies of the impact of transposition on interindividual variation. Here, we describe a global mapping approach for identifying transposon locations in any genome, using a combination of transposon-specific DNA extraction and microarray-based comparative hybridization analysis. We use this approach to map the repertoire of endogenous transposons in different laboratory strains of Saccharomyces cerevisiae and demonstrate that transposons are a source of extensive genomic variation. We also apply this method to mapping bacterial transposon insertion sites in a yeast genomic library. This unique whole genome view of transposon location will facilitate our exploration of transposon dynamics, as well as defining bases for individual differences and adaptive potential.