Proteogenomic review of the changes in primate apoC-I during evolution

Apolipoprotein C-I has evolved more rapidly than any of the other soluble apolipoproteins. During the course of primate evolution, the gene for this apolipoprotein was duplicated. Prompted by our observation that the two resulting genes encode two distinct forms of apoC-I in great apes, we have reviewed both the genomic and proteomic data to examine what changes have occurred during the course of primate evolution. We have found data showing that one of the duplicated genes, known to be a pseudogene in humans, was also a pseudogene in Denisovans and Neandertals. Using genomic and proteomic data for primates, we will provide in this review evidence that the duplication took place after the divergence of New World monkeys from the human lineage and that the formation of the pseudogene took place after the divergence of the bonobos and chimpanzees from the human lineage.     

Characterization of the gorilla carboxyl ester lipase locus, and the appearance of the carboxyl ester lipase pseudogene during primate evolution.

In this study we report on the isolation and characterization of the gorilla carboxyl ester lipase gene, CEL, and the corresponding CEL pseudogene. We also report on the age of the CEL pseudogene.The gorilla CEL gene is 10.5kb long and comprises 11exons intervened by introns similar to the situation in man, mouse and rat. The encoded protein is 998amino acids long and includes a 23amino acid-long leader peptide. Comparison of the coding sequence, excluding exon 11, of CEL from gorilla and man reveals a 97% similarity. Exon 11, which encodes the characteristic proline rich repeats, contains 39 repeated units in gorilla compared to 16 in man. A truncated CEL pseudogene, with the same organization as that found in man, is also shown to be present in the gorilla genome. The gorilla CEL pseudogene is 4.9kb in length and consists of 5exons interrupted by introns. Southern analysis of the gorilla CEL locus shows that the locus is arranged in a similar way as in man with the functional CEL gene being the most 5' one.To bring further insight to the events involved in the rearrangement of the CEL locus, genomic Southern analyses were performed across several primates; Homo sapiens, Pan troglodytes, Gorilla gorilla, Pongo pygmaeus and Macaca arctoides. Results presented show that the CEL gene duplication occurred prior to the separation of Hominidae (man, chimpanzee, gorilla and orangutan) from Old World monkeys (macaque). The deletion of the original CEL gene giving rise to the truncated version of the CEL gene seems, however, to be restricted to man and the great apes only.     

Mass spectral analyses of the two major apolipoproteins of great ape high density lipoproteins

The two major apolipoproteins associated with human and chimpanzee (Pan troglodytes) high density lipoproteins (HDL) are apoA-I and dimeric apoA-II. Although humans are closely related to great apes, apolipoprotein data do not exist for bonobos (Pan paniscus), western lowland gorillas (Gorilla gorilla gorilla) and the Sumatran orangutans (Pongo abelii). In the absence of any data, other great apes simply have been assumed to have dimeric apoA-II while other primates and most other mammals have been shown to have monomeric apoA-II. Using mass spectrometry, we have measured the molecular masses of apoA-I and apoA-II associated with the HDL of these great apes. Each was observed to have dimeric apoA-II. Being phylogenetically related, one would anticipate these apolipoproteins having a high percentage of invariant sequences when compared with human apolipoproteins. However, the orangutan, which diverged from the human lineage between 16 and 21 million years ago, had an apoA-II with the lowest monomeric mass, 8031.3 Da and the highest apoA-I value, 28,311.7 Da, currently reported for various mammals. Interestingly, the gorilla that diverged from the lineage leading to the human–chimpanzee branch after the orangutan had almost identical mass values to those reported for human apoA-I and apoA-II. But chimpanzee and the bonobo that diverged more recently had identical apoA-II mass values that were slightly larger than reported for the human apolipoprotein. The chimpanzee A-I mass values were very close to those of humans; however, the bonobo had values intermediate to the molecular masses of orangutan and the other great apes. With the already existing genomic data for chimpanzee and the recent entries for the orangutan and gorilla, we were able to demonstrate a close agreement between our mass spectral data and the calculated molecular weights determined from the predicted primary sequences of the respective apolipoproteins. Post-translational modification of these apolipoproteins, involving truncation and oxidation of methionine, are also reported.     

Mass spectrometric measurements of the apolipoproteins of bovine (Bos taurus) HDL

As is the case in most mammals, high density lipoproteins (HDL) also comprise the major group of lipid carriers that circulate in bovine (Bos taurus) blood. As a continuation of our proteogenomic studies of mammalian apolipoproteins, we have obtained molecular masses for several of the apolipoproteins associated with bovine HDL. The major apolipoprotein on the HDL surface is apoA-I, but other apolipoproteins were also detected. Using electrospray-ionization mass spectrometry (ESI-MS), we report on values for apolipoproteins, A-I, proA-I and A-II, as well as post-translationally modified apoA-I. Analyses of tryptic fragments did reveal the presence of apoA-IV and apoC-III. However, in contrast to our previous studies of other mammalian HDL, we did not detect apoC-I. Interestingly, examination of the current assembly for the bovine genome does not show any evidence for an apoC-I gene.     

Human type I hair keratin pseudogene ϕ hHaA has functional orthologs in the chimpanzee and gorilla: evidence for recent inactivation of the human gene after the Pan-Homo divergence

In addition to nine functional genes, the human type I hair keratin gene cluster contains a pseudogene, phihHaA (KRTHAP1), which is thought to have been inactivated by a single base-pair substitution that introduced a premature TGA termination codon into exon 4. Large-scale genotyping of human, chimpanzee, and gorilla DNAs revealed the homozygous presence of the phihHaA nonsense mutation in humans of different ethnic backgrounds, but its absence in the functional orthologous chimpanzee (cHaA) and gorilla (gHaA) genes. Expression analyses of the encoded cHaA and gHaA hair keratins served to highlight dramatic differences between the hair keratin phenotypes of contemporary humans and the great apes. The relative numbers of synonymous and non-synonymous substitutions in the phihHaA and cHaA genes, as inferred by using the gHaA gene as an outgroup, suggest that the human hHaA gene was inactivated only recently, viz., less than 240,000 years ago. This implies that the hair keratin phenotype of hominids prior to this date, and after the Pan-Homo divergence some 5.5 million years ago, could have been identical to that of the great apes. In addition, the homozygous presence of the phihHaA exon 4 nonsense mutation in some of the earliest branching lineages among extant human populations lends strong support to the "single African origin" hypothesis of modern humans.     

Human apolipoprotein C-I accounts for the ability of plasma high density lipoproteins to inhibit the cholesteryl ester transfer protein activity

The aim of the present study was to identify the protein that accounts for the cholesteryl ester transfer protein (CETP)-inhibitory activity that is specifically associated with human plasma high density lipoproteins (HDL). To this end, human HDL apolipoproteins were fractionated by preparative polyacrylamide gradient gel electrophoresis, and 30 distinct protein fractions with molecular masses ranging from 80 down to 2 kDa were tested for their ability to inhibit CETP activity. One single apolipoprotein fraction was able to completely inhibit CETP activity. The N-terminal sequence of the 6-kDa protein inhibitor matched the N-terminal sequence of human apoC-I, the inhibition was completely blocked by specific anti-apolipoprotein C-I antibodies, and mass spectrometry analysis confirmed the identity of the isolated inhibitor with full-length human apoC-I. Pure apoC-I was able to abolish CETP activity in a concentration-dependent manner and with a high efficiency (IC(50) = 100 nmol/liter). The inhibitory potency of total delipidated HDL apolipoproteins completely disappeared after a treatment with anti-apolipoprotein C-I antibodies, and the apoC-I deprivation of native plasma HDL by immunoaffinity chromatography produced a mean 43% rise in cholesteryl ester transfer rates. The main localization of apoC-I in HDL and not in low density lipoprotein in normolipidemic plasma provides further support for the specific property of HDL in inhibiting CETP activity.     

Dynamics of DNA Methylation in Recent Human and Great Ape Evolution

DNA methylation is an epigenetic modification involved in regulatory processes such as cell differentiation during development, X-chromosome inactivation, genomic imprinting and susceptibility to complex disease. However, the dynamics of DNA methylation changes between humans and their closest relatives are still poorly understood. We performed a comparative analysis of CpG methylation patterns between 9 humans and 23 primate samples including all species of great apes (chimpanzee, bonobo, gorilla and orangutan) using Illumina Methylation450 bead arrays. Our analysis identified ∼800 genes with significantly altered methylation patterns among the great apes, including ∼170 genes with a methylation pattern unique to human. Some of these are known to be involved in developmental and neurological features, suggesting that epigenetic changes have been frequent during recent human and primate evolution. We identified a significant positive relationship between the rate of coding variation and alterations of methylation at the promoter level, indicative of co-occurrence between evolution of protein sequence and gene regulation. In contrast, and supporting the idea that many phenotypic differences between humans and great apes are not due to amino acid differences, our analysis also identified 184 genes that are perfectly conserved at protein level between human and chimpanzee, yet show significant epigenetic differences between these two species. We conclude that epigenetic alterations are an important force during primate evolution and have been under-explored in evolutionary comparative genomics.     

Linkage analysis of the genetic determinants of high density lipoprotein concentrations and composition: evidence for involvement of the apolipoprotein A-II and cholesteryl ester transfer protein loci

We have tested for evidence of linkage between the genetic loci determining concentrations and composition of plasma high density lipoproteins (HDL) with the genes for the major apolipoproteins and enzymes participating in lipoprotein metabolism. These genes include those encoding various apolipoproteins (apo), including apoA-I, apoA-II, apoA-IV, apoB, apoC-I, apoC-II, apoC-III, apoE, and apo(a), cholesteryl ester transfer protein (CETP), HDL-binding protein, lipoprotein lipase, and the low density lipoprotein (LDL) receptor. Polymorphisms of these genes, and nearby highly polymorphic simple sequence repeat markers, were examined by quantitative sib-pair linkage analysis in 30 coronary artery disease families consisting of a total of 366 individuals. Evidence for linkage was observed between a marker locus D16S313 linked to the CETP locus and a locus determining plasma HDL-cholesterol concentration (P = 0.002), and the genetic locus for apoA-II and a locus determining the levels of the major apolipoproteins of HDL, apoA-I and apoA-II (P = 0.009 and 0.02, respectively). HDL level was also influenced by the variation at the apo(a) locus on chromosome 6 (P = 0.02). Thus, these data indicate the simultaneous involvement of at least two different genetic loci in the determination of the levels of HDL and its associated lipoproteins.     

Two copies of the human apolipoprotein C-I gene are linked closely to the apolipoprotein E gene

The gene for human apolipoprotein (apo) C-I was selected from human genomic cosmid and lambda libraries. Restriction endonuclease analysis showed that the gene for apoC-I is located 5.5 kilobases downstream of the gene for apoE. A copy of the apoC-I gene, apoC-I', is located 7.5 kilobases downstream of the apoC-I gene. Both genes contain four exons and three introns; the apoC-I gene is 4653 base pairs long, the apoC-I' gene 4387 base pairs. In each gene, the first intron is located 20 nucleotides upstream from the translation start signal; the second intron, within the codon of Gly-7 of the signal peptide region; and the third intron, within the codon for Arg39 of the mature plasma protein coding region. The upstream apoC-I gene encodes the known apoC-I plasma protein and differs from the downstream apoC-I' gene in about 9% of the exon nucleotide positions. The most important difference between the exons results in a change in the codon for Gln-2 of the signal peptide region, which introduces a translation stop signal in the downstream gene. Major sequence differences are found in the second and third introns of the apoC-I and apoC-I' genes, which contain 9 and 7.5 copies, respectively, of Alu family sequences. The apoC-I gene is expressed primarily in the liver, and it is activated when monocytes differentiate into macrophages. In contrast, no mRNA product of the apoC-I' gene can be detected in any tissue, suggesting that it may be a pseudogene. The similar structures and the proximity of the apoE and apoC-I genes suggest that they are derived from a common ancestor. Furthermore, they may be considered to be constituents of a family of seven apolipoprotein genes (apoE, -C-I, -C-II, -C-III, -A-I, -A-II, and -A-IV) that have a common evolutionary origin.     

Plasma turnover of HDL apoC-I,apoC-III,and apoE in humans: in vivo evidence for a link between HDL apoC-III and apoA-I metabolism

Numerous factors are known to affect the plasma metabolism of HDL, including lipoprotein receptors, lipid transfer protein, lipolytic enzymes and HDL apolipoproteins. In order to better define the role of HDL apolipoproteins in determining plasma HDL concentrations, the aims of the present study were: a) to compare the in vivo rate of plasma turnover of HDL apolipoproteins [i.e., apolipoprotein A-I (apoA-I), apoC-I, apoC-III, and apoE], and b) to investigate to what extent these metabolic parameters are related to plasma HDL levels. We thus studied 16 individuals with HDL cholesterol levels ranging from 0.56-1.66 mmol/l and HDL apoA-I levels ranging from 89-149 mg/dl. Plasma kinetics of HDL apolipoproteins were investigated using a primed constant (12 h) infusion of deuterated leucine. Plasma HDL apolipoprotein levels were 41.8 +/- 1.5, 9.7 +/- 0.5, 4.9 +/- 0.5, and 0.7 +/- 0.1 micromol/l for apoA-I, apoC-I, apoC-III and apoE. Plasma transport rates (TRs) were 388.6 +/- 24.7, 131.5 +/- 12.5, 66.5 +/- 9.1, and 31.4 +/- 3.3 nmol.kg-1.day-1; and residence times (RTs) were 5.1 +/- 0.4, 3.7 +/- 0.3, 3.6 +/- 0.3, and 1.1 +/- 0.1 days, respectively. HDL cholesterol and apoA-I levels were significantly correlated with HDL apoA-I RT (r = 0.69 and r = 0.56), and were not significantly correlated with HDL apoA-I TR. In contrast, HDL apoC-I, apoC-III, and apoB levels were all positively related to their TRs and not their RTs. HDL apoC-III TR was positively correlated with levels of HDL apoC-III (r = 0.73, P < 0.01), and with those of HDL cholesterol and apoA-I (r = 0.54 and r = 0.53, P < 0.05, respectively). HDL apoC-III TR was in turn related to HDL apoA-I RT (r = 0.51, P < 0.05). Together, these results provide in vivo evidence for a link between the metabolism of HDL apoC-III and apoA-I, and suggest a role for apoC-III in the regulation of plasma HDL levels.     

Novel use of a chimpanzee pseudogene for chromosomal mapping of human cytochrome c oxidase subunit IV

We have isolated a chimpanzee processed pseudogene for subunit IV of cytochrome c oxidase (COX; EC 1.9.3.1) by screening a chimpanzee genomic library in lambda Charon 32 with a bovine liver cDNA encoding COX subunit IV (COX IV), and localized it to a 1.9-kb HindIII fragment. Southern-blot analysis of genomic DNA from five primates showed that DNAs from human, gorilla, and chimpanzee each contained the 1.9-kb pseudogene fragment, whereas orangutan and pigtail macaque monkey DNA did not. This result clearly indicates that the pseudogene arose before the divergence of the chimpanzee and gorilla from the primate lineage. By screening Chinese hamster x human hybrid panels with the human COX4 cDNA, we have mapped COX4 genes to two human chromosomes, 14 and 16. The 1.9-kb HindIII fragment containing the pseudogene, COX4P1, can be assigned to chromosome 14, and by means of rearranged chromosomes in somatic cell hybrids, to 14q21-qter. Similarly, the functional gene, COX4, has been mapped to 16q22-qter.     

Structure and expression of dog apolipoprotein A-I, E, and C-I mRNAs: implications for the evolution and functional constraints of apolipoprotein structure

Dog apolipoprotein (apo) C-I, A-I, and E cDNA clones were identified in a dog liver cDNA library in lambda gt10 by hybridization to synthetic oligonucleotide probes with the corresponding human DNA sequences. The longest clone for each apolipoprotein was completely sequenced. The apoC-I cDNA sequence predicts a protein of 62 residue mature peptide preceded by a 26 amino acid signal peptide. The apoA-I cDNA sequence predicts a 242 residue mature peptide, a 6 residue pro-segment, and an 18 residue signal peptide. The apoE cDNA, which lacks the signal peptide region, predicts a mature peptide of 291 amino acid residues. Slot blot hybridization of total RNA isolated from various dog tissues to dog apoC-I, A-I, and E cDNA probes indicates that apoC-I mRNA is detectable in liver only, apoA-I mRNA is present in liver and small intestine, though the concentration in the latter tissue is only approximately 15% of that in the liver, and apoE mRNA is present in multiple tissues including liver, jejunum, urinary bladder, ileum, colon, brain, kidney, spleen, pancreas, and testis with relative concentrations (%) of 100, 17.5, 7.5, 6.9, 5.9, 5.5, 5.0, 3.3, 1.0, and 1.0, respectively. These tissue distributions indicate that nascent lipoprotein particles produced in the dog small intestine would contain apoA-I and apoE but not apoC-I. The widespread tissue distribution of apoE mRNA indicates that like other mammals, peripheral synthesis of apoE contributes significantly to the total apoE pool in dog. We next compared the cDNA sequences among different vertebrate species for apoC-I (human and dog), A-I (human, rat, dog, rabbit and chicken), and E (human, rat, dog and rabbit) and calculated the rate of nucleotide substitution for each gene. Our results indicate that apoC-I has evolved rather rapidly and that on the whole, apoA-I is more conservative than apoE, contradictory to an earlier suggestion. ApoA-I is also more conservative than a region (residues 4204-4536) at the carboxyl-terminal portion, but less conservative than a region (residues 595-979) at the amino-terminal portion of apoB-100. Some regions in each of the apolipoproteins studied are better conserved than others and the rate of evolution of individual regions seems to be related to the stringency of functional requirements. Finally, we estimate that the human apoC-I pseudogene arose more than 35 million years ago, becoming nonfunctional soon after its formation.     

Chimpanzee Malaria Parasites Related to Plasmodium ovale in Africa

Since the 1970''s, the diversity of Plasmodium parasites in African great apes has been neglected. Surprisingly, P. reichenowi, a chimpanzee parasite, is the only such parasite to have been molecularly characterized. This parasite is closely phylogenetically related to P. falciparum, the principal cause of the greatest malaria burden in humans. Studies of malaria parasites from anthropoid primates may provide relevant phylogenetic information, improving our understanding of the origin and evolutionary history of human malaria species. In this study, we screened 130 DNA samples from chimpanzees (Pan troglodytes) and gorillas (Gorilla gorilla) from Cameroon for Plasmodium infection, using cytochrome b molecular tools. Two chimpanzees from the subspecies Pan t. troglodytes presented single infections with Plasmodium strains molecularly related to the human malaria parasite P. ovale. These chimpanzee parasites and 13 human strains of P. ovale originated from a various sites in Africa and Asia were characterized using cytochrome b and cytochrome c oxidase 1 mitochondrial partial genes and nuclear ldh partial gene. Consistent with previous findings, two genetically distinct types of P. ovale, classical and variant, were observed in the human population from a variety of geographical locations. One chimpanzee Plasmodium strain was genetically identical, on all three markers tested, to variant P. ovale type. The other chimpanzee Plasmodium strain was different from P. ovale strains isolated from humans. This study provides the first evidence of possibility of natural cross-species exchange of P. ovale between humans and chimpanzees of the subspecies Pan t. troglodytes.     

Human-specific amino acid changes found in 103 protein-coding genes

We humans have many characteristics that are different from those of the great apes. These human-specific characters must have arisen through mutations accumulated in the genome of our direct ancestor after the divergence of the last common ancestor with chimpanzee. Gene trees of human and great apes are necessary for extracting these human-specific genetic changes. We conducted a systematic analysis of 103 protein-coding genes for human, chimpanzee, gorilla, and orangutan. Nucleotide sequences for 18 genes were newly determined for this study, and those for the remaining genes were retrieved from the DDBJ/EMBL/GenBank database. The total number of amino acid changes in the human lineage was 147 for 26,199 codons (0.56%). The total number of amino acid changes in the human genome was, thus, estimated to be about 80,000. We applied the acceleration index test and Fisher's synonymous/nonsynonymous exact test for each gene tree to detect any human-specific enhancement of amino acid changes compared with ape branches. Six and two genes were shown to have significantly higher nonsynonymous changes at the human lineage from the acceleration index and exact tests, respectively. We also compared the distribution of the differences of the nonsynonymous substitutions on the human lineage and those on the great ape lineage. Two genes were more conserved in the ape lineage, whereas one gene was more conserved in the human lineage. These results suggest that a small proportion of protein-coding genes started to evolve differently in the human lineage after it diverged from the ape lineage.     

Mass spectral analysis of the apolipoproteins on mouse high density lipoproteins. Detection of post-translational modifications

Using mass spectrometry, we have recently reported on molecular masses of the apolipoproteins associated with porcine and equine HDL. In addition to obtaining accurate masses for the various apolipoproteins, we also were able to detect mass variations due to post-translational modifications. In the present study, we have used these same approaches to characterize the apolipoproteins in two inbred mouse strains, C57BL/6 and BALB/c. Comparing our molecular mass data with calculated values for molecular weight, we were able to identify the correct sequences for several of the major apolipoproteins. Analyses were carried out on the apolipoproteins of ultracentrifugally isolated HDL. Prior to analyses by electrospray ionization mass spectrometry (ESI-MS), the apolipoproteins were separated either by size exclusion or reverse phase chromatography. The molecular masses of apoA-I, proapoA-I, apoA-II, proapoA-II, apoC-I and apoC-III were obtained. Comparing the values obtained for the two strains, differences in the molecular masses of apoA-I, apoA-II and apoC-III were observed. In this study, post-translationally modified apolipoproteins, involving loss of amino acids from both the N- and C-termini, oxidation of methionine residues and possible acylation, were noted following reverse-phase separation. Further analyses by tandem mass spectrometry (MSMS) done on the tryptic digests of apolipoproteins separated by reverse phase chromatography enabled us to confirm sequence differences between the two strains, to verify selected apoA-I sequences that had been entered into the GenBank and to identify which methionines in apoA-I, apoC-III and apoE had been converted to methionine sulfoxides.     

Human apoA-I expression in CETP transgenic rats leads to lower levels of apoC-I in HDL and to magnification of CETP-mediated lipoprotein changes

Plasma cholesteryl ester transfer protein (CETP) has a profound effect on neutral lipid transfers between HDLs and apolipoprotein B (apoB)-containing lipoproteins when it is expressed in combination with human apoA-I in HuAI/CETP transgenic (Tg) rodents. In the present study, human apoA-I-mediated lipoprotein changes in HuAI/CETPTg rats are characterized by 3- to 5-fold increments in the apoB-containing lipoprotein-to-HDL cholesterol ratio, and in the cholesteryl ester-to-triglyceride ratio in apoB-containing lipoproteins. These changes occur despite no change in plasma CETP concentration in HuAI/CETPTg rats, as compared with CETPTg rats. A number of HDL apolipoproteins, including rat apoA-I and rat apoC-I are removed from the HDL surface as a result of human apoA-I overexpression. Rat apoC-I, which is known to constitute a potent inhibitor of CETP, accounts for approximately two-thirds of CETP inhibitory activity in HDL from wild-type rats, and the remainder is carried by other HDL-bound apolipoprotein inhibitors. It is concluded that human apoA-I overexpression modifies HDL particles in a way that suppresses their ability to inhibit CETP. An apoC-I decrease in HDL of HuAI/CETPTg rats contributes chiefly to the loss of the CETP-inhibitory potential that is normally associated with wild-type HDL.     

Helical domains that mediate lipid solubilization and ABCA1-specific cholesterol efflux in apolipoproteins C-I and A-II

Many of the apolipoproteins in HDL can elicit cholesterol efflux via ABCA1, a critical initial step in HDL formation. Recent work has indicated that omnipresent amphipathic helices play a critical role, and these have been studied intensively in the most common HDL protein, apolipoprotein (apo)A-I. However, little information exists about helical domain arrangement in other apolipoproteins. We studied two of the smallest apolipoproteins known to interact with ABCA1, human apoA-II and apoC-I, in terms of ability to reorganize phospholipid (PL) bilayers and to promote ABCA1-mediated cholesterol. We found that both proteins contained helical domains that were fast and slow with respect to solubilizing PL. ABCA1-medated efflux required a minimum of a bihelical polypeptide comprised of at least one each of a slow and fast lipid reorganizing domain. In both proteins, the fast helix was located at the C terminus preceded by a slow helix. Helical placement in apoC-I was not critical for ABCA1 activity, but helix swaps in apoA-II dramatically disrupted cholesterol efflux, indicating that the tertiary structure of the longer apolipoprotein is important for the pathway. This work has implications for a more complete molecular understanding of apolipoprotein-mediated cholesterol efflux.     

The mitochondrial genome of the chimpanzee louse,Pediculus schaeffi: insights into the process of mitochondrial genome fragmentation in the blood-sucking lice of great apes

Background

Blood-sucking lice in the genera Pediculus and Pthirus are obligate ectoparasites of great apes. Unlike most bilateral animals, which have 37 mitochondrial (mt) genes on a single circular chromosome, the sucking lice of humans have extensively fragmented mt genomes. The head louse, Pediculus capitis, and the body louse, Pe. humanus, have their 37 mt genes on 20 minichromosomes. The pubic louse, Pthirus pubis, has its 34 mt genes known on 14 minichromosomes. To understand the process of mt genome fragmentation in the sucking lice of great apes, we sequenced the mt genome of the chimpanzee louse, Pe. schaeffi, and compared it with the three human lice.

Results

We identified all of the 37 mt genes typical of bilateral animals in the chimpanzee louse; these genes are on 18 types of minichromosomes. Seventeen of the 18 minichromosomes of the chimpanzee louse have the same gene content and gene arrangement as their counterparts in the human head louse and the human body louse. However, five genes, cob, trnS₁, trnN, trnE and trnM, which are on three minichromosomes in the human head louse and the human body louse, are together on one minichromosome in the chimpanzee louse.

Conclusions

Using the human pubic louse, Pt. pubis, as an outgroup for comparison, we infer that a single minichromosome has fragmented into three in the lineage leading to the human head louse and the human body louse since this lineage diverged from the chimpanzee louse ~6 million years ago. Our results provide insights into the process of mt genome fragmentation in the sucking lice in a relatively fine evolutionary scale.

Electronic supplementary material

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