Structural rearrangements and insertions of dispersed elements in pericentromeric alpha satellites occur preferably at kinkable DNA sites

Mechanisms of LTR-Retroelement Transposition: Lessons from Drosophila melanogaster

Long terminal repeat (LTR) retrotransposons occupy a special place among all mobile genetic element families. The structure of LTR retrotransposons that have three open reading frames is identical to DNA forms of retroviruses that are integrated into the host genome. Several lines of evidence suggest that LTR retrotransposons share a common ancestry with retroviruses and thus are highly relevant to understanding mechanisms of transposition. Drosophila melanogaster is an exceptionally convenient model for studying the mechanisms of retrotransposon movement because many such elements in its genome are transpositionally active. Moreover, two LTRretrotransposons of D. melanogaster, gypsy and ZAM, have been found to have infectious properties and have been classified as errantiviruses. Despite numerous studies focusing on retroviral integration process, there is still no clear understanding of integration specificity in a target site. Most LTR retrotransposons non-specifically integrate into a target site. Site-specificity of integration at vertebrate retroviruses is rather relative. At the same time, sequence-specific integration is the exclusive property of errantiviruses and their derivatives with two open reading frames. The possible basis for the errantivirus integration specificity is discussed in the present review.

New insights into the interplay between codon bias determinants in plants

Codon bias is the non-random use of synonymous codons, a phenomenon that has been observed in species as diverse as bacteria, plants and mammals. The preferential use of particular synonymous codons may reflect neutral mechanisms (e.g. mutational bias, G|C-biased gene conversion, genetic drift) and/or selection for mRNA stability, translational efficiency and accuracy. The extent to which these different factors influence codon usage is unknown, so we dissected the contribution of mutational bias and selection towards codon bias in genes from 15 eudicots, 4 monocots and 2 mosses. We analysed the frequency of mononucleotides, dinucleotides and trinucleotides and investigated whether the compositional genomic background could account for the observed codon usage profiles. Neutral forces such as mutational pressure and G|C-biased gene conversion appeared to underlie most of the observed codon bias, although there was also evidence for the selection of optimal translational efficiency and mRNA folding. Our data confirmed the compositional differences between monocots and dicots, with the former featuring in general a lower background compositional bias but a higher overall codon bias.

Transposable elements are a significant contributor to tandem repeats in the human genome

Sequence repeats are an important phenomenon in the human genome, playing important roles in genomic alteration often with phenotypic consequences. The two major types of repeat elements in the human genome are tandem repeats (TRs) including microsatellites, minisatellites, and satellites and transposable elements (TEs). So far, very little has been known about the relationship between these two types of repeats. In this study, we identified TRs that are derived from TEs either based on sequence similarity or overlapping genomic positions. We then analyzed the distribution of these TRs among TE families/subfamilies. Our study shows that at least 7,276 TRs or 23% of all minisatellites/satellites is derived from TEs, contributing ∼0.32% of the human genome. TRs seem to be generated more likely from younger/more active TEs, and once initiated they are expanded with time via local duplication of the repeat units. The currently postulated mechanisms for origin of TRs can explain only 6% of all TE-derived TRs, indicating the presence of one or more yet to be identified mechanisms for the initiation of such repeats. Our result suggests that TEs are contributing to genome expansion and alteration not only by transposition but also by generating tandem repeats.

Integration specificity of LTR-retrotransposons and retroviruses in the Drosophila melanogaster genome

Integration of DNA copies in a host genome is a necessary stage in the life cycle of retroviruses and LTR-retrotransposons. There is still no clear understanding of integration specificity of retroelements into a target site. The selection of the target DNA is believed to potentially affect a number of factors such as transcriptional status, association with histones and other DNA-binding proteins, and DNA bending. The authors performed a comprehensive computer analysis of the integration specificity of Drosophila melanogaster LTR-retrotransposons and retroviruses including an analysis of the nucleotide composition of targets, terminal sequences of LTRs, and integrase sequences. A classification of LTR-retrotransposons based on the integration specificity was developed. All the LTR-retrotransposons of the gypsy group with three open frames (errantiviruses) and their derivatives with two open frames demonstrate strict specificity to a target DNA selection. Such specificity correlates with the structural features of the target DNA: bendability, A-philicity, or protein-induced deformability. The remaining LTR-retrotransposons (copia and BEL groups, blastopia and 412 subgroups of the gypsy group) do not show specificity of integration. Chromodomain is present in the integrase structures of blastopia and 412 subgroup LTR-retrotransposons and may facilitate the process of non-specific integration.

Evolution of LEP150 sub-repeat array within the ribosomal IGS of the clam shrimp Leptestheria dahalacensis (Crustacea Branchiopoda Conchostraca)

Leptestheria dahalacensis genome harbours repeats of the LEP150 satellite DNA family linked to 5S gene, within the ribosomal intergenic spacer. In genetically isolated samples, the sequence analysis of the region (5S, flanking region, first satellite monomer: unit A, second satellite monomer: unit B) evidenced three 5S variants. The alpha and gamma variants share a greater homology. They co-occur in the Central European samples, while in the Italian one, the highly divergent alpha and beta variants are present. In phylogenetic analyses, A and B LEP150 monomers show a peculiar clustering; this was further confirmed through the sequencing for the alpha variant of four monomers at the 5' and 3' tails (units A, B, C, D and D', C', B', A', respectively). Horizontal homogenisation was observed only across C, D, C' and D' units. Furthermore, repeat sequence diversity decrease toward terminal repeats, at variance of literature data. The pattern of variation observed is explained taking into account the presence at the LEP150 array borders of two loci under natural selection: the 5S rRNA gene, upstream, and the rDNA transcription promoter, downstream. These elements may drive the dynamics of flanking regions and linked repeats in a process similar to selective sweep. At variance of classical genetic hitchhiking, the selective sweep here scored should be realized and maintained through an interplay of selection and molecular drive.

Consensus higher order repeats and frequency of string distributions in human genome

Key string algorithm (KSA) could be viewed as robust computational generalization of restriction enzyme method. KSA enables robust and effective identification and structural analyzes of any given genomic sequences, like in the case of NCBI assembly for human genome. We have developed a method, using total frequency distribution of all r-bp key strings in dependence on the fragment length l, to determine the exact size of all repeats within the given genomic sequence, both of monomeric and HOR type. Subsequently, for particular fragment lengths equal to each of these repeat sizes we compute the partial frequency distribution of r-bp key strings; the key string with highest frequency is a dominant key string, optimal for segmentation of a given genomic sequence into repeat units. We illustrate how a wide class of 3-bp key strings leads to a key-string-dependent periodic cell which enables a simple identification and consensus length determinations of HORs, or any other highly convergent repeat of monomeric or HOR type, both tandem or dispersed. We illustrated KSA application for HORs in human genome and determined consensus HORs in the Build 35.1 assembly. In the next step we compute suprachromosomal family classification and CENP-B box / pJalpha distributions for HORs. In the case of less convergent repeats, like for example monomeric alpha satellite (20-40% divergence), we searched for optimal compact key string using frequency method and developed a concept of composite key string (GAAAC--CTTTG) or flexible relaxation (28 bp key string) which provides both monomeric alpha satellites as well as alpha monomer segmentation of internal HOR structure. This method is convenient also for study of R-strand (direct) / S-strand (reverse complement) alpha monomer alternations. Using KSA we identified 16 alternating regions of R-strand and S-strand monomers in one contig in choromosome 7. Use of CENP-B box and/or pJalpha motif as key string is suitable both for identification of HORs and monomeric pattern as well as for studies of CENP-B box / pJalpha distribution. As an example of application of KSA to sequences outside of HOR regions we present our finding of a tandem with highly convergent 3434-bp Long monomer in chromosome 5 (divergence less then 0.3%).

The intragenomic polymorphism of a partially inverted repeat (PIR) in Gallus gallus domesticus, potential role of inverted repeats in satellite DNAs evolution

We report here the molecular characterization of the basic repeating unit of a novel repetitive family, partially inverted repeat (PIR), previously identified from chicken genome. This repetitive DNA family shares a close evolutionary relationship with XhoI/EcoRI repeats and chicken nuclear-membrane-associated (CNM) repeat. Sequence analyses reveal the 1430 bp basic repeating unit can be divided into two regions: the central region ( approximately 1000 bp) and the flanking region ( approximately 430 bp). Within the central region, a pair of repeats (86 bp) flanks the central core ( approximately 828 bp) in inversed orientation. Due to the tandem array feature shared by the repeating units, the inverted repeats fall between the central core and flanking region. Southern blot analyses further reveal the intragenomic polymorphism of PIR, and the molecular size of repeating units ranges from 1.1 kb to 1.6 kb. The identified monomer variants may result from multiple crossing-over events, implying the potential roles of inverted repeats in satellite DNAs variation.

Molecular cytogenetic characterization of chromosome 9-derived material in a human thyroid cancer cell line

The incidence of papillary thyroid carcinoma (PTC) increases significantly after exposure of the head and neck region to ionizing radiation, yet we know neither the steps involved in malignant transformation of thyroid epithelium nor the specific carcinogenic mode of action of radiation. Such increased tumor frequency became most evident in children after the 1986 nuclear accident in Chernobyl, Ukraine. In the eight years following the accident, the average incidence of childhood PTCs (chPTC) increased 70-fold in Belarus, 200-fold in Gomel, 10-fold in the Ukraine and 50-fold in Tschnigov, Kiev, Rovno, Shitomyr and Tscherkassy compared to the rate of about 1 tumor incidence per 106 children per year prior to 1986 (Likhtarev et al., 1995; Sobolev et al., 1997; Jacob et al., 1998). To study the etiology of radiation-induced thyroid cancer, we formed an international consortium to investigate chromosomal changes and altered gene expression in cases of post-Chernobyl chPTC. Our approach is based on karyotyping of primary cultures established from chPTC specimens, establishment of cell lines and studies of genotype-phenotype relationships through high resolution chromosome analysis, DNA/cDNA micro-array studies, and mouse xenografts that test for tumorigenicity. Here, we report the application of fluorescence in situ hybridization (FISH)-based techniques for the molecular cytogenetic characterization of a highly tumorigenic chPTC cell line, S48TK, and its subclones. Using chromosome 9 rearrangements as an example, we describe a new approach termed 'BAC-FISH' to rapidly delineate chromosomal breakpoints, an important step towards a better understanding of the formation of translocations and their functional consequences.

Origin and diversification of minisatellites derived from human Alu sequences

We analyze minisatellites derived from Alu fragments corresponding approximately to the first 44 bases of human Alu consensus sequences from different subfamilies. The origin of Alu-derived minisatellites appears to have been mediated by short flanking repeats, as first proposed by Haber and Louis [Haber, J.E., Louis, E.J., 1998. Minisatellite origins in yeast and humans. Genomics 48, 132-135.]. We also present evidence for base substitutions and deletions introduced to minisatellites by gene conversion with partially similar but unrelated flanking regions. Segments flanked by short direct repeats are relatively common in different regions of Alu and other repetitive sequences. Our analysis shows that they can be effectively used in comparative studies of the overall sequence context which may contribute to instability of DNA segments flanked by short direct repeats.

Genomic dynamics of a low-copy-number satellite DNA family in Leptestheria dahalacensis (Crustacea, Branchiopoda, Conchostraca)

The LEP150 satellite DNA (satDNA) family found in Leptestheria dahalacensis (Ruppel, 1837) (Conchostraca) is a low-copy-number satellite with a canonical monomer of 150 bp. Nucleotide variation analyses suggest a 14-bp palindromic region as a possible protein binding site with constraints acting on the whole sequence but a 25-bp variable box. Besides the head-to-tail arrangement of 150 bp monomers, multimers analyses evidenced incomplete monomers, one duplication event, and three inversions. Both observed rearrangements and the higher values of sequence variability scored suggest that rearranged monomers reside in regions with a lower degree of homogenisation efficiency. Sixty-seven percent of the breakpoints occurs at kinkable dinucleotides, thus supporting their role in rearrangements as documented in alphoid satDNA recombination events. Monomers of different lengths may result from crossing over between repeats misaligned through the direct and inverted subrepeats of LEP150 monomers. ANOVA results indicate that the same range of sequence diversity is experienced at the individual and population ranks; therefore, the evolution of the L. dahalacensis satDNA is concerted.

Formation of human hepatocytes by human hematopoietic stem cells in sheep

We took advantage of the proliferative and permissive environment of the developing preimmune fetus to develop a noninjury large animal model in sheep, in which the transplantation of defined populations of human hematopoietic stem cells resulted in the establishment of human hematopoiesis and led to the formation of significant numbers of long-lasting, functional human liver cells, with some animals exhibiting levels as high as 20% of donor (human) hepatocytes 11 months after transplantation. A direct correlation was found between hepatocyte activity and phenotype of transplanted cells, cell dose administered, source of cells used on a cell-per-cell basis (bone marrow, cord blood, mobilized peripheral blood), and time after transplantation. Human hepatocytes generated in this model retained functional properties of normal hepatocytes, constituted hepatic functional units with the presence of human endothelial and biliary duct cells, and secreted human albumin that was detected in circulation. Transplanting populations of hematopoietic stem cells can efficiently generate significant numbers of functional hepatic cells in this noninjury large animal model and thus could be a means of ameliorating or curing genetic diseases in which a deficiency of liver cells or their products threatens the life of the fetus or newborn.

The ribosomal RNA gene promoter and adjacent cis-acting DNA sequences govern plasmid DNA partitioning and stable inheritance in the parasitic protozoan Leishmania

Detailed analysis of the Leishmania donovani ribosomal RNA (rRNA) gene promoter region has allowed the identification of cis-acting sequences involved in plasmid DNA partitioning and stable plasmid inheritance. We report that plasmids bearing the 350 bp rRNA promoter along with the 200 bp region immediately 3' to the promoter exhibited a 6.5-fold increase in transformation frequency and were transmitted to daughter cells as single-copy molecules. This is in contrast to what has been observed for plasmid molecules in this organism so far. Moreover, we show that these low-copy-number plasmids displayed a remarkable mitotic stability in the absence of selective pressure. The region in the vicinity of the RNA pol I transcription initiation site, and also in the adjacent 200 nt, displays a complex structural organization and shares sequence similarity to the yeast autonomously replicating consensus sequence and centromere DNA elements. Deletion analyses indicated that these elements were necessary but not sufficient for plasmid DNA partitioning and stable inheritance, and that the rRNA promoter region was required for optimal function. These results suggest an interplay between RNA pol I transcription, DNA replication, DNA partitioning and mitotic stability in trypanosomatids. This is the first example of defined DNA elements for plasmid partitioning and stable inheritance in the protozoan parasite Leishmania.

Shuffling of genes within low-copy repeats on 22q11 (LCR22) by Alu-mediated recombination events during evolution

Low-copy repeats, or segmental duplications, are highly dynamic regions in the genome. The low-copy repeats on chromosome 22q11.2 (LCR22) are a complex mosaic of genes and pseudogenes formed by duplication processes; they mediate chromosome rearrangements associated with velo-cardio-facial syndrome/DiGeorge syndrome, der(22) syndrome, and cat-eye syndrome. The ability to trace the substrates and products of recombination events provides a unique opportunity to identify the mechanisms responsible for shaping LCR22s. We examined the genomic sequence of known LCR22 genes and their duplicated derivatives. We found Alu (SINE) elements at the breakpoints in the substrates and at the junctions in the truncated products of recombination for USP18, GGT, and GGTLA, consistent with Alu-mediated unequal crossing-over events. In addition, we were able to trace a likely interchromosomal Alu-mediated fusion between IGSF3 on 1p13.1 and GGT on 22q11.2. Breakpoints occurred inside Alu elements as well as in the 5' or 3' ends of them. A possible stimulus for the 5' or 3' terminal rearrangements may be the high sequence similarities between different Alu elements, combined with a potential recombinogenic role of retrotransposon target-site duplications flanking the Alu element, containing potentially kinkable DNA sites. Such sites may represent focal points for recombination. Thus, genome shuffling by Alu-mediated rearrangements has contributed to genome architecture during primate evolution.

Interspersed repeats are found predominantly in the “old” α satellite families

The biased distribution of dispersed repeat insertions in various types of primate specific alpha satellites (AS) is being discussed in the literature in relation to the modes of AS evolution and their possible roles in maintenance and disruption of functional centromeres. However, such a bias has not been properly documented on a genome-wide scale so far. In this work, using a representative sample of about 100 insertions we show that the "old" AS contains at least 10 times more dispersed repeats than the "new" one. In the new arrays insertions accumulate mostly in poorly homogenized areas, presumably in the edges, and in the old AS, throughout the whole array length. Dating of L1 insertions in the old AS revealed that their massive accumulation started at or after the time when the new AS emerged and expanded in the genome and the centromere function had shifted to the new AS arrays.