Evolution of the mammalian mitochondrial control region--comparisons of control region sequences between monotreme and therian mammals

A new type of tandem repeats in Myotis petax (Chiroptera, Vespertilionidae) mitochondrial control region

BACKGROUND

Tandem repeats in mitochondrial DNA control region are known to different animal taxa, including bat species of the family Vespertilionidae. The long R1-repeats in the bat ETAS-domain are often presented in a variable copy number and may exhibit both inter-individual and intra-individual sequence diversity. The function of repeats in the control region is still unclear, but it has been shown that repetitive sequences in some animal groups (shrews, cats and sheep) may include parts of ETAS1 and ETAS2 conservative blocks of mitochondrial DNA.

METHODS AND RESULTS

Analysis of the control region sequences for 31 Myotis petax specimens allowed the identification of the inter-individual variability and clarification of the composition of the R1-repeats. The copy number of the R1-repeats varies from 4 to 7 in individuals. The specimens examined do not exhibit a size heteroplasmy previously described for Myotis species. The unusual short 30 bp R1-repeats have been detected in M. petax for the first time. The ten specimens from Amur Region and Primorsky Territory have one or two copies of these additional repeats.

CONCLUSIONS

It was determined that the R1-repeats in M. petax control region consist of parts of the ETAS1 and ETAS2 blocks. The origin of the additional repeats seems to be related to the 51 bp deletion in the central part of the R1-repeat unit and subsequent duplication. Comparison of repetitive sequences in the control region of closely-related Myotis species identified the occurrence of incomplete repeats also resulting from the short deletions, but distinct from additional repeats of M. petax.

Mitochondrial Genomes Assembled from Non-Invasive eDNA Metagenomic Scat Samples in Critically Endangered Mammals

The abundance of many large-bodied vertebrates, both in marine and terrestrial environments, has declined substantially due to global and regional climate stressors that define the Anthropocene. The development of genetic tools that can serve to monitor population’s health non-intrusively and inform strategies for the recovery of these species is crucial. In this study, we formally evaluate whether whole mitochondrial genomes can be assembled from environmental DNA (eDNA) metagenomics scat samples. Mitogenomes of four different large vertebrates, the panda bear (Ailuropoda melanoleuca), the moon bear (Ursus thibetanus), the Java pangolin (Manis javanica), and the the North Atlantic right whale (Eubalaena glacialis) were assembled and circularized using the pipeline GetOrganelle with a coverage ranging from 12x to 480x in 14 out of 18 different eDNA samples. Partial mitochondrial genomes were retrieved from three other eDNA samples. The complete mitochondrial genomes of the studied species were AT-rich and comprised 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a putative D-loop/control region. Synteny observed in all assembled mitogenomes was identical to that reported for specimens of the same and other closely related species. This study demonstrates that it is possible to assemble accurate whole mitochondrial chromosomes from eDNA samples (scats) using forthright bench and bioinformatics workflows. The retrieval of mitochondrial genomes from eDNA samples represents a tool to support bioprospecting, bio-monitoring, and other non-intrusive conservation strategies in species considered ‘vulnerable’, ‘endangered’, and/or ‘critically endangered’ by the IUCN Red List of Threatened Species.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Complete mitochondrial genome sequence of Lepus yarkandensis Günther, 1875 (Lagomorpha, Leporidae): characterization and phylogenetic analysis

Lepusyarkandensis is a national second-class protected animal endemic to China and distributed only in the hot and arid Tarim Basin in Xinjiang. We sequenced and described the complete mitogenome of L.yarkandensis to analyze its characteristics and phylogeny. The species’ DNA is a 17,047 bp circular molecule that includes 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. The overall base composition was as follows: A, 31.50%; T, 29.40%; G, 13.30% and C, 25.80%, with a high A+T bias of 60.9%. In the PCGs, ND6 had deviation ranges for AT skew (–0.303) and GC skew (0.636). The Ka/Ks values of ND1 (1.067) and ND6 (1.352) genes were >1, indicating positive selection, which might play an important role in the adaptation of L.yarkandensis to arid and hot environments. The conserved sequence block, the central conserved domain, and the extended termination-associated sequences of the control region and their features were identified and described. The phylogenetic tree based on the complete mitogenome showed that L.yarkandensis was closely related to the sympatric Lepustibetanuspamirensis. These novel datasets of L.yarkandensis can supply basic data for phylogenetic studies of Lepus spp., apart from providing essential and important resource for further genetic research and the protection of this species.

Structural Variation of the Turtle Mitochondrial Control Region

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.

Forensic validation of a panel of 12 SNPs for identification of Mongolian wolf and dog

Wolf (Canis lupus) is a species included in appendices of CITES and is often encountered in cases of alleged poaching and trafficking of their products. When such crimes are suspected, those involved may attempt to evade legal action by claiming that the animals involved are domestic dogs (C. l. familiaris). To respond effectively to such claims, law enforcement agencies require reliable and robust methods to distinguish wolves from dogs. Reported molecular genetic methods are either unreliable (mitogenome sequence based), or operationally cumbersome and require much DNA (un-multiplexed microsatellites), or financially expensive (genome wide SNP genotyping). We report on the validation of a panel of 12 ancestral informative single nucleotide polymorphism (SNP) markers for discriminating wolves from dogs. A SNaPshot multiplex genotyping system was developed for the panel, and 97 Mongolian wolves (C. l. chanco) and 108 domestic dogs were used for validation. Results showed this panel had high genotyping success (0.991), reproducibility (1.00) and origin assignment accuracy (0.97 ± 0.05 for dogs and 1.00 ± 0.03 for wolves). Species-specificity testing suggested strong tolerance to DNA contamination across species, except for Canidae. The minimum DNA required for reliable genotyping was 6.25 pg/μl. The method and established gene frequency database are available to support identification of wolves and dogs by law enforcement agencies.

Weird Animals, Sex, and Genome Evolution

Making my career in Australia exposed me to the tyranny of distance, but it gave me opportunities to study our unique native fauna. Distantly related animal species present genetic variation that we can use to explore the most fundamental biological structures and processes. I have compared chromosomes and genomes of kangaroos and platypus, tiger snakes and emus, devils (Tasmanian) and dragons (lizards). I particularly love the challenges posed by sex chromosomes, which, apart from determining sex, provide stunning examples of epigenetic control and break all the evolutionary rules that we currently understand. Here I describe some of those amazing animals and the insights on genome structure, function, and evolution they have afforded us. I also describe my sometimes-random walk in science and the factors and people who influenced my direction. Being a woman in science is still not easy, and I hope others will find encouragement and empathy in my story.

The First Mitogenome of the Cyprus Mouflon (Ovis gmelini ophion): New Insights into the Phylogeny of the Genus Ovis

Sheep are thought to have been one of the first livestock to be domesticated in the Near East, thus playing an important role in human history. The current whole mitochondrial genome phylogeny for the genus Ovis is based on: the five main domestic haplogroups occurring among sheep (O. aries), along with molecular data from two wild European mouflons, three urials, and one argali. With the aim to shed some further light on the phylogenetic relationship within this genus, the first complete mitochondrial genome sequence of a Cypriot mouflon (O. gmelini ophion) is here reported. Phylogenetic analyses were performed using a dataset of whole Ovis mitogenomes as well as D-loop sequences. The concatenated sequence of 28 mitochondrial genes of one Cypriot mouflon, and the D-loop sequence of three Cypriot mouflons were compared to sequences obtained from samples representatives of the five domestic sheep haplogroups along with samples of the extant wild and feral sheep. The sample included also individuals from the Mediterranean islands of Sardinia and Corsica hosting remnants of the first wave of domestication that likely went then back to feral life. The divergence time between branches in the phylogenetic tree has been calculated using seven different calibration points by means of Bayesian and Maximum Likelihood inferences. Results suggest that urial (O. vignei) and argali (O. ammon) diverged from domestic sheep about 0.89 and 1.11 million years ago (MYA), respectively; and dates the earliest radiation of domestic sheep common ancestor at around 0.3 MYA. Additionally, our data suggest that the rise of the modern sheep haplogroups happened in the span of time between six and 32 thousand years ago (KYA). A close phylogenetic relationship between the Cypriot and the Anatolian mouflon carrying the X haplotype was detected. The genetic distance between this group and the other ovine haplogroups supports the hypothesis that it may be a new haplogroup never described before. Furthermore, the updated phylogenetic tree presented in this study determines a finer classification of ovine species and may help to classify more accurately new mitogenomes within the established haplogroups so far identified.

Mitochondrial DNA control region of three mackerels, genus Rastrelliger: structure, molecular diversity and phylogenetic relationship

The complete mitochondrial control regions (CR) of three mackerels (Rastrelliger spp.) were examined and analyzed. The CR contained three domains, in which three termination-associated sequences (TAS-I, TAS-II and TAS-III), two central conserved sequence blocks (CSB-E, CSB-D), three conserved sequence blocks (CSB-I, CSB-II, and CSB-III) and a putative promoter were detected. Molecular indices analyses of the aligned complete CR sequences showed high level of haplotype diversities and genetic divergences among the three species. The intraspecific divergence among species of this genus ranked from 0.25% to 1.62% and interspecific divergence from 1.90% to 4.30%. The phylogenetic tree shows monophyly with R. brachysoma as a basal species of Rastrelliger. Applying the average divergence rate for fish control regions, the results suggest that the time of separation among Rastrelligers could have occurred in the middle Pleistocene era.

Comparative mitochondrial genomics and phylogenetic relationships of the Crossoptilon species (Phasianidae, Galliformes)

BACKGROUND

Phasianidae is a family of Galliformes containing 38 genera and approximately 138 species, which is grouped into two tribes based on their morphological features, the Pheasants and Partridges. Several studies have attempted to reconstruct the phylogenetic relationships of the Phasianidae, but many questions still remain unaddressed, such as the taxonomic status and phylogenetic relationships among Crossoptilon species. The mitochondrial genome (mitogenome) has been extensively used to infer avian genetic diversification with reasonable resolution. Here, we sequenced the entire mitogenomes of three Crossoptilon species (C. harmani, C. mantchuricum and C. crossoptilon) to investigate their evolutionary relationship among Crossoptilon species.

RESULTS

The complete mitogenomes of C. harmani, C. mantchuricum and C. crossoptilon are 16682 bp, 16690 bp and 16680 bp in length, respectively, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a putative control region. C. auritum and C. mantchuricum are more closely related genetically, whereas C. harmani is more closely related to C. crossoptilon. Crossoptilon has a closer relationship with Lophura, and the following phylogenetic relationship was reconstructed: ((Crossoptilon + Lophura) + (Phasianus + Chrysolophus)). The divergence time between the clades C. harmani-C. crossoptilon and C. mantchuricum-C. auritum is consistent with the uplift of the Tibetan Plateau during the Tertiary Pliocene. The Ka/Ks analysis showed that atp8 gene in the Crossoptilon likely experienced a strong selective pressure in adaptation to the plateau environment.

CONCLUSIONS

C. auritum with C. mantchuricum and C. harmani with C. crossoptilon form two pairs of sister groups. The genetic distance between C. harmani and C. crossoptilon is far less than the interspecific distance and is close to the intraspecific distance of Crossoptilon, indicating that C. harmani is much more closely related to C. crossoptilon. Our mito-phylogenomic analysis supports the monophyly of Crossoptilon and its closer relationship with Lophura. The uplift of Tibetan Plateau is suggested to impact the divergence between C. harmani-C. crossoptilon clade and C. mantchuricum-C. auritum clade during the Tertiary Pliocene. Atp8 gene in the Crossoptilon species might have experienced a strong selective pressure for adaptation to the plateau environment.

Molecular characteristics and evolution of the mitochondrial control region in three genera (Hipposideridae: Hipposideros Aselliscus and Coelops) of leaf-nosed bats

The mitochondrial control region (CR) was sequenced for three genera of Hipposideridae to give a detailed overview of its features. The CR of leaf-nosed bats (1288-1560 bp) was divided into three domains like that of other mammals. In addition to the common conserved blocks (ETAS1, ETAS2, F-B boxes, CSB1, CSB2, and CSB3) found in all species, a CSB1-like element was also detected in the conserved sequence blocks (CSB). Repeated motifs were examined in the ETAS of Aselliscus stoliczkanus (26 bp) and Hipposideros bicolor (80 bp) and were present in the CSB of all individuals (6, 8, 16, and 20 bp). Phylogenetic reconstructions using the CR sequences indicated that the phylogenetic relationships among Hipposideros species were consistent with the results of other molecular and phenetic analyses. Aselliscus and Coelops had a closer relationship. But the central domain could not be used for phylogenetic analyses at family and genus levels due to its high conservation.

The strength and timing of the mitochondrial bottleneck in salmon suggests a conserved mechanism in vertebrates

In most species mitochondrial DNA (mtDNA) is inherited maternally in an apparently clonal fashion, although how this is achieved remains uncertain. Population genetic studies show not only that individuals can harbor more than one type of mtDNA (heteroplasmy) but that heteroplasmy is common and widespread across a diversity of taxa. Females harboring a mixture of mtDNAs may transmit varying proportions of each mtDNA type (haplotype) to their offspring. However, mtDNA variants are also observed to segregate rapidly between generations despite the high mtDNA copy number in the oocyte, which suggests a genetic bottleneck acts during mtDNA transmission. Understanding the size and timing of this bottleneck is important for interpreting population genetic relationships and for predicting the inheritance of mtDNA based disease, but despite its importance the underlying mechanisms remain unclear. Empirical studies, restricted to mice, have shown that the mtDNA bottleneck could act either at embryogenesis, oogenesis or both. To investigate whether the size and timing of the mitochondrial bottleneck is conserved between distant vertebrates, we measured the genetic variance in mtDNA heteroplasmy at three developmental stages (female, ova and fry) in chinook salmon and applied a new mathematical model to estimate the number of segregating units (N(e)) of the mitochondrial bottleneck between each stage. Using these data we estimate values for mtDNA Ne of 88.3 for oogenesis, and 80.3 for embryogenesis. Our results confirm the presence of a mitochondrial bottleneck in fish, and show that segregation of mtDNA variation is effectively complete by the end of oogenesis. Considering the extensive differences in reproductive physiology between fish and mammals, our results suggest the mechanism underlying the mtDNA bottleneck is conserved in these distant vertebrates both in terms of it magnitude and timing. This finding may lead to improvements in our understanding of mitochondrial disorders and population interpretations using mtDNA data.

The structure of the Australian and South American marsupial mitochondrial control region

BACKGROUND AND AIMS

The mitochondrial control region (CR) was studied across five marsupialian orders, in order to give a detailed overview of its features.

RESULTS

The CR is organised into three domains similar to the CR of placental mammals. However, the conservation of different features among the marsupial orders is in general more strict. In the first domain, two conserved blocks extended termination-associated sequences (ETAS 1 and ETAS 2) are present in all marsupial orders. In the third domain, the three conserved sequence blocks (CSB 1, CSB 2 and CSB 3) are present and complete, with CSB 1 being duplicated in four of five marsupial orders.

CONCLUSIONS

The nucleotide frequency and secondary structures of the repeats were typical for marsupial species. The repeats are generally AT-rich except in Dasyuridae and Paucituberculata, which show a significant increase in GC content.

ND6 gene "lost" and found: evolution of mitochondrial gene rearrangement in Antarctic notothenioids

Evolution of Antarctic notothenioids in the frigid and oxygen-rich Southern Ocean had led to remarkable genomic changes, most notably the gain of novel antifreeze glycoproteins and the loss of oxygen-binding hemoproteins in the icefish family. Recently, the mitochondrial (mt) NADH dehydrogenase subunit 6 (ND6) gene and the adjacent transfer RNA(Glu) (tRNA(Glu)) were also reportedly lost. ND6 protein is crucial for the assembly and function of Complex I of the mt electron transport chain that produces adenosine triphosphate (ATP) essential for life; thus, ND6 absence would be irreconcilable with Antarctic notothenioids being thriving species. Here we report our discovery that the ND6 gene and tRNA(Glu) were not lost but had been translocated to the control region (CR) from their canonical location between ND5 and cytochrome b genes. We characterized the CR and adjacent sequences of 22 notothenioid species representing all eight families of Notothenioidei to elucidate the mechanism and evolutionary history of this mtDNA rearrangement. Species of the three basal non-Antarctic families have the canonical vertebrate mt gene order, whereas species of all five Antarctic families have a rearranged CR bearing the embedded ND6 (ND6(CR)) and tRNA(Glu), with additional copies of tRNA(Thr), tRNA(Pro), and noncoding region in various lineages. We hypothesized that an initial duplication of the canonical mt region from ND6 through CR occurred in the common ancestor to the Antarctic clade, and we deduced the succession of loss or modification of the duplicated region leading to the extant patterns of mt DNA reorganization that is consistent with notothenioid evolutionary history. We verified that the ND6(CR) gene in Antarctic notothenioids is transcribed and therefore functional. However, ND6(CR)-encoded protein sequences differ substantially from basal non-Antarctic notothenioid ND6, and we detected lineage-specific positive selection on the branch leading to the Antarctic clade of ND6(CR) under the branch-site model. Collectively, the novel mt ND6(CR) genotype of the Antarctic radiation represents another major molecular change in Antarctic notothenioid evolution and may reflect an adaptive change conducive to the functioning of the protein (Complex I) machinery of mt respiration in the polar environment, driven by the advent of freezing, oxygen-rich conditions in the Southern Ocean.

Population genetic structure of long-tailed pygmy rice rats (Oligoryzomys longicaudatus) from Argentina and Chile based on the mitochondrial control region

The rodent Oligoryzomys longicaudatus (Bennett, 1832) (Rodentia, Cricetidae) inhabits southern forests of Argentina and Chile, a region severely affected by glaciations during the Pleistocene–Holocene periods. We evaluate here the diversity of the mitochondrial control region to characterize the genetic structure of this species from forests and bushy areas of seven populations from Argentina and four populations from Chile. Statistical analyses showed shallow haplotype trees and mismatch distributions compatible with recent range expansions. The presence of “private” haplotypes indicates that current levels of gene flow among populations of each country would be low to moderate. Significant differences in haplotype frequencies were detected between eastern and western populations, indicating that the Andes mountains would be an effective geographic barrier for gene flow despite the existing valleys that could act as corridors for dispersion. A single clade containing all the haplotypes was recovered in the phylogenetic trees, suggesting postglacial dispersion from a single refugium during the Last Glacial Maximum. The higher effective size and levels of polymorphism in populations from Chile suggest that the refugium was located in this country. The asymmetric gene flow from Chile to Argentina may reflect a recent colonization of the eastern populations.

Complete nucleotide mtDNA sequence of Barbary sheep (Ammotragus lervia)

In this report we describe the complete sequence of the mtDNA genome of Ammotragus lervia (Barbary sheep or aoudad) as obtained by PCR and sequencing with primer walking using flanking sequences. The molecule is 16,530 base pairs in length, resulting similar to those of goat and sheep. The genome organization matches to those of other mammalian mitochondrial genomes. The phylogenetic position of the Ammotragus relative to twelve other mammalian species was assessed based on protein-coding sequences. A goat-aoudad split of 9.9-8.9 MYBP has been estimated, whereas the divergence ancestor/caprines was dated at 12.1-10.8 MYBP.

Stability and Association with the Cytomatrix of Mitochondrial DNA in Spontaneously Immortalized Mouse Embryo Fibroblasts Containing or Lacking the Intermediate Filament Protein Vimentin

To extend previous observations demonstrating differences in number, morphology, and activity of mitochondria in spontaneously immortalized vim(+) and vim(-) fibroblasts derived from wild-type and vimentin knockout mice, some structural and functional aspects of mitochondrial genome performance and integrity in both types of cells were investigated. Primary Vim(+/+) and Vim(-/-) fibroblasts, which escaped terminal differentiation by immortalization were characterized by an almost twofold lower mtDNA content in comparison to that of their primary precursor cells, whereby the average mtDNA copy number in two clones of vim(+) cells was lower by a factor of 0.6 than that in four clones of vim(-) cells. However, during serial subcultivation up to high passage numbers, the vim(+) and vim() fibroblasts increased their mtDNA copy number 1.5- and 2.5-fold, respectively. While early-passage cells of the vim(+) and vim(-) fibroblast clones differed only slightly in the ratio between mtDNA content and mitochondrial mass represented by mtHSP70 protein, after ca. 300 population doublings the average mtDNA/mtmass ratio in the vim(+) and vim() cells was increased by a factor of 2 and 4.5, respectively. During subcultivation, both types of cells acquired the fully transformed phenotype. These findings suggest that cytoskeletal vimentin filaments exert a strong influence on the mechanisms controlling mtDNA copy number during serial subcultivation of immortalized mouse embryo fibroblasts, and that vimentin deficiency causes a disproportionately enhanced mtDNA content in high-passage vim(-) fibroblasts. Such a role of vimentin filaments was supported by the stronger retention potential for mtDNA and mtDNA polymerase (gamma) detected in vim(+) fibroblasts by Triton X-100 extraction of mitochondria and agaroseembedded cells. Moreover, although the vim(+) and vim(-) fibroblasts were equally active in generating free radicals, the vim(-) cells exhibited higher levels of immunologically detectable 8-oxoG and mismatch repair proteins MSH2 and MLH1 in their mitochondria. Because in vim(-) fibroblasts only one point mutation was detected in the mtDNA D-loop control region, these cells are apparently able to efficiently remove oxidatively damaged nucleobases. On the other hand, a number of large-scale mtDNA deletions were found in high-passage vim(-) fibroblasts, but not in low-passage vim(-) cells and vim(+) cells of both low and high passage. Large mtDNA deletions were also induced in young vim(-) fibroblasts by treatment with the DNA intercalator ethidium bromide, whereas no such deletions were found after treatment of vim(+) cells. These results indicate that in immortalized vim(-) fibroblasts the mitochondrial genome is prone to large-scale rearrangements, probably due to insufficient control of mtDNA repair and recombination processes in the absence of vimentin.

Structure of the mitochondrial control region of the Eurasian otter (Lutra lutra; Carnivora, Mustelidae): patterns of genetic heterogeneity and implications for conservation of the species in Italy

In this study we determined the complete sequence of the mitochondrial DNA (mtDNA) control region of the Eurasian otter (Lutra lutra). We then compared these new sequences with orthologues of nine carnivores belonging to six families (Mustelidae, Mephitidae, Canidae, Hyaenidae, Ursidae, and Felidae). The comparative analyses identified all the conserved regions previously found in mammals. The Eurasian otter and seven other species have a single location with tandem repeats in the right domain, while the spotted hyena (Hyaenidae) and the tiger (Felidae) have repeated sequences in both the right and left domains. To assess the degree of genetic heterogeneity of the Eurasian otter in Italy we sequenced two fragments of the gene and analyzed length polymorphisms of repeated sequences and heteroplasmy in 32 specimens. The study includes 23 museum specimens collected in northern, central, and southern Italy; most of these specimens are from extinct populations, while the southern Italian samples belong to the sole extant Italian population of the Eurasian otter. The study also includes all the captive-reared animals living in the colony "Centro Lontra, Caramanico Terme" (Pescara, central Italy). The colony is maintained for reintroduction of the species. We found a low level of genetic polymorphism; a single haplotype is dominant, but our data indicate the presence in central and southern Italy of two slightly divergent haplotypes. One haplotype belongs to an extinct population, the other is present in the single extant Italian population. Analyses of length polymorphisms and heteroplasmy indicate that the autochthonous Italian samples are characterized by a distinct array of repeated sequences from captive-reared animals.

Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes

Analysis of mitochondrial DNA sequence variation has been used extensively to study the evolutionary relationships of individuals and populations, both within and across species. So ubiquitous and easily acquired are mtDNA data that it has been suggested that such data could serve as a taxonomic 'barcode' for an objective species classification scheme. However, there are technical pitfalls associated with the acquisition of mtDNA data. One problem is the presence of translocated pieces of mtDNA in the nuclear genome of many taxa that may be mistaken for authentic organellar mtDNA. We assessed the extent to which such 'numt' sequences may pose an overlooked problem in analyses of mtDNA from humans and apes. Using long-range polymerase chain reaction (PCR), we generated necessarily authentic mtDNA sequences for comparison with sequences obtained using typical methods for a segment of the mtDNA control region in humans, chimpanzees, bonobos, gorillas and orangutans. Results revealed that gorillas are notable for having such a variety of numt sequences bearing high similarity to authentic mtDNA that any analysis of mtDNA using standard approaches is rendered impossible. Studies on humans, chimpanzees, bonobos or orangutans are apparently less problematic. One implication is that explicit measures need to be taken to authenticate mtDNA sequences in newly studied taxa or when any irregularities arise. Furthermore, some taxa may not be amenable to analysis of mtDNA variation at all.

Are the native giant tortoises from the Seychelles really extinct? A genetic perspective based on mtDNA and microsatellite data

The extinction of the giant tortoises of the Seychelles Archipelago has long been suspected but is not beyond doubt. A recent morphological study of the giant tortoises of the western Indian Ocean concluded that specimens of two native Seychelles species survive in captivity today alongside giant tortoises of Aldabra, which are numerous in zoos as well as in the wild. This claim has been controversial because some of the morphological characters used to identify these species, several measures of carapace morphology, are reputed to be quite sensitive to captive conditions. Nonetheless, the potential survival of giant tortoise species previously thought extinct presents an exciting scenario for conservation. We used mitochondrial DNA sequences and nuclear microsatellites to examine the validity of the rediscovered species of Seychelles giant tortoises. Our results indicate that the morphotypes suspected to represent Seychelles species do not show levels of variation and genetic structuring consistent with long periods of reproductive isolation. We found no variation in the mitochondrial control region among 55 individuals examined and no genetic structuring in eight microsatellite loci, pointing to the survival of just a single lineage of Indian Ocean tortoises.

Patterns of linkage disequilibrium in mitochondrial DNA of 16 ruminant populations

Mitochondrial DNA (mtDNA) is a widely employed molecular tool in phylogeography, in the inference of human evolutionary history, in dating the domestication of livestock and in forensic science. In humans and other vertebrates the popularity of mtDNA can be partially attributed to an assumption of strict maternal inheritance, such that there is no recombination between mitochondrial lineages. The recent demonstration that linkage disequilibrium (LD) declines as a function of distance between polymorphic sites in hominid mitochondrial genomes has been interpreted as evidence of recombination between mtDNA haplotypes, and hence nonclonal inheritance. However, critics of mtDNA recombination have suggested that this association is an artefact of an inappropriate measure of LD or of sequencing error, and subsequent studies of other populations have failed to replicate the initial finding. Here we report the analysis of 16 ruminant populations and present evidence that LD significantly declines with distance in five of them. A meta-analysis of the data indicates a nonsignificant trend of LD declining with distance. Most of the earlier criticisms of patterns between LD and distance in hominid mtDNA are not applicable to this data set. Our results suggest that either ruminant mtDNA is not strictly clonal or that compensatory selection has influenced patterns of variation at closely linked sites within the mitochondrial control region. The potential impact of these processes should be considered when using mtDNA as a tool in vertebrate population genetic, phylogenetic and forensic studies.

The crocodilian mitochondrial control region: general structure, conserved sequences, and evolutionary implications

We present the first comprehensive analysis of the crocodilian control region. We have analyzed sequences from all three families of Crocodylia (Crocodylidae, Gavialidae, Alligatoridae), incorporating all genera except Paleosuchus and Melanosuchus. Within the control region of other vertebrates, several sequence motifs and their order appear to be conserved. Herein, we compare aligned crocodilian D-loop sequences to homologous sequences from other vertebrates ranging from fish to birds. Among other findings, we have discovered that while domain I tends to be shorter than the same region in mammals and birds, it contains sequences similar in structure to both the goose-hairpin and termination associated sequences (TAS). Domain II is highly conservative with regard to size among the taxa examined and contains several of the conserved sequence boxes characterized in other vertebrates. Domain III contains several interesting sequence motifs including tandemly repeated sequences, a long poly-A region in the Crocodylidae, and possible bidirection promoter sequences.

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.

Structure and patterns of sequence variation in the mitochondrial DNA control region of the great cats

Mitochondrial DNA control region structure and variation were determined in the five species of the genus Panthera. Comparative analyses revealed two hypervariable segments, a central conserved region, and the occurrence of size and sequence heteroplasmy. As observed in the domestic cat, but not commonly seen in other animals, two repetitive sequence arrays (RS-2 with an 80-bp motif and RS-3 with a 6-10-bp motif) were identified. The 3' ends of RS-2 and RS-3 were highly conserved among species, suggesting that these motifs have different functional constraints. Control region sequences provided improved phylogenetic resolution grouping the sister taxa lion (Panthera leo) and leopard (Panthera pardus), with the jaguar (Panthera onca).

Phylogeography and population structure of an ecotonal marsupial, Bettongia tropica, determined using mtDNA and microsatellites

The northern bettong, Bettongia tropica, is an endangered species of Potoroidae with a restricted distribution in the wet tropics of north Queensland, Australia. The species is only found within a thin strip of sclerophyll forest along the western margin of rainforest. This tight association with rainforest boundaries is predicted to have resulted in population isolation as rainforest contracted during the Pleistocene, though some have proposed that the northern bettong was not present in the wet tropics until the late Pleistocene. The dispersal ability of the species, and of the family, is not known. This study examined gene flow among populations within areas of continuous habitat complemented by a broader analysis of phylogeography. Individuals trapped at each of the four known regions (one region was subsampled at three different sites), were sequenced for 547 base pairs of the mitochondrial DNA (mtDNA) control region and typed for seven microsatellite loci. The mtDNA phylogeny showed congruence with a biogeographical hypothesis, a relatively deep split suggesting historical isolation in separate northern and southern refugia. The two divergent clades were both present within the Lamb Range, indicating an expansion from these refuges and subsequent admixture at one site. mtDNA allele frequencies indicated relatively limited gene flow within the Lamb Range over distances as short as nine km. Tests of population divergence using microsatellites (FST and assignment tests) strongly supported this result. A molecular signal indicative of a recent bottleneck was unexpectedly detected in one of the Lamb Range subpopulations. This lead us to examine the behaviour of the statistics used in this bottleneck test under a linear stepping-stone model with varying migration rates. We found that it may be more difficult to detect molecular signatures for recent bottlenecks under conditions of very low migration rates than for isolated populations and, conversely, that 'false' bottleneck signatures may be observed at higher migration rates. The Lamb Range FST estimate clearly fell within the category of potentially 'false' bottleneck signals. Despite relatively limited gene flow, evidence for asymmetric dispersal suggests more complicated population dynamics than a simple linear stepping-stone model.

Slow rate of evolution in the mitochondrial control region of gulls (Aves: Laridae)

We sequenced part of the mitochondrial control region and the cytochrome b gene in 72 specimens from 32 gull species (Laridae, Larini) and 2 outgroup representatives (terns: Laridae, Sternini). Our control region segment spanned the conserved central domain II and the usually hypervariable 3' domain III. Apart from some heteroplasmy at the 3' end of the control region, domain III was not more variable than domain II or the cytochrome b gene. Furthermore, variation in the tempo of evolution of domain III was apparent between phyletic species groups. The lack of variation of the gull control region could not be explained by an increase in the proportion of conserved sequences in these birds, and the gull control region showed an organization similar to those of other avian control regions studied to date. A novel invariant direct repeat was identified in domain II of gulls, and in domain III, two to three inverted, sometimes imperfect, repeats are able to form a significantly stable stem-and-loop structure. These putative secondary structures have not been reported before, and a comparison between species groups showed that they are more stable in the group with the more conserved control region. The unusually slow rate of evolution of control region part III of the gulls could thus be partly explained by the existence of secondary structures in domain III of these species.

Phylogeographical population structure of tiger quolls Dasyurus maculatus (Dasyuridae: Marsupialia), an endangered carnivorous marsupial

Tiger quolls, Dasyurus maculatus, are the largest carnivorous marsupials still extant on the mainland of Australia, and occupy an important ecological niche as top predators and scavengers. Two allopatric subspecies are recognized, D.m. gracilis in north Queensland, and D.m. maculatus in the southeast of the mainland and Tasmania. D.m. gracilis is considered endangered while D.m. maculatus is listed as vulnerable to extinction; both subspecies are still in decline. Phylogeographical subdivision was examined to determine evolutionarily significant units (ESUs) and management units (MUs) among populations of tiger quolls to assist in the conservation of these taxa. Ninety-three tiger quolls from nine representative populations were sampled from throughout the species range. Six nuclear microsatellite loci and the mitochondrial DNA (mtDNA) control region (471 bp) were used to examine ESUs and MUs in this species. We demonstrated that Tasmanian tiger quolls are reciprocally monophyletic to those from the mainland using mtDNA analysis, but D.m. gracilis was not monophyletic with respect to mainland D.m. maculatus. Analysis of microsatellite loci also revealed significant differences between the Tasmanian and mainland tiger quolls, and between D.m. gracilis and mainland D.m. maculatus. These results indicate that Tasmanian and mainland tiger quolls form two distinct evolutionary units but that D.m. gracilis and mainland D.m. maculatus are different MUs within the same ESU. The two marker types used in this study revealed different male and female dispersal patterns and indicate that the most appropriate units for short-term management are local populations. A revised classification and management plan are needed for tiger quolls, particularly in relation to conservation of the Tasmanian and Queensland populations.

Phylogeographic differentiation in the mitochondrial control region in the koala, Phascolarctos cinereus (Goldfuss 1817)

The koala, Phascolarctos cinereus, is a geographically widespread species endemic to Australia, with three currently recognized subspecies: P.c. adustus, P.c. cinereus, and P.c. victor. Intraspecific variation in the mitochondrial DNA (mtDNA) control region was examined in over 200 animals from 16 representative populations throughout the species' range. Eighteen different haplotypes were defined in the approximately 860 bp mtDNA control region, as determined by heteroduplex analysis/temperature gradient gel electrophoresis (HDA/TGGE). Any single population typically possessed only one or two haplotypes yielding an average within-population haplotypic diversity of 0.180 +/- 0.003, and nucleotide diversity of 0.16%. Overall, mtDNA control region sequence diversity between populations averaged 0.67%, and ranged from 0% to 1.56%. Nucleotide divergence between populations averaged 0.51%, and ranged from 0% to 1.53%. Neighbour-joining methods revealed limited phylogenetic distinction between geographically distant populations of koalas, and tentative support for a single evolutionarily significant unit (ESU). This is consistent with previous suggestions that the morphological differences formalized by subspecific taxonomy may be interpreted as clinal variation. Significant differentiation in mtDNA-haplotype frequencies between localities suggested that little gene flow currently exists among populations. When combined with microsatellite analysis, which has revealed substantial differentiation among koala populations, we conclude that the appropriate short-term management unit (MU) for koalas is the local population.

A mitochondrial control region and cytochrome b phylogeny of sika deer (Cervus nippon) and report of tandem repeats in the control region

Sika deer (Cervus nippon Temminck) are endemic to mainland and insular Asia. Numerous subspecies have been named, but they are not quantitatively well defined. Portions of the mitochondrial cytochrome b gene (450 bp) and control region (512 bp) were sequenced from 28 individuals belonging to five sika subspecies and two Cervus elaphus subspecies. Phylogenetic trees constructed using these sequences clearly demonstrated that sika are monophyletic with respect to C. elaphus. A survey of variation in the control region showed that approximately half the variation occurred in a 100-base segment between positions 150 and 250 in the left domain of the control region. Within this region there were three tandemly repeated copies of a 39-base motif. In addition, two of the samples (C. n. aplodontus and C. n. hortulorum) contained, respectively, two and four additional copies of the repeated motif.

A high observed substitution rate in the human mitochondrial DNA control region

The rate and pattern of sequence substitutions in the mitochondrial DNA (mtDNA) control region (CR) is of central importance to studies of human evolution and to forensic identity testing. Here, we report a direct measurement of the intergenerational substitution rate in the human CR. We compared DNA sequences of two CR hypervariable segments from close maternal relatives, from 134 independent mtDNA lineages spanning 327 generational events. Ten substitutions were observed, resulting in an empirical rate of 1/33 generations, or 2.5/site/Myr. This is roughly twenty-fold higher than estimates derived from phylogenetic analyses. This disparity cannot be accounted for simply by substitutions at mutational hot spots, suggesting additional factors that produce the discrepancy between very near-term and long-term apparent rates of sequence divergence. The data also indicate that extremely rapid segregation of CR sequence variants between generations is common in humans, with a very small mtDNA bottleneck. These results have implications for forensic applications and studies of human evolution.

The complete mitochondrial genome of the wallaroo (Macropus robustus) and the phylogenetic relationship among Monotremata, Marsupialia, and Eutheria

The complete mitochondrial DNA (mtDNA) (16,896 nt) of the wallaroo (Macropus robustus) was sequenced. The concatenated amino acid sequences of 12 mitochondrial protein-coding genes of the wallaroo plus those of a number of other mammals were included in a phylogenetic study of early mammalian divergences. The analysis joined monotremes and marsupials (the Marsupionta hypothesis) to the exclusion of eutherians. The analysis rejected significantly the commonly acknowledged Theria hypothesis, according to which Marsupialia and Eutheria are grouped together to the exclusion of Monotremata. The region harboring the gene for lysine tRNA (tRNA-Lys) in the mtDNA of other vertebrates is in the wallaroo occupied by a sequence (tRNA-Lys) that lacks both an anticodon loop as well as the anticodon for the amino acid lysine. An alternative tRNA-Lys gene was not identified in any other region of the mtDNA of the wallaroo, suggesting that a tRNA-Lys of nuclear origin is imported into marsupial mitochondria. Previously described RNA editing of tRNA-Asp and rearrangement of some tRNA genes were reconfirmed in the mtDNA of the wallaroo. The divergence between Monotremata/Marsupialia and Eutheria was timed to approximately 130 million years before present (MYBP). The same calculations suggested that Monotremata and Marsupialia diverged approximately 115 MYBP and that Australian and American marsupials separated approximately 75 MYBP. The findings also show that many, probably most, extant eutherian orders had their origin in middle to late Cretaceous times, 115-65 MYBP.