A Pseudolikelihood Approach for Simultaneous Analysis of Array Comparative Genomic Hybridizations (aCGH)

DNA sequence copy number has been shown to be associated with cancer development and progression. Array-based Comparative Genomic Hybridization (aCGH) is a recent development that seeks to identify the copy number ratio at large numbers of markers across the genome. Due to experimental and biological variations across chromosomes and across hybridizations, current methods are limited to analyses of single chromosomes. We propose a more powerful approach that borrows strength across chromosomes and across hybridizations. We assume a Gaussian mixture model, with a hidden Markov dependence structure, and with random effects to allow for intertumoral variation, as well as intratumoral clonal variation. For ease of computation, we base estimation on a pseudolikelihood function. The method produces quantitative assessments of the likelihood of genetic alterations at each clone, along with a graphical display for simple visual interpretation. We assess the characteristics of the method through simulation studies and through analysis of a brain tumor aCGH data set. We show that the pseudolikelihood approach is superior to existing methods both in detecting small regions of copy number alteration and in accurately classifying regions of change when intratumoral clonal variation is present.

EXPLORATION, NORMALIZATION, AND GENOTYPE CALLS OF HIGH DENSITY OLIGONUCLEOTIDE SNP ARRAY DATA

In most microarray technologies, a number of critical steps are required to convert raw intensity measurements into the data relied upon by data analysts, biologists and clinicians. These data manipulations, referred to as preprocessing, can influence the quality of the ultimate measurements. In the last few years, the high-throughput measurement of gene expression is the most popular application of microarray technology. For this application, various groups have demonstrated that the use of modern statistical methodology can substantially improve accuracy and precision of gene expression measurements, relative to ad-hoc procedures introduced by designers and manufacturers of the technology. Currently, other applications of microarrays are becoming more and more popular. In this paper we describe a preprocessing methodology for a technology designed for the identification of DNA sequence variants in specific genes or regions of the human genome that are associated with phenotypes of interest such as disease. In particular we describe methodology useful for preprocessing Affymetrix SNP chips and obtaining genotype calls with the preprocessed data. We demonstrate how our procedure improves existing approaches using data from three relatively large studies including one in which large number independent calls are available. Software implementing these ideas are avialble from the Bioconductor oligo package.

A modular LHC built on the DNA three-way junction

A light-harvesting complex composed of a p-stacked multichromophoric array in a DNA three-way junction is described. The modular design allows for a ready exchange of non-covalently attached energy acceptors

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

From Ribbons to Networks: Hierarchical Organization of DNA-Grafted Supramolecular Polymers

DNA-grafted supramolecular polymers (SPs) allow the programmed organization of DNA in a highly regular, one-dimensional array. Oligonucleotides are arranged along the edges of pyrene-based helical polymers. Addition of complementary oligonucleotides triggers the assembly of individual nanoribbons resulting in the development of extended supramolecular networks. Network formation is enabled by cooperative coaxial stacking interactions of terminal GC base pairs. The process is accompanied by structural changes in the pyrene polymer core that can be followed spectroscopically. Network formation is reversible, and disassembly into individual ribbons is realized either via thermal denaturation or by addition of a DNA separator strand.

Supramolecular Polymers in DNA Nanotechnology

Creation of biocompatible functional materials is an important task in supramolecular chemistry. In this contribution, we report on noncovalent synthesis of DNA-grafted supramolecular polymers (SPs). DNA-grafted SPs enable programmed arrangement of oligonucleotides in a regular, tightly packed one-dimensional array. Further interactions of DNA-grafted SPs with complementary DNA strands leads to the formation of networks through highly cooperative G-C blunt-end stacking interactions. The structural changes in the polymeric core enable to monitor spectroscopically the stepwise formation of networks. Such stimuli-responsive supramolecular networks may lead to the development of DNA-based smart materials.

Clinical relevance of polymorphic DNA copy number variation (CNV) in African American women with stage I-II breast cancer

Objectives. The chief goal of this study was to analyze copy number variation (CNV) in breast cancer tumors from 25 African American women with early stage breast cancer (BC) using molecular inversion probes (MIP) in order to: (1) compare the degree of CNV in tumors compared to normal lymph nodes, and (2) determine whether gains and/or losses of genes in specific chromosomes differ between pathologic subtypes of breast cancer defined by known prognostic markers, (3) determine whether gains/losses in CN are associated with known oncogenes or tumor suppressor genes, and (4) determine whether increased gains/losses in CN for specific chromosomes were associated with differences in breast cancer recurrence. ^ Methods. Twenty to 37 nanograms of DNA extracted from 25 formalin-fixed paraffin embedded (FFPE) tumor samples and matched normal lymph nodes were added to individual tubes. Oligonucleotide probes with recognition sequences at each terminus were hybridized with a genomic target sequence to form a circular structure. Probes are released from genomic DNA obtained from FFPE samples, and those which have been correctly "circularized" in the proper allele/nucleotide reaction combination are amplified using polymerase chain reaction (PCR) primers. Amplicons were fluorescently labeled and the tag sequences released from the genome homology regions by treatment with uracil-N-glycosylase to cleave the probe at the site where uracils are present, and detected using a complementary tag array developed by Affymetrix. ^ Results. Analysis of CN gains and losses from tumors and normal tissues showed marked differences in tumors with numerous chromosomes affected. Similar changes were not observed in normal lymph nodes. When tumors were stratified into four groups based on expression or lack of expression of the estrogen receptor and HER2/neu, distinct patterns of CNV for different chromosomes were observed. Gains or losses in CN for specific chromosomes correlated with amplifications/deletions of particular oncogenes or tumor suppressor genes (i.e. such as found on chromosome 17) known to be associated with aggressive tumor phenotype and poor prognosis. There was a trend for increases in CN observed for chromosome 17 to correlate inversely with time to recurrence of BC (p=0.14 for trend). CNV was also observed for chromosomes 5, 8, 10, 11, and 16, which are known sites for several breast cancer susceptibility alleles. ^ Conclusions. This study is the first to validate the MIP technique, to correlate differences in gene expression with known prognostic tumor markers, and to correlate significant increases/decreases in CN with known tumor markers associated with prognosis. The results of this study may have far reaching public health implications towards identifying new high-risk groups based on genomic differences in CNP, both with respect to prognosis and response to therapy, and to eventually identify new therapeutic targets for prevention and treatment of this disease. ^

Automated parallel DNA sequencing on multiple channel microchips

We report automated DNA sequencing in 16-channel microchips. A microchip prefilled with sieving matrix is aligned on a heating plate affixed to a movable platform. Samples are loaded into sample reservoirs by using an eight-tip pipetting device, and the chip is docked with an array of electrodes in the focal plane of a four-color scanning detection system. Under computer control, high voltage is applied to the appropriate reservoirs in a programmed sequence that injects and separates the DNA samples. An integrated four-color confocal fluorescent detector automatically scans all 16 channels. The system routinely yields more than 450 bases in 15 min in all 16 channels. In the best case using an automated base-calling program, 543 bases have been called at an accuracy of >99%. Separations, including automated chip loading and sample injection, normally are completed in less than 18 min. The advantages of DNA sequencing on capillary electrophoresis chips include uniform signal intensity and tolerance of high DNA template concentration. To understand the fundamentals of these unique features we developed a theoretical treatment of cross-channel chip injection that we call the differential concentration effect. We present experimental evidence consistent with the predictions of the theory.

p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32–36°, whereas wt p53 bends it by 51–57°. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by ≈35° and ≈70°, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53–DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

A novel application of gene arrays: Escherichia coli array provides insight into the biology of the obligate endosymbiont of tsetse flies

Symbiotic associations with microorganisms are pivotal in many insects. Yet, the functional roles of obligate symbionts have been difficult to study because it has not been possible to cultivate these organisms in vitro. The medically important tsetse fly (Diptera: Glossinidae) relies on its obligate endosymbiont, Wigglesworthia glossinidia, a member of the Enterobacteriaceae, closely related to Escherichia coli, for fertility and possibly nutrition. We show here that the intracellular Wigglesworthia has a reduced genome size smaller than 770 kb. In an attempt to understand the composition of its genome, we used the gene arrays developed for E. coli. We were able to identify 650 orthologous genes in Wigglesworthia corresponding to ≈85% of its genome. The arrays were also applied for expression analysis using Wigglesworthia cDNA and 61 gene products were detected, presumably coding for some of its most abundant products. Overall, genes involved in cell processes, DNA replication, transcription, and translation were found largely retained in the small genome of Wigglesworthia. In addition, genes coding for transport proteins, chaperones, biosynthesis of cofactors, and some amino acids were found to comprise a significant portion, suggesting an important role for these proteins in its symbiotic life. Based on its expression profile, we predict that Wigglesworthia may be a facultative anaerobic organism that utilizes ammonia as its major source of nitrogen. We present an application of E. coli gene arrays to obtain broad genome information for a closely related organism in the absence of complete genome sequence data.

DNA analysis and diagnostics on oligonucleotide microchips.

We present a further development in the technology of sequencing by hybridization to oligonucleotide microchips (SHOM) and its application to diagnostics for genetic diseases. A robot has been constructed to manufacture sequencing "microchips." The microchip is an array of oligonucleotides immobilized into gel elements fixed on a glass plate. Hybridization of the microchip with fluorescently labeled DNA was monitored in real time simultaneously for all microchip elements with a two-wavelength fluorescent microscope equipped with a charge-coupled device camera. SHOM has been used to detect beta-thalassemia mutations in patients by hybridizing PCR-amplified DNA with the microchips. A contiguous stacking hybridization technique has been applied for the detection of mutations; it can simplify medical diagnostics and enhance its reliability. The use of multicolor monitoring of contiguous stacking hybridization is suggested for large-scale diagnostics and gene polymorphism studies. Other applications of the SHOM technology are discussed.

DNA binding specificity of two homeodomain proteins in vitro and in Drosophila embryos.

In previous experiments, the homeodomain proteins even-skipped and fushi-tarazu were found to UV cross-link to a surprisingly wide array of DNA sites in living Drosophila embryos. We now show that UV cross-linking gives a highly accurate measure of DNA binding by these proteins. In addition, the binding of even-skipped and fushi-tarazu proteins has been measured in vitro to the same DNA fragments that were examined in vivo. This analysis shows that these proteins have broad DNA recognition properties in vitro that are likely to be important determinants of their distribution on DNA in vivo, but it also shows that in vitro DNA binding specificity alone is not sufficient to explain the distribution of these proteins in embryos.

Equine severe combined immunodeficiency: a defect in V(D)J recombination and DNA-dependent protein kinase activity.

V(D)J rearrangement is the molecular mechanism by which an almost infinite array of specific immune receptors are generated. Defects in this process result in profound immunodeficiency as is the case in the C.B-17 SCID mouse or in RAG-1 (recombination-activating gene 1) or RAG-2 deficient mice. It has recently become clear that the V(D)J recombinase most likely consists of both lymphoid-specific factors and ubiquitously expressed components of the DNA double-strand break repair pathway. The deficit in SCID mice is in a factor that is required for both of these pathways. In this report, we show that the factor defective in the autosomal recessive severe combined immunodeficiency of Arabian foals is required for (i) V(D)J recombination, (ii) resistance to ionizing radiation, and (iii) DNA-dependent protein kinase activity.

Trapping of branched DNA in microfabricated structures.

We have observed electrostatic trapping of tribranched DNA molecules undergoing electrophoresis in a microfabricated pseudo-two-dimensional array of posts. Trapping occurs in a unique transport regimen in which the electrophoretic mobility is extremely sensitive to polymer topology. The arrest of branched polymers is explained by considering their center-of-mass motion; in certain conformations, owing to the constraints imposed by the obstacles a molecule cannot advance without the center of mass first moving a short distance backwards. The depth of the resulting local potential well can be much greater than the thermal energy so that escape of an immobilized molecule can be extremely slow. We summarize the expected behavior of the mobility as a function of field strength and topology and point out that the microfabricated arrays are highly suitable for detecting an extremely small number of branched molecules in a very large population of linear molecules.