LZ complexity distance of DNA sequences and its application in phylogenetic tree reconstruction

DNA sequences can be treated as finite-length symbol strings over a four-letter alphabet （A, C, T, G）. As a universal and computable complexity measure, LZ complexity is valid to describe the complexity of DNA sequences. In this study, a concept of conditional LZ complexity between two sequences is proposed according to the principle of LZ complexity measure. An LZ complexity distance metric between two nonnull sequences is defined by utilizing conditional LZ complexity. Based on LZ complexity distance, a phylogenetic tree of 26 species of placental mammals （Eutheria） with three outgroup species was reconstructed from their complete mitochondrial genomes. On the debate that which two of the three main groups of placental mammals, namely Primates, Ferungulates, and Rodents, are more closely related, the phylogenetic tree reconstructed based on LZ complexity distance supports the suggestion that Primates and Ferungulates are more closely related.     

Sequence complexity profiles of prokaryotic genomic sequences: a fast algorithm for calculating linguistic complexity

MOTIVATION: One of the major features of genomic DNA sequences, distinguishing them from texts in most spoken or artificial languages, is their high repetitiveness. Variation in the repetitiveness of genomic texts reflects the presence and density of different biologically important messages. Thus, deviation from an expected number of repeats in both directions indicates a possible presence of a biological signal. Linguistic complexity corresponds to repetitiveness of a genomic text, and potential regulatory sites may be discovered through construction of typical patterns of complexity distribution. RESULTS: We developed software for fast calculation of linguistic sequence complexity of DNA sequences. Our program utilizes suffix trees to compute the number of subwords present in genomic sequences, thereby allowing calculation of linguistic complexity in time linear in genome size. The measure of linguistic complexity was applied to the complete genome of Haemophilus influenzae. Maps of complexity along the entire genome were obtained using sliding windows of 40, 100, and 2000 nucleotides. This approach provided an efficient way to detect simple sequence repeats in this genome. In addition, local profiles of complexity distribution around the starts of translation were constructed for 21 complete prokaryotic genomes. We hypothesize that complexity profiles correspond to evolutionary relationships between organisms. We found principal differences in profiles of the GC-rich and other (non-GC-rich) genomes. We also found characteristic differences in profiles of AT genomes, which probably reflect individual species variations in translational regulation. AVAILABILITY: The program is available upon request from Alexander Bolshoy or at http://csweb.haifa.ac.il/library/#complex.     

Comparative analysis of polyadenylated RNA complexity in soybean hypocotyl tissue and cultured suspension cells

Growth parameters of suspension culture cells of soybean (Glycine max L.) were compared between cells grown in medium with (+) auxin and without (−) auxin. Growth rates were greater for (+) auxin cells. Cells transferred to (−) auxin medium primarily expanded in size while (+) auxin cells initially divided and then expanded. Two methods were used to estimate polyadenylated RNA sequence complexity. Kinetic analysis gave a sum of component complexity values of 36,000 and 64,000 diverse poly(A) RNA sequences of about 1,400 nucleotides in (+) and (−) auxin grown cells, respectively. The most striking difference between these cell populations was the increase in the poly(A) RNA sequence complexity in cells grown in medium without auxin. RNA complexities were also determined by the saturation of `single' copy DNA by poly(A) RNAs from (+) and (−) auxin suspension cells. These saturation studies estimated the total complexity of (+) and (−) auxin suspension cells as 41,000 and 57,000 diverse sequences, respectively. Suspension cells in auxin-depleted medium produced about 20,000 more diverse sequences than (+) auxin cells. Comparisons of poly(A) complexities were also made among auxin-treated and untreated hypocotyl cells from the intact plant relative to suspension culture cells. Mixed populations of poly(A) RNA from these tissues and cells allowed the determination of shared sequences among them. When all combinations of poly(A) RNA were mixed, the percentage of `single' copy DNA that saturated was equivalent to diverse sequence complexity estimates of about 60,000. When mixed poly(A) RNA from suspension cells from (+) and (−) auxin medium were compared, they shared about 40,000 sequences and (−) auxin cells contained an additional 20,000. Both (+) and (−) tissue culture cells shared a subset of about 20,000 sequences with cells from (+) and (−) auxin treated hypocotyl. A third subset of about 20,000 sequences was shared by (−) auxin suspension cells and hypocotyl treated with or without auxin, a subset most of which were not shared by (+) auxin suspension cells. Kinetic and saturation data estimates of poly(A) RNA complexity compared favorably and indicated that exogenous auxin treatment can dramatically alter the complexity of all classes of poly(A) RNAs in cultured cells.     

Lack of Major Genome Instability in Tumors of p53 Null Rats

Tumorigenesis is often associated with loss of tumor suppressor genes (such as TP53), genomic instability and telomere lengthening. Previously, we generated and characterized a rat p53 knockout model in which the homozygous rats predominantly develop hemangiosarcomas whereas the heterozygous rats mainly develop osteosarcomas. Using genome-wide analyses, we find that the tumors that arise in the heterozygous and homozygous Tp53^C273X mutant animals are also different in their genomic instability profiles. While p53 was fully inactivated in both heterozygous and homozygous knockout rats, tumors from homozygous animals show very limited aneuploidy and low degrees of somatic copy number variation as compared to the tumors from heterozygous animals. In addition, complex structural rearrangements such as chromothripsis and breakage-fusion-bridge cycles were never found in tumors from homozygous animals, while these were readily detectable in tumors from heterozygous animals. Finally, we measured telomere length and telomere lengthening pathway activity and found that tumors of homozygous animals have longer telomeres but do not show clear telomerase or alternative lengthening of telomeres (ALT) activity differences as compared to the tumors from heterozygous animals. Taken together, our results demonstrate that host p53 status in this rat p53 knockout model has a large effect on both tumor type and genomic instability characteristics, where full loss of functional p53 is not the main driver of large-scale structural variations. Our results also suggest that chromothripsis primarily occurs under p53 heterozygous rather than p53 null conditions.     

DNA strand break-sensing molecule poly(ADP-Ribose) polymerase cooperates with p53 in telomere function, chromosome stability, and tumor suppression

Genomic instability is often caused by mutations in genes that are involved in DNA repair and/or cell cycle checkpoints, and it plays an important role in tumorigenesis. Poly(ADP-ribose) polymerase (PARP) is a DNA strand break-sensing molecule that is involved in the response to DNA damage and the maintenance of telomere function and genomic stability. We report here that, compared to single-mutant cells, PARP and p53 double-mutant cells exhibit many severe chromosome aberrations, including a high degree of aneuploidy, fragmentations, and end-to-end fusions, which may be attributable to telomere dysfunction. While PARP(-/-) cells showed telomere shortening and p53(-/-) cells showed normal telomere length, inactivation of PARP in p53(-/-) cells surprisingly resulted in very long and heterogeneous telomeres, suggesting a functional interplay between PARP and p53 at the telomeres. Strikingly, PARP deficiency widens the tumor spectrum in mice deficient in p53, resulting in a high frequency of carcinomas in the mammary gland, lung, prostate, and skin, as well as brain tumors, reminiscent of Li-Fraumeni syndrome in humans. The enhanced tumorigenesis is likely to be caused by PARP deficiency, which facilitates the loss of function of tumor suppressor genes as demonstrated by a high rate of loss of heterozygosity at the p53 locus in these tumors. These results indicate that PARP and p53 interact to maintain genome integrity and identify PARP as a cofactor for suppressing tumorigenesis.     

Combined effects of leukocyte telomere length,p53 polymorphism and human papillomavirus infection on esophageal squamous cell carcinoma in a Han Chinese population

Telomere shortening has been suggested to be a genetic predictor for various cancers. However, evidences about this point with respect to esophageal squamous cell carcinoma (ESCC) in Han Chinese populations remain limited. Our previous study demonstrated that p53 Arg72Pro polymorphism was associated with the risk of human papillomavirus (HPV)-related ESCC. Telomeres and p53 play important roles in maintaining genomic stability and regulating the cell cycle. HPV impacts both telomere length stabilization and p53 degradation. Given the roles of the three factors, we evaluated leukocyte telomere length, p53 variants and HPV-16 serology to examine the potential associations between them and ESCC risk in a case–control study with 308 patients and 309 cancer-free controls matched by age and sex. Compared with long telomere length, short telomere length was significantly associated with an increased risk of ESCC (adjusted OR 2.01; 95% CI 1.41–2.80). Moreover, this association was enhanced when combined with HPV-16 seropositivity and p53 Arg/Arg or Arg/Pro genotypes. Notably, individuals with short telomere length, Arg/Pro or Arg/Arg genotypes and HPV-16 seropositivity had a 12.08-fold (95% CI 5.49–26.56) increased risk of ESCC compared to those with none of the three investigated risk factors. Taken together, these results indicate that short telomere length in peripheral blood leukocytes is a biomarker for ESCC risk, and has statistically additive effects with p53 variants and HPV seropositivity with regard to the risk of ESCC in a Han Chinese population.     

The role of p53-mediated apoptosis as a crucial anti-tumor response to genomic instability: lessons from mouse models

Genomic instability is a major force driving human cancer development. A cellular safeguard against such genetic destabilization, which can ensue from defects in telomere maintenance, DNA repair, and checkpoint function, is activation of the p53 tumor suppressor protein, which commonly responds to these DNA damage signals by inducing apoptosis. If, however, p53 becomes inactivated, as is typical of many tumors and pre-cancerous lesions, then cells with compromised genome integrity pathways survive inappropriately, and the accrual of oncogenic lesions can fuel the carcinogenic process. Studies of mouse models have been instrumental in providing support for this idea. Mouse knockouts in genes important for telomere function, DNA damage checkpoint activation and DNA repair - both non-homologous end joining and homologous recombination - are prone to the development of genomic instability. As a consequence of these DNA damage signals, p53 becomes activated in cells of these mutant mice, leading to the induction of apoptosis, sometimes at the expense of organismal viability. This apoptotic response can be rescued through crosses to p53-deficient mice, but has dire consequences: mice predisposed to genomic instability and lacking p53 are susceptible to tumorigenesis. Thus p53-mediated apoptosis provides a crucial tumor suppressive mechanism to eliminate cells succumbing to genomic instability.     

Oxidative Stress Induces Persistent Telomeric DNA Damage Responsible for Nuclear Morphology Change in Mammalian Cells

One main function of telomeres is to maintain chromosome and genome stability. The rate of telomere shortening can be accelerated significantly by chemical and physical environmental agents. Reactive oxygen species are a source of oxidative stress and can produce modified bases (mainly 8-oxoG) and single strand breaks anywhere in the genome. The high incidence of guanine residues in telomeric DNA sequences makes the telomere a preferred target for oxidative damage. Our aim in this work is to evaluate whether chromosome instability induced by oxidative stress is related specifically to telomeric damage. We treated human primary fibroblasts (MRC-5) in vitro with hydrogen peroxide (100 and 200 µM) for 1 hr and collected data at several time points. To evaluate the persistence of oxidative stress-induced DNA damage up to 24 hrs after treatment, we analysed telomeric and genomic oxidative damage by qPCR and a modified comet assay, respectively. The results demonstrate that the genomic damage is completely repaired, while the telomeric oxidative damage persists. The analysis of telomere length reveals a significant telomere shortening 48 hrs after treatment, leading us to hypothesise that residual telomere damage could be responsible for the telomere shortening observed. Considering the influence of telomere length modulation on genomic stability, we quantified abnormal nuclear morphologies (Nucleoplasmic Bridges, Nuclear Buds and Micronuclei) and observed an increase of chromosome instability in the same time frame as telomere shortening. At subsequent times (72 and 96 hrs), we observed a restoration of telomere length and a reduction of chromosome instability, leaving us to conjecture a correlation between telomere shortening/dysfunction and chromosome instability. We can conclude that oxidative base damage leads to abnormal nuclear morphologies and that telomere dysfunction is an important contributor to this effect.     

Sequence complexity of cDNA transcribed from a diverse mRNA population

Mouse liver poly(A)+mRNA was reverse transcribed using oligo-p(dT) or random oligonucleotides as primers to yield cDNA about equal to the mass of the template RNA. The size profile of the oligo-p(dT)-primedd cDNA was similar to that of the template RNA. RNA or cDNA driven saturation annealing of labeled single copy genomic DNA (scDNA) showed that 2% of the scDNA was complementary in either case indicating the sequence complexity of cDNA was equivalent to that of the template mRNA. These results establish for the first time that cDNA represents essentially all of the sequence complexity of a diverse template RNA population in which individual mRNA species are present in vastly different concentrations. RNA driven hydridization of the cDNA showed that about 40% of the cDNA mass represents most of the sequence complexity of the template RNA. Also, kinetics of this hybridization indicate a complexity of 58,000 kb for the template RNA, a value similar to that obtained by scDNA hybridization. We conclude that appropriately characterized cDNA probes can be used to make valid qualitative and quantitative comparisons of the complex, infrequent class mRNAs of different cells and tissues.     

Re-usable DNA template for the polymerase chain reaction.

DNA covalently bound to an uncharged nylon membrane was used for consecutive amplifications of several different genes by PCR. Successful PCR amplifications were obtained for membrane-bound genomic and plasmid DNA. Membrane-bound genomic DNA templates were re-used at least 15 times for PCR with specific amplification of the desired gene each time. PCR amplifications of specific sequences of p53, p16, CYP1A1, CYP2D6, GSTM1 and GSTM3 were performed independently on the same strips of uncharged nylon membrane containing genomic DNA. PCR products were subjected to restriction fragment length polymorphism analysis, single-strand conformational polymorphism analysis and/or dideoxy sequencing to confirm PCR-amplified gene sequences. We found that PCR fragments obtained by amplification from bound genomic DNA as template were identical in sequence to those of PCR products obtained from free genomic DNA in solution. PCR was performed using as little as 5 ng genomic or 4 fg plasmid DNA bound to membrane. These results suggest that DNA covalently bound to membrane can be re-used for sample-specific PCR amplifications, providing a potentially unlimited source of DNA for PCR.     

Signaling pathway requirements for induction of senescence by telomere homolog oligonucleotides

Cellular senescence is a major defense against cancer. In human fibroblasts, suppressing both the p53 and pRb pathways is necessary to bypass replicative senescence as well as senescence induced by ectopic expression of a dominant negative form of the telomere repeat binding factor 2, TRF2(DN). We recently reported that exposure to oligonucleotides homologous to the telomere 3' overhang (T-oligos) activates both the p53 and pRb pathways and leads to senescence in primary human fibroblasts. To further characterize T-oligo-induced senescence, we compared established isogenic fibroblast cell lines lacking functional p53 and/or pRb pathways to the normal parental line. Here, we report that, as in physiologic senescence, inactivation of both the p53 and pRb pathways is necessary to suppress T-oligo-induced senescence. Moreover, T-oligo rapidly induces senescence in a malignant fibroblast-derived cell line, demonstrating the potential of using T-oligo as a novel anticancer therapeutic. Our data support the hypothesis that exposure of the TTAGGG tandem repeat telomere 3' overhang sequence is the event that initiates signaling through DNA damage response pathways after experimental telomere disruption, serial passage, or acute genomic damage of normal cells.     

Cloning and genetic mapping of wheat telomere-associated sequences

Wheat telomere-associated sequences (TASs) were cloned using a Vectorette approach and sequenced. Reverse primers specific to the TASs were combined with labelled degenerate telomere primers in PCR reactions containing total genomic DNA as template. Amplification products were separated on sequencing gels. In total, seventeen primer combinations provided 47 polymorphic fragments. Nine of these mapped beyond the most distal RFLP markers and defined the ends of seven chromosome arms. Seven of the nine terminal fragments were derived from a 118-bp tandem repeat, indicating that subtelomeric tandem repeat sequences provide an efficient means to target chromosome ends. A telomere cloning strategy and the terminal and interstitial location of TASs are discussed. Received: 13 September 1996 / Accepted: 22 January 1997     

Kinetic complexity of chloroplastal deoxyribonucleic acid and mitochondrial deoxyribonucleic acid from higher plants

1. Chloroplasts and mitochondria were isolated by aqueous and non-aqueous cell-fractionation techniques. In a variety of higher plants the mitochondrial DNA bands in a caesium chloride gradient at 1.706g.cm.(-3), whereas chloroplastal DNA has a buoyant density of 1.697g.cm.(-3). 2. In total cellular DNA of moderate molecular weight, the chloroplastal DNA is found within the Gaussian distribution of the nuclear DNA and is not resolved as a satellite. 3. Both chloroplastal DNA and mitochondrial DNA from lettuce renature rapidly. 4. The kinetic complexity of mitochondrial DNA is > 10(8) daltons. 5. Chloroplastal DNA is made up from fast and slow renaturing sequences with kinetic complexities of 3x10(6) and 1.2x10(8) daltons respectively. 6. From the discrepancy between analytical and kinetic complexity it is concluded that chloroplastal DNA is extensively reiterated.     

Dual roles of telomere dysfunction in initiation and suppression of tumorigenesis

Human carcinomas arise through the acquisition of genetic changes that endow precursor cancer cells with a critical threshold of cancer-relevant genetic lesions. This complex genomic alterations confer upon precursor cancer cells the ability to grow indefinitely and to metastasize to distant sites. One important mechanism underlying a cell's tumorigenic potential is the status of its telomere. Telomeres are G-rich simple repeat sequences that serve to prevent chromosomal ends from being recognized as DNA double-strand breaks (DSBs). Dysfunctional telomeres resemble DSBs, leading to the formation of dicentric chromosomes that fuel high degrees of genomic instability. In the setting of an intact p53 pathway, this instability promotes cellular senescence, a potent tumor suppressor mechanism. However, rare cells that stochastically lose p53 function emerge from this sea of genomic instability and progress towards cancer. In this review, we describe the use of mouse models to probe the impact of dysfunctional telomeres on tumor initiation and suppression.