Insertion of L1 elements into sites that can form non-B DNA

Differences in Alu vs L1-rich chromosome bands underpin architectural reorganization of the inactive-X chromosome and SAHFs

The linear DNA sequence of mammalian chromosomes is organized in large blocks of DNA with similar sequence properties, producing a pattern of dark and light staining bands on mitotic chromosomes. Cytogenetic banding is essentially invariant between people and cell-types and thus may be assumed unrelated to genome regulation. We investigate whether large blocks of Alu-rich R-bands and L1-rich G-bands provide a framework upon which functional genome architecture is built. We examine two models of large-scale chromatin condensation: X-chromosome inactivation and formation of senescence-associated heterochromatin foci (SAHFs). XIST RNA triggers gene silencing but also formation of the condensed Barr Body (BB), thought to reflect cumulative gene silencing. However, we find Alu-rich regions are depleted from the L1-rich BB, supporting it is a dense core but not the entire chromosome. Alu-rich bands are also gene-rich, affirming our earlier findings that genes localize at the outer periphery of the BB. SAHFs similarly form within each territory by coalescence of syntenic L1 regions depleted for highly Alu-rich DNA. Analysis of senescent cell Hi-C data also shows large contiguous blocks of G-band and R-band DNA remodel as a segmental unit. Entire dark-bands gain distal intrachromosomal interactions as L1-rich regions form the SAHF. Most striking is that sharp Alu peaks within R-bands resist these changes in condensation. We further show that Chr19, which is exceptionally Alu rich, fails to form a SAHF. Collective results show regulation of genome architecture corresponding to large blocks of DNA and demonstrate resistance of segments with high Alu to chromosome condensation.

Mitochondrial DNA-like sequences in the nucleus (NUMTs): insights into our African origins and the mechanism of foreign DNA integration

Nuclear mitochondrial DNA sequences (NUMTs) are common in eukaryotes. However, the mechanism by which they integrate into the nuclear genome remains a riddle. We analyzed 247 NUMTs in the human nuclear DNA (nDNA), along with their flanking regions. This analysis revealed that some NUMTs have accumulated many changes, and thus have resided in the nucleus a long time, while others are >94% similar to the reference human mitochondrial DNA (mtDNA), and thus must be recent. Among the latter, two NUMTs, encompassing the COI gene, carry a set of transitions characteristic of the extant African-specific L macrohaplogroup mtDNAs and are more homologous to human mtDNA than to chimp. Screening for one of these NUMTs revealed its presence in all human samples tested, confirming that the African macrohaplogroup L mtDNAs were present in the earliest modern humans and thus were the first human mtDNAs. An analysis of flanking sequences of the NUMTs revealed that 59% were within 150 bp of repetitive elements, with 26% being within 15 bp of and 33% being within 15-150 bp of repetitive elements. Only 14% were integrated into a repetitive element. This association of NUMTs with repetitive elements is highly nonrandom (p<0.001). These data suggest that the vicinity of transposable elements influences the ongoing integration of mtDNA sequences and their subsequent duplication within the nDNA. Finally, NUMTs appear to preferentially integrate into DNA with different GC content than the surrounding chromosomal band. Our results suggest that chromosomal structure might influence integration of NUMTs.

L1 (LINE-1) retrotransposon evolution and amplification in recent human history

L1 (LINE-1) elements constitute a large family of mammalian retrotransposons that have been replicating and evolving in mammals for more than 100 Myr and now compose 20% or more of the DNA of some mammals. Here, we investigated the evolutionary dynamics of the active human Ta L1 family and found that it arose approximately 4 MYA and subsequently differentiated into two major subfamilies, Ta-0 and Ta-1, each of which contain additional subsets. Ta-1, which has not heretofore been described, is younger than Ta-0 and now accounts for at least 50% of the Ta family. Although Ta-0 contains some active elements, the Ta-1 subfamily has replaced it as the replicatively dominant subfamily in humans; 69% of the loci that contain Ta-1 inserts are polymorphic for the presence or absence of the insert in human populations, as compared with 29% of the loci that contain Ta-0 inserts. This value is 90% for loci that contain Ta-1d inserts, which are the youngest subset of Ta-1 and now account for about two thirds of the Ta-1 subfamily. The successive emergence and amplification of distinct Ta L1 subfamilies shows that L1 evolution has been as active in recent human history as it has been found to be for rodent L1 families. In addition, Ta-1 elements have been accumulating in humans at about the same rate per generation as recently evolved active rodent L1 subfamilies.

The biological properties and evolutionary dynamics of mammalian LINE-1 retrotransposons

Mammalian LINE-1 (L1) elements belong to the superfamily of autonomously replicating retrotransposable elements that lack the long terminal repeated (LTR) sequences typical of retroviruses and retroviral-like retrotransposons. The non-LTR superfamily is very ancient and L1-like elements are ubiquitous in nature, having been found in plants, fungi, invertebrates, and various vertebrate classes from fish to mammals. L1 elements have been replicating and evolving in mammals for at least the past 100 million years and now constitute 20% or more of some mammalian genomes. Therefore, L1 elements presumably have had a profound, perhaps defining, effect on the evolution, structure, and function of mammalian genomes. L1 elements contain regulatory signals and encode two proteins: one is an RNA-binding protein and the second one presumably functions as an integrase-replicase, because it has both endonuclease and reverse transcriptase activities. This work reviews the structure and biological properties of L1 elements, including their regulation, replication, evolution, and interaction with their mammalian hosts. Although each of these processes is incompletely understood, what is known indicates that they represent challenging and fascinating biological phenomena, the resolution of which will be essential for fully understanding the biology of mammals.

Tdr2, a new zebrafish transposon of the Tc1 family

We describe here Tdr2, a new class of Tc1-like transposons in zebrafish. Tdr2 was identified from the genomic sequence of a zebrafish PAC (P1 artificial chromosome) clone, and fragments of Tdr2 were found in several zebrafish EST (expressed sequence tag) sequences. Predicted translation of the Tdr2 transposase gene showed that it was most closely related to Caenorhabditis elegans Tc3A, suggesting an ancient origin of the Tdr2 transposon. Tdr2 spans 1. 1kb and is flanked by inverted repeats of approx. 100bp. The 5' repeat is itself composed of an inverted repeat, raising the possibility of the formation of a cruciform DNA structure. Tdr2 transposons may facilitate the development of novel transposon-based tools for the genetic analysis of zebrafish.

AT-rich sequences flanking the 5'-end breakpoint of the 4977-bp deletion of human mitochondrial DNA are located between two bent-inducing DNA sequences that assume distorted structure in organello

The 4977-bp deletion is the most common deletion among more than 90 large-scale deletions of human mitochondrial DNA (mtDNA) that are associated with aging and mitochondrial myopathies. The reason why the frequency of occurrence of this common deletion is so high in aged and myopathic human tissues is not clear. Since several studies proved that unusual DNA structures play very important roles in a number of recombination events, we hypothesized that some kind of unusual DNA structure may flank the breakpoints of the 4977-bp mtDNA deletion. We used two-dimensional (2-D) gel electrophoresis to assess the mobility abnormalities of the PCR-amplified DNA fragments encompassing the sequences of nucleotide position (np) 7901 to 9058 of human mtDNA. The results showed that the sequences of np 7901-8732 and np 8251-9058 exhibited retarded and increased mobilities, respectively, and that the sequence of np 8285-8676 showed normal mobility in the 2-D gel. This indicates that the 5'-end breakpoint of the 4977-bp deletion is located within the junction site of two flanking bent-inducing DNA sequences. We confirmed this notion by using osmium tetroxide (OsO4) to probe mtDNA in organello. The results showed that the two AT-rich sequences flanking the 5'-end breakpoint of the 4977-bp deletion are susceptible to OsO4 modification. These findings suggest that the DNA sequences of the 5'-end breakpoint of the common mtDNA deletion are rendered to assume a more distorted structure than B-DNA by these two flanking bent-inducing DNA sequences in organello and thereby render this region to be more vulnerable to attack by reactive oxygen species and free radicals.

De novo generation of simple sequence during gene amplification

Mammalian cells that have undergone gene amplification and/or gene rearrangement have been used as resources to gain insight into the questions of chromosome structure and dynamics. The multidrug resistant murine cell line J7.V2-1 has been shown previously to contain two distinct forms of the highly amplified mdr2 gene, a member of the mouse gene family responsible for the multidrug resistant (MDR) phenotype [Kirschner, L. S. (1995) DNA Cell Biol. 14, 47-59]. Characterization of both forms of the gene revealed that one form corresponded to the wild-type structure of the gene, whereas the other represented a rearrangement. Investigation of this altered gene demonstrated a deletion of 1.6 kb of the wild-type sequence, and replacement of this region with a poly(AT) tract that appears to have been generated de novo. Analysis of the native sequence in this region demonstrated the absence of repetitive elements, but was notable for the presence of two long stretches of polypurine: polypyrimidine strand asymmetry. Analysis of mdr2 transcripts in this cell line revealed that nearly all of the mRNA is transcribed from the rearranged form of the gene. This message is unable to code for a functional mdr2 gene product, owing to a deletion of the fourth exon during this event. Mechanisms of the rearrangement, as well as the significance of this curious effect on transcription, are discussed.

Formation of a combined H-DNA/open TATA box structure in the promoter sequence of the human Na,K-ATPase alpha2 gene

Structural variation of DNA within the promoter of the human Na, K-ATPase alpha2 gene, which contains a 35-base pair (bp) homopyrimidine.homopurine (Py.Pu) tract adjacent to a TATA box has been studied. The Py.Pu tract contains a 26-bp quasi-mirror repeat sequence with a potential for intramolecular triplex formation. As analyzed by two-dimensional agarose gel electrophoresis, a plasmid containing 151 bp of the promoter sequence including the 35-bp Py.Pu tract undergoes structural transitions under moderately acidic pH. Chemical probing with chloroacetaldehyde, dimethyl sulfate, and potassium permanganate is consistent with the formation of triplex DNA within the Py.Pu tract at native superhelical density as isolated from Escherichia coli. Chemical probing was used to determine a supercoil dependence for the formation of this combined unwound structure. At the superhelical density sufficient to locally unwind DNA, an H-y3 isomer of intermolecular triplex likely forms. However, at higher superhelical tension an H-y5 structure forms in the Py.Pu tract, and with increasing supercoiling the local DNA unwinding extends into the abutting TATA box. The H-y5/open TATA box combination structure might be favorable at higher superhelical densities since it relaxes more supercoils. The possible involvement of the H-y5/open TATA box structure in transcription is discussed.

The role of constrained self-organization in genome structural evolution

A hypothesis of genome structural evolution is explored. Rapid and cohesive alterations in genome organization are viewed as resulting from the dynamic and constrained interactions of chromosomal subsystem components. A combination of macromolecular boundary conditions and DNA element involvement in far-from-equilibrium reactions is proposed to increase the complexity of genomic subsystems via the channelling of genome turnover; interactions between subsystems create higher-order subsystems expanding the phase space for further genetic evolution. The operation of generic constraints on structuration in genome evolution is suggested by i) universal, homoplasic features of chromosome organization and ii) the metastable nature of genome structures where lower-level flux is constrained by higher-order structures. Phenomena such as 'genomic shock', bursts of transposable element activity, concerted evolution, etc., are hypothesized to result from constrained systemic responses to endogenous/exogenous, micro/macro perturbations. The constraints operating on genome turnover are expected to increase with chromosomal structural complexity, the number of interacting subsystems, and the degree to which interactions between genomic components are tightly ordered.

The chicken beta-globin gene promoter forms a novel "cinched" tetrahelical structure

We have previously shown that the G-rich sequence G16CG(GGT)2GG in the promoter region of the chicken beta-globin gene poses a formidable barrier to DNA synthesis in vitro (Woodford et al., 1994, J. Biol. Chem. 269, 27029-27035). The K+ requirement, template-strand specificity, template concentration independence, and involvement of Hoogsteen bonding suggested that the underlying basis of this new type of DNA synthesis arrest site might be an intrastrand tetrahelical structure. However, the arrest site lacks the four G-rich repeats that are a hallmark of previously described intramolecular tetraplexes and contains a number of noncanonical bases that would be expected to greatly destabilize such a structure. Here we report evidence for an unusual K+-dependent intrastrand "cinched" tetraplex. This structure has several unique features including the incorporation of bases other than guanine into the stem of the tetraplex, interaction between loop bases and bases in the flanking region, and base pairing between bases 3 and 5 of the tetrahelix-forming region to form a molecular "cinch." This finding extends the range of sequences capable of tetraplex formation as well as our appreciation of the conformational complexity of the chicken beta-globin promoter.

Rat karyotyping by fluorescence in situ hybridization (FISH): localization of oncogene c-raf to 4q42, retinoblastoma antioncogene to 15q12, and mitochondrial D-loop-like sequences to the Y chromosome

Fluorescence in situ hybridization (FISH) has been used to karyotype the rat genome with long interspersed repetitive elements (LINEs). Two-color FISH experiments were used to localize the rat c-raf oncogene to 4q42 and the rat retinoblastoma anti-oncogene to 15q12. In addition, sequences similar to the rat mitochondrial origin of replication (D-loop-like sequences) have been found to be concentrated among repetitive element sequences on the Y chromosome. This report extends hybridization-based karyotype technology to the rat, thus facilitating the application of FISH to rat gene mapping.

The evolution of long interspersed repeated DNA (L1, LINE 1) as revealed by the analysis of an ancient rodent L1 DNA family

All modern mammals contain a distinctive, highly repeated (> or = 50,000 members) family of long interspersed repeated DNA called the L1 (LINE 1) family. While the modern L1 families were derived from a common ancestor that predated the mammalian radiation approximately 80 million years ago, most of the members of these families were generated within the last 5 million years. However, recently we demonstrated that modern murine (Old World rats and mice) genomes share an older long interspersed repeated DNA family that we called Lx. Here we report our analysis of the DNA sequence of Lx family members and the relationship of this family to the modern L1 families in mouse and rat. The extent of DNA sequence divergence between Lx members indicates that the Lx amplification occurred about 12 million years ago, around the time of the murine radiation. Parsimony analysis revealed that Lx elements were ancestral to both the modern rat and mouse L1 families. However, we found that few if any of the evolutionary intermediates between the Lx and the modern L1 families were extensively amplified. Because the modern L1 families have evolved under selective pressure, the evolutionary intermediates must have been capable of replication. Therefore, replication-competent L1 elements can reside in genomes without undergoing extensive amplification. We discuss the bearing of our findings on the evolution of L1 DNA elements and the mammalian genome.

A survey of protein-DNA interaction sites within the murine immunoglobulin heavy chain locus reveals a particularly complex pattern around the DQ52 element

The expression of immunoglobulin (Ig) genes is regulated at two levels: rearrangement of individual gene segments and transcription of continuous genes. To find transacting factors involved in mediating locus- and segment specific gene activation and expression, we surveyed a 3600-bp genomic region of the murine Ig heavy chain locus, spanning from the DQ52 element to the Ig heavy chain intron enhancer. We discovered nine, previously undescribed, protein-DNA complexes and estimated their individual binding-affinity preferences (Kr) by quantitative gel shift measurements. We observed one novel protein DNA interaction at the enhancer, two in the JH1 region and six within a 300-bp region immediately 5' to the DQ52 locus. The latter show a complex and specific binding pattern when comparing nuclear extracts derived from pre-B cells and fibroblasts. Further characterization of the interactions at the DQ52 locus by electron microscopy revealed the preferential formation of a protein complex binding to the DQ52 locus with pre B cell extracts. This behavior and the clustering of interaction sites 5' to the DQ52 element suggest that this region is involved in the regulation of heavy chain gene assembly.

Genome canalization: the coevolution of transposable and interspersed repetitive elements with single copy DNA

Transposable and interspersed repetitive elements (TIREs) are ubiquitous features of both prokaryotic and eukaryotic genomes. However, controversy has arisen as to whether these sequences represent useless 'selfish' DNA elements, with no cellular function, as opposed to useful genetic units. In this review, we selected two insect species, the Dipteran Drosophila and the Lepidopteran Bombyx mori (the silkmoth), in an attempt to resolve this debate. These two species were selected on the basis of the special interest that our laboratory has had over the years in Bombyx with its well known molecular and developmental biology, and the wealth of genetic data that exist for Drosophila. In addition, these two species represent contrasting repetitive element types and patterns of distribution. On one hand, Bombyx exhibits the short interspersion pattern in which Alu-like TIREs predominate while Drosophila possesses the long interspersion pattern in which retroviral-like TIREs are prevalent. In Bombyx, the main TIRE family is Bm-1 while the Drosophila group contains predominantly copia-like elements, non-LTR retroposons, bacterial-type retroposons and fold-back transposable elements sequences. Our analysis of the information revealed highly non-random patterns of both TIRE biology and evolution, more indicative of these sequences acting as genomic symbionts under cellular regulation rather than useless or selfish junk DNA. In addition, we extended our analysis of potential TIRE functionality to what is known from other eukaryotic systems. From this study, it became apparent that these DNA elements may have originated as innocuous or selfish sequences and then adopted functions. The mechanism for this conversion from non-functionality to specific roles is a process of coevolution between the repetitive element and other cellular DNA often times in close physical proximity. The resulting interdependence between repetitive elements and other cellular sequences restrict the number of evolutionarily successful mutational changes for a given function or cistron. This mutual limitation is what we call genome canalization. Well documented examples are discussed to support this hypothesis and a mechanistic model is presented for how such genomic canalization can occur. Also proposed are empirical studies which would support or invalidate aspects of this hypothesis.

Base-pairing shift in the major groove of (CA)n tracts by B-DNA crystal structures

the crystal packing of the B-DNA dodecamer d(ACCG-GCGCCACA).d(TGTGGCGCCGGT) is characterized by the reciprocal fit of double helices with specific base-backbone interactions in the major groove. Cooling the crystals below -10 degrees C stabilizes a new conformational state with a long-range sequence-dependent one-step shift in the major-groove base pairing. The tilt of the bases leads to the disruption of the Watson-Crick pairing in the major groove and to the formation of interactions with the 5' neighbour of their complement. This alteration propagates along the helical axis over more than half a turn. As a result, the molecular structure is normal when seen from the minor groove side and mismatched in the major groove. Comparison with a parent isomorphous dodecamer structure corresponding to the codon 10-13 of the c-Ha-ras proto-oncogene show that this new structural feature is sequence dependent and clearly favoured by (CA)n tracts. As(CA)n tracts of DNA are involved both in recombination and in transcription, this new recognition pattern should be considered in the analysis of the various processes involving the reading of the genetic information.