Nucleotide substitution rate of mammalian mitochondrial genomes

Mitonuclear match: optimizing fitness and fertility over generations drives ageing within generations

Many conserved eukaryotic traits, including apoptosis, two sexes, speciation and ageing, can be causally linked to a bioenergetic requirement for mitochondrial genes. Mitochondrial genes encode proteins involved in cell respiration, which interact closely with proteins encoded by nuclear genes. Functional respiration requires the coadaptation of mitochondrial and nuclear genes, despite divergent tempi and modes of evolution. Free-radical signals emerge directly from the biophysics of mosaic respiratory chains encoded by two genomes prone to mismatch, with apoptosis being the default penalty for compromised respiration. Selection for genomic matching is facilitated by two sexes, and optimizes fitness, adaptability and fertility in youth. Mismatches cause infertility, low fitness, hybrid breakdown, and potentially speciation. The dynamics of selection for mitonuclear function optimize fitness over generations, but the same selective processes also operate within generations, driving ageing and age-related diseases. This coherent view of eukaryotic energetics offers striking insights into infertility and age-related diseases.

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

Apoptosis induced by persistent single-strand breaks in mitochondrial genome: critical role of EXOG (5'-EXO/endonuclease) in their repair

Reactive oxygen species (ROS), continuously generated as by-products of respiration, inflict more damage on the mitochondrial (mt) than on the nuclear genome because of the nonchromatinized nature and proximity to the ROS source of the mitochondrial genome. Such damage, particularly single-strand breaks (SSBs) with 5'-blocking deoxyribose products generated directly or as repair intermediates for oxidized bases, is repaired via the base excision/SSB repair pathway in both nuclear and mt genomes. Here, we show that EXOG, a 5'-exo/endonuclease and unique to the mitochondria unlike FEN1 or DNA2, which, like EXOG, has been implicated in the removal of the 5'-blocking residue, is required for repairing endogenous SSBs in the mt genome. EXOG depletion induces persistent SSBs in the mtDNA, enhances ROS levels, and causes apoptosis in normal cells but not in mt genome-deficient rho0 cells. Thus, these data show for the first time that persistent SSBs in the mt genome alone could provide the initial trigger for apoptotic signaling in mammalian cells.

DNA marker technology for wildlife conservation

Use of molecular markers for identification of protected species offers a greater promise in the field of conservation biology. The information on genetic diversity of wildlife is necessary to ascertain the genetically deteriorated populations so that better management plans can be established for their conservation. Accurate classification of these threatened species allows understanding of the species biology and identification of distinct populations that should be managed with utmost care. Molecular markers are versatile tools for identification of populations with genetic crisis by comparing genetic diversities that in turn helps to resolve taxonomic uncertainties and to establish management units within species. The genetic marker analysis also provides sensitive and useful tools for prevention of illegal hunting and poaching and for more effective implementation of the laws for protection of the endangered species. This review summarizes various tools of DNA markers technology for application in molecular diversity analysis with special emphasis on wildlife conservation.

Evolutionary patterns of the mitochondrial genome in Metazoa: exploring the role of mutation and selection in mitochondrial protein coding genes

The mitochondrial genome is a fundamental component of the eukaryotic domain of life, encoding for several important subunits of the respiratory chain, the main energy production system in cells. The processes by means of which mitochondrial DNA (mtDNA) replicates, expresses itself and evolves have been explored over the years, although various aspects are still debated. In this review, we present several key points in modern research on the role of evolutionary forces in affecting mitochondrial genomes in Metazoa. In particular, we assemble the main data on their evolution, describing the contributions of mutational pressure, purifying, and adaptive selection, and how they are related. We also provide data on the evolutionary fate of the mitochondrial synonymous variation, related to the nonsynonymous variation, in comparison with the pattern detected in the nucleus.

Elevated mutational pressure characterizes the evolution of the mitochondrial synonymous variation, whereas purging selection, physiologically due to phenomena such as cell atresia and intracellular mtDNA selection, guarantees coding sequence functionality. This enables mitochondrial adaptive mutations to emerge and fix in the population, promoting mitonuclear coevolution.

"Patchy-tachy" leads to false positives for recombination

Indirect tests have detected recombination in mitochondrial DNA (mtDNA) from many animal lineages, including mammals. However, it is possible that features of the molecular evolutionary process without recombination could be incorrectly inferred by indirect tests as being due to recombination. We have identified one such example, which we call "patchy-tachy" (PT), where different partitions of sequences evolve at different rates, that leads to an excess of false positives for recombination inferred by indirect tests. To explore this phenomena, we characterized the false positive rates of six widely used indirect tests for recombination using simulations of general models for mtDNA evolution with PT but without recombination. All tests produced 30-99% false positives for recombination, although the conditions that produced the maximal level of false positives differed between the tests. To evaluate the degree to which conditions that exacerbate false positives are found in published sequence data, we turned to 20 animal mtDNA data sets in which recombination is suggested by indirect tests. Using a model where different regions of the sequences were free to evolve at different rates in different lineages, we demonstrated that PT is prevalent in many data sets in which recombination was previously inferred using indirect tests. Taken together, our results argue that PT without recombination is a viable alternative explanation for detection of widespread recombination in animal mtDNA using indirect tests.

Association of mtDNA D-loop polymorphisms with risk of gastric cancer in Chinese population

The aim of present study was to evaluate the association of common polymorphisms detected in mitochondrial DNA (mtDNA) D-loop region (mononucleotide repetitive D310, single nucleotide polymorphism (SNP) D16521) with susceptibility to gastric cancer (GC) in northwestern Chinese population. A total of 180 GC patients and 218 healthy controls were investigated by using PCR- denaturing high performance liquid chromatography (DHPLC) assay. Genotype and allele distributions and haplotype construction were analyzed in case-control study. We found D310 and D16521 heteroplasmy were significantly different between GC cases and controls (p < 0.05), and D16521 homoplasmy showed association with histological grade of GC (p < 0.05). Haplotype 7C/T, 8C/C and 9C/C had significant association with GC risk implied from analysis of D310 and D16521. Taken together, these findings suggested that mtDNA D-Loop polymorphisms and haplotypes may contribute to genetic susceptibility to GC in Chinese population.

The large-scale evolution by generating new genes from gene duplication; similarity and difference between monoploid and diploid organisms

On the basis of the concept of biological activity, the large-scale evolution by generating new genes from gene duplication is theoretically compared between the monoploid organism and the diploid organism. The comparison is carried out not only for the process of generating one new gene but also for the process of generating two or more kinds of new genes from successive gene duplication. This comparison reveals the following difference in evolutionary pattern between the monoploids and diploids. The monoploid organism is more suitable to generate one or two new genes step by step but its successive gene duplication is obliged to generate smaller sizes of genes by the severer lowering of biological activity or self-reproducing rate. This is consistent with the evolutionary pattern of prokaryotes having steadily developed chemical syntheses, O₂-releasing photosynthesis and O₂-respiration in the respective lineages. On the other hand, the diploid organism with the plural number of homologous chromosome pairs has a chance to get together many kinds of new genes by the hybridization of variants having experienced different origins of gene duplication. Although this strategy of hybridization avoids the severe lowering of biological activity, it takes the longer time to establish the homozygotes of the more kinds of new genes. During this long period, furthermore different types of variants are accumulated in the population, and their successive hybridization sometimes yields various styles of new organisms. This evolutionary pattern explains the explosive divergence of body plans that has occasionally occurred in the diploid organisms, because the cell differentiation is a representative character exhibited by many kinds of genes and its evolution to the higher hierarchy constructs body plans.

Impairment of mitochondrial tRNAIle processing by a novel mutation associated with chronic progressive external ophthalmoplegia

We report a sporadic case of chronic progressive external ophthalmoplegia associated with ragged red fibers. The patient presented with enlarged mitochondria with deranged internal architecture and crystalline inclusions. Biochemical studies showed reduced activities of complex I, III and IV in skeletal muscle. Molecular genetic analysis of all mitochondrial tRNAs revealed a G to A transition at nt 4308; the G is a highly conserved nucleotide that participates in a GC base-pair in the T-stem of mammalian mitochondrial tRNA(Ile). The mutation was detected at a high level (approx. 50%) in muscle but not in blood. The mutation co-segregated with the phenotype, as the mutation was absent from blood and muscle in the patient's healthy mother. Functional characterization of the mutation revealed a six-fold reduced rate of tRNA(Ile) precursor 3' end maturation in vitro by tRNAse Z. Furthermore, the mutated tRNA(Ile) displays local structural differences from wild-type. These results suggest that structural perturbations reduce efficiency of tRNA(Ile) precursor 3' end processing and contribute to the molecular pathomechanism of this mutation.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.

The complete mitochondrial genome of dhole Cuon alpinus: phylogenetic analysis and dating evolutionary divergence within Canidae

The dhole (Cuon alpinus) is the only existent species in the genus Cuon (Carnivora: Canidae). In the present study, the complete mitochondrial genome of the dhole was sequenced. The total length is 16672 base pairs which is the shortest in Canidae. Sequence analysis revealed that most mitochondrial genomic functional regions were highly consistent among canid animals except the CSB domain of the control region. The difference in length among the Canidae mitochondrial genome sequences is mainly due to the number of short segments of tandem repeated in the CSB domain. Phylogenetic analysis was progressed based on the concatenated data set of 14 mitochondrial genes of 8 canid animals by using maximum parsimony (MP), maximum likelihood (ML) and Bayesian (BI) inference methods. The genera Vulpes and Nyctereutes formed a sister group and split first within Canidae, followed by that in the Cuon. The divergence in the genus Canis was the latest. The divarication of domestic dogs after that of the Canis lupus laniger is completely supported by all the three topologies. Pairwise sequence divergence data of different mitochondrial genes among canid animals were also determined. Except for the synonymous substitutions in protein-coding genes, the control region exhibits the highest sequence divergences. The synonymous rates are approximately two to six times higher than those of the non-synonymous sites except for a slightly higher rate in the non-synonymous substitution between Cuon alpinus and Vulpes vulpes. 16S rRNA genes have a slightly faster sequence divergence than 12S rRNA and tRNA genes. Based on nucleotide substitutions of tRNA genes and rRNA genes, the times since divergence between dhole and other canid animals, and between domestic dogs and three subspecies of wolves were evaluated. The result indicates that Vulpes and Nyctereutes have a close phylogenetic relationship and the divergence of Nyctereutes is a little earlier. The Tibetan wolf may be an archaic pedigree within wolf subspecies. The genetic distance between wolves and domestic dogs is less than that among different subspecies of wolves. The domestication of dogs was about 1.56-1.92 million years ago or even earlier.

Phylogeny and biogeography of western Indian Ocean Rousettus (Chiroptera: Pteropodidae)

We examined patterns of genetic variation in Rousettus madagascariensis from Madagascar and R. obliviosus from the Comoros (Grande Comore, Anjouan, and Mohéli). Genetic distances among individuals on the basis of 1,130 base pairs of the mitochondrial cytochrome b (Cytb) locus were estimated from specimens collected from 17 sites on Madagascar, 3 sites on Grande Comore, 3 sites on Anjouan, and 2 sites on Mohéli. We observed little variation in Madagascar and nearshore island samples (maximum 1.1%) and interisland Comoros samples (maximum 1.8%). In contrast, pairwise distances between different sampled sites on Madagascar and the Comoros varied from 8.5% to 13.2%. For 131 Malagasy animals, 69 unique haplotypes were recovered with 86 variable sites, and for 44 Comorian individuals, 17 unique haplotypes were found with 30 variable sites. No haplotype was shared between Madagascar and the Comoros, adding to previous morphological evidence that these 2 populations should be considered separate species. Cytb data showed that Rousettus populations of Madagascar (including nearshore islands) and the Comoros are respectively monophyletic and display no geographic structure in haplotype diversity, and that R. madagascariensis and R. obliviosus are strongly supported as sister to each other relative to other Rousettus species. Genotypic data from 6 microsatellite loci confirm lack of geographic structure in either of the 2 species. In pairwise tests of population differentiation, the only significant values were between samples from the Comoro Islands and Madagascar (including nearshore islands). Estimates of current and historical demographic parameters support population expansion in both the Comoros and Madagascar. These data suggest a more recent and rapid demographic expansion in Madagascar in comparison with greater population stability on the Comoros. On the basis of available evidence, open-water crossings approaching 300 km seem rarely traversed by Rousettus, and, if successful, can result in genetic isolation and subsequent differentiation.

Conservation genetics of a rare Gerbil species: a comparison of the population genetic structures and demographic histories of the locally rare Pygmy Gerbil and the common Anderson's Gerbil

Background

One of the major challenges in evolutionary biology is identifying rare species and devising management plans to protect them while also sustaining their genetic diversity. However, in attempting a broad understanding of rarity, single-species studies provide limited insights because they do not reveal whether the factors that affect rare species differ from those that affect more common species. To illustrate this important concept and to arrive at a better understanding of the form of rarity characterizing the rare Gerbillus henleyi, we explored its population genetic structure alongside that of the locally common Gerbillus andersoni allenbyi. We trapped gerbils in several locations in Israel's western and inner Negev sand dunes. We then extracted DNA from ear samples, and amplified two mitochondrial sequences: the control region (CR) and the cytochrome oxidase 2 gene (CO2).

Results

Nucleotide diversity was low for all sequences, especially for the CR of G. a. allenbyi, which showed no diversity. We could not detect any significant population genetic structure in G. henleyi. In contrast, G. a. allenbyi's CO2 sequence showed significant population genetic structure. Pairwise PhiPT comparisons showed low values for G. henleyi but high values for G. a. allenbyi. Analysis of the species' demographic history indicated that G. henleyi's population size has not changed recently, and is under the influence of an ongoing bottleneck. The same analysis for G. a. allenbyi showed that this species has undergone a recent population expansion.

Conclusions

Comparing the two species, the populations of G. a. allenbyi are more isolated from each other, likely due to the high habitat specificity characterizing this species. The bottleneck pattern found in G. henleyi may be the result of competition with larger gerbil species. This result, together with the broad habitat use and high turnover rate characterizing G. henleyi, may explain the low level of differentiation among its populations. The evidence for a recent population expansion of G. a. allenbyi fits well with known geomorphological data about the formation of the Negev sand dunes and paleontological data about this species' expansion throughout the Levant. In conclusion, we suggest that adopting a comparative approach as presented here can markedly improve our understanding of the causes and effects of rarity, which in turn can allow us to better protect biodiversity patterns.

Comparative analysis of mitochondrial control region in polyploid hybrids of red crucian carp (Carassius auratus) x blunt snout bream (Megalobrama amblycephala)

The entire sequences of the mitochondrial (mt)DNA control region (CR) and portions of its flanking genes in the red crucian carp (RC) and blunt snout bream (BSB) as well as their polyploid hybrids (3nRB, 4nRB and 5nRB) were determined and subjected to a comparative analysis. The mtDNA-CRs of these five fish species ranged from 923 to 937 bp in length, they had the same flanking gene arrangement as other vertebrates and the pattern of nucleotide substitution bias was also similar to that in other vertebrates. Our data are consistent with the viewpoint of three domains [extended terminal associated sequence (ETAS domain), central conserved sequence block domain and conserved sequence block (CSB) domain] within the mtDNA-CR of mammals. On the basis our comparative analysis of the mtDNA-CRs of these five fish species, we were able to identify the consensus sequences of functional conserved units, including the ETAS, CSB-F, CSB-D, CSB-E, CSB1, CSB2 and CSB3 and putative promoter. The percentage of variable nucleotide positions (41.98%) in the central domain was lower than those in the ETAS and conserved domain (71.70 and 47.12%, respectively), suggesting that the central domain was the most conserved part of the mtDNA-CR. These results provide useful and important information for the further study of mtDNA-CR structure in fish. The sequence similarities of mtDNA-CR among the 3nRB, 4nRB, 5nRB hybrids and their respective female parents were higher than those among the 3nRB, 4nRB, 5nRB hybrids and their respective male parents, providing the direct evidence of stringent maternal inheritance of mtDNA-CR in the 3nRB, 4nRB and 5nRB hybrids.

Evaluating nested clade phylogeographic analysis under models of restricted gene flow

Nested clade phylogeographic analysis (NCPA) is a widely used method that aims to identify past demographic events that have shaped the history of a population. In an earlier study, NCPA has been fully automated, allowing it to be tested with simulated data sets generated under a null model in which samples simulated from a panmictic population are geographically distributed. It was noted that NCPA was prone to inferring false positives, corroborating earlier findings. The present study aims to evaluate both single-locus and multilocus NCPA under the scenario of restricted gene flow among spatially distributed populations. We have developed a new program, ANeCA-ML, which implements multilocus NCPA. Data were simulated under 3 models of gene flow: a stepping stone model, an island model, and a stepping stone model with some long-distance dispersal. Results indicate that single-locus NCPA tends to give a high frequency of false positives, but, unlike the random-mating scenario presented previously, inferences are not limited to restricted gene flow with isolation by distance or contiguous range expansion. The proportion of single-locus data sets that contained false inferences was 76% for the panmictic case, 87% for the stepping stone model, 79% for the stepping stone model with long-distance dispersal, and more than 99% for the island model. The frequency of inferences is inversely related to the amount of gene flow between demes. We performed multilocus NCPA by grouping the simulated loci into data sets of 5 loci. The false-positive rate was reduced in multilocus NCPA for some inferences but remained high for others. The proportion of multilocus data sets that contained false inferences was 17% for the panmictic case, 30% for the stepping stone model, 4% for the stepping stone model with long-distance dispersal, and 54% for the island model. Multilocus NCPA reduces the false-positive rate by restricting the sensitivity of the method but does not appear to increase the accuracy of the approach. Three classical tests-the analysis of molecular variance method, Fu's Fs, and the Mantel test-show that there is information in the data that gives rise to explicable results using these standard approaches. In conclusion, for the scenarios that we have examined, our simulation study suggests that the NCPA method is unreliable and its inferences may be misleading. We suggest that the NCPA method should not be used without objective simulation-based testing by independent researchers.

Mitochondrial phylogenomics of the Bivalvia (Mollusca): searching for the origin and mitogenomic correlates of doubly uniparental inheritance of mtDNA

Background

Doubly uniparental inheritance (DUI) is an atypical system of animal mtDNA inheritance found only in some bivalves. Under DUI, maternally (F genome) and paternally (M genome) transmitted mtDNAs yield two distinct gender-associated mtDNA lineages. The oldest distinct M and F genomes are found in freshwater mussels (order Unionoida). Comparative analyses of unionoid mitochondrial genomes and a robust phylogenetic framework are necessary to elucidate the origin, function and molecular evolutionary consequences of DUI. Herein, F and M genomes from three unionoid species, Venustaconcha ellipsiformis, Pyganodon grandis and Quadrula quadrula have been sequenced. Comparative genomic analyses were carried out on these six genomes along with two F and one M unionoid genomes from GenBank (F and M genomes of Inversidens japanensis and F genome of Lampsilis ornata).

Results

Compared to their unionoid F counterparts, the M genomes contain some unique features including a novel localization of the trnH gene, an inversion of the atp8-trnD genes and a unique 3'coding extension of the cytochrome c oxidase subunit II gene. One or more of these unique M genome features could be causally associated with paternal transmission. Unionoid bivalves are characterized by extreme intraspecific sequence divergences between gender-associated mtDNAs with an average of 50% for V. ellipsiformis, 50% for I. japanensis, 51% for P. grandis and 52% for Q. quadrula (uncorrected amino acid p-distances). Phylogenetic analyses of 12 protein-coding genes from 29 bivalve and five outgroup mt genomes robustly indicate bivalve monophyly and the following branching order within the autolamellibranch bivalves: ((Pteriomorphia, Veneroida) Unionoida).

Conclusion

The basal nature of the Unionoida within the autolamellibranch bivalves and the previously hypothesized single origin of DUI suggest that (1) DUI arose in the ancestral autolamellibranch bivalve lineage and was subsequently lost in multiple descendant lineages and (2) the mitochondrial genome characteristics observed in unionoid bivalves could more closely resemble the DUI ancestral condition. Descriptions and comparisons presented in this paper are fundamental to a more complete understanding regarding the origins and consequences of DUI.

Molecular systematics and global phylogeography of angel sharks (genus Squatina)

Angel sharks of the genus Squatina represent a group comprising 22 extant benthic species inhabiting continental shelves and upper slopes. In the present study, a comprehensive phylogenetic reconstruction of 17 Squatina species based on two mitochondrial markers (COI and 16S rRNA) is provided. The phylogenetic reconstructions are used to test biogeographic patterns. In addition, a molecular clock analysis is conducted to estimate divergence times of the emerged clades. All analyses show Squatina to be monophyletic. Four geographic clades are recognized, of which the Europe-North Africa-Asia clade is probably a result of the Tethys Sea closure. A second sister group relationship emerged in the analyses, including S. californica (eastern North Pacific) and S. dumeril (western North Atlantic), probably related to the rise of the Panamanian isthmus. The molecular clock analysis show that both lineage divergences coincide with the estimated time of these two geological events.

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.

Mitochondrial haplotype does not affect sperm velocity in the zebra finch Taeniopygia guttata

Mitochondrial DNA (mtDNA) variation has been suggested as a possible cause of variation in male fertility because sperm activity is tightly coupled to mitochondrial oxidative phosphorylation and ATP production, both of which are sensitive to mtDNA mutations. Since male-specific phenotypes such as sperm have no fitness consequences for mitochondria due to maternal mitochondrial (and mtDNA) inheritance, mtDNA mutations that are deleterious in males but which have negligible or no fitness effect in females can persist in populations. How often such mutations arise and persist is virtually unknown. To test whether there were associations between mtDNA variation and sperm performance, we haplotyped 250 zebra finches Taeniopygia guttata from a large pedigreed-population and measured sperm velocity using computer-assisted sperm analysis. Using quantitative genetic 'animal' models, we found no effect of mtDNA haplotype on sperm velocity. Therefore, there is no evidence that in this system mitochondrial mutations have asymmetric fitness effects on males and females, leading to genetic variation in male fertility that is blind to natural selection.

Using species-specific repeat and PCR-RFLP in typing of DNA derived from blood of human and animal species

Species determination of tissue specimens, including blood, is an important component of forensic analysis to distinguish human from animal remains. DNA markers based on a method of species-specific PCR and amplifying the 359-base pair (bp) fragment of the mitochondrially encoded cytochrome-b gene and then digestion with the TaqI restriction enzyme were developed for detection and discrimination of human, cattle, buffalo, horse, sheep, pig, dog, cat and chicken blood samples. The results reveal that PCR-amplification of the gene encoding the species-specific repeat (SSR) region generated 603 bp in cattle and buffalo, 221 bp in horse, 374 bp in sheep, <or=100 bp in pig, 808 bp in dog, 672 bp in cat and 50 bp in chicken. Restriction analysis of the amplified 359-bp portion of the cytochrome-b gene using the TaqI restriction enzyme results in species-specific restriction fragment length polymorphism (RFLP) between buffalo, cattle and human. Two different bands were generated in buffalo (191 and 168 bp) and human (209 and 150 bp), with no digestion in cattle (359 bp). Cytochrome-b is a highly conserved region and consequently a good molecular marker for diagnostic studies. Therefore, the two complementary techniques, SSR-PCR and PCR-RFLP, could be used successfully as routine methods in forensics for sensitive, rapid, simple and inexpensive identification of the species in bloodstains.

Identification of forensically informative SNPs in the domestic dog mitochondrial control region

Dog hair is often found at crime scenes either due to the dog's involvement in the crime or secondary transfer. As little nuclear DNA is present in shed hair, a 1000 base pair fragment of the mitochondrial control region (mtCR) from 552 dogs was assessed for forensically useful sequence variation. Through pairwise alignment to a standard reference sequence, existing haplotypes were further described and 36 new haplotypes and 24 new single nucleotide polymorphisms were identified. The probability of exclusion was found to be 0.957. Breeds were found to have similar sequences, although not identical. No genetic basis was found for grouping dogs by either purebred or mixed or geographic location within the continental United States. Our research demonstrates that the domestic dog mtCR has not been thoroughly surveyed for sequence variation and that a single database comprised of purebred and mixed breed dogs is sufficient for the continental United States.

Mitochondrial genome DNA analysis of the domestic dog: identifying informative SNPs outside of the control region

While the mitochondrial control region has proven successful for human forensic evaluations by indicating ethnic origin, domestic dogs (Canis lupus familiaris) of seemingly unrelated breeds often form large groups based on identical control region sequences. In an attempt to break up these large haplotype groups, we have analyzed the remaining c. 15,484 base pairs of the canine mitochondrial genome for 79 dogs and used phylogenetic and population genetic methods to search for additional variability in the form of single nucleotide polymorphisms (SNPs). We have identified 356 SNPs and 65 haplotypes in the remainder of the mitochondrial genome excluding the control region. The exclusion capacity was found to be 0.018. The mitochondrial control region was also evaluated for the same 79 dogs. The signals from the different fragments do not conflict, but instead support one another and provide a larger fragment of DNA that can be analyzed as forensic evidence.

Comparative genomics of Drosophila mtDNA: Novel features of conservation and change across functional domains and lineages

To gain insight on mitochondrial DNA (mtDNA) evolution, we assembled and analyzed the mitochondrial genomes of Drosophila erecta, D. ananassae, D. persimilis, D. willistoni, D. mojavensis, D. virilis and D. grimshawi together with the sequenced mtDNAs of the melanogaster subgroup. Genomic comparisons across the well-defined Drosophila phylogeny impart power for detecting conserved mtDNA regions that maintain metabolic function and regions that evolve uniquely on lineages. Evolutionary rate varies across intergenic regions of the mtDNA. Rapidly evolving intergenic regions harbor the majority of mitochondrial indel divergence. In contrast, patterns of nearly perfect conservation within intergenic regions reveal a refined set of nucleotides underlying the binding of transcription termination factors. Sequencing of 5' cDNA ends indicates that cytochrome C oxidase I (CoI) has a novel (T/C)CG start codon and that perfectly conserved regions upstream of two NADH dehydrogenase (ND) genes are transcribed and likely extend these protein sequences. Substitutions at synonymous sites in the Drosophila mitochondrial proteomes reflect a mutation process that is biased toward A and T nucleotides and differs between mtDNA strands. Differences in codon usage bias across genes reveal that weak selection at silent sites may offset the mutation bias. The mutation-selection balance at synonymous sites has also diverged between the Drosophila and Sophophora lineages. Rates of evolution are highly heterogeneous across the mitochondrial proteome, with ND accumulating many more amino acid substitutions than CO. These oxidative phosphorylation complex-specific rates of evolution vary across lineages and may reflect physiological and ecological change across the Drosophila phylogeny.

Identification of complete mitochondrial genome of the tufted deer

The tufted deer Elaphodus cephalophus are endangered animals in the world and little is understood about their mitochondrial (mt) genome. In our study, the mt genome of the tufted deer is identified--which is about 16 kb in length and contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a non-coding sequence (control region). The distinguishing feature is that GTG is the start codon of the NADH4L gene and the cyt b gene has a subterminal AAA followed by the stop codon TAG. According to 12 H strand protein-coding genes and phylogenetic analysis, Elaphodus may have a sister relationship with another deer group Muntiacus.

Do geological or climatic processes drive speciation in dynamic archipelagos? The tempo and mode of diversification in Southeast Asian shrews

Geological and climatic processes potentially alter speciation rates by generating and modifying barriers to dispersal. In Southeast Asia, two processes have substantially altered the distribution of land. Volcanic uplift produced many new islands during the Miocene-Pliocene and repeated sea level fluctuations during the Pleistocene resulted in intermittent land connections among islands. Each process represents a potential driver of diversification. We use a phylogenetic analysis of a group of Southeast Asian shrews (Crocidura) to examine geographic and temporal processes of diversification. In general, diversification has taken place in allopatry following the colonization of new areas. Sulawesi provides an exception, where we cannot reject within-island speciation for a clade of eight sympatric and syntopic species. We find only weak support for temporally declining diversification rates, implying that neither volcanic uplift nor sea level fluctuations had a strong effect on diversification rates. We suggest that dynamic archipelagos continually offer new opportunities for allopatric diversification, thereby sustaining high speciation rates over long periods of time, or Southeast Asian shrews represent an immature radiation on a density-dependent trajectory that has yet to fill geographic and ecological space.

Using classical population genetics tools with heterochroneous data: time matters!

BACKGROUND

New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations.

METHODOLOGY/PRINCIPAL FINDINGS

We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)).

CONCLUSIONS/SIGNIFICANCE

Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale.

Comparing evolutionary patterns and variability in the mitochondrial control region and cytochrome B in three species of baleen whales

The rapidly evolving mitochondrial control region remains an important source of information on phylogeography and demographic history for cetaceans and other vertebrates, despite great uncertainty in the rate of nucleotide substitution across both nucleotide positions and lineages. Patterns of variation in linked markers with slower rates of evolution can potentially be used to calibrate the rate of nucleotide substitution in the control region and to better understand the interplay of evolutionary and demographic forces across the mitochondrial genome above and below the species level. We have examined patterns of diversity within and between three baleen whale species (gray, humpback, and Antarctic minke whales) in order to determine how patterns of molecular evolution differ between cytochrome b and the control region. Our results show that cytochrome b is less variable than expected given the diversity in the control region for gray and humpback whales, even after functional differences are taken into account, but more variable than expected for minke whales. Differences in the frequency distributions of polymorphic sites and in best-fit models of nucleotide substitution indicate that these patterns may be the result of hypervariability in the control region in gray and humpback whales but, in minke whales, may result from a large, stable or expanding population size coupled with saturation at the control region. Using paired cytochrome b and control region data across individuals, we show that the average rate of nucleotide substitution in the control region may be on average 2.6 times higher than phylogenetically derived estimates in cetaceans. These results highlight the complexity of making inferences from control region data alone and suggest that applying simple rules of DNA sequence analyses across species may be difficult.

Reconsidering the generation time hypothesis based on nuclear ribosomal ITS sequence comparisons in annual and perennial angiosperms

Background

Differences in plant annual/perennial habit are hypothesized to cause a generation time effect on divergence rates. Previous studies that compared rates of divergence for internal transcribed spacer (ITS1 and ITS2) sequences of nuclear ribosomal DNA (nrDNA) in angiosperms have reached contradictory conclusions about whether differences in generation times (or other life history features) are associated with divergence rate heterogeneity. We compared annual/perennial ITS divergence rates using published sequence data, employing sampling criteria to control for possible artifacts that might obscure any actual rate variation caused by annual/perennial differences.

Results

Relative rate tests employing ITS sequences from 16 phylogenetically-independent annual/perennial species pairs rejected rate homogeneity in only a few comparisons, with annuals more frequently exhibiting faster substitution rates. Treating branch length differences categorically (annual faster or perennial faster regardless of magnitude) with a sign test often indicated an excess of annuals with faster substitution rates. Annuals showed an approximately 1.6-fold rate acceleration in nucleotide substitution models for ITS. Relative rates of three nuclear loci and two chloroplast regions for the annual Arabidopsis thaliana compared with two closely related Arabidopsis perennials indicated that divergence was faster for the annual. In contrast, A. thaliana ITS divergence rates were sometimes faster and sometimes slower than the perennial. In simulations, divergence rate differences of at least 3.5-fold were required to reject rate constancy in > 80 % of replicates using a nucleotide substitution model observed for the combination of ITS1 and ITS2. Simulations also showed that categorical treatment of branch length differences detected rate heterogeneity > 80% of the time with a 1.5-fold or greater rate difference.

Conclusion

Although rate homogeneity was not rejected in many comparisons, in cases of significant rate heterogeneity annuals frequently exhibited faster substitution rates. Our results suggest that annual taxa may exhibit a less than 2-fold rate acceleration at ITS. Since the rate difference is small and ITS lacks statistical power to reject rate homogeneity, further studies with greater power will be required to adequately test the hypothesis that annual and perennial plants have heterogeneous substitution rates. Arabidopsis sequence data suggest that relative rate tests based on multiple loci may be able to distinguish a weak acceleration in annual plants. The failure to detect rate heterogeneity with ITS in past studies may be largely a product of low statistical power.

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L. crocea individuals and two L. polyactis individuals. Interestingly, all seven CR sequences from L. crocea were identical while the three sequences from L. polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L. crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L. crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys.