Mitochondrial DNA reveals hidden diversity for tardigrades from across the Antarctic realm

Antarctica contains some of the most challenging environmental conditions on the planet due to freezing temperatures, prolonged winters and lack of liquid water. Whereas 99.7% of Antarctica is permanently covered by ice and snow, some coastal areas and mountain ridges have remained ice-free and are able to sustain populations of microinvertebrates. Tardigrades are one of the more dominant groups of microfauna in soil and limno-terrestrial habitats, but little is known of their diversity and distribution across Antarctica. Here, we examine tardigrades sampled from across an extensive region of continental Antarctica, and analyse and compare their partial mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences with those from the Antarctic Peninsula, maritime and sub-Antarctica, Tierra del Fuego and other worldwide locations in order to recognise operational taxonomic units (OTUs). From 439 new tardigrade COI sequences, we identified 98 unique haplotypes (85 from Antarctica) belonging to Acutuncus, Diphascon, Echiniscus, Macrobiotus, Milnesium and unidentified Parachela. Operational taxonomic units were delimited by Poisson tree processes and general mixed Yule coalescent methods, resulting in 58 and 55 putative species, respectively. Most tardigrades appear to be locally endemic (i.e. restricted to a single geographic region), but some (e.g. Acutuncus antarcticus (Richters, 1904)) are widespread across continental Antarctica. Our molecular results reveal: (i) greater diversity than has previously been appreciated with distinct OTUs that potentially represent undescribed species, and (ii) a lack of connectivity between most OTUs from continental Antarctica and those from other Antarctic geographical zones.

The high mobility group protein Abf2p influences the level of yeast mitochondrial DNA recombination intermediates in vivo

Abf2p is a high mobility group (HMG) protein found in yeast mitochondria that is required for the maintenance of wild-type (ρ+) mtDNA in cells grown on fermentable carbon sources, and for efficient recombination of mtDNA markers in crosses. Here, we show by two-dimensional gel electrophoresis that Abf2p promotes or stabilizes Holliday recombination junction intermediates in ρ+ mtDNA in vivo but does not influence the high levels of recombination intermediates readily detected in the mtDNA of petite mutants (ρ−). mtDNA recombination junctions are not observed in ρ+ mtDNA of wild-type cells but are elevated to detectable levels in cells with a null allele of the MGT1 gene (Δmgt1), which codes for a mitochondrial cruciform-cutting endonuclease. The level of recombination intermediates in ρ+ mtDNA of Δmgt1 cells is decreased about 10-fold if those cells contain a null allele of the ABF2 gene. Overproduction of Abf2p by ≥ 10-fold in wild-type ρ+ cells, which leads to mtDNA instability, results in a dramatic increase in mtDNA recombination intermediates. Specific mutations in the two Abf2p HMG boxes required for DNA binding diminishes these responses. We conclude that Abf2p functions in the recombination of ρ+ mtDNA.

Philopatry of male marine turtles inferred from mitochondrial DNA markers

Recent studies of mitochondrial DNA (mtDNA) variation among marine turtle populations are consistent with the hypothesis that females return to beaches in their natal region to nest as adults. In contrast, less is known about breeding migrations of male marine turtles and whether they too are philopatric to natal regions. Studies of geographic structuring of restriction fragment and microsatellite polymorphisms at anonymous nuclear loci in green turtle (Chelonia mydas) populations indicate that nuclear gene flow is higher than estimates from mtDNA analyses. Regional populations from the northern and southern Great Barrier Reef were distinct for mtDNA but indistinguishable at nuclear loci, whereas the Gulf of Carpentaria (northern Australia) population was distinct for both types of marker. To assess whether this result was due to reduced philopatry of males across the Great Barrier Reef, we determined the mtDNA haplotypes of breeding males at courtship areas for comparison with breeding females from the same three locations. We used a PCR-restriction fragment length polymorphism approach to determine control region haplotypes and designed mismatch primers for the identification of specific haplotypes. The mtDNA haplotype frequencies were not significantly different between males and females at any of the three areas and estimates of Fst among the regions were similar for males and females (Fst = 0.78 and 0.73, respectively). We conclude that breeding males, like females, are philopatric to courtship areas within their natal region. Nuclear gene flow between populations is most likely occurring through matings during migrations of both males and females through nonnatal courtship areas.

Mitochondrial DNA rearrangements of Podospora anserina are under the control of the nuclear gene grisea

Podospora anserina is a filamentous fungus with a limited life span. Life span is controlled by nuclear and extranuclear genetic traits. Herein we report the nature of four alterations in the nuclear gene grisea that lead to an altered morphology, a defect in the formation of female gametangia, and an increased life span. Three sequence changes are located in the 5′ upstream region of the grisea ORF. One mutation is a G → A transition at the 5′ splice site of the single intron of the gene, leading to a RNA splicing defect. This loss-of-function affects the amplification of the first intron of the mitochondrial cytochrome c oxidase subunit I gene (COI) and the specific mitochondrial DNA rearrangements that occur during senescence of wild-type strains. Our results indicate that the nuclear gene grisea is part of a molecular machinery involved in the control of mitochondrial DNA reorganizations. These DNA instabilities accelerate but are not a prerequisite for the aging of P. anserina cultures.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last ≈1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.

Analysis of microsatellite mutations in the mitochondrial DNA of Saccharomyces cerevisiae

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

Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in Myoclonus Epilepsy and Ragged-Red Fibers (MERRF) syndrome by a multiplex Molecular Beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNALys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.

Mitochondrial DNA ligase function in Saccharomyces cerevisiae

The Saccharomyces cerevisiae CDC9 gene encodes a DNA ligase protein that is targeted to both the nucleus and the mitochondria. While nuclear Cdc9p is known to play an essential role in nuclear DNA replication and repair, its role in mitochondrial DNA dynamics has not been defined. It is also unclear whether additional DNA ligase proteins are present in yeast mitochondria. To address these issues, mitochondrial DNA ligase function in S.cerevisiae was analyzed. Biochemical analysis of mitochondrial protein extracts supported the conclusion that Cdc9p was the sole DNA ligase protein present in this organelle. Inactivation of mitochondrial Cdc9p function led to a rapid decline in cellular mitochondrial DNA content in both dividing and stationary yeast cultures. In contrast, there was no apparent defect in mitochondrial DNA dynamics in a yeast strain deficient in Dnl4p (Δdnl4). The Escherichia coli EcoRI endonuclease was targeted to yeast mitochondria. Transient expression of this recombinant EcoRI endonuclease led to the formation of mitochondrial DNA double-strand breaks. While wild-type and Δdnl4 yeast were able to rapidly recover from this mitochondrial DNA damage, clones deficient in mitochondrial Cdc9p were not. These results support the conclusion that yeast rely upon a single DNA ligase, Cdc9p, to carry out mitochondrial DNA replication and recovery from both spontaneous and induced mitochondrial DNA damage.

A novel human hexameric DNA helicase: expression, purification and characterization

We have cloned, expressed and purified a hexameric human DNA helicase (hHcsA) from HeLa cells. Sequence analysis demonstrated that the hHcsA has strong sequence homology with DNA helicase genes from Saccharomyces cerevisiae and Caenorhabditis elegans, indicating that this gene appears to be well conserved from yeast to human. The hHcsA gene was cloned and expressed in Escherichia coli and purified to homogeneity. The expressed protein had a subunit molecular mass of 116 kDa and analysis of its native molecular mass by size exclusion chromatography suggested that hHcsA is a hexameric protein. The hHcsA protein had a strong DNA-dependent ATPase activity that was stimulated ≥5-fold by single-stranded DNA (ssDNA). Human hHcsA unwinds duplex DNA and analysis of the polarity of translocation demonstrated that the polarity of DNA unwinding was in a 5′→3′ direction. The helicase activity was stimulated by human and yeast replication protein A, but not significantly by E.coli ssDNA-binding protein. We have analyzed expression levels of the hHcsA gene in HeLa cells during various phases of the cell cycle using in situ hybridization analysis. Our results indicated that the expression of the hHcsA gene, as evidenced from the mRNA levels, is cell cycle-dependent. The maximal level of hHcsA expression was observed in late G1/early S phase, suggesting a possible role for this protein during S phase and in DNA synthesis.

Escherichia coli Nth and human hNTH1 DNA glycosylases are involved in removal of 8-oxoguanine from 8-oxoguanine/guanine mispairs in DNA

The spectrum of DNA damage caused by reactive oxygen species includes a wide variety of modifications of purine and pyrimidine bases. Among these modified bases, 7,8-dihydro-8-oxoguanine (8-oxoG) is an important mutagenic lesion. Base excision repair is a critical mechanism for preventing mutations by removing the oxidative lesion from the DNA. That the spontaneous mutation frequency of the Escherichia coli mutT mutant is much higher than that of the mutM or mutY mutant indicates a significant potential for mutation due to 8-oxoG incorporation opposite A and G during DNA replication. In fact, the removal of A and G in such a situation by MutY protein would fix rather than prevent mutation. This suggests the need for differential removal of 8-oxoG when incorporated into DNA, versus being generated in situ. In this study we demonstrate that E.coli Nth protein (endonuclease III) has an 8-oxoG DNA glycosylase/AP lyase activity which removes 8-oxoG preferentially from 8-oxoG/G mispairs. The MutM and Nei proteins are also capable of removing 8-oxoG from mispairs. The frequency of spontaneous G:C→C:G transversions was significantly increased in E.coli CC103mutMnthnei mutants compared with wild-type, mutM, nth, nei, mutMnei, mutMnth and nthnei strains. From these results it is concluded that Nth protein, together with the MutM and Nei proteins, is involved in the repair of 8-oxoG when it is incorporated opposite G. Furthermore, we found that human hNTH1 protein, a homolog of E.coli Nth protein, has similar DNA glycosylase/AP lyase activity that removes 8-oxoG from 8-oxoG/G mispairs.

Evidence for aminoacylation-induced conformational changes in human mitochondrial tRNAs.

Analysis by acid polyacrylamide/urea gel electrophoresis of 14 individual mitochondrial tRNAs (mt-tRNAs) from human cells has revealed a variable decrease in mobility of the aminoacylated relative to the nonacylated form, with the degree of separation of the two forms not being correlated with the mass, polar character, or charge of the amino acid. Separation of the charged and uncharged species has been found to be independent of tRNA denaturation, being observed also in the absence of urea. In another approach, electrophoresis through a perpendicular denaturing gradient gel of several individual mt-tRNAs has shown a progressive unfolding of the tRNA with increasing denaturant concentration, which is consistent with an initial disruption of tertiary interactions, followed by the sequential melting of the four stems of the cloverleaf structure. A detailed analysis of the unfolding process of charged and uncharged tRNALys and tRNALeu(UUR) has revealed that the separation of the two forms of these tRNAs persisted throughout the almost entire range of denaturant concentrations used and was lost upon denaturation of the last helical domain(s), which most likely included the amino acid acceptor stem. These observations strongly suggest that the electrophoretic retardation of the charged species reflects an aminoacylation-induced conformational change of the 3'-end of these mt-tRNAs, with possible significant implications in connection with the known role of the acceptor end in tRNA interactions with the ribosomal peptidyl transferase center and the elongation factor Tu.

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

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