Sequences flanking the repeat arrays of human minisatellites: association with tandem and dispersed repeat elements

Advances in the discovery and analyses of human tandem repeats

Long-read sequencing platforms provide unparalleled access to the structure and composition of all classes of tandemly repeated DNA from STRs to satellite arrays. This review summarizes our current understanding of their organization within the human genome, their importance with respect to disease, as well as the advances and challenges in understanding their genetic diversity and functional effects. Novel computational methods are being developed to visualize and associate these complex patterns of human variation with disease, expression, and epigenetic differences. We predict accurate characterization of this repeat-rich form of human variation will become increasingly relevant to both basic and clinical human genetics.

Transposable Elements as a Source of Novel Repetitive DNA in the Eukaryote Genome

The impact of transposable elements (TEs) on the evolution of the eukaryote genome has been observed in a number of biological processes, such as the recruitment of the host’s gene expression network or the rearrangement of genome structure. However, TEs may also provide a substrate for the emergence of novel repetitive elements, which contribute to the generation of new genomic components during the course of the evolutionary process. In this review, we examine published descriptions of TEs that give rise to tandem sequences in an attempt to comprehend the relationship between TEs and the emergence of de novo satellite DNA families in eukaryotic organisms. We evaluated the intragenomic behavior of the TEs, the role of their molecular structure, and the chromosomal distribution of the paralogous copies that generate arrays of repeats as a substrate for the emergence of new repetitive elements in the genome. We highlight the involvement and importance of TEs in the eukaryote genome and its remodeling processes.

Conversion of DNA Sequences: From a Transposable Element to a Tandem Repeat or to a Gene

Eukaryotic genomes are rich in repetitive DNA sequences grouped in two classes regarding their genomic organization: tandem repeats and dispersed repeats. In tandem repeats, copies of a short DNA sequence are positioned one after another within the genome, while in dispersed repeats, these copies are randomly distributed. In this review we provide evidence that both tandem and dispersed repeats can have a similar organization, which leads us to suggest an update to their classification based on the sequence features, concretely regarding the presence or absence of retrotransposons/transposon specific domains. In addition, we analyze several studies that show that a repetitive element can be remodeled into repetitive non-coding or coding sequences, suggesting (1) an evolutionary relationship among DNA sequences, and (2) that the evolution of the genomes involved frequent repetitive sequence reshuffling, a process that we have designated as a "DNA remodeling mechanism". The alternative classification of the repetitive DNA sequences here proposed will provide a novel theoretical framework that recognizes the importance of DNA remodeling for the evolution and plasticity of eukaryotic genomes.

Analysis of a c0t-1 library enables the targeted identification of minisatellite and satellite families in Beta vulgaris

BACKGROUND

Repetitive DNA is a major fraction of eukaryotic genomes and occurs particularly often in plants. Currently, the sequencing of the sugar beet (Beta vulgaris) genome is under way and knowledge of repetitive DNA sequences is critical for the genome annotation. We generated a c0t-1 library, representing highly to moderately repetitive sequences, for the characterization of the major B. vulgaris repeat families. While highly abundant satellites are well-described, minisatellites are only poorly investigated in plants. Therefore, we focused on the identification and characterization of these tandemly repeated sequences.

RESULTS

Analysis of 1763 c0t-1 DNA fragments, providing 442 kb sequence data, shows that the satellites pBV and pEV are the most abundant repeat families in the B. vulgaris genome while other previously described repeats show lower copy numbers. We isolated 517 novel repetitive sequences and used this fraction for the identification of minisatellite and novel satellite families. Bioinformatic analysis and Southern hybridization revealed that minisatellites are moderately to highly amplified in B. vulgaris. FISH showed a dispersed localization along most chromosomes clustering in arrays of variable size and number with exclusion and depletion in distinct regions.

CONCLUSION

The c0t-1 library represents major repeat families of the B. vulgaris genome, and analysis of the c0t-1 DNA was proven to be an efficient method for identification of minisatellites. We established, so far, the broadest analysis of minisatellites in plants and observed their chromosomal localization providing a background for the annotation of the sugar beet genome and for the understanding of the evolution of minisatellites in plant genomes.

Evolutionary conservation of orthoretroviral long terminal repeats (LTRs) and ab initio detection of single LTRs in genomic data

Background

Retroviral LTRs, paired or single, influence the transcription of both retroviral and non-retroviral genomic sequences. Vertebrate genomes contain many thousand endogenous retroviruses (ERVs) and their LTRs. Single LTRs are difficult to detect from genomic sequences without recourse to repetitiveness or presence in a proviral structure. Understanding of LTR structure increases understanding of LTR function, and of functional genomics. Here we develop models of orthoretroviral LTRs useful for detection in genomes and for structural analysis.

Principal Findings

Although mutated, ERV LTRs are more numerous and diverse than exogenous retroviral (XRV) LTRs. Hidden Markov models (HMMs), and alignments based on them, were created for HML- (human MMTV-like), general-beta-, gamma- and lentiretroviruslike LTRs, plus a general-vertebrate LTR model. Training sets were XRV LTRs and RepBase LTR consensuses. The HML HMM was most sensitive and detected 87% of the HML LTRs in human chromosome 19 at 96% specificity. By combining all HMMs with a low cutoff, for screening, 71% of all LTRs found by RepeatMasker in chromosome 19 were found. HMM consensus sequences had a conserved modular LTR structure. Target site duplications (TG-CA), TATA (occasionally absent), an AATAAA box and a T-rich region were prominent features. Most of the conservation was located in, or adjacent to, R and U5, with evidence for stem loops. Several of the long HML LTRs contained long ORFs inserted after the second A rich module. HMM consensus alignment allowed comparison of functional features like transcriptional start sites (sense and antisense) between XRVs and ERVs.

Conclusion

The modular conserved and redundant orthoretroviral LTR structure with three A-rich regions is reminiscent of structurally relaxed Giardia promoters. The five HMMs provided a novel broad range, repeat-independent, ab initio LTR detection, with prospects for greater generalisation, and insight into LTR structure, which may aid development of LTR-targeted pharmaceuticals.

Population structure within and between subspecies of the Mediterranean triplefin fish Tripterygion delaisi revealed by highly polymorphic microsatellite loci

Although F(ST) values are widely used to elucidate population relationships, in some cases, when employing highly polymorphic loci, they should be regarded with caution, particularly when subspecies are under consideration. Tripterygion delaisi presents two subspecies that were investigated here, using 10 microsatellite loci. A Bayesian approach allowed us to clearly identify both subspecies as two different evolutionary significant units. However, low F(ST) values were found between subspecies as a consequence of the large number of alleles per locus, while homoplasy could be disregarded as indicated by the standardized genetic distance G'(ST). Heterozygosity saturation was observed in highly polymorphic loci containing more than 15 alleles, and this threshold was used to define two loci pools. The less variable loci pool revealed higher genetic variance between subspecies, while the more variable pool showed higher genetic variance between populations. Furthermore, higher differentiation was also observed between populations using G'(ST) with the more variable loci. Nonetheless, a more reliable population structure within subspecies was obtained when all loci were included in the analyses. In T. d. xanthosoma, isolation by distance was detected between the eight analysed populations, and six genetically homogeneous clusters were inferred by Bayesian analyses that are in accordance with F(ST) values. The neighbourhood-size method also indicated rather small dispersal capabilities. In conclusion, in fish with limited adult and larval dispersal capabilities, continuous rocky habitat seems to allow contact between populations and prevent genetic differentiation, while large discontinuities of sand or deep-water channels seems to reduce gene flow.

Crystal structure of d(gcGXGAgc) with X=G: a mutation at X is possible to occur in a base-intercalated duplex for multiplex formation

DNA fragments with the sequences d(gcGX[Y]n Agc) (n=1, X=A, and Y=A, T, or G)form base-intercalated duplexes, which is a basic unit for formation of multiplexes such as octaplex and hexaplex. To examine the stability of multiplexes, a DNA with X=Y=G and n=1 was crystallized under conditions different from those of the previously determined sequences, and its crystal structure has been determined. The two strands are coupled in an anti-parallel fashion to form a base-intercalated duplex, in which the first and second residues form Watson-Crick type G:C pairs and the third and sixth residues form a sheared G:A pairs at both ends of the duplex. The G4 and G5 bases are stacked alternatively on those of the counter strand to form a long G column of G3-G4-G5*-G5-G4*-G3*, the central four Gs being protruded. In addition, the three duplexes are associated to form a hexaplex around a mixture of calcium and sodium cations on the crystallographic threefold axis. These structural features are similar to those of the previous crystals, though slightly different in detail. The present study indicates that mutation at the 4th position is possible to occur in a base-intercalated duplex for multiplex formations, suggesting that DNA fragments with any sequence sandwiched between the two triplets gcG and Agc can form a multiplex.

Local mutagenic impact of insertions of LTR retrotransposons on the mouse genome

Solitary LTR loci are the predominant form of LTR retrotransposons in most eukaryotic genomes. They originate from recombination between the two LTRs of an ancestral retrovirus and are therefore incapable of transposition. Despite this inactivity, they appear to have a substantial impact on the host genome. Here we use the murine RMER10 LTR family as an example to describe how such elements can reshape regions of the genome through multiple mutations on an evolutionary time scale. Specifically, we use phylogenetic analysis of multiple copies of RMER10 in rodent species, as well as comparisons of orthologous pairs in mouse and rat, to argue that insertions of members of this family have locally induced the emergence of tandem repeat loci as well as many indels. Analysis of structural aspects of these sequences (secondary structures and transcription factors signals) may explain why RMER10 can become endogenous "mutagenic" factors through induction of replication fork blockages and/or error-prone repair of aberrant DNA structures. This hypothesis is also consistent with features of other interspersed repeated elements.

Crystal structures of a DNA octaplex with I-motif of G-quartets and its splitting into two quadruplexes suggest a folding mechanism of eight tandem repeats

Recent genomic analyses revealed many kinds of tandem repeats of specific sequences. Some of them are related to genetic diseases, but their biological functions and structures are still unknown. Two X-ray structures of a short DNA fragment d(gcGA[G]1Agc) show that four base-intercalated duplexes are assembled to form an octaplex at a low K+ concentration, in which the eight G5 residues form a stacked double G-quartet in the central part. At a higher K+ concentration, however, the octaplex is split into just two halves. These structural features suggest a folding process of eight tandem repeats of d(ccGA[G]4Agg), according to a double Greek-key motif. Such a packaging of the repeats could facilitate slippage of a certain sequence during DNA replication, to induce increase or decrease of the repeats.

Role of transposable elements in the propagation of minisatellites in the rice genome

A survey of minisatellites (MSs) in 5.3 Mb of randomly selected rice DNA sequences from public databases was carried out to clarify the role of transposable elements (TEs) in the dispersal of MSs in the rice genome. The estimated frequency of MSs in this sample was one per 23.4 kb, and this frequency is approximately equivalent to that of Class I microsatellites in the rice genome. Of the MSs in the 5.3-Mb sequence sample, 82% were found to be present in multiple copies in the rice genome, and all of these were a part of TE sequences. In this study at least 61 TE groups were identified as MS carriers. It was also shown that the GC-rich MS pOs6.2H, which was previously reported to be one of the interspersed MSs in the rice genome, is a component of an En/Spm-like element. These results indicate that the majority of MSs in the rice genome are maintained in TEs, and amplified and dispersed as components of the TEs. The G+C content of the multi-locus MS sequences reflected that of the TE sequences containing those MSs, but no obvious bias towards the high G+C content of DNA was observed. Single locus MSs also did not show any obvious bias towards the high G+C content of DNA in the rice genome. In this respect, the MSs in the rice genome are quite different from those in the human genome: in the latter, the majority of MSs show an obvious bias towards the high G+C content of DNA.

Extremely complex repeat shuffling during germline mutation at human minisatellite B6.7

Human minisatellite B6.7 is a highly variable locus showing extensive heterozygosity with alleles ranging from six to >500 repeat units. Paternal and maternal mutation rates to new length alleles were estimated from pedigrees at 7.0 and 3.9% per gamete, respectively, indicating that B6.7 is one of the most unstable minisatellites isolated to date. Mutation at this locus was also analysed by small pool PCR of sperm and blood DNA. Male germline instability varied from <0.8 to 14% per allele and increased with tandem array size. In contrast, the frequency of mutants in somatic (blood) DNA was far lower (<0.5%), consistent with a meiotic origin of germline mutants. Sperm mutants were further characterized by minisatellite variant repeat mapping using four major polymorphic sites within the B6.7 repeats. This highly informative system revealed a wide variety of changes in allele structure, including simple intra-allelic duplications and deletions and more complicated inter- and intra-allelic transfers of repeat blocks, as seen at other human minisatellites. The main mode of sperm mutation, however, resulted in extremely complex allele reorganization with evidence of inter-allelic transfer plus the generation of novel repeats by rearrangement at the sub-repeat level, suggesting that recombinational instability at B6.7 is a complex multistep process.

Comparative Sequence Analysis of Human Minisatellites Showing Meiotic Repeat Instability

The highly variable human minisatellites MS32 (D1S8), MS31A (D7S21), and CEB1 (D2S90) all show recombination-based repeat instability restricted to the germline. Mutation usually results in polar interallelic conversion or occasionally in crossovers, which, at MS32 at least, extend into DNA flanking the repeat array, defining a localized recombination hotspot and suggesting that cis-acting elements in flanking DNA can influence repeat instability. Therefore, comparative sequence analysis was performed to search for common flanking elements associated with these unstable loci. All three minisatellites are located in GC-rich DNA abundant in dispersed and tandem repetitive elements. There were no significant sequence similarities between different loci upstream of the unstable end of the repeat array. Only one of the three loci showed clear evidence for putative coding sequences near the minisatellite. No consistent patterns of thermal stability or DNA secondary structure were shared by DNA flanking these loci. This work extends previous data on the genomic environment of minisatellites. In addition, this work suggests that recombinational activity is not controlled by primary or secondary characteristics of the DNA sequence flanking the repeat array and is not obviously associated with gene promoters as seen in yeast.

[The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF048727(CEB1), AF048728 (MS31A), and AF048729 (MS32).]

Comparative sequence analysis of human minisatellites showing meiotic repeat instability

The highly variable human minisatellites MS32 (D1S8), MS31A (D7S21), and CEB1 (D2S90) all show recombination-based repeat instability restricted to the germline. Mutation usually results in polar interallelic conversion or occasionally in crossovers, which, at MS32 at least, extend into DNA flanking the repeat array, defining a localized recombination hotspot and suggesting that cis-acting elements in flanking DNA can influence repeat instability. Therefore, comparative sequence analysis was performed to search for common flanking elements associated with these unstable loci. All three minisatellites are located in GC-rich DNA abundant in dispersed and tandem repetitive elements. There were no significant sequence similarities between different loci upstream of the unstable end of the repeat array. Only one of the three loci showed clear evidence for putative coding sequences near the minisatellite. No consistent patterns of thermal stability or DNA secondary structure were shared by DNA flanking these loci. This work extends previous data on the genomic environment of minisatellites. In addition, this work suggests that recombinational activity is not controlled by primary or secondary characteristics of the DNA sequence flanking the repeat array and is not obviously associated with gene promoters as seen in yeast.

Distribution of the mammalian-wide interspersed repeats (MIRs) in the isochores of the human genome

The distribution of MIRs (mammalian-wide interspersed repeats) was investigated in 164 human sequences (> or = 100 kb), which were assigned, according to their GC level, to isochore families L, H1, H2 and H3. MIR elements, whose total number in the genome was estimated to be about 3.3 x 10(5), were found to be unevenly distributed in human isochores. The majority of MIRs (55%) were found in the L isochore family. In contrast, MIR density was highest in H2, closely followed by H1, whereas densities in L and H3 were 2- and 3-fold lower than in H2, respectively. For this reason, the assessment of MIR distribution by inter-repeat PCR led to an overestimation of MIR numbers in H2 isochore and an underestimation in L isochores.

Linkage of two distinct AT-rich minisatellites at multiple loci in the genome of Theileria parva

Minisatellite tandem repeat elements are well known components of vertebrate genomes, but have not yet been extensively characterized in lower eukaryotes. We describe two unusual, AT-rich minisatellites of the protozoan parasite Theileria parva whose sequences are unrelated to the G/C-rich i minisatellite superfamily' of vertebrate and plant genomes. The T. parva tandem repeats, one with a conserved sequence T2-5ACACA (6-17 copies), and the other with a 6-bp core sequence of either ACTATA or TATACT associated with additional variable sequences in repeats of 10-17bp (3-7 copies), were closely linked at more than 20 sites in the T. parva genome, separated by 390, 510 and 660bp at three loci analysed in detail. Such linkage is without precedent in minisatellites so far analysed in other organisms. The minisatellite loci were widely dispersed on 13 out of 33 genomic SfiI fragments, on all four T. parva chromosomes and did not exhibit a telomeric bias in their distribution. Analysis of flanking sequences revealed no obvious conserved sequences between the five loci, or other multicopy repeat sequences outside the minisatellite regions. The T2-5 ACACA minisatellite was highly effective as a multilocus fingerprinting probe for discrimination of T. parva isolates. Analysis of two individual minisatellite loci revealed variation between the genomic DNAs of two T. parva isolates in the copy number of the constituent repeats within the array, similar to that typical of vertebrate minisatellites. 1998 Elsevier Science B.V.

High-resolution mapping of crossovers in human sperm defines a minisatellite-associated recombination hotspot

Little is known about the fine-scale distribution of meiotic crossovers in human chromosomes. Methods have therefore been developed for detecting and mapping recombination products directly in human sperm DNA. Analysis of crossovers adjacent to the GC-rich minisatellite MS32, which is known to mutate by conversion and crossover within the repeat array, revealed an intense and highly localized recombination hotspot centered upstream of the locus and extending into the beginning of the minisatellite. Allele-specific cosuppression of crossovers and repeat instability suggests that the hotspot is responsible for driving repeat turnover at MS32 and thus that minisatellites might evolve as by-products of localized meiotic recombination in the human genome.

A novel unstable mouse VNTR family expanded from SINE B1 elements

Human hypervariable minisatellites provide highly informative loci for analyzing processes of tandem repeat turnover, but little is known about mechanisms of instability in other species. We have therefore screened the mouse genome for analogous VNTR loci. One of the probes we isolated, MMS10, detects a highly variable rodent-specific family of multiple loci derived by expansion of a common GGCAGA repeat unit from within a subset of B1 short interspersed elements. The mean germline mutation rate for loci detected by MMS10 was estimated at 1.7% per offspring band, though analysis of individual loci showed substantial variation in germinal instability apparently related to repeat array size. This MMS10 family of expanded hexanucleotide repeat loci provides a novel resource for investigating mechanisms of tandem repeat turnover in the mouse and an efficient means for monitoring germline mutations induced by external agents such as ionizing radiation.

Minisatellite origins in yeast and humans

Variable numbers of tandem repeats are valuable markers in genetic studies. The arrays of interest are simple microsatellites, containing repetitions of 1-5 nucleotides, and minisatellites, with multiple iterations of approximately 10 to 100 bp. Microsatellite origins can be explained by replication errors in regions fortuitously containing two or more adjacent short repeats. Microsatellite variation arises by replication errors in the absence of mismatch correction (R. Parsons et al., 1993, Cell 75: 1227-1236; M. Strand et al., 1993, Nature 365: 274-276). Variation in the size of minisatellites is thought to involve homologous recombination processes, including gene conversions (J. Buard and G. Vergnaud, 1994, EMBO J. 13: 3203-3210; A. J. Jeffreys et al., 1994, Nature Genet. 6: 136-145) and possibly unequal exchanges among repeats. However, the origins of minisatellites are less obvious. The probability of finding a direct tandem repeat of minisatellite size by chance alone is very low [<4-(10 to 100)]. Here we report the finding of short direct repeats of 5 to 10 bp flanking many yeast and human minisatellites that may be involved in their origins through replication slippage or unequal crossings over.

Linkage analysis of bronchial hyperreactivity and atopy with chromosome 11q13

There are two key clinical features of asthma: allergy and bronchial hyperreactivity (BHR). Some pedigree studies of atopy have indicated linkage with the high affinity IgE receptor (Fc epsilon RI-beta) gene on chromosome 11q13, but others failed to confirm this linkage. We examined the genetic linkage of three polymorphic microsatellite markers to atopy and BHR in 120 affected sibling pairs recruited from the general community. We found no linkage to atopy at any of the three 11q13 loci studied. Our findings also do not favour linkage between BHR and loci approximately 8-9 cM either side of the Fc epsilon RI-beta gene.

The Australian Electrophoresis Society annual conference at Chromatography '96: Separation Sciences. Rosehill, NSW, Australia, July 9-11, 1996

A set of eleven biallelic and three multiallelic molecular markers have been developed to analyze populations of Histoplasma capsulatum. All markers are amplified by polymerase chain reaction (PCR) and can be readily scored using minimal amounts of template DNA. The 11 biallelic loci have polymorphic restriction endonuclease sites or small insertions or deletions which may be assessed by agarose gel electrophoresis. These markers are inherited in an unambiguous manner and are ideal for assessing structure and gene flow within US populations of H. capsulatum, but are monomorphic in non-US populations. Both length and sequence variation are present in the multiallelic loci, which can be scored by direct sequencing, polyacrylamide gel electrophoresis, or single-strand conformation polymorphism (SSCP): As they are hypervariable, the multiallelic loci can be used to type isolates and to assess the level of genetic variation within populations. Preliminary results indicate that the three multiallelic markers presented are sufficient to distinguish isolates at the individual level and are polymorphic in both US and non-US populations. This collection of molecular markers will be a useful tool in population and epidemiology studies of H. capsulatum.

DNA secondary structures and the evolution of hypervariable tandem arrays

Tandem repeats are ubiquitous in nature and constitute a major source of genetic variability in populations. This variability is associated with a number of genetic disorders in humans including triplet expansion diseases such as Fragile X syndrome and Huntington's disease. The mechanism responsible for the variability/instability of these tandem arrays remains contentious. We show here that formation of secondary structures, in particular intrastrand tetraplexes, is an intrinsic property of some of the more unstable arrays. Tetraplexes block DNA polymerase progression and may promote instability of tandem arrays by increasing the likelihood of reiterative strand slippage. In the course of doing this work we have shown that some of these tetraplexes involve unusual base interactions. These interactions not only generate tetraplexes with novel properties but also lead us to conclude that the number of sequences that can form stable tetraplexes might be much larger than previously thought.

Framework YAC contig anchored into a 3.2-Mb high-resolution physical map in proximal 11q13

Despite the presence on band q13 of chromosome 11 of a number of genes predisposing individuals to various human diseases, most of this genomic region remains loosely mapped. Moreover, there is a relative dearth of yeast artificial chromosome (YAC) contigs from genome-wide studies: YACs are irregularly distributed over this chromosomal region and have not been arranged into contigs. We have thus undertaken fine-scale mapping of a 3.2-Mb region flanked by ACTN3 and FGF3. Since this region has demonstrated a high degree of YAC instability, we have established a framework contig by anchoring YACs and cosmids into a high-resolution physical map based on fluorescence in situ hybridization and long-range restriction mapping. The 3.2-Mb area studied includes the boundaries of regions thought to contain genes predisposing individuals to osteoporosis-pseudoglioma syndrome and insulin-dependent diabetes mellitus, as well as genes driving amplification events in human carcinomas. Another feature of this genomic area is that it cross-hybridizes to nonsyntenic regions of the genome. In addition, it spans the region where syntenic conservation with mouse chromosome 19 ends, making clones that we have anchored there valuable tools in understanding genome evolution.

Short alleles revealed by PCR demonstrate no heterozygote deficiency at minisatellite loci D1S7, D7S21, and D12S11

Short VNTR alleles that go undetected after conventional Southern blot hybridization may constitute an alternative explanation for the heterozygosity deficiency observed at some minisatellite loci. To examine this hypothesis, we have employed a screening procedure based on PCR amplification of those individuals classified as homozygotes in our databases for the loci D1S7, D7S21, and D12S11. The results obtained indicate that the frequency of these short alleles is related to the heterozygosity deficiency observed. For the most polymorphic locus, D1S7, approximately 60% of those individuals previously classified as homozygotes were in fact heterozygotes for a short allele. After the inclusion of these new alleles, the agreement between observed and expected heterozygosity, along with other statistical tests employed, provide additional evidence for lack of population substructuring. Comparisons of allele frequency distributions reveal greater differences between racial groups than between closely related populations.

Isolation of telomere junction fragments by anchored polymerase chain reaction

We describe a simple polymerase chain reaction (PGR)-based method for isolating short stretches of nontelomeric DNA adjacent to arrays of telomere repeat units, in principle applicable to any species for which the telomere repeat sequence is known. Application of this approach to human DNA resulted in the isolation of many candidate telomere junction clones, at least some of which were shown to be derived from telomere-adjacent regions. Most of the isolated clones detect multiple sequences in the human genome which represent one or a few sequence families present at the ends of most or all autosomes and variably truncated before the start of the telomere repeat array. Substantial sequence divergence between different members of these sequence families suggests a low rate of sequence homogenization by telomere exchange processes. The pseudoautosomal telomere junction has also been isolated and contains a shortened version of a recently described family of short interspersed repetitive elements (SINEs), only 14 base pairs (b.p.) from the start of the telomere.

A recombination hot spot in the Rh genes revealed by analysis of unrelated donors with the rare D-- phenotype

We have studied the arrangement of Rh (rhesus) genes in donors who are completely null for the products of one of them, RHCE. We show that five of six homozygous individuals with the so-called Rh D-- phenotype, who express no red-cell antigens of the C/c and E/e series, have rearranged RHCE genes in which internal sequences have been replaced by the corresponding sequences from RHD. Moreover, although there is heterogeneity at the 3' end, the 5' boundary of this chimerism is within the same small interval around exon 2. This interval is characterized by an exceptionally high degree of sequence homology between RHCE and RHD, a high density of dispersed repetitive elements, and the presence of an alternating purine-pyrimidine copolymer tract. We suggest that these features may explain the mechanistic basis for the origin of the rearrangement.

A panel of VNTR markers in pigs

By cloning tandemly repeated sequences from the pig genome by use of non-porcine minisatellite probes for library screening, five novel polymorphic VNTR loci were isolated: three minisatellites and two satellite-like loci. Four of them could be mapped onto chromosomes by linkage analysis and/or in situ hybridization. They were assigned to Chromosomes (Chrs) 5, 6, 14, and 16. Physical mapping on both presumed satellites and on one of the minisatellites revealed that the former resided near or at the centromere and the latter towards the chromosome ends. The location of the minisatellite is of particular interest since, together with data on three other minisatellites previously isolated, it supports the idea that, as in humans, minisatellites may preferentially be subtelomeric also in pigs.

The specificity of association of the IgD molecule with the accessory proteins BAP31/BAP29 lies in the IgD transmembrane sequence

Mature B cells co-express on their cell surface two classes of antigen receptor, the IgM and IgD immunoglobulins. The structural and functional differences between the two receptor classes are poorly understood. Recently two proteins of 29 and 31 kDa (BAP29 and BAP31) have been described that are preferentially associated with membrane IgD but only weakly with membrane IgM. We describe here the cloning of full-length murine and human BAP31 cDNAs encoding proteins of 245 and 246 amino acids respectively. The two BAP31 proteins are 95% identical. The BAP31 gene is ubiquitously expressed in murine tissues and is located on the X chromosome in both mouse and man. The murine BAP31 protein has 43% sequence identity to murine BAP29. Both proteins have a hydrophobic N-terminus and an alpha-helical C-terminus which ends with a KKXX motif implicated in vesicular transport. By a mutational analysis we have identified amino acids in the transmembrane sequence of the delta m chain that are critical for binding to BAP31/BAP29. A structural model of the BAPs and their potential functions are discussed.

Does monozygotic twinning occur in mice?

Published reports suggest that the incidence of monozygotic twinning in women is increased after hormonally induced ovulation. Since some statistical evidence exists to indicate that monozygotic twinning may also occur in mice, we attempted to devise a mouse system in which the incidence of such twinning could be compared after spontaneous versus hormonally induced ovulation, in order to analyse the developmental basis of such an effect. We used phenotypic identity in litters segregating for ten genetic loci (not all independent) to indicate possible twin pairs. DNA finger-printing using three human minisatellite probes was then performed blind on these pairs and on sibling controls. From a total of over 2000 mice born, 40 apparently identical pairs were identified, on which DNA finger-printing was successfully carried out on 35 pairs. All proved to be derived from different zygotes. We conclude that monozygotic twin pairs are either extremely rare in the stock of mice that we studied, or have such reduced viability that their chance of surviving to weaning is low.

Use of a mammalian interspersed repetitive (MIR) element in the coding and processing sequences of mammalian genes

The mammalian interspersed repetitive (MIR) element was amplified in mammals 130 million years ago. The MIR element is at least 260 bp in length and is found in approximately 105 copies in the mammalian genome. We analyzed copies of the MIR element in the DNA of various mammals to determine its relationship to the structure and function of genes, in an attempt to identify specific uses of the MIR element within the mammalian genome. We found that alternative splicing within the acetylcholine receptor gene in humans takes place within the MIR element and results in the incorporation of part of the MIR element into the coding sequence of this gene. Furthermore, the polyadenylation signal (AATAAA) at the 3' end of four different mammalian genes is derived from the MIR element. These uses of the MIR element suggest that other regulatory sequences found within the mammalian genome originated from ancient transposable elements, many of which may no longer be recognizable.

Chromosome 11q13 markers and D-type cyclins in breast cancer

One in six primary human breast cancers has DNA amplification centered on the cyclin D1 gene (CCND1) on chromosome 11q13. This genetic abnormality is preferentially associated with estrogen-receptor positive tumors and may define a sub-class of patients with an adverse prognosis. Although CCND1 has the credentials of a cellular oncogene, being a target for chromosomal translocation and retroviral integration, the 11q13 amplicon encompasses several other markers and CCND1 is not the only candidate for the key gene on the amplified DNA. To assess their relative importance, we have constructed a physical map of the amplified DNA and compared the extent and frequency of amplification across the region. Since it is likely that the gene providing the selective force for amplification will be expressed at elevated levels, we have also examined expression of both RNA and protein. By these criteria, cyclin D1 remains the strongest candidate for the key oncogene on the amplicon and we are currently investigating the functional consequences of its over-expression.

The detection of subtelomeric chromosomal rearrangements in idiopathic mental retardation

A major challenge for human genetics is to identify new causes of mental retardation, which, although present in about 3% of individuals, is unexplained in more than half of all cases. We have developed a strategy to screen for the abnormal inheritance of subtelomeric DNA polymorphisms in individuals with mental retardation and have detected three abnormalities in 99 patients with normal routine karyotypes. Pulsed-field gel electrophoresis and reverse chromosome painting showed that one case arose from an interstitial or terminal deletion and two from the de novo inheritance of derivative translocation chromosomes. At least 6% of unexplained mental retardation is accounted for by these relatively small chromosomal abnormalities, which will be an important resource in the characterization of the genetic basis of neurodevelopment.

Ubiquitous mammalian-wide interspersed repeats (MIRs) are molecular fossils from the mesozoic era

Short interspersed elements (SINEs) are ubiquitous in mammalian genomes. Remarkable variety of these repeats among placental orders indicates that most of them amplified in each lineage independently, following mammalian radiation. Here, we present an ancient family of repeats, whose sequence divergence and common occurrence among placental mammals, marsupials and monotremes indicate their amplification during the Mesozoic era. They are called MIRs for abundant Mammalian-wide Interspersed Repeats. With approximately 120,000 copies still detectable in the human genome (0.2-0.3% DNA), MIRs represent a 'fossilized' record of a major genetic event preceding the radiation of placental orders.

Minisatellite mutation rate variation associated with a flanking DNA sequence polymorphism

Human minisatellite mutation in the male germline frequently involves complex interallelic gene conversion events restricted to one end of the tandem repeat array. Some alleles at minisatellite MS32 show reduced variability in human populations and are associated with a G to C transversion upstream of the array. Analysis of single sperm demonstrated a frequently profound reduction in mutation rate at alleles carrying the C variant. This mutation suppression acts in cis, but does not affect the ability of an allele to act as sequence donor during gene conversion. This mutation rate polymorphism provides strong evidence for elements near the minisatellite that regulate tandem repeat instability.

Identification, isolation and characterization of canine minisatellite sequences

A Charomid ordered-array library containing a 2-16 Kb size fraction of MboI-digested canine genomic DNA has been screened with the Jeffreys multilocus probes, 33.6 and 33.15, to identify and isolate canine minisatellite sequences. Of the 48 positive clones identified, 7 were found to contain polymorphic minisatellites with heterozygosities in the range 20-88%. The majority of the remainder were either monomorphic or dimorphic in the animals tested. Analysis of intrabreed variation in Bedlington Terriers using two polymorphic minisatellites has shown that a significant reduction occurs in the number of alleles seen compared to an agglomerated population sample, correlating with the high level of inbreeding within this breed. Flanking DNA sequence and partial repeat sequence is presented for the most polymorphic minisatellite thus far identified, cCfaMP5. The variable region in this minisatellite is similar to human minisatellites which show a distinct purine or pyrimidine strand bias.

Familial atopy in Australian pedigrees: adventitious linkage to chromosome 8 is not confirmed nor is there evidence of linkage to the high affinity IgE receptor

Atopy frequently displays autosomal dominant inheritance and recent studies have favoured genetic linkage between atopy and the human chromosome 11q13. We have studied 12 extended families with aggregation of atopy consistent with autosomal dominant inheritance. The families have been studied for linkage of asthma and atopy to loci on chromosome 8p following the observation that one family suggested preliminary evidence of linkage to an anonymous hypervariable locus cloned from a DNA fingerprint and mapped to 8pter-p22. Subsequent analysis shows this putative linkage to be adventitious as the remaining 11 families do not support linkage between atopy and 8p. We have analysed the same families for evidence of linkage of atopy to loci on 11q13. In these families there is no evidence of association between atopy and the 11q loci stronger than that expected by chance alone; furthermore there is no suggestion that a subpopulation of these families display linkage between atopy and the loci. In addition neither the 8p loci nor the 11q loci exhibit evidence of linkage to atopy by affected sib-pair analysis. This also conflicts with previously published data for 11q.

How is it that microsatellites and random oligonucleotides uncover DNA fingerprint patterns?

Minisatellites, microsatellites, and short random oligonucleotides all uncover highly polymorphic DNA fingerprint patterns in Southern analysis of genomic DNA that has been digested with a restriction enzyme having a 4-bp specificity. The polymorphic nature of the fragments is attributed to tandem repeat number variation of embedded minisatellite sequences. This explains why DNA fingerprint fragments are uncovered by minisatellite probes, but does not explain how it is that they are also uncovered by microsatellite and random oligonucleotide probes. To clarify this phenomenon, we sequenced a large bovine genomic BamHI restriction fragment hybridizing to the Jeffreys 33.6 minisatellite probe and consisting of small and large Sau3A-resistant subfragments. The large Sau3A subfragment was found to have a complex architecture, consisting of two different minisatellites, flanked and separated by stretches of unique DNA. The three unique sequences were characterized by sequence simplicity, that is, a higher than chance occurrence of tandem or dispersed repetition of simple sequence motifs. This complex repetitive structure explains the absence of Sau3A restriction sites in the large Sau3A subfragment, yet provides this subfragment with the ability to hybridize to a variety of probe sequences. It is proposed that a large class of interspersed tracts sharing this complex yet simplified sequence structure is found in the genome. Each such tract would have a broad ability to hybridize to a variety of probes, yet would exhibit a dearth of restriction sites. For each restriction enzyme having 4-bp specificity, a subclass of such tracts, completely lacking the corresponding restriction sites, will be present. On digestion with the given restriction enzyme, each such tract would form a large fragment.(ABSTRACT TRUNCATED AT 250 WORDS)

Characterization of minisatellites in Arabidopsis thaliana with sequence similarity to the human minisatellite core sequence

A strategy based on random PCR amplification was used to isolate new repetitive elements of Arabidopsis thaliana. One of the random PCR product analyzed by this approach contained a tandem repetitive minisatellite sequence composed of 33 bp repeated units. The genomic locus corresponding to this PCR product was isolated by screening a lambda genomic library. New related loci were also isolated from the genomic library by screening with a 14 mer oligonucleotide representing a region conserved among the different repeated units. Alignment of the consensus sequence for each minisatellite locus allowed the definition of an Arabidopsis thaliana core sequence that shows strong sequence similarities with the human core sequence and with the generalized recombination signal Chi of Escherichia coli. The minisatellites were tested for their ability to detect polymorphism, and their chromosomal position was established.

DNA variability and parentage testing in captive Waldrapp ibises

The Waldrapp ibis (Geronticus eremita) is one of the most endangered birds in the wild. The only real hope for its survival seems to rest with captive propagation, which has fortunately proved straightforward. Preservation of genetic diversity for maintaining viable populations is a major goal in any species conservation programme. Therefore sensitive tools have to be developed for the detection and regular monitoring of population variability. Individual identity has to be verified and each individual's pedigree defined in order to enable effective selective breeding. In this study, extensive DNA fingerprinting technology was applied to a group of 39 captive Waldrapp ibises. Five minisatellite probes were used. Despite the fact that the colony was founded by only six birds and that some of the founders were apparently related, highly polymorphic banding patterns were obtained. The human- and pig-derived probes MS1, MS31, p lambda g3 and pCMS12 each detected multiple distinct loci suitable for individual identification and resolution of parentage. Probe MS51 was locus-specific and revealed six different alleles. In addition, MS31 and p lambda g3 co-detected a band present in all females but none of the males and thus could be used for sex determination. For all but four individuals both parents could be identified unambiguously. The exceptions were offspring derived from a full sib mating in which DNA results allowed no distinction between the parents and the grandparents. Comparison of the parentage assigned by DNA analysis and by observation of breeding behaviour, respectively, showed agreement in all but five instances. Here, DNA data uncovered interchange of individuals in the stud book. Application of these minisatellite probes therefore offers an extremely sensitive means for individual identification and verification and completion of pedigree records.

Structural organisation of microsatellite families in the Leishmania genome and polymorphisms at two (CA)n loci

In the present study, we have analysed the frequency and distribution of several microsatellite DNAs [(CA)n, (GGT)n and (GCA)n] in the genome of Leishmania. Hybridisation analysis on the molecular karyotypes of different Leishmania strains showed the presence of these three microsatellites on all chromosomes of the parasite. The number of microsatellite clusters appeared grossly similar among strains from different Old World complexes. However, these three microsatellite families showed an uneven distribution among heterologous chromosomes of the same strain. Moreover, restriction analysis of chromosome I in various strains of Leishmania infantum showed a strong clustering of these microsatellites in the same chromosomal region. A partial genomic library was screened with a (CA)n probe, and 21 positive clones were isolated. The sequencing of these clones confirmed the association of various microsatellites such as (CA)n, (CT)n, and (GCA)n. Finally, specific polymerase chain reaction amplification of two cloned (CA)n loci demonstrated allelic size polymorphisms among strains within L. infantum and Leishmania donovani. Most of the 34 strains analysed were found to be monoallelic, while two alleles were found in a small number of strains. The interest of these sequences for studies on ploidy and population genetics of the parasite is discussed.

Chromosome allele loss in colorectal liver metastases and its association with clinical features

Molecular genetic changes are better characterized in colorectal carcinoma than in other common adult tumours. Consistent allele losses, or loss of heterozygosity (LOH), on chromosomes 5q, 17p and 18q have been well established. These changes are associated with the prognosis of the disease. Little is known of such changes in liver metastases of colorectal origin. The extent of allelic loss and its association with clinical features were investigated in 19 patients with colorectal liver metastases by using 24 probes to detect restriction fragment length polymorphism. A high frequency of LOH on chromosomes 5q, 17p and 18q was found in these secondary tumours. No consistent loss has so far been shown in any other chromosome. The frequency of allele loss correlated significantly with prognostic features such as the number and size of liver secondaries (P < 0.005), metastasis to the lymph nodes (P < 0.01) and curative or palliative operation (P < 0.02).

Complex gene conversion events in germline mutation at human minisatellites

Mutation at the human minisatellites MS32, MS205 and MS31A has been investigated by characterizing mutant alleles in pedigrees and in the case of MS32 by direct analysis of mutant molecules in single sperm. Most mutations at all three loci are polar, involving the preferential gain of a few repeat units at one end of the tandem repeat array. Incoming repeats can be derived from the same allele or the homologous chromosome, through they are frequently rearranged during mutation. Lack of exchange of flanking markers suggests the involvement of complex conversion-like events in the generation of mutant alleles. At MS32, high frequency mutation processes in sperm appear to be largely germline specific and to occur at a constant rate irrespective of allele size. Together with mutational polarity, this implies that germline instability is controlled by elements outside the tandem repeat array.

Chemically induced changes in the spectrum of amplifications of the human minisatellite MS1 integrated in chromosome III of a haploid yeast strain

To study chemically induced DNA amplifications we used the haploid Saccharomyces cerevisiae strain TR(MS1)-1 carrying an integrated chromosomal copy of the human minisatellite. MS1. Chemicals with different mechanisms of action were tested in this strain: methyl methanesulphonate, ethylene oxide (EO), propylene oxide (PO), camptothecin, 2,3,7,8-tetrachlorodibenso-p-dioxin (TCDD) and reserpine. No increase in frequency of new MS1 length alleles was seen with any of the tested chemicals relative to the spontaneous frequency of approximately 30%. EO and TCDD induced changes in the amplification spectrum, i.e., the frequency distribution of MS1 length alleles longer than the original 1.42 kb allele. PO and camptothecin increased the frequency of plasmid "pop-out" events. It seems likely that several mechanisms e.g. unequal exchanges, replication slippage and loop formation leading to deletion of a ring of tandem repeats, are involved in the generation of new MS1 length alleles. A loop-forming deletion mechanism is supported by the tendency to multimodality shown in the deamplification (loss of repeat units) spectra, i.e. the frequency distribution of new MS1 length alleles shorter than the original allele. EO and TCDD induced "longer" MS1 length alleles as compared to the control. The frequent generation of new MS1 length alleles in this haploid yeast strain further demonstrates the instability of such sequences and their possible relevance to genetic toxicology and the mechanisms of induction of cancer as well as other diseases. This study is a first step towards the development of an assay for DNA amplification without the use of a selective agent.

Possible role of natural selection in the formation of tandem-repetitive noncoding DNA

A simulation model of sequence-dependent amplification, unequal crossing over and mutation is analyzed. This model predicts the spontaneous formation of tandem-repetitive patterns of noncoding DNA from arbitrary sequences for a wide range of parameter values. Natural selection is found to play an essential role in this self-organizing process. Natural selection which is modeled as a mechanism for controlling the length of a nucleotide string but not the sequence itself favors the formation of tandem-repetitive structures. Two measures of sequence heterogeneity, inter-repeat variability and repeat length, are analyzed in detail. For fixed mutation rate, both inter-repeat variability and repeat length are found to increase with decreasing rates of (unequal) crossing over. The results are compared with data on micro-, mini- and satellite DNAs. The properties of minisatellites and satellite DNAs resemble the simulated structures very closely. This suggests that unequal crossing over is a dominant long-range ordering force which keeps these arrays homogeneous even in regions of very low recombination rates, such as at satellite DNA loci. Our analysis also indicates that in regions of low rates of (unequal) crossing over, inter-repeat variability is maintained at a low level at the expense of much larger repeat units (multimeric repeats), which are characteristic of satellite DNA. In contrast, the microsatellite data do not fit the proposed model well, suggesting that unequal crossing over does not act on these very short tandem arrays.

Abundance, polymorphism and genetic mapping of microsatellites in rice

Dinucleotide microsatellites have been characterized and used as genetic markers in rice. Screening of a rice genomic library with poly(dG-dA).(dC-dT) and poly(dG-dT).(dC-dA) probes indicated that (GA)n repeats occurred, on average, once every 225 kb and (GT)n repeats once every 480 kb. DNA sequencing of ten randomly selected microsatellites indicated that the numbers of repeats ranged from 12 to 34 and that the patterns of microsatellites in rice were similar to those of humans and other mammals. Primers to these microsatellite loci as well as to four published microsatellite-containing sequences have been designed and degrees of polymorphism has been examined with 20 rice accessions. Multiple alleles, ranging from 5 to 11, have been observed at all the microsatellite loci in 20 rice accessions. Alleles specific to two cultivated subspecies, indica and japonica, were found in some microsatellite loci. Heterozygosity values of all the microsatellite markers were significantly higher than those of RFLP markers, based upon a parallel comparison. Ten microsatellite loci have been genetically mapped to four rice chromosomes. The genomic distribution of microsatellites appears to be random in rice.

Characterization of tomato DNA clones with sequence similarity to human minisatellites 33.6 and 33.15

A tomato lambda genomic library was screened with the human minisatellites 33.6 and 33.15. Similar tomato sequences are estimated to occur on average every 4000 kb. In thirteen hybridizing clones characterized, the size of minisatellite arrays varied between 100 bp and 3 kb. The structure of the repetitive elements is complex as the human core sequence is interspersed with other elements. In three cases, sequences similar to the human minisatellites were part of a higher-order tandem repeat. The chromosomal position of these sequences was established by ascertaining linkage to previously mapped RFLP markers. In contrast to the human genome, no clustering of minisatellite loci was observed in tomato. The fingerprints generated by hybridizing tomato minisatellites to genomic DNA of a set of cultivars were, in two cases, more variable than those obtained with 33.6 or 33.15. Two of the characterized probes detected 4-8 alleles of a single locus, which displayed 10-15 times more polymorphism than random RFLP clones. Some minisatellites contain di- and tri-nucleotide microsatellite repeated motifs which may account for the high level of polymorphism detected with these clones.

VNTR and microsatellite polymorphisms within the subtelomeric region of 7q

The molecular basis of a highly polymorphic RFLP marker, HTY146c3 (D7S591), within the subtelomeric region of human chromosome 7q was determined by restriction-fragment and DNA sequence analysis. Two polymorphic systems were found--a simple base-substitution polymorphism and a GC-rich VNTR element with a core structure of C3AG2C2. In addition, a compound-imperfect CA dinucleotide-repeat element was identified approximately 10-20 kb from the telomeric sequence repeat (T2AG3), demonstrating that microsatellites can extend essentially to the ends of human chromosomes. The microsatellite marker, sAVH-6 (D7S594), is highly polymorphic, with 10 alleles and an observed heterozygosity of 84% found with the CEPH (Centre d'Etude du Polymorphisme Humain) reference pedigree collection. In combination with the RFLPs, the informativeness of the markers contained within 240 kb at the telomere approaches 100%. A unique genetic and physical STS marker, sAVH-6, defines the endpoint of the long arm of human chromosome 7.