Primate phylogenomics uncovers multiple rapid radiations and ancient interspecific introgression.

Our understanding of the evolutionary history of primates is undergoing continual revision due to ongoing genome sequencing efforts. Bolstered by growing fossil evidence, these data have led to increased acceptance of once controversial hypotheses regarding phylogenetic relationships, hybridization and introgression, and the biogeographical history of primate groups. Among these findings is a pattern of recent introgression between species within all major primate groups examined to date, though little is known about introgression deeper in time. To address this and other phylogenetic questions, here, we present new reference genome assemblies for 3 Old World monkey (OWM) species: Colobus angolensis ssp. palliatus (the black and white colobus), Macaca nemestrina (southern pig-tailed macaque), and Mandrillus leucophaeus (the drill). We combine these data with 23 additional primate genomes to estimate both the species tree and individual gene trees using thousands of loci. While our species tree is largely consistent with previous phylogenetic hypotheses, the gene trees reveal high levels of genealogical discordance associated with multiple primate radiations. We use strongly asymmetric patterns of gene tree discordance around specific branches to identify multiple instances of introgression between ancestral primate lineages. In addition, we exploit recent fossil evidence to perform fossil-calibrated molecular dating analyses across the tree. Taken together, our genome-wide data help to resolve multiple contentious sets of relationships among primates, while also providing insight into the biological processes and technical artifacts that led to the disagreements in the first place.

The Origin and Population History of the Endangered Golden Snub-Nosed Monkey (Rhinopithecus roxellana).

The origin and population history of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) remain largely unavailable and/or controversial. We here integrate analyses of multiple genomic markers, including mitochondrial (mt) genomes, Y-chromosomes, and autosomes of 54 golden monkey individuals from all three geographic populations (SG, QL, and SNJ). Our results reveal contrasting population structures. Mt analyses suggest a division of golden monkeys into five lineages: one in SNJ, two in SG, and two in QL. One of the SG lineages (a mixed SG/QL lineage) is basal to all other lineages. In contrast, autosomal analyses place SNJ as the most basal lineage and identify one QL and three SG lineages. Notably, Y-chromosome analyses bear features similar to mt analyses in placing the SG/QL-mixed lineage as the first diverging lineage and dividing SG into two lineages, while resembling autosomal analyses in identifying one QL lineage. We further find bidirectional gene flow among all three populations at autosomal loci, while asymmetric gene flow is suggested at mt genomes and Y-chromosomes. We propose that different population structures and gene flow scenarios are the result of sex-linked differences in the dispersal pattern of R. roxellana. Moreover, our demographic simulation analyses support an origin hypothesis suggesting that the ancestral R. roxellana population was once widespread and then divided into SNJ and non-SNJ (SG and QL) populations. This differs from previous mt-based "mono-origin (SG is the source population)" and "multiorigin (SG is a fusion of QL and SNJ)" hypotheses. We provide a detailed and refined scenario for the origin and population history of this endangered primate species, which has a broader significance for Chinese biogeography. In addition, this study highlights the importance to investigate multiple genomic markers with different modes of inheritance to trace the complete evolutionary history of a species, especially for those exhibiting differential or mixed patterns of sex dispersal.

Nomenclature report 2019: major histocompatibility complex genes and alleles of Great and Small Ape and Old and New World monkey species.

The major histocompatibility complex (MHC) is central to the innate and adaptive immune responses of jawed vertebrates. Characteristic of the MHC are high gene density, gene copy number variation, and allelic polymorphism. Because apes and monkeys are the closest living relatives of humans, the MHCs of these non-human primates (NHP) are studied in depth in the context of evolution, biomedicine, and conservation biology. The Immuno Polymorphism Database (IPD)-MHC NHP Database (IPD-MHC NHP), which curates MHC data of great and small apes, as well as Old and New World monkeys, has been upgraded. The curators of the database are responsible for providing official designations for newly discovered alleles. This nomenclature report updates the 2012 report, and summarizes important nomenclature issues and relevant novel features of the IPD-MHC NHP Database.

Alu master copies serve as the drivers of differential SINE transposition in recent primate genomes.

Alu elements, averaging ~300bp in length, are a family of primate-specific short intersperse nuclear elements (SINEs) with more than one million copies and contributing to ~11% of primate genomes. Despite mostly being shared among primates, our recent study revealed highly differential recent Alu transposition among the genomes of primates from Hominidae and Cercopithecidae families. To understand the underlying mechanism, we analyzed six primate genomes and revealed species- and lineage-specific Alu profile exclusively defined by AluY composition. Among all Alus from the 6 genomes, we identified 5401 Alu master copies with 99% being from the AluY subfamily. The numbers of Alu master copies are positively correlated to the number of AluY elements in the genomes with the baboon genome having the largest number of most recent Alu master copies at high activities, while the crab-eating macaque genome having a low number of Alu master copies with low activity. Furthermore, the expression level of Alu master copies is positively correlated with their transposition activity. Our results support the concept that Alu transposition in primate genomes is driven by a small number of master copies, the number and relative activity of which contribute to the differential Alu transposition in recent primate genomes.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene.

BACKGROUND: The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates. RESULTS: The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination. CONCLUSIONS: We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

The evolution of S100A7 in primates: a model of concerted and birth-and-death evolution.

The human S100A7 resides in the epidermal differentiation complex (EDC) and has been described as a key effector of innate immunity. In humans, there are five S100A7 genes located in tandem-S100A7A, S100A7P1, S100AL2, S100A7, and S100AP2. The presence of several retroelements in the S100A7A/S100A7P1 and S100A7/S100A7P2 clusters suggests that these genes were originated from a duplication around ~ 35 million years ago, during or after the divergence of Platyrrhini and Catarrhini primates. To test this hypothesis, and taking advantage of the high number of genomic sequences available in the public databases, we retrieved S100A7 gene sequences of 12 primates belonging to the Cercopithecoidea and Hominoidea (Catarrhini species). Our results support the duplication theory, with at least one gene of each cluster being identified in both Cercopithecoidea and Hominoidea species. Moreover, given the presence of an ongoing gene conversion event between S100A7 and S100A7A, a high rate of mutation in S100A7L2 and the presence of pseudogenes, we proposed a model of concerted and birth-and-death evolution to explain the evolution of S100A7 gene family. Indeed, our results suggest that S100A7L2 most likely suffered a neofunctionalization in the Catarrhini group. Being S100A7 a major protein in innate defense, we believe that our findings could open new doors in the study of this gene family in immune system.

Superoxide dismutase 1 is positively selected to minimize protein aggregation in great apes.

Positive (adaptive) selection has recently been implied in human superoxide dismutase 1 (SOD1), a highly abundant antioxidant protein with energy signaling and antiaging functions, one of very few examples of direct selection on a human protein product (exon); the molecular drivers of this selection are unknown. We mapped 30 extant SOD1 sequences to the recently established mammalian species tree and inferred ancestors, key substitutions, and signatures of selection during the protein's evolution. We detected elevated substitution rates leading to great apes (Hominidae) at ~1 per 2 million years, significantly higher than in other primates and rodents, although these paradoxically generally evolve much faster. The high evolutionary rate was partly due to relaxation of some selection pressures and partly to distinct positive selection of SOD1 in great apes. We then show that higher stability and net charge and changes at the dimer interface were selectively introduced upon separation from old world monkeys and lesser apes (gibbons). Consequently, human, chimpanzee and gorilla SOD1s have a net charge of -6 at physiological pH, whereas the closely related gibbons and macaques have -3. These features consistently point towards selection against the malicious aggregation effects of elevated SOD1 levels in long-living great apes. The findings mirror the impact of human SOD1 mutations that reduce net charge and/or stability and cause ALS, a motor neuron disease characterized by oxidative stress and SOD1 aggregates and triggered by aging. Our study thus marks an example of direct selection for a particular chemical phenotype (high net charge and stability) in a single human protein with possible implications for the evolution of aging.

The Relationship between the (In-)Stability of NORs and Their Chromosomal Location: The Example of Cercopithecidae and a Short Review of Other Primates.

Amongst Cercopithecidae, the species of the Cercopithecini tribe underwent a very active chromosome evolution, principally by fissions, which increased their chromosome number up to 72. In contrast, all the species of Papionini have fairly similar karyotypes with 42 chromosomes. In animals, nucleolus organizer regions (NORs) are generally considered as instable structures, which frequently vary in size, number, and location at both infra- and interspecific levels. Although in Cercopithecinae the NORs, involved in breaks, exchanges, and translocations, behave like fragile sites in somatic cells, their number and location appear to be very stable between species. Fluorescence in situ hybridization of a 28S rDNA probe on metaphase chromosomes displayed a unique interstitial location in either an acrocentric pair (in 12 species of Cercopithecini) or a metacentric pair (in 6 species of Papionini). A non-exhaustive survey of literature data on NOR location in other primates shows that numerical variations of the NORs principally depend on their location: most multiple NORs are in terminal positions, while almost all unique NORs are in interstitial positions. We propose that this correlation is the consequence of the selection against gametic imbalances involving the chromosomal material distal to the NORs, which is effective when they are interstitially, but not terminally, located. Thus, the consequences of the interstitial NOR instability for reproduction are essentially limited to their size variations, as observed in Cercopithecidae.

Differential evolutionary fate of an ancestral primate endogenous retrovirus envelope gene, the EnvV syncytin, captured for a function in placentation.

Syncytins are envelope genes of retroviral origin that have been co-opted for a role in placentation. They promote cell-cell fusion and are involved in the formation of a syncytium layer--the syncytiotrophoblast--at the materno-fetal interface. They were captured independently in eutherian mammals, and knockout mice demonstrated that they are absolutely required for placenta formation and embryo survival. Here we provide evidence that these "necessary" genes acquired "by chance" have a definite lifetime with diverse fates depending on the animal lineage, being both gained and lost in the course of evolution. Analysis of a retroviral envelope gene, the envV gene, present in primate genomes and belonging to the endogenous retrovirus type V (ERV-V) provirus, shows that this captured gene, which entered the primate lineage >45 million years ago, behaves as a syncytin in Old World monkeys, but lost its canonical fusogenic activity in other primate lineages, including humans. In the Old World monkeys, we show--by in situ analyses and ex vivo assays--that envV is both specifically expressed at the level of the placental syncytiotrophoblast and fusogenic, and that it further displays signs of purifying selection based on analysis of non-synonymous to synonymous substitution rates. We further show that purifying selection still operates in the primate lineages where the gene is no longer fusogenic, indicating that degeneracy of this ancestral syncytin is a slow, lineage-dependent, and multi-step process, in which the fusogenic activity would be the first canonical property of this retroviral envelope gene to be lost.

Significance of the evolutionary α1,3-galactosyltransferase (GGTA1) gene inactivation in preventing extinction of apes and old world monkeys.

The α1,3-galactosyltransferase (α1,3GT or GGTA1) gene displays unique evolutionary characteristics. This gene appeared early in mammalian evolution and is absent in other vertebrates. The α1,3GT gene is active in marsupials, nonprimate placental mammals, lemurs (prosimians) and New World monkeys, encoding the α1,3GT enzyme that synthesizes a carbohydrate antigen called "α-gal epitope." The α-gal epitope is present in large numbers on cell membrane glycolipids and glycoproteins. The α1,3GT gene was inactivated in ancestral Old World monkeys and apes by frameshift single-base deletions forming premature stop codons. Because of this gene inactivation, humans, apes, and Old World monkeys lack α-gal epitopes and naturally produce an antibody called the "anti-Gal antibody" which binds specifically to α-gal epitopes and which is the most abundant antibody in humans. The evolutionary event that resulted in the inactivation of the α1,3GT gene in ancestral Old World primates could have been mediated by a pathogen endemic to Eurasia-Africa landmass that exerted pressure for selection of primate populations lacking the α-gal epitope. Once the α-gal epitope was eliminated, primates could produce the anti-Gal antibody, possibly as means of defense against pathogens expressing this epitope. This assumption is supported by the fossil record demonstrating an almost complete extinction of apes in the late Miocene and failure of Old World monkeys to radiate into multiple species before that period. A present outcome of this evolutionary event is the anti-Gal-mediated rejection of mammalian xenografts expressing α-gal epitopes in humans, apes, and Old World monkeys.

The ABO blood group is a trans-species polymorphism in primates.

The ABO histo-blood group, the critical determinant of transfusion incompatibility, was the first genetic polymorphism discovered in humans. Remarkably, ABO antigens are also polymorphic in many other primates, with the same two amino acid changes responsible for A and B specificity in all species sequenced to date. Whether this recurrence of A and B antigens is the result of an ancient polymorphism maintained across species or due to numerous, more recent instances of convergent evolution has been debated for decades, with a current consensus in support of convergent evolution. We show instead that genetic variation data in humans and gibbons as well as in Old World monkeys are inconsistent with a model of convergent evolution and support the hypothesis of an ancient, multiallelic polymorphism of which some alleles are shared by descent among species. These results demonstrate that the A and B blood groups result from a trans-species polymorphism among distantly related species and has remained under balancing selection for tens of millions of years-to date, the only such example in hominoids and Old World monkeys outside of the major histocompatibility complex.

Divergence and diversity of ULBP2 genes in rhesus and cynomolgus macaques.

Non-human primates such as rhesus macaque and cynomolgus macaque are important animals for medical research fields and they are classified as Old World monkey, in which genome structure is characterized by gene duplications. In the present study, we investigated polymorphisms in two genes for ULBP2 molecules that are ligands for NKG2D. A total of 15 and 11 ULBP2.1 alleles and 11 and 10 ULBP2.2 alleles were identified in rhesus macaques and cynomolgus macaques, respectively. Nucleotide sequences of exons for extra cellular domain were highly polymorphic and more than 70 % were non-synonymous variations in both ULBP2.1 and ULBP2.2. In addition, phylogenetic analyses revealed that the ULBP2.2 was diverged from a branch of ULBP2.1 along with ULBP2s of higher primates. Moreover, when 3D structural models were constructed for the rhesus ULBP2 molecules, residues at presumed contact sites with NKG2D were polymorphic in ULBP2.1 and ULBP2.2 in the rhesus macaque and cynomolgus macaque, respectively. These observations suggest that amino acid replacements at the interaction sites with NKG2D might shape a specific nature of ULBP2 molecules in the Old World monkeys.

Evolution of primate α and θ defensins revealed by analysis of genomes.

Defensins are endogenous peptides with cysteine-rich antimicrobial ability that contribute to host defence against bacterial, fungal and viral infections. There are three subfamilies of defensins in primates: α, ß and θ-defensins. α-defensins are most present in neutrophils and Paneth cells; ß-defensins are involved in protecting the skin and the mucous membranes of the respiratory, genitourinary and gastrointestinal tracts; and θ-defensins are physically distinguished as the only known fully-cyclic peptides of animal origin, which are first isolated from rhesus macaques. All three kinds of defensins have six conserved cysteines, three intramolecular disulfide bonds, a net positive charge, and ß-sheet regions. α and θ-defensins are closely related, comparative amino acid sequences showed that the difference between them is that θ-defensins have an additional stop codon limits the initial defensin domain peptides to 12 residues. Humans, chimpanzees and gorillas do not produce θ-defensin peptides due to a premature stop codon present in the signal sequence of all θ-defensin pseudogenes. By using comprehensive computational searches, here we report the discovery of complete repertoires of the α and θ-defensin gene family in ten primate species. Consistent with previous studies, our phylogenetic analyses showed all primate θ-defensins evident formed one distinct clusters evolved from α-defensins. ß-defensins are ancestors of both α and θ-defensins. Human has two copies of DEFA1 and DEFT1P, and two extra DEFA3 and DEFA10P genes compared with gorilla. As different primates inhabit in quite different ecological niches, the production of species-specific α and θ-defensins and these highly evolved θ-defensins in old world monkeys would presumably allow them to better respond to the specific microbial challenges that they face.

Sequencing of rhesus macaque Y chromosome clarifies origins and evolution of the DAZ (Deleted in AZoospermia) genes.

Studies of Y chromosome evolution often emphasize gene loss, but this loss has been counterbalanced by addition of new genes. The DAZ genes, which are critical to human spermatogenesis, were acquired by the Y chromosome in the ancestor of Old World monkeys and apes. We and our colleagues recently sequenced the rhesus macaque Y chromosome, and comparison of this sequence to human and chimpanzee enables us to reconstruct much of the evolutionary history of DAZ. We report that DAZ arrived on the Y chromosome about 38 million years ago via the transposition of at least 1.1 megabases of autosomal DNA. This transposition also brought five additional genes to the Y chromosome, but all five genes were subsequently lost through mutation or deletion. As the only surviving gene, DAZ experienced extensive restructuring, including intragenic amplification and gene duplication, and has been the target of positive selection in the chimpanzee lineage. Editor's suggested further reading in BioEssays Should Y stay or should Y go: The evolution of non-recombining sex chromosomes Abstract.

Reduction of bitter taste receptor gene family in folivorous colobine primates relative to omnivorous cercopithecine primates.

Bitter taste perception is important in preventing animals from ingesting potentially toxic compounds. Whole-genome assembly (WGA) data have revealed that bitter taste receptor genes (TAS2Rs) comprise a multigene family with dozens of intact and disrupted genes in primates. However, publicly available WGA data are often incomplete, especially for multigene families. In this study, we employed a targeted capture (TC) approach specifically probing TAS2Rs for ten species of cercopithecid primates with diverse diets, including eight omnivorous cercopithecine species and two folivorous colobine species. We designed RNA probes for all TAS2Rs that we modeled to be intact in the common ancestor of cercopithecids ("ancestral-cercopithecid TAS2R gene set"). The TC was followed by short-read and high-depth massive-parallel sequencing. TC retrieved more intact TAS2R genes than found in WGA databases. We confirmed a large number of gene "births" at the common ancestor of cercopithecids and found that the colobine common ancestor and the cercopithecine common ancestor had contrasting trajectories: four gene "deaths" and three gene births, respectively. The number of intact TAS2R genes was markedly reduced in colobines (25-28 detected via TC and 20-26 detected via WGA analysis) as compared with cercopithecines (27-36 via TC and 19-30 via WGA). Birth or death events occurred at almost every phylogenetic-tree branch, making the composition of intact genes variable among species. These results show that evolutionary change in intact TAS2R genes is a complex process, refute a simple general prediction that herbivory favors more TAS2R genes, and have implications for understanding dietary adaptations and the evolution of detoxification abilities.

Functional divergence of the rapidly evolving miR-513 subfamily in primates.

BACKGROUND: The miR-513 subfamily belongs to an X-linked primate-specific miR506-514 cluster. Across primate species, there have been several duplication events and different species each possess a variety of miR-513 copies, indicating it underwent rapid evolution. Evidence suggests that this subfamily is preferentially expressed in the testis, but otherwise, to date, the evolutionary history and functional significance of this miRNA subfamily has remained largely unexplored. RESULTS: We analyzed the evolutionary pattern of gene duplications and their functional consequence for the miR-513 subfamily in primates. Sequence comparisons showed that the duplicated copies of miR-513 were derived from transposable element (MER91C). Moreover, duplication events of the miR-513 subfamily seem to have occurred independently in Platyrrhini (New World monkeys) and Catarrhini (Old World monkeys, apes and humans) after they diverged. Different copies of the miR-513 subfamily (miR-513a/b/c) have different seed sequences, due to after-duplication sequence divergences, which eventually led to functional divergences. The results of functional assays indicated that miR-513b could inhibit the expression of its target gene, the down-regulator of transcription 1 (DR1) at both the mRNA and protein levels. In the developing testis of rhesus macaques, we observed a temporal coupling of expression levels between miR-513b and DR1, suggesting that miR-513b could affect male sexual maturation by negatively regulating the development-stage related functioning of DR1. CONCLUSIONS: The miR-513 subfamily underwent multiple independent gene duplications among five different lineages of primates. The after-duplication sequence divergences among the different copies of miR-513 led to functional divergence of these copies in primates.

The evolution and expression of the snaR family of small non-coding RNAs.

We recently identified the snaR family of small non-coding RNAs that associate in vivo with the nuclear factor 90 (NF90/ILF3) protein. The major human species, snaR-A, is an RNA polymerase III transcript with restricted tissue distribution and orthologs in chimpanzee but not rhesus macaque or mouse. We report their expression in human tissues and their evolution in primates. snaR genes are exclusively in African Great Apes and some are unique to humans. Two novel families of snaR-related genetic elements were found in primates: CAS (catarrhine ancestor of snaR), limited to Old World Monkeys and apes; and ASR (Alu/snaR-related), present in all monkeys and apes. ASR and CAS appear to have spread by retrotransposition, whereas most snaR genes have spread by segmental duplication. snaR-A and snaR-G2 are differentially expressed in discrete regions of the human brain and other tissues, notably including testis. snaR-A is up-regulated in transformed and immortalized human cells, and is stably bound to ribosomes in HeLa cells. We infer that snaR evolved from the left monomer of the primate-specific Alu SINE family via ASR and CAS in conjunction with major primate speciation events, and suggest that snaRs participate in tissue- and species-specific regulation of cell growth and translation.

Evolutionary analysis of the short-type peptidoglycan-recognition protein gene (PGLYRP1) in primates.

Short-type peptidoglycan (PGN)-recognition protein 1 (PGLYRP1), an innate immunity protein that directly breaks down the structure of microbial cell wall PGNs, plays an important role both in antibacterial defenses and several inflammatory diseases. To explore the adaptive evolution of the PGLYRP1 gene in primates and provide insight into the function of this antibacterial protein, we sequenced the entire PGLYRP1 gene from Macaca thibetana and Rhinopithecus roxellana, identified the corresponding sequences from the draft genome of 8 other primates, including humans, and conducted related statistical analyses. Homology analysis showed that the identity of nucleotide and deduced amino acid sequences of PGLYRP1 among 10 primates ranged from 82.0 to 99.0% and 74.5 to 98.5%, respectively. The R value (transition/transversion) and disparity index per site also presented relatively low-base composition biases. Selective pressure analysis for the PGLYRP1 sequences among major primates revealed that both the whole gene and the substructure of PGLYRP1 are under strong purifying selection at similar levels of selective pressure among 6 major primate lineages (human, great ape, lesser ape, Old World monkey, New World monkey, and prosimian monkey). Using the Bayes empirical Bayes procedure, we also detected 2 positively selected codons (121L and 141T sites) that are independent of PGN-binding and PGLYRP-specific regions, implying 2 potential key sites for the functional effect of the PGLYRP1 protein. These results demonstrated that PGLYRP1 was highly conserved at the molecular level and subjected to strong functional constraints during primate evolution.

Cytogenetic studies of small ape (Hylobatidae) chromosomes.

Each genus of small apes has a highly distinctive karyotype (karyomorph) at every level of cytogenetic analysis. Early workers using classical staining and banding had problems integrating the karyolocial data with that of other primates. Chromosome painting allowed syntenic homology maps to be constructed for each of the four karyomorphs (2n = 38, 44, 50 and 52). They revealed that the great apes and Old World monkeys had strongly conserved karyotypes while those of small apes were highly rearranged. However, they provided contradictory phylogenetic results to other bio-molecular tree of small ape evolution. More recently BAC-FISH investigations using a panel of about 900 BACs defined each breakpoint by spanning or flanking BAC clones The syntenic map was refined and now includes small segments of homology which had previously gone undected, marker order (synteny block orientation) and the location of ancestral and Evolutionarily New Centromeres. However, the BAC-FISH data similar to other biomolecular methods used up to now could not resolve the phylogenetic tree of hylobatids. These difficulties may be explained by the rapid divergence of crown hylobatids, reticulate evolution and incomplete lineage sorting. The lack of significant cytogenetic landmarks at the nodes of the gibbon tree could indicate that chromosomal rearrangements did not play a primary role in hylobatid speciation.

Framework of the Alu Subfamily Evolution in the Platyrrhine Three-Family Clade of Cebidae, Callithrichidae, and Aotidae.

The history of Alu retroposons has been choreographed by the systematic accumulation of inherited diagnostic nucleotide substitutions to form discrete subfamilies, each having a distinct nucleotide consensus sequence. The oldest subfamily, AluJ, gave rise to AluS after the split between Strepsirrhini and what would become Catarrhini and Platyrrhini. The AluS lineage gave rise to AluY in catarrhines and to AluTa in platyrrhines. Platyrrhine Alu subfamilies Ta7, Ta10, and Ta15 were assigned names based on a standardized nomenclature. However, with the subsequent intensification of whole genome sequencing (WGS), large scale analyses to characterize Alu subfamilies using the program COSEG identified entire lineages of subfamilies simultaneously. The first platyrrhine genome with WGS, the common marmoset (Callithrix jacchus; [caljac3]), resulted in Alu subfamily names sf0 to sf94 in an arbitrary order. Although easily resolved by alignment of the consensus sequences, this naming convention can become increasingly confusing as more genomes are independently analyzed. In this study, we reported Alu subfamily characterization for the platyrrhine three-family clade of Cebidae, Callithrichidae, and Aotidae. We investigated one species/genome from each recognized family of Callithrichidae and Aotidae and of both subfamilies (Cebinae and Saimiriinae) of the family Cebidae. Furthermore, we constructed a comprehensive network of Alu subfamily evolution within the three-family clade of platyrrhines to provide a working framework for future research. Alu expansion in the three-family clade has been dominated by AluTa15 and its derivatives.