Statistical methods for the analysis of DNA sequences: application to dinucleotide distribution in the human genome

Questa tesi si inserisce nell'ambito delle analisi statistiche e dei metodi stocastici applicati all'analisi delle sequenze di DNA. Nello specifico il nostro lavoro è incentrato sullo studio del dinucleotide CG (CpG) all'interno del genoma umano, che si trova raggruppato in zone specifiche denominate CpG islands. Queste sono legate alla metilazione del DNA, un processo che riveste un ruolo fondamentale nella regolazione genica. La prima parte dello studio è dedicata a una caratterizzazione globale del contenuto e della distribuzione dei 16 diversi dinucleotidi all'interno del genoma umano: in particolare viene studiata la distribuzione delle distanze tra occorrenze successive dello stesso dinucleotide lungo la sequenza. I risultati vengono confrontati con diversi modelli nulli: sequenze random generate con catene di Markov di ordine zero (basate sulle frequenze relative dei nucleotidi) e uno (basate sulle probabilità di transizione tra diversi nucleotidi) e la distribuzione geometrica per le distanze. Da questa analisi le proprietà caratteristiche del dinucleotide CpG emergono chiaramente, sia dal confronto con gli altri dinucleotidi che con i modelli random. A seguito di questa prima parte abbiamo scelto di concentrare le successive analisi in zone di interesse biologico, studiando l’abbondanza e la distribuzione di CpG al loro interno (CpG islands, promotori e Lamina Associated Domains). Nei primi due casi si osserva un forte arricchimento nel contenuto di CpG, e la distribuzione delle distanze è spostata verso valori inferiori, indicando che questo dinucleotide è clusterizzato. All’interno delle LADs si trovano mediamente meno CpG e questi presentano distanze maggiori. Infine abbiamo adottato una rappresentazione a random walk del DNA, costruita in base al posizionamento dei dinucleotidi: il walk ottenuto presenta caratteristiche drasticamente diverse all’interno e all’esterno di zone annotate come CpG island. Riteniamo pertanto che metodi basati su questo approccio potrebbero essere sfruttati per migliorare l’individuazione di queste aree di interesse nel genoma umano e di altri organismi.

PCR identification of HIV-1 DNA sequences in brain tissue of patients with AIDS encephalopathy.

We analyzed brain tissue from 39 patients for the presence of proviral HIV-1 sequences, using the polymerase chain reaction (PCR) for the amplification of segments of the viral LTR and gag genes. A novel primer extension procedure allowed the detection of a single HIV-1 copy in 1 micrograms DNA. We detected proviral HIV-1 DNA in 16 of 25 brain samples from AIDS patients. Semiquantitative evaluation of the amplified DNAs indicated considerable variation in viral load. Highest levels of proviral DNA were present in brain samples from six patients with clinical evidence of HIV-associated cognitive/motor complex and the histopathologic correlate of HIV leukoencephalopathy or HIV encephalitis. An additional 11 brain samples contained smaller amounts of proviral DNA. In these patients, clinical data were inconclusive regarding the diagnosis of HIV-1 encephalopathy and histopathologically there was no evidence of HIV-1-induced tissue lesions. In nine of 25 seropositive patients with AIDS (36%), brain samples scored negative or did not contain an unequivocal signal indicating the presence of proviral DNA. HIV-1 sequences were not detected in any of 14 control brain samples from HIV-1 seronegative patients. Our data indicate that HIV-1 is present in the central nervous system of the majority (two thirds) of AIDS patients and that the highest levels of proviral DNA in brain tissue are associated with HIV encephalopathy.

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.

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 ^

Minimal DNA sequences that control the cell lineage-specific expression of the pro$\alpha$2(I) collagen promoter in transgenic mice

The pattern of expression of the pro$\alpha$2(I) collagen gene is highly tissue-specific in adult mice and shows its strongest expression in bones, tendons, and skin. Transgenic mice were generated harboring promoter fragments of the mouse pro$\alpha$2(I) collagen gene linked to the Escherichia coli $\beta$-galactosidase or firefly luciferase genes to examine the activity of these promoters during development. A region of the mouse pro$\alpha$2(I) collagen promoter between $-$2000 and +54 exhibited a pattern of $\beta$-galactosidase activity during embryonic development that corresponded to the expression pattern of the endogenous pro$\alpha$2(I) collagen gene as determined by in situ hybridization. A similar pattern of activity was also observed with much smaller promoter fragments containing either 500 or 350 bp of upstream sequence relative to the start of transcription. Embryonic regions expressing high levels of $\beta$-galactosidase activity included the valves of the developing heart, sclerotomes, meninges, limb buds, connective tissue fascia between muscle fibers, osteoblasts, tendon, periosteum, dermis, and peritoneal membranes. The pattern of $\beta$-galactosidase activity was similar to the extracellular immunohistochemical localization of transforming growth factor-$\beta$1 (TGF-$\beta$1). The $-$315 to $-$284 region of the pro$\alpha$2(I) collagen promoter was previously shown to mediate the stimulatory effects of TGF-$\beta$1 on the pro$\alpha$2(I) collagen promoter in DNA transfection experiments with cultured fibroblasts. A construct containing this sequence tandemly repeated 5$\sp\prime$ to both a very short $\alpha$2(I) collagen promoter ($-$40 to +54) and a heterologous minimal promoter showed preferential activity in tail and skin of 4-week old transgenic mice. The pattern of expression mimics that of the $-$350 to +54 pro$\alpha$2(I) collagen promoter linked to a luciferase reporter gene in transgenic mice. ^

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Repetitive-DNA elements are similarly distributed on Caenorhabditis elegans autosomes

The positions of ≈4,800 individual miniature inverted-repeat transposable element (MITE)-like repeats from four families were mapped on the Caenorhabditis elegans chromosomes. These families represent 1–2% of the total sequence of the organism. The four MITE families (Cele1, Cele2, Cele14, and Cele42) displayed distinct chromosomal distribution profiles. For example, the Cele14 MITEs were observed clustering near the ends of the autosomes. In contrast, the Cele2 MITEs displayed an even distribution through the central autosome domains, with no evidence for clustering at the ends. Both the number of elements and the distribution patterns of each family were conserved on all five C. elegans autosomes. The distribution profiles indicate chromosomal polarity and suggest that the current genetic and physical maps of chromosomes II, III, and X are inverted with respect to the other chromosomes. The degree of conservation of both the number and distribution of these elements on the five autosomes suggests a role in defining specific chromosomal domains.

Analysis of microsatellite mutations in the mitochondrial DNA of Saccharomyces cerevisiae

In the nuclear genome of Saccharomyces cerevisiae, simple, repetitive DNA sequences (microsatellites) mutate at rates much higher than nonrepetitive sequences. Most of these mutations are deletions or additions of repeat units. The yeast mitochondrial genome also contains many microsatellites. To examine the stability of these sequences, we constructed a reporter gene (arg8m) containing out-of-frame insertions of either poly(AT) or poly(GT) tracts within the coding sequence. Yeast strains with this reporter gene inserted within the mitochondrial genome were constructed. Using these strains, we showed that poly(GT) tracts were considerably less stable than poly(AT) tracts and that alterations usually involved deletions rather than additions of repeat units. In contrast, in the nuclear genome, poly(GT) and poly(AT) tracts had similar stabilities, and alterations usually involved additions rather than deletions. Poly(GT) tracts were more stable in the mitochondria of diploid cells than in haploids. In addition, an msh1 mutation destabilized poly(GT) tracts in the mitochondrial genome.

Differences in water release for the binding of EcoRI to specific and nonspecific DNA sequences.

The free energy difference between complexes of the restriction nuclease EcoRI with nonspecific DNA and with the enzyme's recognition sequence is linearly dependent on the water chemical potential of the solution, set using several very different solutes, ranging from glycine and glycerol to triethylene glycol and sucrose. This osmotic dependence indicates that the nonspecific complex sequesters some 110 waters more than the specific complex with the recognition sequence. The insensitivity of the difference in number of waters released to the solute identity further indicates that this water is sequestered in a space that is sterically inaccessible to solutes, most likely at the protein-DNA interface of the nonspecific complex. Calculations based on the structure of the specific complex suggest that the apposing DNA and protein surfaces in the nonspecific complex retain approximately a full hydration layer of water.

Does behavior reflect phylogeny in swiftlets (Aves: Apodidae)? A test using cytochrome b mitochondrial DNA sequences.

Swiftlets are small insectivorous birds, many of which nest in caves and are known to echolocate. Due to a lack of distinguishing morphological characters, the taxonomy of swiftlets is primarily based on the presence or absence of echolocating ability, together with nest characters. To test the reliability of these behavioral characters, we constructed an independent phylogeny using cytochrome b mitochondrial DNA sequences from swiftlets and their relatives. This phylogeny is broadly consistent with the higher classification of swifts but does not support the monophyly of swiftlets. Echolocating swiftlets (Aerodramus) and the nonecholocating "giant swiftlet" (Hydrochous gigas) group together, but the remaining nonecholocating swiftlets belonging to Collocalia are not sister taxa to these swiftlets. While echolocation may be a synapomorphy of Aerodramus (perhaps secondarily lost in Hydrochous), no character of Aerodramus nests showed a statistically significant fit to the molecular phylogeny, indicating that nest characters are not phylogenetically reliable in this group.

The effects of sequence context on DNA curvature.

Recent experiments have exposed significant discrepancies between experimental data and predictive models for DNA structure. These results strongly suggest that DNA structural parameters incorporated in the models are not always sufficient to account for the influence of sequence context and of specific ion effects. In an attempt to evaluate these two effects, we have investigated repetitive DNA sequences with the sequence motif GAGAG.CTCTC located in different helical phasing arrangements with respect to poly(A) tracts and GGGCCC.GGGCCC sequence motifs. Methods used are ligase-mediated cyclization and gel mobility experiments along with DNase I cutting and chemical probe studies. The results provide new evidence for curvature in poly(A) tracts. They also show that the sequence context in which bending and flexible sequence elements are found is an important aspect of sequence-dependent DNA conformation. Although dinucleotide models generally have good predictive power, this work demonstrates that in some instances sequence elements larger than the dinucleotide must be taken into account, and hence it provides a starting point for the appropriate modification and refinement of existing structural models for DNA.

DNA sequences of Alu elements indicate a recent replacement of the human autosomal genetic complement.

DNA sequences of neutral nuclear autosomal loci, compared across diverse human populations, provide a previously untapped perspective into the mode and tempo of the emergence of modern humans and a critical comparison with published clonally inherited mitochondrial DNA and Y chromosome measurements of human diversity. We obtained over 55 kilobases of sequence from three autosomal loci encompassing Alu repeats for representatives of diverse human populations as well as orthologous sequences for other hominoid species at one of these loci. Nucleotide diversity was exceedingly low. Most individuals and populations were identical. Only a single nucleotide difference distinguished presumed ancestral alleles from descendants. These results differ from those expected if alleles from divergent archaic populations were maintained through multiregional continuity. The observed virtual lack of sequence polymorphism is the signature of a recent single origin for modern humans, with general replacement of archaic populations.

Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

Coiled bodies contain U7 small nuclear RNA and associate with specific DNA sequences in interphase human cells.

Coiled bodies (CBs) are nuclear organelles whose structures appear to be highly conserved in evolution. In rapidly cycling cells, they are typically located in the nucleoplasm but are often found in contact with the nucleolus. The CBs in human cells contain a unique protein, called p80-coilin. Studies on amphibian oocyte nuclei have revealed a protein within the "sphere" organelle that shares significant structural similarity to p80-coilin. Spheres and CBs are also highly enriched in small nuclear ribonucleoproteins and other RNA-processing components. We present evidence that, like spheres, CBs contain U7 small nuclear RNA (snRNA) and associate with specific chromosomal loci. Using biotinylated 2'-O-methyl oligonucleotides complementary to the 5' end of U7 snRNA and fluorescence in situ hybridization, we show that U7 is distributed throughout the nucleoplasm, excluding nucleoli, and is concentrated in CBs. Interestingly, we found that CBs often associate with subsets of the histone, U1, and U2 snRNA gene loci in interphase HeLa-ATCC and HEp-2 monolayer cells. However, in a strain of suspension-grown HeLa cells, called HeLa-JS1000, we found a much lower rate of association between CBs and snRNA genes. Possible roles for CBs in the metabolism of these various histone and snRNAs are discussed.