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. ^

Mapping and cloning of genes from portions of the human genome using somatic cell hybrids

Human x rodent somatic cell hybrids have played an important role in human genetics research. They have been especially useful for assigning genes to chromosomes and isolating DNA markers from specific regions of the human genome.^ By employing a combination of somatic cell genetic, recombinant DNA, and cytogenetic techniques, human DNA excision repair gene ERCC4 was mapped regionally to human 16p13.13-13.2, even though the gene has not been cloned. Human x Chinese hamster ovary (CHO) cell hybrids selected for human ERCC4 activity and containing 16p13.1-p13.3 as the only human genetic material were identified. These hybrids were used to order DNA markers located in 16p13.1-p13.3. New DNA markers physically close to ERCC4 were isolated from such hybrids. Using amplified human DNA from the hybrids as probe in fluorescent in situ hybridization, the short arm breakpoint in the chromosome 16 inversion associated with acute myelomonocytic leukemia (AMML) was found to be physically close to the ERCC4 gene. The physical mapping and eventually, the cloning of the ERCC4 gene, will benefit the understanding of the DNA repair system and the study of other important biomedical problems such as tumorigenesis.^ To facilitate the cloning of ERCC4 gene and, in general, the cloning of genes from any defined regions of the human genome, a method was developed for the direct isolation of human transcribed genes ffom somatic cell hybrids. cDNA was prepared from human x rodent hybrid by using consensus 5$\sp\prime$ splice site sequences as primers. These primers were designed to select immature, unspliced messenger RNA (still retaining species specific repeat sequences) as templates. Screening of a derived cDNA library for human repeat sequences resulted in the isolation of human clones at the anticipated frequency with characteristics expected of exons of transcribed human genes. The usefulness of the splice site specific primers was analyzed and the cDNA synthesis conditions with these primers were optimized. The procedure was shown to be sensitive enough to clone weakly expressed genes. Studying the expression of the represented genes with the isolated clones was shown to be feasible. Such regional specific human gene fragments will be very valuable for many human genetic studies such as the search of inherited disease genes and the construction of a cDNA map of the human genome. ^

Identification of the Causative Bacteria in Musculoskeletal Infections Using 16s rDNA - Denaturing Gradient Gel Electrophoresis Analysis

Musculoskeletal infections are infections of the bone and surrounding tissues. They are currently diagnosed based on culture analysis, which is the gold standard for pathogen identification. However, these clinical laboratory methods are frequently inadequate for the identification of the causative agents, because a large percentage (25-50%) of confirmed musculoskeletal infections are false negatives in which no pathogen is identified in culture. My data supports these results. The goal of this project was to use PCR amplification of a portion of the 16S rRNA gene to test an alternative approach for the identification of these pathogens and to assess the diversity of the bacteria involved. The advantages of this alternative method are that it should increase sample sensitivity and the speed of detection. In addition, bacteria that are non-culturable or in low abundance can be detected using this molecular technique. However, a complication of this approach is that the majority of musculoskeletal infections are polymicrobial, which prohibits direct identification from the infected tissue by DNA sequencing of the initial 16S rDNA amplification products. One way to solve this problem is to use denaturing gradient gel electrophoresis (DGGE) to separate the PCR products before DNA sequencing. Denaturing gradient gel electrophoresis (DGGE) separates DNA molecules based on their melting point, which is determined by their DNA sequence. This analytical technique allows a mixture of PCR products of the same length that electrophoreses through agarose gels as one band, to be separated into different bands and then used for DNA sequence analysis. In this way, the DGGE allows for the identification of individual bacterial species in polymicrobial-infected tissue, which is critical for improving clinical outcomes. By combining the 16S rDNA amplification and the DGGE techniques together, an alternative approach for identification has been used. The 16S rRNA gene PCR-DGGE method includes several critical steps: DNA extraction from tissue biopsies, amplification of the bacterial DNA, PCR product separation by DGGE, amplification of the gel-extracted DNA, and DNA sequencing and analysis. Each step of the method was optimized to increase its sensitivity and for rapid detection of the bacteria present in human tissue samples. The limit of detection for the DNA extraction from tissue was at least 20 Staphylococcus aureus cells and the limit of detection for PCR was at least 0.05 pg of template DNA. The conditions for DGGE electrophoreses were optimized by using a double gradient of acrylamide (6 – 10%) and denaturant (30-70%), which increased the separation between distinct PCR products. The use of GelRed (Biotium) improved the DNA visualization in the DGGE gel. To recover the DNA from the DGGE gels the gel slices were excised, shredded in a bead beater, and the DNA was allowed to diffuse into sterile water overnight. The use of primers containing specific linkers allowed the entire amplified PCR product to be sequenced and then analyzed. The optimized 16S rRNA gene PCR-DGGE method was used to analyze 50 tissue biopsy samples chosen randomly from our collection. The results were compared to those of the Memorial Hermann Hospital Clinical Microbiology Laboratory for the same samples. The molecular method was congruent for 10 of the 17 (59%) culture negative tissue samples. In 7 of the 17 (41%) culture negative the molecular method identified a bacterium. The molecular method was congruent with the culture identification for 7 of the 33 (21%) positive cultured tissue samples. However, in 8 of the 33 (24%) the molecular method identified more organisms. In 13 of the 15 (87%) polymicrobial cultured tissue samples the molecular method identified at least one organism that was also identified by culture techniques. Overall, the DGGE analysis of 16S rDNA is an effective method to identify bacteria not identified by culture analysis.

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

A molecular study of the rhodopsin locus in patients with retinitis pigmentosa

DNA for this study was collected from a sample of 133 retinitis pigmentosa (RP) patients and the rhodopsin locus molecularly analyzed by linkage and for disease specific mutations. The cohort of patients consisted of 85 individuals diagnosed with autosomal dominant RP (adRP), and 48 patients representing other forms of retinitis pigmentosa or retinal dystrophy related disease. In three large families with adRP rhodopsin was excluded from linkage to the disease locus. A search for subtle mutations in the rhodopsin coding region using single strand conformational polymorphisms (SSCP) and sequencing detected a total of 14 unique sequence variants in 24 unrelated patients. These variants included one splicing variant, 5168 -1G-A, one deletion variant of 17 base pairs causing a frame shift at codon 332, and 12 misense variants: Pro23His, Leu46Arg, Gly106Trp, Arg135Pro, Pro171Glu, Pro180Ala, Glu181Lys, Asp190Asn, His211Arg, Ser270Arg, Leu328Pro and Pro347Thr. All but three of the missense variants change amino acids that are evolutionarily conserved. The Pro23His mutation was found in 10 unrelated individuals with family histories of adRP and not in any normal controls (over 80 chromosomes tested). The Pro180Ala mutation was present in a patient with simplex RP and probably represents a new mutation. Three normal polymorphic nucleotide substitutions, A-269-G, T-3982-C, and G-5145-A, were also identified. We conclude, based on this study, that 25% of adRP cases are attributable to rhodopsin mutations.^ Clinical data, including ERG results and visual field testing, was available for patients with eleven different mutations. The eleven patients were all diagnosed with RP, however the severity of the disease varied with five patients mildly affected and diagnosed with type II adRP and 5 patients severely affected and diagnosed with type I adRP. The patient with simplex RP was mildly affected. The location of the mutations within the rhodopsin protein was randomly associated with the severity of the disease in those patients evaluated. However, four mutations, Pro23His, Leu46Arg, Pro347Thr, and 5168 -1G-A, are particularly interesting. The Pro23His mutation appears to have radiated from a recent common ancestor of the affected patients as all of them share a common haplotype at the rhodopsin locus. The Leu46Arg mutation causes an unusually severe form of RP. Hydropathy analysis of the mutated sequence revealed a marked change in the hydrophobicity of this first transmembrane spanning region. Codon 347 has been the target of multiple mutations with at least six documented changes at the position, significantly more than expected by a random distribution of mutations. Finally the splice-site variant is extremely variable in its expression in the family studied. Similar mutations have been reported in other cases of adRP and postulated to be involved in autosomal recessive RP (arRP). Mechanisms to account for the variable expression of rhodopsin mutations in relation to RP heterogeneity are discussed. (Abstract shortened by UMI.) ^

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

Non B-DNA structures and genomic instability associated with myotonic dystrophy type 2 (CCTG).(CAGG) repeats

There are many diseases associated with the expansion of DNA repeats in humans. Myotonic dystrophy type 2 is one of such diseases, characterized by expansions of a (CCTG)•(CAGG) repeat tract in intron 1 of zinc finger protein 9 (ZNF9) in chromosome 3q21.3. The DM2 repeat tract contains a flanking region 5' to the tract that consists of a polymorphic repetitive sequence (TG)14-25(TCTG)4-11(CCTG) n. The (CCTG)•(CAGG) repeat is typically 11-26 repeats in persons without the disease, but can expand up to 11,000 repeats in affected individuals, which is the largest expansion seen in DNA repeat diseases to date. This DNA tract remains one of the least characterized disease-associated DNA repeats, and mechanisms causing the repeat expansion in humans have yet to be elucidated. Alternative, non B-DNA structures formed by the expanded repeats are typical in DNA repeat expansion diseases. These sequences may promote instability of the repeat tracts. I determined that slipped strand structure formation occurs for (CCTG)•(CAGG) repeats at a length of 42 or more. In addition, Z-DNA structure forms in the flanking human sequence adjacent to the (CCTG)•(CAGG) repeat tract. I have also performed genetic assays in E. coli cells and results indicate that the (CCTG)•(CAGG) repeats are more similar to the highly unstable (CTG)•(CAG) repeat tracts seen in Huntington's disease and myotonic dystrophy type 1, than to those of the more stable (ATTCT)•(AGAAT) repeat tracts of spinocerebellar ataxia type 10. This instability, however, is RecA-independent in the (CCTG)•(CAGG) and (ATTCT)•(AGAAT) repeats, whereas the instability is RecA-dependent in the (CTG)•(CAG) repeats. Structural studies of the (CCTG)•(CAGG) repeat tract and the flanking sequence, as well as genetic selection assays may reveal the mechanisms responsible for the repeat instability in E. coli, and this may lead to a better understanding of the mechanisms contributing to the human disease state. ^