Examining Phylogenetic Relationships Among Gibbon Genera Using Whole Genome Sequence Data Using an Approximate Bayesian Computation Approach

Gibbons are believed to have diverged from the larger great apes ∼16.8 MYA and today reside in the rainforests of Southeast Asia. Based on their diploid chromosome number, the family Hylobatidae is divided into four genera, Nomascus, Symphalangus, Hoolock, and Hylobates. Genetic studies attempting to elucidate the phylogenetic relationships among gibbons using karyotypes, mitochondrial DNA (mtDNA), the Y chromosome, and short autosomal sequences have been inconclusive . To examine the relationships among gibbon genera in more depth, we performed second-generation whole genome sequencing (WGS) to a mean of ∼15× coverage in two individuals from each genus. We developed a coalescent-based approximate Bayesian computation (ABC) method incorporating a model of sequencing error generated by high coverage exome validation to infer the branching order, divergence times, and effective population sizes of gibbon taxa. Although Hoolock and Symphalangus are likely sister taxa, we could not confidently resolve a single bifurcating tree despite the large amount of data analyzed. Instead, our results support the hypothesis that all four gibbon genera diverged at approximately the same time. Assuming an autosomal mutation rate of 1 × 10⁻⁹/site/year this speciation process occurred ∼5 MYA during a period in the Early Pliocene characterized by climatic shifts and fragmentation of the Sunda shelf forests. Whole genome sequencing of additional individuals will be vital for inferring the extent of gene flow among species after the separation of the gibbon genera.     

Molecular phylogeny of gibbons inferred from mitochondrial DNA sequences: Preliminary report

We analyzed the 896 base-pair (bp) mitochondrial DNA (mtDNA) sequences for seven gibbons, representative of three out of four subgenera. The result from our molecular analysis is consistent with previous studies as to the monophyly of subgenus Hylobates species, yet the relationship among subgenera remains slightly ambiguous. A striking result of the analysis is the phylogenetic location of Kloss's gibbon (H. klossii). Kloss's gibbon has been considered to be an initial off-shoot of the subgenus Hylobates because of its morphological primitiveness. However, our molecular data strongly suggest that Kloss's gibbon speciated most recently within the subgenus Hylobates. Correspondence to: S. Horai     

Concordance between vocal and genetic diversity in crested gibbons

Background  

Gibbons or small apes are, next to great apes, our closest living relatives, and form the most diverse group of contemporary hominoids. A characteristic trait of gibbons is their species-specific song structure, which, however, exhibits a certain amount of inter- and intra-individual variation. Although differences in gibbon song structure are routinely applied as taxonomic tool to identify subspecies and species, it remains unclear to which degree acoustic and phylogenetic differences are correlated. To trace this issue, we comparatively analyse song recordings and mitochondrial cytochrome b gene sequence data from 22 gibbon populations representing six of the seven crested gibbon species (genus Nomascus). In addition, we address whether song similarity and geographic distribution can support a recent hypothesis about the biogeographic history of crested gibbons.     

Mitochondrial genomes reveal slow rates of molecular evolution and the timing of speciation in beavers (Castor), one of the largest rodent species

Background

Beavers are one of the largest and ecologically most distinct rodent species. Little is known about their evolution and even their closest phylogenetic relatives have not yet been identified with certainty. Similarly, little is known about the timing of divergence events within the genus Castor.

Methodology/Principal Findings

We sequenced complete mitochondrial genomes from both extant beaver species and used these sequences to place beavers in the phylogenetic tree of rodents and date their divergence from other rodents as well as the divergence events within the genus Castor. Our analyses support the phylogenetic position of beavers as a sister lineage to the scaly tailed squirrel Anomalurus within the mouse related clade. Molecular dating places the divergence time of the lineages leading to beavers and Anomalurus as early as around 54 million years ago (mya). The living beaver species, Castor canadensis from North America and Castor fiber from Eurasia, although similar in appearance, appear to have diverged from a common ancestor more than seven mya. This result is consistent with the hypothesis that a migration of Castor from Eurasia to North America as early as 7.5 mya could have initiated their speciation. We date the common ancestor of the extant Eurasian beaver relict populations to around 210,000 years ago, much earlier than previously thought. Finally, the substitution rate of Castor mitochondrial DNA is considerably lower than that of other rodents. We found evidence that this is correlated with the longer life span of beavers compared to other rodents.

Conclusions/Significance

A phylogenetic analysis of mitochondrial genome sequences suggests a sister-group relationship between Castor and Anomalurus, and allows molecular dating of species divergence in congruence with paleontological data. The implementation of a relaxed molecular clock enabled us to estimate mitochondrial substitution rates and to evaluate the effect of life history traits on it.     

Niche Divergence versus Neutral Processes: Combined Environmental and Genetic Analyses Identify Contrasting Patterns of Differentiation in Recently Diverged Pine Species

Background and Aims

Solving relationships of recently diverged taxa, poses a challenge due to shared polymorphism and weak reproductive barriers. Multiple lines of evidence are needed to identify independently evolving lineages. This is especially true of long-lived species with large effective population sizes, and slow rates of lineage sorting. North American pines are an interesting group to test this multiple approach. Our aim is to combine cytoplasmic genetic markers with environmental information to clarify species boundaries and relationships of the species complex of Pinus flexilis, Pinus ayacahuite, and Pinus strobiformis.

Methods

Mitochondrial and chloroplast sequences were combined with previously obtained microsatellite data and contrasted with environmental information to reconstruct phylogenetic relationships of the species complex. Ecological niche models were compared to test if ecological divergence is significant among species.

Key Results and Conclusion

Separately, both genetic and ecological evidence support a clear differentiation of all three species but with different topology, but also reveal an ancestral contact zone between P. strobiformis and P. ayacahuite. The marked ecological differentiation of P. flexilis suggests that ecological speciation has occurred in this lineage, but this is not reflected in neutral markers. The inclusion of environmental traits in phylogenetic reconstruction improved the resolution of internal branches. We suggest that combining environmental and genetic information would be useful for species delimitation and phylogenetic studies in other recently diverged species complexes.     

Population genetic evidence for speciation pattern and gene flow between Picea wilsonii,P. morrisonicola and P. neoveitchii

Background and Aims

Genetic drift due to geographical isolation, gene flow and mutation rates together make it difficult to determine the evolutionary relationships of present-day species. In this study, population genetic data were used to model and decipher interspecific relationships, speciation patterns and gene flow between three species of spruce with similar morphology, Picea wilsonii, P. neoveitchii and P. morrisonicola. Picea wilsonii and P. neoveitchii occur from central to north-west China, where they have overlapping distributions. Picea morrisonicola, however, is restricted solely to the island of Taiwan and is isolated from the other two species by a long distance.

Methods

Sequence variations were examined in 18 DNA fragments for 22 populations, including three fragments from the chloroplast (cp) genome, two from the mitochondrial (mt) genome and 13 from the nuclear genome.

Key Results

In both the cpDNA and the mtDNA, P. morrisonicola accumulated more species-specific mutations than the other two species. However, most nuclear haplotypes of P. morrisonicola were shared by P. wilsonii, or derived from the dominant haplotypes found in that species. Modelling of population genetic data supported the hypothesis that P. morrisonicola derived from P. wilsonii within the more recent past, most probably indicating progenitor–derivative speciation with a distinct bottleneck, although further gene flow from the progenitor to the derivative continued. In addition, the occurrence was detected of an obvious mtDNA introgression from P. neoveitchii to P. wilsonii despite their early divergence.

Conclusions

The extent of mutation, introgression and lineage sorting taking place during interspecific divergence and demographic changes in the three species had varied greatly between the three genomes. The findings highlight the complex evolutionary histories of these three Asian spruce species.     

Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1664-4) contains supplementary material, which is available to authorized users.     

On the value of nuclear and mitochondrial gene sequences for reconstructing the phylogeny of vanilloid orchids (Vanilloideae, Orchidaceae)

Background and Aims

Most molecular phylogenetic studies of Orchidaceae have relied heavily on DNA sequences from the plastid genome. Nuclear and mitochondrial loci have only been superficially examined for their systematic value. Since 40% of the genera within Vanilloideae are achlorophyllous mycoheterotrophs, this is an ideal group of orchids in which to evaluate non-plastid gene sequences.

Methods

Phylogenetic reconstructions for Vanilloideae were produced using independent and combined data from the nuclear 18S, 5·8S and 26S rDNA genes and the mitochondrial atpA gene and nad1b-c intron.

Key Results

These new data indicate placements for genera such as Lecanorchis and Galeola, for which plastid gene sequences have been mostly unavailable. Nuclear and mitochondrial parsimony jackknife trees are congruent with each other and previously published trees based solely on plastid data. Because of high rates of sequence divergence among vanilloid orchids, even the short 5·8S rDNA gene provides impressive levels of resolution and support.

Conclusions

Orchid systematists are encouraged to sequence nuclear and mitochondrial gene regions along with the growing number of plastid loci available.Key words: 26S rDNA, 18S rDNA, 5·8S rDNA, atpA, nad1, orchids, plastid, Vanilla, vanilloid orchids, Vanilloideae     

The origin of malarial parasites in orangutans

Background

Recent findings of Plasmodium in African apes have changed our perspectives on the evolution of malarial parasites in hominids. However, phylogenetic analyses of primate malarias are still missing information from Southeast Asian apes. In this study, we report molecular data for a malaria parasite lineage found in orangutans.

Methodology/Principal Findings

We screened twenty-four blood samples from Pongo pygmaeus (Kalimantan, Indonesia) for Plasmodium parasites by PCR. For all the malaria positive orangutan samples, parasite mitochondrial genomes (mtDNA) and two antigens: merozoite surface protein 1 42 kDa (MSP-1₄₂) and circumsporozoite protein gene (CSP) were amplified, cloned, and sequenced. Fifteen orangutans tested positive and yielded 5 distinct mitochondrial haplotypes not previously found. The haplotypes detected exhibited low genetic divergence among them, indicating that they belong to one species. We report phylogenetic analyses using mitochondrial genomes, MSP-1₄₂ and CSP. We found that the orangutan malaria parasite lineage was part of a monophyletic group that includes all the known non-human primate malaria parasites found in Southeast Asia; specifically, it shares a recent common ancestor with P. inui (a macaque parasite) and P. hylobati (a gibbon parasite) suggesting that this lineage originated as a result of a host switch. The genetic diversity of MSP-1₄₂ in orangutans seems to be under negative selection. This result is similar to previous findings in non-human primate malarias closely related to P. vivax. As has been previously observed in the other Plasmodium species found in non-human primates, the CSP shows high polymorphism in the number of repeats. However, it has clearly distinctive motifs from those previously found in other malarial parasites.

Conclusion

The evidence available from Asian apes indicates that these parasites originated independently from those found in Africa, likely as the result of host switches from other non-human primates.     

Pollinator shifts as triggers of speciation in painted petal irises (Lapeirousia: Iridaceae)

Background and Aims

Adaptation to different pollinators has been hypothesized as one of the main factors promoting the formation of new species in the Cape region of South Africa. Other researchers favour alternative causes such as shifts in edaphic preferences. Using a phylogenetic framework and taking into consideration the biogeographical scenario explaining the distribution of the group as well as the distribution of pollinators, this study compares pollination strategies with substrate adaptations to develop hypotheses of the primary factors leading to speciation in Lapeirousia (Iridaceae), a genus of corm-bearing geophytes well represented in the Cape and presenting an important diversity of pollination syndromes and edaphic preferences.

Methods

Phylogenetic relationships are reconstructed within Lapeirousia using nuclear and plastid DNA sequence data. State-of-the-art methods in biogeography, divergence time estimation, character optimization and diversification rate assessments are used to examine the evolution of pollination syndromes and substrate shifts in the history of the group. Based on the phylogenetic results, ecological factors are compared for nine sister species pairs in Lapeirousia.

Key Results

Seventeen pollinator shifts and ten changes in substrate types were inferred during the evolution of the genus Lapeirousia. Of the nine species pairs examined, all show divergence in pollination syndromes, while only four pairs present different substrate types.

Conclusions

The available evidence points to a predominant influence of pollinator shifts over substrate types on the speciation process within Lapeirousia, contrary to previous studies that favoured a more important role for edaphic factors in these processes. This work also highlights the importance of biogeographical patterns in the study of pollination syndromes.     

Genome digging: insight into the mitochondrial genome of Homo

Background

A fraction of the Neanderthal mitochondrial genome sequence has a similarity with a 5,839-bp nuclear DNA sequence of mitochondrial origin (numt) on the human chromosome 1. This fact has never been interpreted. Although this phenomenon may be attributed to contamination and mosaic assembly of Neanderthal mtDNA from short sequencing reads, we explain the mysterious similarity by integration of this numt (mtAncestor-1) into the nuclear genome of the common ancestor of Neanderthals and modern humans not long before their reproductive split.

Principal Findings

Exploiting bioinformatics, we uncovered an additional numt (mtAncestor-2) with a high similarity to the Neanderthal mtDNA and indicated that both numts represent almost identical replicas of the mtDNA sequences ancestral to the mitochondrial genomes of Neanderthals and modern humans. In the proteins, encoded by mtDNA, the majority of amino acids distinguishing chimpanzees from humans and Neanderthals were acquired by the ancestral hominins. The overall rate of nonsynonymous evolution in Neanderthal mitochondrial protein-coding genes is not higher than in other lineages. The model incorporating the ancestral hominin mtDNA sequences estimates the average divergence age of the mtDNAs of Neanderthals and modern humans to be 450,000–485,000 years. The mtAncestor-1 and mtAncestor-2 sequences were incorporated into the nuclear genome approximately 620,000 years and 2,885,000 years ago, respectively.

Conclusions

This study provides the first insight into the evolution of the mitochondrial DNA in hominins ancestral to Neanderthals and humans. We hypothesize that mtAncestor-1 and mtAncestor-2 are likely to be molecular fossils of the mtDNAs of Homo heidelbergensis and a stem Homo lineage. The d_N/d_S dynamics suggests that the effective population size of extinct hominins was low. However, the hominin lineage ancestral to humans, Neanderthals and H. heidelbergensis, had a larger effective population size and possessed genetic diversity comparable with those of chimpanzee and gorilla.     

Phylogenetic Separation in Limb Use in Captive Gibbons (Hylobatidae): A Comparison Across the Primate Order

Although there have been few studies of self‐scratching in primates, some have reported distinct differences in whether hands or feet are used, and these variations seem to reflect the evolutionary history of the Order. Monkeys and prosimians use both hands and feet to self‐scratch while African great apes use hands almost exclusively. Gibbons represent an evolutionary divergence between monkeys and great apes and incidental observations at the Gibbon Conservation Center pointed to a difference in self‐scratching among the four extant gibbon genera (Hoolock, Nomascus, Symphalangus, and Hylobates). To validate and further explore these preliminary observations, we collected systematic data on self‐scratching from 32 gibbons, including nine species and all four genera. To supplement gibbon data, we also collected self‐scratching information from 18 great apes (four species), five prosimians (two species), 26 New World Monkeys (nine species) and 20 Old World Monkeys (seven species). All monkeys and some prosimians used both hands and feet to self‐scratch, whereas one prosimian species used only feet. All African great apes used hands exclusively (orangutans were an exception displaying occasional foot‐use). This appears to represent a fundamental difference between monkeys and great apes in limb use. Interestingly, there was a clear difference in self‐scratching between the four gibbon genera. Hylobates and Symphalangus self‐scratched only with hands (like all African great apes), while Hoolock and Nomascus self‐scratched with both hands and feet (like monkeys and prosimians). This difference in gibbon behavior may reflect the evolutionary history of gibbons as Hoolock and Nomascus are thought to have evolved before both Hylobates and Symphalangus. What evolutionary pressures led to this divergent pattern is currently opaque; however, this shift in limb preference may result from niche separation across the order facilitating differences in the behavioral repertoire associated with hind and forelimbs. Am. J. Primatol. 74:1035‐1043, 2012. © 2012 Wiley Periodicals, Inc.     

Incomplete Lineage Sorting Is Common in Extant Gibbon Genera

We sequenced reduced representation libraries by means of Illumina technology to generate over 1.5 Mb of orthologous sequence from a representative of each of the four extant gibbon genera (Nomascus, Hylobates, Symphalangus, and Hoolock). We used these data to assess the evolutionary relationships between the genera by evaluating the likelihoods of all possible bifurcating trees involving the four taxa. Our analyses provide weak support for a tree with Nomascus and Hylobates as sister taxa and with Hoolock and Symphalangus as sister taxa, though bootstrap resampling suggests that other phylogenetic scenarios are also possible. This uncertainty is due to short internal branch lengths and extensive incomplete lineage sorting across taxa. The true phylogenetic relationships among gibbon genera will likely require a more extensive whole-genome sequence analysis.     

Association of G-quadruplex forming sequences with human mtDNA deletion breakpoints

Background

Mitochondrial DNA (mtDNA) deletions cause disease and accumulate during aging, yet our understanding of the molecular mechanisms underlying their formation remains rudimentary. Guanine-quadruplex (GQ) DNA structures are associated with nuclear DNA instability in cancer; recent evidence indicates they can also form in mitochondrial nucleic acids, suggesting that these non-B DNA structures could be associated with mtDNA deletions. Currently, the multiple types of GQ sequences and their association with human mtDNA stability are unknown.

Results

Here, we show an association between human mtDNA deletion breakpoint locations (sites where DNA ends rejoin after deletion of a section) and sequences with G-quadruplex forming potential (QFP), and establish the ability of selected sequences to form GQ in vitro. QFP contain four runs of either two or three consecutive guanines (2G and 3G, respectively), and we identified four types of QFP for subsequent analysis: intrastrand 2G, intrastrand 3G, duplex derived interstrand (ddi) 2G, and ddi 3G QFP sequences. We analyzed the position of each motif set relative to either 5'' or 3'' unique mtDNA deletion breakpoints, and found that intrastrand QFP sequences, but not ddi QFP sequences, showed significant association with mtDNA deletion breakpoint locations. Moreover, a large proportion of these QFP sequences occur at smaller distances to breakpoints relative to distribution-matched controls. The positive association of 2G QFP sequences persisted when breakpoints were divided into clinical subgroups. We tested in vitro GQ formation of representative mtDNA sequences containing these 2G QFP sequences and detected robust GQ structures by UV–VIS and CD spectroscopy. Notably, the most frequent deletion breakpoints, including those of the "common deletion", are bounded by 2G QFP sequence motifs.

Conclusions

The potential for GQ to influence mitochondrial genome stability supports a high-priority investigation of these structures and their regulation in normal and pathological mitochondrial biology. These findings emphasize the potential importance of helicases that subsequently resolve GQ to maintain the stability of the mitochondrial genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-677) contains supplementary material, which is available to authorized users.     

Crested gibbon (Hylobates [Nomascus]) identification using noninvasively obtained DNA

The development of captive breeding programs for the crested gibbons (Hylobates [Nomascus]) of Indochina is hindered by the difficulty of sorting individuals into their correct species and subspecies groups. We describe taxon-specific DNA sequence variation in a 252-base pair region of the mitochondrial cytochrome b gene, which can be used to identify individuals of the three taxa of crested gibbon commonly found in Western zoos. These molecular genetic markers can be amplified from plucked hair and will permit the identification of morphologically similar females and crested gibbons of either sex of unknown origin. Noninvasive genotyping of captive animals will facilitate the genetic management of these critically endangered primates. © 1994 Wiley-Liss, Inc.     

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

Background

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal “barcoding gene”. An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

Methods and Findings

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a “barcode gap” using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

Conclusions

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a “species”.     

An Ancient Origin for the Enigmatic Flat-Headed Frogs (Bombinatoridae: Barbourula) from the Islands of Southeast Asia

Background

The complex history of Southeast Asian islands has long been of interest to biogeographers. Dispersal and vicariance events in the Pleistocene have received the most attention, though recent studies suggest a potentially more ancient history to components of the terrestrial fauna. Among this fauna is the enigmatic archaeobatrachian frog genus Barbourula, which only occurs on the islands of Borneo and Palawan. We utilize this lineage to gain unique insight into the temporal history of lineage diversification in Southeast Asian islands.

Methodology/Principal Findings

Using mitochondrial and nuclear genetic data, multiple fossil calibration points, and likelihood and Bayesian methods, we estimate phylogenetic relationships and divergence times for Barbourula. We determine the sensitivity of focal divergence times to specific calibration points by jackknife approach in which each calibration point is excluded from analysis. We find that relevant divergence time estimates are robust to the exclusion of specific calibration points. Barbourula is recovered as a monophyletic lineage nested within a monophyletic Costata. Barbourula diverged from its sister taxon Bombina in the Paleogene and the two species of Barbourula diverged in the Late Miocene.

Conclusions/Significance

The divergences within Barbourula and between it and Bombina are surprisingly old and represent the oldest estimates for a cladogenetic event resulting in living taxa endemic to Southeast Asian islands. Moreover, these divergence time estimates are consistent with a new biogeographic scenario: the Palawan Ark Hypothesis. We suggest that components of Palawan''s terrestrial fauna might have “rafted” on emergent portions of the North Palawan Block during its migration from the Asian mainland to its present-day position near Borneo. Further, dispersal from Palawan to Borneo (rather than Borneo to Palawan) may explain the current day disjunct distribution of this ancient lineage.     

Genetic Differentiation of the Mitochondrial Cytochrome Oxidase c Subunit I Gene in Genus Paramecium (Protista,Ciliophora)

Background

The mitochondrial cytochrome c oxidase subunit I (COI) gene is being used increasingly for evaluating inter- and intra-specific genetic diversity of ciliated protists. However, very few studies focus on assessing genetic divergence of the COI gene within individuals and how its presence might affect species identification and population structure analyses.

Methodology/Principal findings

We evaluated the genetic variation of the COI gene in five Paramecium species for a total of 147 clones derived from 21 individuals and 7 populations. We identified a total of 90 haplotypes with several individuals carrying more than one haplotype. Parsimony network and phylogenetic tree analyses revealed that intra-individual diversity had no effect in species identification and only a minor effect on population structure.

Conclusions

Our results suggest that the COI gene is a suitable marker for resolving inter- and intra-specific relationships of Paramecium spp.     

The Relationship between Leukocyte Mitochondrial DNA Copy Number and Telomere Length in Community-Dwelling Elderly Women

Purpose

Both telomere length and mitochondrial function are accepted as reflective indices of aging. Recent studies have shown that telomere dysfunction may influence impaired mitochondrial biogenesis and function. However, there has been no study regarding the possible association between telomere and mitochondrial function in humans. Therefore, the purpose of the study was to identify any relationships between mitochondrial and telomere function.

Methods

The present study included 129 community-dwelling, elderly women. The leukocyte mitochondrial DNA copy number and telomere length were measured using a quantitative real-time polymerase chain reaction method. Anthropometric measurement, biochemical blood testing, a depression screening questionnaire using a 15-question geriatric depression scale (GDS-15), and a cognitive function test using the Korean version of the mini mental state examination (K-MMSE) were performed.

Results

Leukocyte mtDNA copy number was positively associated with telomere length (r=0.39, p=<0.0001) and K-MMSE score (r=0.06, p=0.02). Additionally, leukocyte mtDNA copy number was negatively correlated with GDS-15 score (r=-0.17, p=0.04). Age (r=-0.15, p=0.09), waist circumference (r=-0.16, p=0.07), and serum ferritin level (r=-0.13, p=0.07) tended to be inversely correlated with leukocyte mtDNA copy number. With a stepwise multiple regression analysis, telomere length was found to be an independent factor associated with leukocyte mtDNA copy number after adjustment for confounding variables including age, body mass index, waist circumference, total cholesterol, HDL-cholesterol, LDL-cholesterol, triglycerides, hs-CRP, serum ferritin, HOMA-IR, K-MMSE, GDS-15, hypertension, diabetes, dyslipidemia, currently smoking, alcohol drinking, and regular exercise.

Conclusions

This study showed that leukocyte mtDNA copy number was positively correlated with leukocyte telomere length in community-dwelling elderly women. Our findings suggest that telomere function may influence mitochondrial function in humans.