Evolution of Repeated Sequence Arrays in the D-Loop Region of Bat Mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, θ and π, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.     

Detection of species-specific long VNTRs in mitochondrial control region and their application to identifying sympatric Hong Kong grouper (Epinephelus akaara) and yellow grouper (Epinephelus awoara)

Length variation and heteroplasmy were observed in PCR products of the first half of mtDNA control region of both Hong Kong grouper (Epinephelus akaara) and yellow grouper (Epinephelus awoara). DNA sequencing unveiled the phenomena were caused by the presence of species-specific long variable number tandem repeats (VNTRs). This is the first report on the mtDNA VNTRs and their heteroplasmy in groupers. Moreover, these VNTRs are also the longest such structure found in teleost fish. Thereafter, we designed two species-specific PCR reverse primers according to the 3' end sequences of the VNTRs and successfully established assays for the identification of these two sympatric grouper species.     

Heteroplasmy in the mtDNA control region of sturgeon (Acipenser, Huso and Scaphirhynchus)

Data from 1238 fishes from 19 sturgeon species and 1 paddlefish were used to analyze heteroplasmy in sturgeon. Lengths of central repeat units ranged from 74 to 83 bp among sturgeon species. No repeat sequence was found in the paddlefish, Polyodon spathula. A general feature of the repeat units was the presence of termination associated sequence (TAS) motifs. About 50% of 138 interspecific mutations observed among the D-loop sequences are located 10 bp down- and upstream from these TAS motifs. Interestingly, most homoplasmic species showed deletions upstream to the TAS motifs, whereas deletions downstream to the TAS motifs observed in two species do not seem to preclude heteroplasmy. Calculations of secondary structures and thermal stabilities of repeat units showed DeltaG values for all heteroplasmic species to be <-8 and for most homoplasmic species DeltaG value to be >-8. Most heteroplasmic fishes had two and/or three repeat units. No homoplasmic sturgeon with >2 repeat units were observed. Molecular phylogeny based on the entire cytochrome b showed that heteroplasmy probably resulted from a single evolutionary event. Our data demonstrate that heteroplasmy is present in most sturgeon species and suggest that the thermal stability of the secondary structure of the repeat unit in combination with mutations downstream of the TAS sequences influences heteroplasmy.     

Heteroplasmy of Short Tandem Repeats in Mitochondrial DNA of Atlantic Cod, Gadus Morhua

The mitochondrial DNA of the Atlantic cod (Gadus morhua) contains a tandem array of 40-bp repeats in the D-loop region of the molecule. Variation among molecules in the copy number of these repeats results in mtDNA length variation and heteroplasmy (the presence of more than one form of mtDNA in an individual). In a sample of fish collected from different localities around Iceland and off George's Bank, each individual was heteroplasmic for two or more mtDNAs ranging in repeat copy number from two (common) to six (rare). An earlier report on mtDNA heteroplasmy in sturgeon (Acipenser transmontanus) presented a competitive displacement model for length mutations in mtDNAs containing tandem arrays and the cod data deviate from this model. Depending on the nature of putative secondary structures and the location of D-loop strand termination, additional mechanisms of length mutation may be needed to explain the range of mtDNA length variants maintained in these populations. The balance between genetic drift and mutation in maintaining this length polymorphism is estimated through a hierarchical analysis of diversity of mtDNA length variation in the Iceland samples. Eighty percent of the diversity lies within individuals, 8% among individuals and 12% among localities. An estimate of theta = 2N(eo) mu greater than 1 indicates that this system is characterized by a high mutation rate and is governed primarily by deterministic dynamics. The sequences of repeat arrays from fish collected in Norway, Iceland and George's Bank show no nucleotide variation suggesting that there is very little substructuring to the North Atlantic cod population.     

Polymorphism and structural features of two noncoding regions of the liver fluke <Emphasis Type="Italic">Fasciola hepatica</Emphasis> (Plathelminthes: Trematoda) mitochondrial genome

Structural characteristics and polymorphism of the long (LNR) and short (SNR) mitochondrial noncoding regions were studied in the liver fluke Fasciola hepatica. Flukes were sampled from several populations of Russia and Belarus. LNR amplification yielded a set of nine fragments, neighboring ones differing in length by one tandem repeat (85 bp), published for Australian flukes. The LNR amplification products of different lengths were cloned and sequenced. A comparison of the LNR sequences of Australian and Belarussian flukes revealed three nucleotide substitutions and one point heteroplasmy in the first positions of the imperfect repeat and four adjacent perfect repeats. The positions of the three mutations coincided in the perfect and imperfect repeats. The frequency of mutations was 4.0–4.7 %, while the frequency of heteroplasmic sites varied from 0.1 to 1.2%. It was shown that the mutations and the heteroplasmy of one site could change the structure and stability of the putative secondary structures of the perfect and imperfect repeats. SNR amplification in F. hepatica from several populations yielded fragments that differed from the published SNR sequence of Australian F. hepatica by one transversion (T → G in position 21). Both noncoding regions had several conserved and potential regulatory sequences. The possible causes of heteroplasmy and a concerted origin of substitutions in different repeats are discussed.     

Tandemly repeated sequences in mtDNA control region of whitefish, Coregonus lavaretus.

Length variation of the mitochondrial DNA control region was observed with PCR amplification of a sample of 138 whitefish (Coregonus lavaretus). Nucleotide sequences of representative PCR products showed that the variation was due to the presence of an approximately 100-bp motif tandemly repeated two, three, or five times in the region between the conserved sequence block-3 (CSB-3) and the gene for phenylalanine tRNA. This is the first report on the tandem array composed of long repeat units in mitochondrial DNA of salmonids.     

Trinucleotide repeats associated with human disease.

Triplet repeat expansion diseases (TREDs) are characterized by the coincidence of disease manifestation with amplification of d(CAG. CTG), d(CGG.CCG) or d(GAA.TTC) repeats contained within specific genes. Amplification of triplet repeats continues in offspring of affected individuals, which generally results in progressive severity of the disease and/or an earlier age of onset, phenomena clinically referred to as 'anticipation'. Recent biophysical and biochemical studies reveal that five of the six [d(CGG)n, d(CCG)n, (CAG)n, d(CTG)n and d(GAA)n] complementary sequences that are associated with human disease form stable hairpin structures. Although the triplet repeat sequences d(GAC)n and d(GTC)n also form hairpins, repeats of the double-stranded forms of these sequences are conspicuously absent from DNA sequence databases and are not anticipated to be associated with human disease. With the exception of d(GAG)n and d(GTG)n, the remaining triplet repeat sequences are unlikely to form hairpin structures at physiological salt and temperature. The details of hairpin structures containing trinucleotide repeats are summarized and discussed with respect to potential mechanisms of triplet repeat expansion and d(CGG.CCG) n methylation/demethylation.     

Mitochondrial inheritance in Schistosoma mansoni: mitochondrial variable number tandem repeat mutation produces noise on top of the signal.

The Schistosoma mansoni mitochondrial genome contains tandemly arrayed copies of a 62-base repeat motif. The tandem array is highly polymorphic with respect to number of repeats and commonly exhibits heteroplasmy. This study shows that a very high rate of mutation rapidly produces new repeat lengths (new haplotypes) for this mitochondrial variable number tandem repeat. A maternal inheritance pattern is also demonstrated for this repeat sequence, while the high mutation rate causes some offspring to exhibit nonmaternal haplotypes. Frequent generation of new haplotypes can be observed within samples of clonal cohorts taken from monomiracidial snail infections. These same clonal cercarial groups, when crossed, produce F1 generations that exhibit the maternal set of haplotypes, across all individuals, with the frequent addition of new mutant haplotypes. In each of 2 crosses, a subset of the recently arisen haplotypes match paternal haplotypes by chance (30.4% and 18.8%), thus giving the false appearance of partial paternal inheritance of mitochondria.     

Instability of CAG and CTG trinucleotide repeats in Saccharomyces cerevisiae.

A quantitative genetic assay was developed to monitor alterations in tract lengths of trinucleotide repeat sequences in Saccharomyces cerevisiae. Insertion of (CAG)50 or (CTG)50 repeats into a promoter that drives expression of the reporter gene ADE8 results in loss of expression and white colony color. Contractions within the trinucleotide sequences to repeat lengths of 8 to 38 restore functional expression of the reporter, leading to red colony color. Reporter constructs including (CAG)50 or (CTG)50 repeat sequences were integrated into the yeast genome, and the rate of red colony formation was measured. Both orientations yielded high rates of instability (4 x 10(-4) to 18 x 10(-4) per cell generation). Instability depended on repeat sequences, as a control harboring a randomized (C,A,G)50 sequence was at least 100-fold more stable. PCR analysis of the trinucleotide repeat region indicated an excellent correlation between change in color phenotype and reduction in length of the repeat tracts. No preferential product sizes were observed. Strains containing disruptions of the mismatch repair gene MSH2, MSH3, or PMS1 or the recombination gene RAD52 showed little or no difference in rates of instability or distributions of products, suggesting that neither mismatch repair nor recombination plays an important role in large contractions of trinucleotide repeats in yeast.     

Direct PCR sequencing of dystrophin polymorphic CACA alleles after purification to remove shadow bands.

A method is described that allows the sequencing of polymerase chain reaction (PCR) products containing CACA repeats. The method was tested using a DNA polymorphism that exists at the 3' end of the dystrophin gene. This polymorphism consists of a variation in the length of a CACA dinucleotide repeat. Four alleles from a total of 16 individuals were sequenced at this locus after the DNA sequence had been amplified by the PCR. Five examples of each of the common alleles were sequenced. For each allele all five sequences were the same. The only example of a rare allele was also sequenced. The PCR products of DNA sequences containing dinucleotide repeats consist of a number of bands differing by 2 bp below the most intense main band. Previously, direct sequencing of the PCR products lead to ambiguities and smearing at and above the CACA repeat. In this paper, the main PCR band was cut out of a sequencing gel and directly sequenced to give a clear DNA sequence. Our results indicate that for a particular allele, all individuals had exactly the same DNA sequence. This implies that with the appropriate choice of oligonucleotide primers, polymorphisms could be detected without electrophoresis.     

Detection of male and female fetal DNA in maternal plasma by multiplex fluorescent polymerase chain reaction amplification of short tandem repeats

The purpose of this study was to develop a fluorescent polymerase chain reaction (PCR) assay for the detection of circulating fetal DNA in maternal plasma. Maternal DNA extracted from plasma samples of pregnant women at term and newborn DNA isolated from cord blood were used to genotype 12 mother/child pairs at nine different polymorphic short tandem repeat loci. Multiplex fluorescent PCR was used to detect fetus-specific alleles in the corresponding maternal plasma samples. Fetus-specific alleles were found in all maternal plasma samples studied. Using these polymorphic repeat sequences, every mother/child pair was informative in at least four of nine loci. Paternally inherited fetal alleles were detected in 84% of informative short tandem repeats. This approach may have implications for non-invasive prenatal diagnosis. Compared with other fetal DNA detection systems that use fetus-derived Y sequences to detect only male fetal DNA in maternal plasma, our proposed technique can be applied to both female and male fetuses.     

To slip or skip, visualizing frameshift mutation dynamics for error-prone DNA polymerases

Three models describing frameshift mutations are "classical" Streisinger slippage, proposed for repetitive DNA, and "misincorporatation misalignment" and "dNTP-stabilized misalignment," proposed for non-repetitive DNA. We distinguish between models using pre-steady state fluorescence kinetics to visualize transiently misaligned DNA intermediates and nucleotide incorporation products formed by DNA polymerases adept at making small frameshift mutations in vivo. Human polymerase (pol) mu catalyzes Streisinger slippage exclusively in repetitive DNA, requiring as little as a dinucleotide repeat. Escherichia coli pol IV uses dNTP-stabilized misalignment in identical repetitive DNA sequences, revealing that pol mu and pol IV use different mechanisms in repetitive DNA to achieve the same mutational end point. In non-repeat sequences, pol mu switches to dNTP-stabilized misalignment. pol beta generates -1 frameshifts in "long" repeats and base substitutions in "short" repeats. Thus, two polymerases can use two different frameshift mechanisms on identical sequences, whereas one polymerase can alternate between frameshift mechanisms to process different sequences.     

FMR1 gene deletion/reversion: a pitfall of fragile X carrier testing

The Fragile X syndrome is, in the majority of cases, caused by CGG trinucleotide amplification within the FMR1 gene. The syndrome is rarely caused by point mutations or deletions. Here we describe a family with 2 sons and 1 daughter affected by Fragile X syndrome and 2 unaffected daughters whose carrier status was unknown prior to this study. Analysis of DNA from each of the 2 daughters revealed two alleles in the normal size range. However, 1 daughter carried one allele of 10 CGG repeats that was not present in either the mother or the father. No evidence for mosaicism could be detected. Haplotype analysis of flanking polymorphic markers revealed that the 10 CGG allele was derived from the mutated allele inherited from the mother. Thus, this case most likely represents an additional case of a reverse mutation from a premutation allele in a female to a normal-sized allele in the offspring. It remains unclear how frequently such reversion events occur. The observation has important consequences for genetic testing, because many laboratories prescreen for the Fragile X syndrome by determining the length of the CGG repeat using PCR. If this shows alleles in the normal size range, a diagnosis of Fragile X syndrome is considered to be excluded. Because the routine PCR and/or Southern blot analyses alone may yield false-negative results in cases of a regression of the number of CGG repeats, we strongly recommend the inclusion of fragment length or haplotype analysis when determining the carrier status within Fragile X syndrome families.     

The use of next generation sequencing technology to study the effect of radiation therapy on mitochondrial DNA mutation

The human mitochondrial genome has an exclusively maternal mode of inheritance. Mitochondrial DNA (mtDNA) is particularly vulnerable to environmental insults due in part to an underdeveloped DNA repair system, limited to base excision and homologous recombination repair. Radiation exposure to the ovaries may cause mtDNA mutations in oocytes, which may in turn be transmitted to offspring. We hypothesized that the children of female cancer survivors who received radiation therapy may have an increased rate of mtDNA heteroplasmy mutations, which conceivably could increase their risk of developing cancer and other diseases. We evaluated 44 DNA blood samples from 17 Danish and 1 Finnish families (18 mothers and 26 children). All mothers had been treated for cancer as children and radiation doses to their ovaries were determined based on medical records and computational models. DNA samples were sequenced for the entire mitochondrial genome using the Illumina GAII system. Mother's age at sample collection was positively correlated with mtDNA heteroplasmy mutations. There was evidence of heteroplasmy inheritance in that 9 of the 18 families had at least one child who inherited at least one heteroplasmy site from his or her mother. No significant difference in single nucleotide polymorphisms between mother and offspring, however, was observed. Radiation therapy dose to ovaries also was not significantly associated with the heteroplasmy mutation rate among mothers and children. No evidence was found that radiotherapy for pediatric cancer is associated with the mitochondrial genome mutation rate in female cancer survivors and their children.     

Molecular instability in the COII-tRNA(Lys) intergenic region of the human mitochondrial genome: multiple origins of the 9-bp deletion and heteroplasmy for expanded repeats.

We have identified two individuals from Glasgow in Scotland who have a deletion of one of two copies of the intergenic 9-bp sequence motif CCCCCTCTA, located between the cytochrome oxidase II (COII) and lysine tRNA (tRNA(Lys)) genes of the human mitochondrial genome. Although this polymorphism is common in Africa and Asia, it has not been reported in Northern Europe. Analysis of the mitochondrial DNA control region sequences of these two individuals suggests that they belong to a lineage that originated independently of the previously characterized African and Asian 9-bp deleted lineages. Among the Scottish population we have also identified a maternal lineage of three generations exhibiting heteroplasmy for two, three and four copies of the CCCCCTCTA motif. Polymerase chain reaction amplification across the COII-tRNA(Lys) intergenic region of these individuals gives different ratios of the three product lengths that are dependent on the concentration of the DNA-binding dye crystal violet. To investigate whether changes in repeat number were generated de novo, we constructed clones containing known numbers of the CCCCCTCTA motif. In the presence of high concentrations of crystal violet we obtained two, three and four copies of this motif when the amplification template contained only four copies. Various DNA-binding drugs are known to stabilize bulged structures in DNA and contribute to the process of slipped-strand mispairing during DNA replication. These results suggest that the COII-tRNA(Lys) intergenic region is unstable owing to slipped-strand mispairing. Although sequences containing four copies of the CCCCCTCTA motif are less stable in vitro, we observed an increase in the proportion of mitochondrial genomes with four repeats between-a mother and a daughter in the heteroplasmic lineage. From this we conclude that drift in the germ-line lineage is a main factor in the maintenance or loss of heteroplasmy.     

Characterization of an unusual DNA length polymorphism 5'' to the human antithrombin III gene

Nucleotide sequence analysis revealed that a DNA length polymorphism 5' to the human antithrombin III gene is due to the presence of 32bp or 108bp nonhomologous nucleotide sequences (variable segments) 345bp upstream from the translation initiation codon. Sequences at the 3' borders of both variable segments can form intrastrand inverted repeat structures with sequences further downstream. An inverted repeat is also found immediately 5' to the site where the variable segments are located. Thus, cruciform structures may form flanking the variable segments of both alleles of this DNA length polymorphism. DNA secondary structure may be detected with single strand specific nucleases. S1 nuclease sensitive sites were mapped in recombinant plasmids containing the cloned alleles of the ATIII length polymorphism. The site most sensitive to S1 is located upstream from the variable segments in an AT-rich segment flanked by 6bp direct repeats. A region of lesser nuclease sensitivity was also observed in the AT-rich loops formed between the inverted repeats 5' to the variable segments.     

Junctions between repetitive DNAs on the PSR chromosome of Nasonia vitripennis: Association of palindromes with recombination

The Paternal-Sex-Ratio (PSR) chromosome of Nasonia vitripennis contains several families of repetitive DNAs that show significant sequence divergence but share two palindromic regions. This study reports on the analysis of junctions between two of these repetitive DNA families (psr2 and psr18). Three lambda clones that hybridized to both repeat families were isolated from PSR-genomic DNA libraries through multiple screenings and analyzed by Southern blots. Analysis of clones showed a region in which the two repeat types are interspersed, flanked by uniform blocks of each repeat type. PCR amplification of genomic DNA confirmed the contiguous arrangement of psr2 and psr18 on PSR and identified an additional junction region between these repeats that was not present in the lambda inserts. We isolated and sequenced 41 clones from the lambda inserts and genomic PCR products containing junction sequences. Sequence analysis showed that all transitions between psr2 and psr18 repeats occurred near one of the two palindromes. Based on the inheritance pattern of PSR, recombination between repeats on this chromosome must be mitotic (rather than meiotic) in origin. The occurrence of exchanges near the palindromes suggests that these sequences enhance recombination between repeat units. Rapid amplification of repetitive DNA may have been an important factor in the evolution of the PSR chromosome. Correspondence to: John H. Werren