Polydnavirus DNA of the braconid wasp Chelonus inanitus is integrated in the wasp's genome and excised only in later pupal and adult stages of the female

Many endoparasitic wasps inject, along with the egg, polydnavirus into their insect hosts, the virus being a prerequisite for successful parasitoid development. The genome of polydnaviruses consists of multiple circular dsDNA molecules of variable size. We show for a 12 kbp segment of the braconid Chelonus inanitus (CiV12) that it is integrated into the wasp genome. This is the first direct demonstration of integration for a bracovirus. PCR data indicated that the integrated form of CiV12 was present in all male and female stages investigated while the excised circular virus DNA only appeared in females after a specific stage in pupal-adult development. The data also indicated that after excision of virus DNA the genomic DNA was rejoined. This has not yet been reported for any polydnavirus. Sequence analyses in the junction regions revealed the presence of an imperfect consensus sequence of 15 nucleotides in CiV12, in each terminus of the integrated virus DNA and in the rejoined genomic DNA. Within these repeats two sequence types (ATA, TAC) were observed in the various virus clones and in the clones encompassing the rejoined genomic DNA; they corresponded to the sequence type in the right and left junction, respectively. To explain this, we propose a model of virus DNA replication in which the genomic DNA is folded to juxtapose the direct repeat of the left with that of the right junction; recombination at specific sites would then yield the two types of virus and rejoined genomic DNA.

A DNA probe for the detection of Dicrocoelium dendriticum in ants of Formica spp. and Lasius spp

Repetitive DNA sequences present in the genome of Dicrocoelium dendriticum were identified by hybridization of genomic DNA that had been digested with different restriction enzymes with 32P-labeled genomic D. dendriticum DNA. DNA fragments containing repetitive sequences were isolated from PstI-digested D. dendriticum DNA and were subcloned into a plasmid vector. Plasmids containing repetitive sequences were identified by colony hybridization. One of these plasmids, designated Ddr-IV, was isolated and used as a probe in further studies. Ddr-IV is specific for D. dendriticum since it does not hybridize to DNA isolated from other trematodes. In addition, Ddr-IV was capable of detecting D. dendriticum metacercariae in ants (Formica cunicularia, F. rufibarbis, and Lasius sp.), which act as second intermediate hosts in the parasite's life cycle. Since metacercariae constitute the infectious stage of the parasite for grazing animals, Ddr-IV will provide a useful tool for epidemiology studies of dicrocoeliosis.

Genetic predictors of bladder cancer with an emphasis on methylation-related genes

Although tobacco exposure remains the prevailing risk factor for bladder cancer (BC), only a small percentage of exposed individuals develop cancer, suggesting that tobacco-related carcinogenesis is modulated by genetic susceptibility and possibly by DNA methylation-related events. Methylation patterns established by DNA methyltransferases (DNMTs) are influenced by dietary folate and genetic polymorphisms in the methylene-tetrahydrofolate reductase gene (MTHFR). Therefore, we hypothesized that DNA methylation-related genes, such as DNMT3B and MTHFR, might modulate BC risk. ^ In a study of 514 Caucasian BC cases and 498 healthy Caucasian controls examining the DNMT3B C46359T polymorphism, CC genotype was found to be a risk factor in women (Odds Ratio (OR) = 1.79), but not in men. This risk was further increased among women who were never smokers, consumed low dietary folate, and had adverse variants of MTHFR. In addition, higher DNMT3B expression among smokers was a risk factor (OR = 4.27) and correlated with genetic variants of the DNMT3B C46359T polymorphism, providing salient evidence for the risk associated with the CC variant. This suggests that the DNMT3B CC variant may confer a predisposition toward aberrant de novo methylation of CpG islands in critical tumor suppressor genes. ^ The convergence of alterations in DNMT3B, associated with promoter methylation, and reduced dietary folate consumption, accompanying global hypomethylation and genetic instability, may act synergistically to promote bladder carcinogenesis, especially in women. The results of this study unveiled new gender-specific paradigms of BC risk for women and demonstrated that this risk can be modified by folate consumption as well as polymorphisms in the folate pathway. ^

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

Mechanism of methylation gene reactivation in human cancer cells

Epigenetic silencing of tumor suppressor genes by DNA hypermethylation at promoter regions is a common event in carcinogenesis and tumor progression. Abrogation of methylation and reversal of epigenetic silencing is a very potent way in cancer treatment. However, the reactivation mechanisms are poorly understood. In this study, we first developed a cell line model system named YB5, derived from SW48 cancer cell line, which bears one copy of stably integrated EGFP gene on Chromosome 1p31.1 region. The GFP gene expression is transcriptionally silenced due to the hypermethylated promoter CMV. However, the GFP expression can be restored using demethylating agent 5-aza-2' deoxycytidine (DAC), and detected by FACS and fluorescent microscopy. Using this system, we observed the heterogeneous reactivation induced by DAC treatment. After flow sorting, GFP negative cells exhibited similar level of incomplete demethylation compared to GFP positive cells on repetitive LINE1 element, tumor suppressor genes such as P16, CDH13, and RASSF1a, and CMV promoter as well. However, the local chromatin of CMV-GFP locus altered to an open structure marked by high H3 lysine 9 acetylation and low H3 lysine 27 tri-methylation in GFP positive cells, while the GFP negative cells retained mostly the original repressive marks. Thus, we concluded that DAC induced DNA hypomethylation alone does not directly determine the level of re-expression, and the resetting of the local chromatin structure under hypomethylation environment is required for gene reactivation. Besides, a lentivirus vector-based shRNA screening was performed using the YB5 system. Although it is the rare chance that vector lands in the neighboring region of GFP, we found that the exogenous vector DNA inserted into the upstream region of GFP gene locus led to the promoter demethylation and reactivated the silenced GFP gene. Thus, epigenetic state can be affected by changing of the adjacent nucleic acid sequences. Further, this hypermethylation silenced system was utilized for epigenetic drug screening. We have found that DAC combined with carboplatin would enhance the GFP% yield and increase expression of other tumor suppressor genes than DAC alone, and this synergistic effect may be related to DNA repair process. In summary, these studies reveal that reversing of methylation silencing requires coordinated alterations of DNA methylation, chromatin structure, and local microenvironment. ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

A cluster of oppositely imprinted transcripts at the Gnas locus in the distal imprinting region of mouse chromosome 2

Imprinted genes tend to occur in clusters. We have identified a cluster in distal mouse chromosome (Chr) 2, known from early genetic studies to contain both maternally and paternally imprinted, but unspecified, genes. Subsequently, one was identified as Gnas, which encodes a G protein α subunit, and there is clinical and biochemical evidence that the human homologue GNAS1, mutated in patients with Albright hereditary osteodystrophy, is also imprinted. We have used representational difference analysis, based on parent-of-origin methylation differences, to isolate candidate imprinted genes in distal Chr 2 and found two oppositely imprinted genes, Gnasxl and Nesp. Gnasxl determines a variant G protein α subunit associated with the trans-Golgi network and Nesp encodes a secreted protein of neuroendocrine tissues. Gnasxl is maternally methylated in genomic DNA and encodes a paternal-specific transcript, whereas Nesp is paternally methylated with maternal-specific expression. Their reciprocal imprinting may offer insight into the distal Chr 2 imprinting phenotypes. Remarkably, Gnasxl, Nesp, and Gnas are all part of the same transcription unit; transcripts for Gnasxl and Nesp are alternatively spliced onto exon 2 of Gnas. This demonstrates an imprinting mechanism in which two oppositely imprinted genes share the same downstream exons.

Transgenic DNA integrated into the oat genome is frequently interspersed by host DNA

Integration of transgenic DNA into the plant genome was investigated in 13 transgenic oat (Avena sativa L.) lines produced using microprojectile bombardment with one or two cotransformed plasmids. In all transformation events, the transgenic DNA integrated into the plant genome consisted of intact transgene copies that were accompanied by multiple, rearranged, and/or truncated transgene fragments. All fragments of transgenic DNA cosegregated, indicating that they were integrated at single gene loci. Analysis of the structure of the transgenic loci indicated that the transgenic DNA was interspersed by the host genomic DNA. The number of insertions of transgenic DNA within the transgene loci varied from 2 to 12 among the 13 lines. Restriction endonucleases that do not cleave the introduced plasmids produced restriction fragments ranging from 3.6 to about 60 kb in length hybridizing to a probe comprising the introduced plasmids. Although the size of the interspersing host DNA within the transgene locus is unknown, the sizes of the transgene-hybridizing restriction fragments indicated that the entire transgene locus must be at least from 35–280 kb. The observation that all transgenic lines analyzed exhibited genomic interspersion of multiple clustered transgenes suggests a predominating integration mechanism. We propose that transgene integration at multiple clustered DNA replication forks could account for the observed interspersion of transgenic DNA with host genomic DNA within transgenic loci.

Association of arsenic-induced malignant transformation with DNA hypomethylation and aberrant gene expression

Inorganic arsenic, a human carcinogen, is enzymatically methylated for detoxication, consuming S-adenosyl-methionine (SAM) in the process. The fact that DNA methyltransferases (MeTases) require this same methyl donor suggests a role for methylation in arsenic carcinogenesis. Here we test the hypothesis that arsenic-induced initiation results from DNA hypomethylation caused by continuous methyl depletion. The hypothesis was tested by first inducing transformation in a rat liver epithelial cell line by chronic exposure to low levels of arsenic, as confirmed by the development of highly aggressive, malignant tumors after inoculation of cells into Nude mice. Global DNA hypomethylation occurred concurrently with malignant transformation and in the presence of depressed levels of S-adenosyl-methionine. Arsenic-induced DNA hypomethylation was a function of dose and exposure duration, and remained constant even after withdrawal of arsenic. Hyperexpressibility of the MT gene, a gene for which expression is clearly controlled by DNA methylation, was also detected in transformed cells. Acute arsenic or arsenic at nontransforming levels did not induce global hypomethylation of DNA. Whereas transcription of DNA MeTase was elevated, the MeTase enzymatic activity was reduced with arsenic transformation. Taken together, these results indicate arsenic can act as a carcinogen by inducing DNA hypomethylation, which in turn facilitates aberrant gene expression, and they constitute a tenable theory of mechanism in arsenic carcinogenesis.

Panhandle PCR strategy to amplify MLL genomic breakpoints in treatment-related leukemias

Panhandle PCR amplifies genomic DNA with known 5′ and unknown 3′ sequences from a template with an intrastrand loop schematically shaped like a pan with a handle. We used panhandle PCR to clone MLL genomic breakpoints in two pediatric treatment-related leukemias. The karyotype in a case of treatment-related acute lymphoblastic leukemia showed the t(4;11)(q21;q23). Panhandle PCR amplified the translocation breakpoint at position 2158 in intron 6 in the 5′ MLL breakpoint cluster region (bcr). The karyotype in a case of treatment-related acute myeloid leukemia was normal, but Southern blot analysis showed a single MLL gene rearrangement. Panhandle PCR amplified the breakpoint at position 1493 in MLL intron 6. Screening of somatic cell hybrid and radiation hybrid DNAs by PCR and reverse transcriptase-PCR analysis of the leukemic cells indicated that panhandle PCR identified a fusion of MLL intron 6 with a previously uncharacterized sequence in MLL intron 1, consistent with a partial duplication. In both cases, the breakpoints in the MLL bcr were in Alu repeats, and there were Alu repeats in proximity to the breakpoints in the partner DNAs, suggesting that Alu sequences were relevant to these rearrangements. This study shows that panhandle PCR is an effective method for cloning MLL genomic breakpoints in treatment-related leukemias. Analysis of additional pediatric cases will determine whether breakpoint distribution deviates from the predilection for 3′ distribution in the bcr that has been found in adult cases.

Inhibition of DNA methyltransferase stimulates the expression of signal transducer and activator of transcription 1, 2, and 3 genes in colon tumor cells

Inhibitors of DNA methyltransferase, typified by 5-aza-2′-deoxycytidine (5-Aza-CdR), induce the expression of genes transcriptionally down-regulated by de novo methylation in tumor cells. We utilized gene expression microarrays to examine the effects of 5-Aza-CdR treatment in HT29 colon adenocarcinoma cells. This analysis revealed the induction of a set of genes that implicated IFN signaling in the HT29 cellular response to 5-Aza-CdR. Subsequent investigations revealed that the induction of this gene set correlates with the induction of signal transducer and activator of transcription (STAT) 1, 2, and 3 genes and their activation by endogenous IFN-α. These observations implicate the induction of the IFN-response pathway as a major cellular response to 5-Aza-CdR and suggests that the expression of STATs 1, 2, and 3 can be regulated by DNA methylation. Consistent with STAT’s limiting cell responsiveness to IFN, we found that 5-Aza-CdR treatment sensitized HT29 cells to growth inhibition by exogenous IFN-α2a, indicating that 5-Aza-CdR should be investigated as a potentiator of IFN responsiveness in certain IFN-resistant tumors.

Checkpoint Defects Leading to Premature Mitosis Also Cause Endoreplication of DNA in Aspergillus nidulans

The G2 DNA damage and slowing of S-phase checkpoints over mitosis function through tyrosine phosphorylation of NIMXcdc2 in Aspergillus nidulans. We demonstrate that breaking these checkpoints leads to a defective premature mitosis followed by dramatic rereplication of genomic DNA. Two additional checkpoint functions, uvsB and uvsD, also cause the rereplication phenotype after their mutation allows premature mitosis in the presence of low concentrations of hydroxyurea. uvsB is shown to encode a rad3/ATR homologue, whereas uvsD displays homology to rad26, which has only previously been identified in Schizosaccharomyces pombe. uvsBrad3 and uvsDrad26 have G2 checkpoint functions over mitosis and another function essential for surviving DNA damage. The rereplication phenotype is accompanied by lack of NIMEcyclinB, but ectopic expression of active nondegradable NIMEcyclinB does not arrest DNA rereplication. DNA rereplication can also be induced in cells that enter mitosis prematurely because of lack of tyrosine phosphorylation of NIMXcdc2 and impaired anaphase-promoting complex function. The data demonstrate that lack of checkpoint control over mitosis can secondarily cause defects in the checkpoint system that prevents DNA rereplication in the absence of mitosis. This defines a new mechanism by which endoreplication of DNA can be triggered and maintained in eukaryotic cells.

Rescue of dystrophin expression in mdx mouse muscle by RNA/DNA oligonucleotides

Chimeric RNA/DNA oligonucleotides (“chimeraplasts”) have been shown to induce single base alterations in genomic DNA both in vitro and in vivo. The mdx mouse strain has a point mutation in the dystrophin gene, the consequence of which is a muscular dystrophy resulting from deficiency of the dystrophin protein in skeletal muscle. To test the feasibility of chimeraplast-mediated gene therapy for muscular dystrophies, we used a chimeraplast (designated “MDX1”) designed to correct the point mutation in the dystrophin gene in mdx mice. After direct injection of MDX1 into muscles of mdx mice, immunohistochemical analysis revealed dystrophin-positive fibers clustered around the injection site. Two weeks after single injections into tibialis anterior muscles, the maximum number of dystrophin-positive fibers (approximately 30) in any muscle represented 1–2% of the total number of fibers in that muscle. Ten weeks after single injections, the range of the number of dystrophin-positive fibers was similar to that seen after 2 wk, suggesting that the expression was stable, as would be predicted for a gene-conversion event. Staining with exon-specific antibodies showed that none of these were “revertant fibers.” Furthermore, dystrophin from MDX1-injected muscles was full length by immunoblot analysis. No dystrophin was detectable by immunohistochemical or immunoblot analysis after control chimeraplast injections. Finally, reverse transcription–PCR analysis demonstrated the presence of transcripts with the wild-type dystrophin sequence only in mdx muscles injected with MDX1 chimeraplasts. These results provide the foundation for further studies of chimeraplast-mediated gene therapy as a therapeutic approach to muscular dystrophies and other genetic disorders of muscle.

Clustered DNA damages induced in isolated DNA and in human cells by low doses of ionizing radiation

Clustered DNA damages—two or more closely spaced damages (strand breaks, abasic sites, or oxidized bases) on opposing strands—are suspects as critical lesions producing lethal and mutagenic effects of ionizing radiation. However, as a result of the lack of methods for measuring damage clusters induced by ionizing radiation in genomic DNA, neither the frequencies of their production by physiological doses of radiation, nor their repairability, nor their biological effects are known. On the basis of methods that we developed for quantitating damages in large DNAs, we have devised and validated a way of measuring ionizing radiation-induced clustered lesions in genomic DNA, including DNA from human cells. DNA is treated with an endonuclease that induces a single-strand cleavage at an oxidized base or abasic site. If there are two closely spaced damages on opposing strands, such cleavage will reduce the size of the DNA on a nondenaturing gel. We show that ionizing radiation does induce clustered DNA damages containing abasic sites, oxidized purines, or oxidized pyrimidines. Further, the frequency of each of these cluster classes is comparable to that of frank double-strand breaks; among all complex damages induced by ionizing radiation, double-strand breaks are only about 20%, with other clustered damage constituting some 80%. We also show that even low doses (0.1–1 Gy) of high linear energy transfer ionizing radiation induce clustered damages in human cells.

Sensitive detection of DNA polymorphisms by the serial invasive signal amplification reaction

The invasive signal amplification reaction has been previously developed for quantitative detection of nucleic acids and discrimination of single-nucleotide polymorphisms. Here we describe a method that couples two invasive reactions into a serial isothermal homogeneous assay using fluorescence resonance energy transfer detection. The serial version of the assay generates more than 107 reporter molecules for each molecule of target DNA in a 4-h reaction; this sensitivity, coupled with the exquisite specificity of the reaction, is sufficient for direct detection of less than 1,000 target molecules with no prior target amplification. Here we present a kinetic analysis of the parameters affecting signal and background generation in the serial invasive signal amplification reaction and describe a simple kinetic model of the assay. We demonstrate the ability of the assay to detect as few as 600 copies of the methylene tetrahydrofolate reductase gene in samples of human genomic DNA. We also demonstrate the ability of the assay to discriminate single base differences in this gene by using 20 ng of human genomic DNA.