Nuclear counterparts of the cytoplasmic mitochondrial 12S rRNA gene: A problem of ancient DNA and molecular phylogenies

Monkey mummy bones and teeth originating from the North Saqqara Baboon Galleries (Egypt), soft tissue from a mummified baboon in a museum collection, and nineteenth/twentieth-century skin fragments from mangabeys were used for DNA extraction and PCR amplification of part of the mitochondrial 12S rRNA gene. Sequences aligning with the 12S rRNA gene were recovered but were only distantly related to contemporary monkey mitochondrial 12S rRNA sequences. However, many of these sequences were identical or closely related to human nuclear DNA sequences resembling mitochondrial 12S rRNA (isolated from a cell line depleted in mitochondria) and therefore have to be considered contamination. Subsequently in a separate study we were able to recover genuine mitochondrial 12S rRNA sequences from many extant species of nonhuman Old World primates and sequences closely resembling the human nuclear integrations. Analysis of all sequences by the neighbor-joining (NJ) method indicated that mitochondrial DNA sequences and their nuclear counterparts can be divided into two distinct clusters. One cluster contained all temporary cytoplasmic mitochondrial DNA sequences and approximately half of the monkey nuclear mitochondriallike sequences. A second cluster contained most human nuclear sequences and the other half of monkey nuclear sequences with a separate branch leading to human and gorilla mitochondrial and nuclear sequences. Sequences recovered from ancient materials were equally divided between the two clusters. These results constitute a warning for when working with ancient DNA or performing phylogenetic analysis using mitochondrial DNA as a target sequence: Nuclear counterparts of mitochondrial genes may lead to faulty interpretation of results.Correspondence to: A.C. van der Kuyl     

Evolutionary hypervariability in the hinge region of the immunoglobulin alpha gene

The hinge region of the immunoglobulin molecule is responsible for antigen-binding and cross-linking reactions, varying the distance between the two antigen-binding sites. As the amino acid sequence of the hinge region is identical among immunoglobulin molecules of the same (sub)class, it has been regarded as a constant region. By comparison of the nucleotide sequences among primate C alpha genes, it is clear that there is a wide variety of length among the hinge regions of hominoid C alpha genes, which basically consist of tandem repeats of a 15 base-pair sequence. This reiterated structure probably facilitates rapid evolutionary changes in the length of the hinge region. The hinge region of the Old World monkey C alpha gene has a non-reiterated structure whose nucleotide sequence is quite different from those of the hominoid C alpha genes, although its surrounding region is conserved during evolution. This unusual hypervariability reveals that the hinge region has evolved as a semi-variable region in contrast to its constant character from an ontogenic viewpoint.     

Evolution of the pseudoautosomal boundary in Old World monkeys and great apes

Mammalian sex chromosomes are divided into sex-specific and pseudoautosomal regions. Sequences in the pseudoautosomal region recombine between the sex chromosomes; the sex-specific sequences normally do not. The interface between sex-specific and pseudoautosomal sequences is the pseudoautosomal boundary. The boundary is the centromeric limit to recombination in the pseudoautosomal region. In man, an Alu repeat element is found inserted at the boundary on the Y chromosome. In the evolutionary comparison conducted here, the Alu repeat element is found at the Y boundary in great apes, but it is not found there in two Old World monkeys. During the evolution of the Old World monkey and great ape lineages, homology between the sex chromosomes was maintained by recombination in the sequences telomeric to the Alu insertion site. The Alu repeat element did not create the present-day boundary; instead, it inserted at the preexisting boundary after the Old World monkey and great ape lineages diverged.     

HERV-K(OLD): ancestor sequences of the human endogenous retrovirus family HERV-K(HML-2)

Sequences homologous to the human endogenous retrovirus (HERV) family HERV-K(HML-2) are present in all Old World primate species. A previous study showed that a central region of the HERV-K(HML-2) gag genes in Hominoidea species displays a 96-bp deletion compared to the gag genes in lower Old World primates. The more ancient HERV-K(HML-2) sequences present in lower Old World primates were apparently not conserved during hominoid evolution, as opposed to the deletion variants. To further clarify the evolutionary origin of the HERV-K(HML-2) family, we screened GenBank with the 96-bp gag-sequence characteristic of lower Old World primates and identified, to date, 10 human sequence entries harboring either full-length or partially deleted proviral structures, probably representing remnants of a more ancient HERV-K(HML-2) variant. The high degree of mutations demonstrates the long-time presence of these HERV-K(OLD) proviruses in the genome. Nevertheless, they still belong to the HML-2 family as deduced from dot matrix and phylogenetic analyses. We estimate, based on the family ages of integrated Alu elements and on long terminal repeat (LTR) divergence data, that the average age of HERV-K(OLD) proviruses is ca. 28 million years, supporting an integration time before the evolutionary split of Hominoidea from lower Old World primates. Analysis of HERV-K(OLD) LTR sequences led to the distinction of two subgroups, both of which cluster with LTRs belonging to an evolutionarily older cluster. Taken together, our data give further insight into the evolutionary history of the HERV-K(HML-2) family during primate evolution.     

Molecular evolution of the psi eta-globin gene locus: gibbon phylogeny and the hominoid slowdown

An 8.4-kb genomic region spanning both the psi eta-globin gene locus and flanking DNA was sequenced from the common gibbon (Hylobates lar). In addition, sequencing of the entire orthologous region from galago (Galago crassicaudatus) was completed. The gibbon and galago sequences, along with published orthologous sequences from 10 other species, were aligned. These noncoding nucleotide sequences represented four human alleles, four apes (chimpanzee, gorilla, organgutan, and gibbon), an Old World monkey (rhesus monkey), two New World monkeys (spider and owl monkeys), tarsier, two strepsirhines (galago and lemur), and goat. Divergence and maximum parsimony analyses of the psi eta genomic region first groups humans and chimpanzees and then, at progressively more ancient branch points, successively joins gorillas, orangutans, gibbons, Old World monkeys, New World monkeys, tarsiers, and strepsirhines (the lemuriform-lorisiform branch of primates). This cladistic pattern supports the taxonomic grouping of all extant hominoids into family Hominidae, the division of Hominidae into subfamilies Hylobatinae (gibbons) and Homininae, the division of Homininae into tribes Pongini (orangutans) and Hominini, and the division of Hominini into subtribes Gorillina (gorillas) and Hominina (chimpanzees and humans). The additional gibbon and galago sequence data provide further support for the occurrence of a graded evolutionary-rate slowdown in the descent of simian primates, with the slowing rate being more pronounced in the great-ape and human lineages than in the gibbon or monkey lineages. A comparison of global versus local molecular clocks reveals that local clock predictions, when focused on a specific number of species within a narrow time frame, provide a more accurate estimate of divergence dates than do those of global clocks.     

Evolutionary rate of immunoglobulin alpha noncoding region is greater in hominoids than in Old World monkeys.

Recent studies on the molecular evolution of primates show that the evolutionary rate among hominoids is considerably slower than that among nonhominoid primates. However, this observation at the nucleotide-sequence level is restricted to the beta-globin family region. In this study, we sequenced orthologous immunoglobulin alpha (C alpha) genes of chimpanzee, gorilla, orangutan, and crab-eating macaque (an Old World monkey) and compared them with that of the human by using noncoding regions for analysis. Since significant differences in rates among hominoids were not found by using the relative rate test, we evaluated the ratio (R) of the evolutionary distance between Old World monkey and human to the distance between orangutan and human. The R value (1.12) for the C alpha gene was much smaller than the expected value (1.38-2.33), showing that the nucleotide substitution rate (= mutation rate per year under selective neutrality) of the C alpha gene is greater in the human lineage than in the Old World monkey lineage. We also did a similar analysis for the gamma 1-, gamma 2-, psi eta-, and delta-globin genes and found a considerable heterogeneity (1.12-2.37) among the R values, including that for the C alpha gene. This indicates that the hominoid slowdown of the evolutionary rate is not a universal phenomenon in primate evolution.     

A Cascade of Complex Subtelomeric Duplications during the Evolution of the Hominoid and Old World Monkey Genomes

Subtelomeric duplications of an obscure tubulin "genic" segment located near the telomere of human chromosome 4q35 have occurred at different evolutionary time points within the last 25 million years of the catarrhine (i.e., hominoid and Old World monkey) evolution. The analyses of these segments reported here indicate an exceptional level of evolutionary instability. Substantial intra- and interspecific differences in copy number and distribution are observed among cercopithecoid (Old World monkey) and hominoid genomes. Characterization of the hominoid duplicated segments reveals a strong positional bias within pericentromeric and subtelomeric regions of the genome. On the basis of phylogenetic analysis from predicted proteins and comparisons of nucleotide-substitution rates, we present evidence of a conserved b-tubulin gene among the duplications. Remarkably, the evolutionary conservation has occurred in a nonorthologous fashion, such that the functional copy has shifted its positional context between hominoids and cercopithecoids. We propose that, in a chimpanzee-human common ancestor, one of the paralogous copies assumed the original function, whereas the ancestral copy acquired mutations and eventually became silenced. Our analysis emphasizes the dynamic nature of duplication-mediated genome evolution and the delicate balance between gene acquisition and silencing.     

Plasticity of human chromosome 3 during primate evolution

Comparative mapping of more than 100 region-specific clones from human chromosome 3 in Bornean and Sumatran orangutans, siamang gibbon, and Old and New World monkeys allowed us to reconstruct ancestral simian and hominoid chromosomes. A single paracentric inversion derives chromosome 1 of the Old World monkey Presbytis cristata from the simian ancestor. In the New World monkey Callithrix geoffroyi and siamang, the ancestor diverged on multiple chromosomes, through utilizing different breakpoints. One shared and two independent inversions derive Bornean orangutan 2 and human 3, implying that neither Bornean orangutans nor humans have conserved the ancestral chromosome form. The inversions, fissions, and translocations in the five species analyzed involve at least 14 different evolutionary breakpoints along the entire length of human 3; however, particular regions appear to be more susceptible to chromosome reshuffling. The ancestral pericentromeric region has promoted both large-scale and micro-rearrangements. Small segments homologous to human 3q11.2 and 3q21.2 were repositioned intrachromosomally independent of the surrounding markers in the orangutan lineage. Breakage and rearrangement of the human 3p12.3 region were associated with extensive intragenomic duplications at multiple orangutan and gibbon subtelomeric sites. We propose that new chromosomes and genomes arise through large-scale rearrangements of evolutionarily conserved genomic building blocks and additional duplication, amplification, and/or repositioning of inherently unstable smaller DNA segments contained within them.     

The mitochondrial control region of Cervidae: evolutionary patterns and phylogenetic content

The mitochondrial control region (CR) sequence, also known as the D- loop, has been determined for six Cervidae (Artiodactyla, Ruminantia): the red and fallow deers (subfamily Cervinae), the brocket deer and two roe deers (subfamily Odocoileinae), and the Chinese water deer (Hydropotinae). These new sequences have been aligned with available cervid and bovid orthologues. Comparative analyses indicate that the 5'- peripheral domain exhibits a 75-bp length polymorphism near sequences associated with the termination of the H-strand replication. The New World Odocoileinae possess the longest cervid CR due to the presence of an additional 47-bp tandem repeat, located in the 3'-peripheral domain, downstream of the initiation site for H-strand replication (OH) and the first conserved sequence block (CSB-1). This insertion represents a duplication spanning the OH to CSB-1 region and constitutes an exclusive synapomorphy for New World Odocoileinae. Phylogenetic analyses of the complete CR support the paraphyly of antlered deers due to the nesting of the antlerless Hydropotes within Odocoileinae. Capreolus is the closest relative of Hydropotes, and the divergence of this Old World Odocoileinae clade may have occurred between 8.7 and 10.4 MYA. The conserved central domain of CR can be aligned across ungulates and indicates the Pecora monophyly, their close association with cetaceans, and the earlier emergence of suiformes.      

The complete mitochondrial genome of Tupaia belangeri and the phylogenetic affiliation of scandentia to other eutherian orders

The complete mitochondrial genome of Tupaia belangeri, a representative of the eutherian order Scandentia, was determined and compared with full-length mitochondrial sequences of other eutherian orders described to date. The complete mitochondrial genome is 16, 754 nt in length, with no obvious deviation from the general organization of the mammalian mitochondrial genome. Thus, features such as start codon usage, incomplete stop codons, and overlapping coding regions, as well as the presence of tandem repeats in the control region, are within the range of mammalian mitochondrial (mt) DNA variation. To address the question of a possible close phylogenetic relationship between primates and Tupaia, the evolutionary affinities among primates, Tupaia and bats as representatives of the Archonta superorder, ferungulates, guinea pigs, armadillos, rats, mice, and hedgehogs were examined on the basis of the complete mitochondrial DNA sequences. The opossum sequence was used as an outgroup. The trees, estimated from 12 concatenated genes encoded on the mitochondrial H-strand, add further molecular evidence against an Archonta monophyly. With the new data described in this paper, most of both the mitochondrial and the nuclear data point away from Scandentia as the closest extant relatives to primates. Instead, the complete mitochondrial data support a clustering of Scandentia with Lagomorpha connecting to the branch leading to ferungulates. This closer phylogenetic relationship of Tupaia to rabbits than to primates first received support from several analyses of nuclear and partial mitochondrial DNA data sets. Given that short sequences are of limited use in determining deep mammalian relationships, the partial mitochondrial data available to date support this hypothesis only tentatively. Our complete mitochondrial genome data therefore add considerably more evidence in support of this hypothesis.     

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.     

Mitochondrial DNA phylogeny of the Old-World monkey tribe Papionini.

The evolution of the Old World monkey tribe Papionini, composed of macaques, baboons, mandrills, drills, and mangabeys, was examined using mitochondrial DNA (mtDNA) sequence data on the cytochrome oxidase subunit II gene. When analyzed cladistically, these data support a baboon clade of savannah (Papio) plus gelada (Theropithecus) baboons, as well as a clade containing drill (Mandrillus) plus mangabey (Cerocebus) genera. This result stands in opposition to most morphological phylogenies, which break up the baboon clade by placing Papio and Mandrillus as sister taxa and Theropithecus as a more distantly related lineage. Analyses of COII gene sequences also suggest that the papionin ancestral stock divided into two lineages, one leading to macaques and the other to the purely African genera. From a molecular evolutionary perspective, the papionin COII gene sequences reveal a pattern of amino acid replacements concentrated in the regions spanning the mitochondrial membrane.     

Parallel evolution in eight-barbel loaches of the genus Lefua (Balitoridae, Cypriniformes) revealed by mitochondrial and nuclear DNA phylogenies

The evolutionary history of eight-barbel loaches of the genus Lefua contains important phylogenetic information that will aid in resolution of the faunal formations and evolutionary histories of Japanese and East Asian freshwater fishes. Our sequencing of the mitochondrial D-loop region in a large number of samples allowed construction of the most comprehensive phylogeny of these loaches to date; we demonstrated monophyly of five Lefua species and identified populations of Lufua. sp. and Lefua echigonia. Loaches inhabiting the Tokai region in Japan were morphologically and ecologically indistinguishable from Lefua sp. However, they were included in the L. echigonia lineage. We determined a novel phylogeny by sequencing the nuclear ribosomal S7 subunit and showed that nuclear DNA phylogeny essentially matched the mitochondrial DNA phylogeny. Loaches from the Tokai region were part of the L. echigonia lineage, indicating parallel evolution between Tokai loaches and Lefua sp. in western Japan. We presented the most robust phylogeny to date using concatenated mitochondrial and nuclear sequences. The wealth of molecular information allowed us to speculate on evolutionary processes in the genus Lefua.     

Concerted evolution of the tandem array encoding primate U2 snRNA occurs in situ, without changing the cytological context of the RNU2 locus.

In primates, the tandemly repeated genes encoding U2 small nuclear RNA evolve concertedly, i.e. the sequence of the U2 repeat unit is essentially homogeneous within each species but differs somewhat between species. Using chromosome painting and the NGFR gene as an outside marker, we show that the U2 tandem array (RNU2) has remained at the same chromosomal locus (equivalent to human 17q21) through multiple speciation events over > 35 million years leading to the Old World monkey and hominoid lineages. The data suggest that the U2 tandem repeat, once established in the primate lineage, contained sequence elements favoring perpetuation and concerted evolution of the array in situ, despite a pericentric inversion in chimpanzee, a reciprocal translocation in gorilla and a paracentric inversion in orang utan. Comparison of the 11 kb U2 repeat unit found in baboon and other Old World monkeys with the 6 kb U2 repeat unit in humans and other hominids revealed that an ancestral U2 repeat unit was expanded by insertion of a 5 kb retrovirus bearing 1 kb long terminal repeats (LTRs). Subsequent excision of the provirus by homologous recombination between the LTRs generated a 6 kb U2 repeat unit containing a solo LTR. Remarkably, both junctions between the human U2 tandem array and flanking chromosomal DNA at 17q21 fall within the solo LTR sequence, suggesting a role for the LTR in the origin or maintenance of the primate U2 array.     

Molecular phylogeny of ateline new world monkeys (Platyrrhini, atelinae) based on gamma-globin gene sequences: evidence that brachyteles is the sister group of lagothrix

Nucleotide sequences, each spanning approximately 7 kb of the contiguous gamma1 and gamma2 globin genomic loci, were determined for seven species representing all extant genera (Ateles, Lagothrix, Brachyteles, and Alouatta) of the New World monkey subfamily Atelinae. After aligning these seven ateline sequences with outgroup sequences from several other primate (non-ateline) genera, they were analyzed by maximum parsimony, maximum likelihood, and neighbor-joining algorithms. All three analyzes estimated the same phylogenetic relationships: [Alouatta [Ateles (Brachyteles, Lagothrix)]]. Brachyteles and Lagothrix are sister-groups supported by 100% of bootstrap replications in the parsimony analyses. Ateles joins this clade, followed by the basal genus Alouatta; these joinings were strongly supported, again with 100% bootstrap values. This cladistic pattern for the four ateline genera is congruent with that obtained in previous studies utilizing epsilon-globin, IRBP, and G6PD nuclear genomic sequences as well as mitochondrial COII sequences. Because the number of aligned nucleotide positions is much larger in the present datasetoff than in any of these other datasets, much stronger support was obtained for the cladistic classification that divides subfamily Atelinae into tribes Alouattini (Alouatta) and Atelini, while the latter divides into subtribes Atelina (Ateles) and Brachytelina (Brachyteles and Lagothrix).     

Phylogeny and zoogeography of six squirrel species of the genus sciurus (mammalia, rodentia), inferred from cytochrome B gene sequences

To investigate the phylogenetic relationships between the New World Sciurus and the Old World Sciurus and their biogeographic history, the partial mitochondrial cytochrome b gene sequences (1,040 base pairs) were analyzed on six Sciurus species: S. aberti, S. carolinensis, S. lis, S. niger, S. stramineus, and S. vulgaris. Phylogenetic trees (maximum parsimony, neighbor-joining, and maximum likelihood methods) commonly showed two groups with high bootstrap values (73-100%): one consisting of the New World Sciurus and the other consisting of the Old World Sciurus. Genetic distances among the New World Sciurus species were remarkably larger than that between two Sciurus species of the Old World, suggesting the earlier radiation of the New World Sciurus than the Old World Sciurus.     

An unusual source of apparent mitochondrial heteroplasmy: duplicate mitochondrial control regions in Thalassarche albatrosses

Molecular ecologists, in search of suitable molecular markers, frequently PCR-amplify regions of mitochondrial DNA from total DNA extracts. This approach, although common, is prone to the co-amplification of nuclear copies of transposed DNA sequences (numts), which can then generate apparent mitochondrial sequence heteroplasmy. In this study we describe the discovery of apparent mitochondrial sequence heteroplasmy in Thalassarche albatrosses but eliminate the possibility of true sequence heteroplasmy and numts and instead reveal the source of the apparent heteroplasmy to be a duplicated control region. The two control regions align easily but are not identical in sequence or in length. Comparisons of functionally significant conserved sequence blocks do not provide evidence of degeneration in either duplicate. Phylogenetic analyses of domain I of both control region copies in five Thalassarche species indicate that they are largely evolving in concert; however, a short section within them is clearly evolving independently. To our knowledge this is the first time contrasting evolutionary patterns have been reported for duplicate control regions. Available evidence suggests that this duplication may be taxonomically widespread, so the results presented here should be considered in future evolutionary studies targeting the control region of all Procellariiformes and potentially other closely related avian groups.