The involucrin genes of the white-fronted capuchin and cottontop tamarin: the platyrrhine middle region.

In all anthropoid species, the coding region of the involucrin gene contains a segment of short tandem repeats that were added sequentially, beginning in a common anthropoid ancestor. The involucrin coding region of each of two platyrrhine species, the white-fronted capuchin (Cebus albifrons) and the cottontop tamarin (Saguinus oedipus), has now been cloned and sequenced. These genes share with the genes of the catarrhines the repeats added in the common anthropoid lineage (the early region). After their divergence, the platyrrhines, like the catarrhines, continued to add repeats vectorially 5' of the early region, to form a middle region. The mechanism that was established in the common anthropoid lineage for the addition of repeats at a definite site in the coding region was transmitted to both platyrrhines and catarrhines, enabling each to generate its middle region independently. The process of vectorial repeat addition continued in two platyrrhine sublineages after their divergence from each other.     

The involucrin gene of Old-World monkeys and other higher primates: synapomorphies and parallelisms resulting from the same gene-altering mechanism.

The involucrin gene of platyrrhines and hominoids contains a segment of 10-codon repeats which were added vectorially at the same site in the coding region. We have now cloned and sequenced the involucrin gene of four cercopithecoid monkeys--two macaques (mulatta and fascicularis) and two Cercopithecus monkeys (aethiops and hamlyni). Each gene contains a similar segment of short repeats; some of these were added in a common anthropoid lineage, others were added in a common catarrhine lineage, and still others were added in a common macaque or Cercopithecus lineage. Repeats added before a lineage diverges become synapomorphies in the sister taxa resulting from the divergence. Repeats added independently in different diverged lineages become parallelisms. The synapomorphies are the result of the action of a targeted duplication mechanism acting in a common ancestral lineage, but the parallelisms are the result of the same duplication mechanism transmitted to successively divergent sublineages and acting independently in each.     

The involucrin gene of the owl monkey: origin of the early region

A large part of the coding region of the hominoid involucrin gene is of recent origin. This part of the gene, which we have called the modern segment, contains numerous repeats of a sequence of 10 codons, created by multiple duplications some of which consist of 3-12 repeats. We have sequenced two alleles of the involucrin gene in the owl monkey and found that the involucrin gene of this species also possesses a modern segment. By comparing the modern segment of the owl monkey with that of the hominoids, we find that only a part of this segment is shared by the two species. We call this part the early region because it must have originated in a common ancestor of the anthropoids. The rest of the hominoid modern segment does not correspond to any groups of repeats in the owl monkey and was therefore created after divergence of the two lineages. As in the hominoids, the latest additions to the modern segment of the owl monkey have been in its 5' half, which possesses different duplication patterns in the two alleles. Lineage divergences within the anthropoids can be detected at different sites within the modern segment.      

The involucrin gene of the orangutan: generation of the late region as an evolutionary trend in the hominoids

In the evolutionary line leading to the higher primates, the coding region of the involucrin gene evolved a segment consisting of numerous repeats of a 10-codon sequence. Additions to this segment of repeats have been made successively, thus generating regions that can be defined as early, middle, and late. The involucrin gene of the orangutan (Pongo pygmaeus abelii) possesses a segment of repeats whose early region has the same repeat structure as that in other anthropoids. The middle region is not similar in repeat structure to that of all anthropoids but is similar to that of other hominoids. The late region is unique to the species; it does not correspond at all in its repeat structure to that of the human or gorilla and is much larger. The late region of the orangutan was generated by duplications of blocks of older repeats clearly belonging to the middle region. Continued duplications extending the late region are an evolutionary trend in the hominoids. The process of addition of repeats at a particular location is a more significant aspect of the evolution of involucrin than are random nucleotide substitutions; in addition, it has proceeded more rapidly.      

Polymorphism due to variable number of repeats in the human involucrin gene.

The coding region of the involucrin gene in higher primates contains a segment consisting of numerous tandem repeats of a 10-codon sequence. The process of repeat addition began in a common ancestor of all higher primates and subsequent repeats were added vectorially. As a result, the principal site of repeat addition has moved in the 3' to 5' direction and the most recently generated repeats (the late region) are close to the 5' end of the segment of repeats. In the human, most of the late region is made up of two different blocks, each consisting of nearly identical repeats. We describe here five polymorphic forms resulting from the addition of differing numbers of repeats to each block. As the variety and nature of the polymorphic alleles are different in different human populations, we postulate that the process of repeat addition is genetically determined.     

Remodeling of the involucrin gene during primate evolution

The protein involucrin is a product of terminal differentiation in the epidermal cell and related cell types. By comparing the nucleotide sequence of the involucrin gene of the lemur with that of the human, it is clear that the gene has undergone unusual evolution in the primates. The coding region of the gene contains an ancestral segment, most of which is common to the lemur and the human, and a species-specific segment of repeats derived from the ancestral segment. Instead of the modern segment of repeats found in the human gene, the lemur gene possesses repeats derived from another sequence at a different location in the ancestral segment. The two kinds of segments of repeats probably represent alternative ways of creating a repeat structure in the involucrin molecule. The modern segment of repeats must have been created after divergence of the higher primates from the prosimians.     

Evolution of beta satellite DNA sequences: evidence for duplication-mediated repeat amplification and spreading

In this article, we report studies on the evolutionary history of beta satellite repeats (BSR) in primates. In the orangutan genome, the bulk of BSR sequences was found organized as very short stretches of approximately 100 to 170 bp, embedded in a 60-kb to 80-kb duplicated DNA segment. The estimated copy number of the duplicon that carries BSR sequences ranges from 70 to 100 per orangutan haploid genome. In both macaque and gibbon, the duplicon mapped to a single chromosomal region at the boundary of the rDNA on the marker chromosome (chromosome 13 and 12, respectively). However, only in the gibbon, the duplicon comprised 100 bp of beta satellite. Thus, the ancestral copy of the duplicon appeared in Old World monkeys ( approximately 25 to approximately 35 MYA), whereas the prototype of beta satellite repeats took place in a gibbon ancestor, after apes/Old World monkeys divergence ( approximately 25 MYA). Subsequently, a burst in spreading of the duplicon that carries the beta satellite was observed in the orangutan, after lesser apes divergence from the great apes-humans lineage ( approximately 18 MYA). The analysis of the orangutan genome also indicated the existence of two variants of the duplication that differ for the length (100 or 170 bp) of beta satellite repeats. The latter organization was probably generated by nonhomologous recombination between two 100-bp repeated regions, and it likely led to the duplication of the single Sau3A site present in the 100-bp variant, which generated the prototype of Sau3A 68-bp beta satellite tandem organization. The two variants of the duplication, although with a different ratios, characterize the hominoid genomes from the orangutan to humans, preferentially involving acrocentric chromosomes. At variance to alpha satellite, which appeared before the divergence of New World and Old World monkeys, the beta satellite evolutionary history began in apes ancestor, where we have first documented a low-copy, nonduplicated BSR sequence. The first step of BSR amplification and spreading occurred, most likely, because the BSR was part of a large duplicon, which underwent a burst dispersal in great apes' ancestor after the lesser apes' branching. Then, after orangutan divergence, BSR acquired the clustered structural organization typical of satellite DNA.     

A high-resolution map of synteny disruptions in gibbon and human genomes

Gibbons are part of the same superfamily (Hominoidea) as humans and great apes, but their karyotype has diverged faster from the common hominoid ancestor. At least 24 major chromosome rearrangements are required to convert the presumed ancestral karyotype of gibbons into that of the hominoid ancestor. Up to 28 additional rearrangements distinguish the various living species from the common gibbon ancestor. Using the northern white-cheeked gibbon (2n = 52) (Nomascus leucogenys leucogenys) as a model, we created a high-resolution map of the homologous regions between the gibbon and human. The positions of 100 synteny breakpoints relative to the assembled human genome were determined at a resolution of about 200 kb. Interestingly, 46% of the gibbon–human synteny breakpoints occur in regions that correspond to segmental duplications in the human lineage, indicating a common source of plasticity leading to a different outcome in the two species. Additionally, the full sequences of 11 gibbon BACs spanning evolutionary breakpoints reveal either segmental duplications or interspersed repeats at the exact breakpoint locations. No specific sequence element appears to be common among independent rearrangements. We speculate that the extraordinarily high level of rearrangements seen in gibbons may be due to factors that increase the incidence of chromosome breakage or fixation of the derivative chromosomes in a homozygous state.     

Origin of the polymorphism of the involucrin gene in Asians.

The involucrin gene, encoding a protein of the terminally differentiated keratinocyte, is polymorphic in the human. There is polymorphism of marker nucleotides a two positions in the coding region, and there are over eight polymorphic forms based on the number and kind of 10-codon tandem repeats in that part of the coding region most recently added in the human lineage. The involucrin alleles of Caucasians and Africans differ in both nucleotides and repeat patterns. We show that the involucrin alleles of East Asians (Chinese and Japanese) can be divided into two populations according to whether they possess the two marker nucleotides typical of Africans or Caucasians. The Asian population bearing Caucasian-type marker nucleotides has repeat patterns similar to those of Caucasians, whereas Asians bearing African-type marker nucleotides have repeat patterns that resemble those of Africans more than those of Caucasians. The existence of two populations of East Asian involucrin alleles gives support for the existence of a Eurasian stem lineage from which Caucasians and a part of the Asian population originated.     

The involucrin gene of the tree shrew: recent repeat additions and the relocation of cysteine codons

The coding region of the involucrin gene of Tupaia glis has been cloned and sequenced. It resembles the involucrin coding region of other non-anthropoid mammals in possessing a segment of related, short tandem repeats at a defined location, but in Tupaia, there has been recent serial duplication of a repeat into which a cysteine codon had earlier been introduced. As a result of the duplication, there is a total of as many as six cysteine codons in the segment of repeats, a number larger than for any other species yet examined. In Rattus there has been a comparable but independent addition of cysteine codons, and both Tupaia and Rattus have eliminated an otherwise conserved cysteine codon 75 located close to but outside the segment of repeats. In Tupaia, this elimination probably occurred by gene conversion. Also independently, the gene of Canis has added cysteine codons to the segment of repeats but has not yet lost cysteine 75. It is proposed that the gain and the loss of cysteine codons are parts of a multi-stage program of cysteine relocation.     

Expansion of mouse involucrin by intra-allelic repeat addition

Involucrin, loricrin and the small proline-rich proteins (SPRRs) are precursors of the cornified envelope of terminally differentiated keratinocytes. The genes for these proteins are closely linked on mouse chromosome 3. Each of the proteins is encoded by a single exon and is largely composed of a segment of short tandem repeats. No size polymorphism of either loricrin or the SPRRs was observed. In contrast, involucrin was found in at least eight polymorphic forms of different size with molecular weights ranging from 51 to 82kDa. Two classes of involucrin alleles were identified. Size polymorphism of involucrin has resulted from the recent expansion of the segment of repeats in one class of alleles, but not in the other. In expanding alleles, repeats were added at a precise location within the segment of repeats, in a 5'-to-3' direction. A study of a large number of allele-specific markers, located on both sides of the site of repeat addition, revealed no evidence for recombination between any of the alleles examined. Expansion of the segment of repeats of the gene for mouse involucrin must result from an intra-allelic process controlled by a cis-acting element, active in one class of alleles, and inactive in the other.     

Molecular evolutionary processes and conflicting gene trees: The hominoid case

Molecular evolutionary processes modify DNA over time, creating both newly derived substitutions shared by related descendant lineages (phylogenetic signal) and “false” similarities which confound phylogenetic reconstruction (homoplasy). However, some types of DNA regions, for example those containing tandem duplicate repeats, are preferentially subject to homoplasy-inducing processes such as sporadically occurring concerted evolution and DNA insertion/deletion. This added level of homoplasic “noise” can make DNA regions with repeats less reliable in phylogenetic reconstruction than those without repeats. Most molecular datasets which distinguish among African hominoids support a human-chimpanzee clade; the most notable exception is from the involucrin gene. However, phylogenetic resolution supporting a chimpanzee-gorilla clade is based entirely on involucrin DNA repeat regions. This is problematic because (1) involucrin repeats are difficult to align, and published alignments are contradictory; (2) involucrin repeats are subject to DNA insertion/deletion; (3) gorillas are polymorphic in that some do not have repeats reported to be synapomorphies linking chimpanzees and gorillas. Gene tree/species tree conflicts can occur due to the sorting of ancestrally polymorphic alleles during speciation. Because hominoid females transfer between groups, mitochondrial and nuclear gene flow occur to the same extent, and the probability of conflict between mitochondrial and nuclear gene trees is theoretically low. When hominoid intraspecific mitochondrial variability is taken into account [based on cytochrome oxidase subunit II (COII) gene sequences], humans and chimpanzees are most closely related, showing the same relative degree of separation from gorillas as when single individuals representing species are analyzed. Conflicting molecular phylogenies can be explained in terms of molecular evolutionary processes and sorting of ancient polymorphisms. This perspective can enhance our understanding of hominoid molecular phylogenies. © 1994 Wiley-Liss, Inc.     

The phylogeny of human globin genes investigated by the maximum parsimony method

Gene phylogenetic trees were constructed by the maximum parsimony method for various sets of ninety six globin chain amino acid sequences spanning plant and animal kingdoms. The method, executed by several computer programs, constructed ancestor and descendant globin messengers on tree topologies which required the least number of nucleotide replacements to account for the evolution of the globins. The human myoglobin-hemoglobin divergence was traced to a gene duplication which occurred either in the first vertebrates or earlier yet in the common ancestor of chordates and annelids, the alpha-beta divergence to a gene duplication in the common ancestor of teleosts and tetrapods, the gamma divergence from typical beta chains to a gene duplication in basal therian mammals, and the delta separation from beta to a duplication in the basal catarrhine primates. Evidence was provided by the globin phylogenies for the hominoid affinities of the gibbon and the close phyletic relationship of the African apes to man. Over the period of teleos-tetrapod divergence the globin messengers evolved at an average rate of 18.5 nucleotide replacements per 100 codons per 10⁸ years, a faster rate than most previous estimates. Very fast and very slow rates were encountered in different globin lineages and at different stages of descent, reducing the effectiveness of globins as molecular clocks. Rates increased with gene duplication and decreased after selection discovered useful specializations in the products of genes which had previously been freer to accept mutations. The early eutherian radiation was characterized by rapid rates of globin evolution, but the later hominoid radiation by extremely slow rates. This pattern was related to more complicated grades of internal organization evolving in human ancestors. The types of nucleotide replacements in the globin messengers over the long course of globin evolution did not seem indicative of any special mutational mechanisms.     

Variation in mitochondrial DNA and population structure of the Taipei treefrog Rhacophorus taipeianus in Taiwan

Mitochondrial DNA (mtDNA) from 40 samples of the Taipei treefrog Rhacophorus taipeianus collected from seven populations in Taiwan were sequenced to document the DNA sequence variation in anuran mtDNA and to elucidate the phylogeographic population structure in the Taipei treefrog. Sequences of 722–764 bases in length, including a 108-bp segment of the cytochrome b gene and a 614–656-bp D-loop segment, were obtained by direct sequencing using polymerase chain reaction (PCR). The variation in length was due to a 40-bp region that tandemly repeated four to five times in the D-loop region. The first repeat is the most conserved one among the five repeats because there are no variable sites in this repeat. Besides the 40-bp length variation, 28 positions in the 764-bp sequences are variable and distributed evenly in the cytochrome b gene fragment and D-loop region. Variation in the D-loop of the Taipei treefrog is comparable to those of other vertebrates. Two well-differentiated lineages (northern and central) differing by mean sequence divergence of 1.7% are identified and concordant with their geographic distributions. The two lineages are inferred to have split from a common ancestral population in the early Pleistocene. However, the interpopulation divergence of the northern lineage (< 0.33%) is apparently lower than that of the central lineage (1.11%), implying that the two lineages evolved independently and had different demographic histories after divergence. This study reveals that anuran D-loop has potential as a genetic marker in phylogenetic and population genetic analyses of anurans.     

The involucrin genes of the mouse and the rat: study of their shared repeats

The involucrin genes of the mouse (Mus musculus) and the rat (Rattus norvegicus) have been cloned and sequenced. The coding region of each gene contains, at site P, a segment of repeats homologous to that of other nonanthropoid mammals. In contrast to the repeats of species belonging to different mammalian orders, many individual repeats of the mouse and the rat can be matched. Both before and after the divergence of the two species, these repeats have been the site of systematic alterations in nucleotide sequence. One of the alterations is the correction of nucleotides of one repeat by those of another. Corrected nucleotides may be closely linked to flanking nucleotides that are uncorrected; the systematic correction process therefore appears to be due to gene conversion. There is a stretch of 18 reiterated CAGs in the segment of repeats of the Mus gene; most of these reiterations were introduced recently, supporting the idea that the gene was generated originally from poly CAG. An antiserum to a synthetic peptide encoded by the segment of repeats of the Mus gene reveals differentiation- specific expression of the gene in the epidermis.      

Systematic repeat addition at a precise location in the coding region of the involucrin gene of wild mice reveals their phylogeny

The involucrin gene encodes a protein of terminally differentiated keratinocytes. Its segment of repeats, which represents up to 80% of the coding region, is highly polymorphic in mouse strains derived from wild progenitors. Polymorphism includes nucleotide substitutions, but is most strikingly due to the recent addition of a variable number of repeats at a precise location within the segment of repeats. Each mouse taxon examined showed consistent and distinctive patterns of evolution of its variable region: very rapid changes in most M. m. domesticus alleles, slow changes in M. m. musculus, and complete arrest in M. spretus. We conclude that changes in the variable region are controlled by the genetic background. One of the M. m. domesticus alleles (DIK-L), which is of M. m. musculus origin, has undergone a recent repeat duplication typical of M. m. domesticus. This suggests that the genetic background controls repeat duplications through trans-acting factors. Because the repeat pattern differs in closely related murine taxa, involucrin reveals with greater sensitivity than random nucleotide substitutions the evolutionary relations of the mouse and probably of all murids.     

Estimation of hominoid ancestral population sizes under bayesian coalescent models incorporating mutation rate variation and sequencing errors

Estimation of population parameters for the common ancestors of humans and the great apes is important in understanding our evolutionary history. In particular, inference of population size for the human-chimpanzee common ancestor may shed light on the process by which the 2 species separated and on whether the human population experienced a severe size reduction in its early evolutionary history. In this study, the Bayesian method of ancestral inference of Rannala and Yang (2003. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics. 164:1645-1656) was extended to accommodate variable mutation rates among loci and random species-specific sequencing errors. The model was applied to analyze a genome-wide data set of approximately 15,000 neutral loci (7.4 Mb) aligned for human, chimpanzee, gorilla, orangutan, and macaque. We obtained robust and precise estimates for effective population sizes along the hominoid lineage extending back approximately 30 Myr to the cercopithecoid divergence. The results showed that ancestral populations were 5-10 times larger than modern humans along the entire hominoid lineage. The estimates were robust to the priors used and to model assumptions about recombination. The unusually low X chromosome divergence between human and chimpanzee could not be explained by variation in the male mutation bias or by current models of hybridization and introgression. Instead, our parameter estimates were consistent with a simple instantaneous process for human-chimpanzee speciation but showed a major reduction in X chromosome effective population size peculiar to the human-chimpanzee common ancestor, possibly due to selective sweeps on the X prior to separation of the 2 species.     

Evolutionary Breakpoints in the Gibbon Suggest Association between Cytosine Methylation and Karyotype Evolution

Gibbon species have accumulated an unusually high number of chromosomal changes since diverging from the common hominoid ancestor 15–18 million years ago. The cause of this increased rate of chromosomal rearrangements is not known, nor is it known if genome architecture has a role. To address this question, we analyzed sequences spanning 57 breaks of synteny between northern white-cheeked gibbons (Nomascus l. leucogenys) and humans. We find that the breakpoint regions are enriched in segmental duplications and repeats, with Alu elements being the most abundant. Alus located near the gibbon breakpoints (<150 bp) have a higher CpG content than other Alus. Bisulphite allelic sequencing reveals that these gibbon Alus have a lower average density of methylated cytosine that their human orthologues. The finding of higher CpG content and lower average CpG methylation suggests that the gibbon Alu elements are epigenetically distinct from their human orthologues. The association between undermethylation and chromosomal rearrangement in gibbons suggests a correlation between epigenetic state and structural genome variation in evolution.     

The involucrin genes of pig and dog: comparison of their segments of repeats with those of prosimians and higher primates

The involucrin genes of the dog and the pig have been cloned and sequenced. Like the corresponding genes of the prosimians, each contains a homologous segment of short tandem repeats at the same position in the coding region. However, the codon sequence of the repeats in the prosimians differs significantly from that of the nonprimate mammals. This evolution has been brought about by a combination of genetic modifications (selective deletions, mutations, and gene conversions). In the anthropoids, this segment of repeats was replaced by a modern one differing in location, sequence, and repeat length. In several of its properties the modern segment has continued the prosimian trend away from the nonprimates. The overall direction of the evolution of this segment has therefore been maintained even though there have been sudden changes in the evolutionary processes acting on the gene.     

A hominoid-specific nuclear insertion of the mitochondrial D-loop: implications for reconstructing ancestral mitochondrial sequences

A nuclear integration of a mitochondrial control region sequence on human chromosome 9 has been isolated. PCR analyses with primers specific for the respective insertion-flanking nuclear regions showed that the insertion took place on the lineage leading to Hominoidea (gibbon, orangutan, gorilla, chimpanzee, and human) after the Old World monkey-Hominoidea split. The sequences of the control region integrations were determined for humans, chimpanzees, gorillas, orangutans, and siamangs. These sequences were then used to construct phylogenetic trees with different methods, relating them with several hominoid, Old Work monkey, and New World monkey mitochondrial control region sequences. Applying maximum-likelihood, neighbor-joining, and parsimony algorithms, the insertion clade was attached to the branch leading to the hominoid mitochondrial sequences as expected from the PCR-determined presence/absence of this integration. An unexpected long branch leading to the internal node that connects all insertion sequences was observed for the different phylogeny reconstruction procedures. This finding is not totally compatible with the lower evolutionary rate in the nucleus than in the mitochondrial compartment. We determined the unambiguous substitutions on the branch leading to the most recent common ancestor (MRCA) of the mitochondrial inserts according to the parsimony criterium. We propose that they are unlikely to have been caused by damage of the transposing nucleic acid and that they are probably due to a change in the evolutionary mode after the transposition.