Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal.

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale-multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8-15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades.

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

Mutations in CERS3 cause autosomal recessive congenital ichthyosis in humans.

Autosomal recessive congenital ichthyosis (ARCI) is a rare genetic disorder of the skin characterized by abnormal desquamation over the whole body. In this study we report four patients from three consanguineous Tunisian families with skin, eye, heart, and skeletal anomalies, who harbor a homozygous contiguous gene deletion syndrome on chromosome 15q26.3. Genome-wide SNP-genotyping revealed a homozygous region in all affected individuals, including the same microdeletion that partially affects two coding genes (ADAMTS17, CERS3) and abolishes a sequence for a long non-coding RNA (FLJ42289). Whereas mutations in ADAMTS17 have recently been identified in autosomal recessive Weill-Marchesani-like syndrome in humans and dogs presenting with ophthalmologic, cardiac, and skeletal abnormalities, no disease associations have been described for CERS3 (ceramide synthase 3) and FLJ42289 so far. However, analysis of additional patients with non-syndromic ARCI revealed a splice site mutation in CERS3 indicating that a defect in ceramide synthesis is causative for the present skin phenotype of our patients. Functional analysis of patient skin and in vitro differentiated keratinocytes demonstrated that mutations in CERS3 lead to a disturbed sphingolipid profile with reduced levels of epidermis-specific very long-chain ceramides that interferes with epidermal differentiation. Taken together, these data present a novel pathway involved in ARCI development and, moreover, provide the first evidence that CERS3 plays an essential role in human sphingolipid metabolism for the maintenance of epidermal lipid homeostasis.

A novel system of polymorphic and diverse NK cell receptors in primates.

There are two main classes of natural killer (NK) cell receptors in mammals, the killer cell immunoglobulin-like receptors (KIR) and the structurally unrelated killer cell lectin-like receptors (KLR). While KIR represent the most diverse group of NK receptors in all primates studied to date, including humans, apes, and Old and New World monkeys, KLR represent the functional equivalent in rodents. Here, we report a first digression from this rule in lemurs, where the KLR (CD94/NKG2) rather than KIR constitute the most diverse group of NK cell receptors. We demonstrate that natural selection contributed to such diversification in lemurs and particularly targeted KLR residues interacting with the peptide presented by MHC class I ligands. We further show that lemurs lack a strict ortholog or functional equivalent of MHC-E, the ligands of non-polymorphic KLR in "higher" primates. Our data support the existence of a hitherto unknown system of polymorphic and diverse NK cell receptors in primates and of combinatorial diversity as a novel mechanism to increase NK cell receptor repertoire.

A gene conversion hotspot in the human growth hormone (GH1) gene promoter.

To assess the evolutionary importance of nonallelic (or interlocus) gene conversion for the highly polymorphic human growth hormone (GH1) gene promoter, sequence variation in this region was studied in four different ethnic groups. For 14 SNPs in the proximal GH1 promoter (535 bp), 60 different haplotypes were observed in 577 individuals (156 Britons, 116 Spaniards, 163 West-Africans, 142 Asians). Using a novel coalescence-based statistical test, significant evidence was found in the British, Spanish, and African groups for GH1 having acted as an acceptor of gene conversion, with at least one of the four paralogous GH gene promoters serving as the donor (and specifically GH2 in the Britons and Spaniards). The average gene conversion tract length was estimated to be 84 bp. A gene conversion hotspot was identified, spanning the GH1 transcriptional initiation site (positions -6 to +25). Although these findings serve to highlight the importance of gene conversion for the recent evolution of the human GH1 promoter, its relative frequency does not appear to be related simply to the presence of specific DNA sequence motifs or secondary structures, the degree of homology between GH paralogs, the distance between them, or their transcriptional orientation. The GH1 promoter was also found to be highly polymorphic in chimpanzee but not in macaque. This may reflect the lower degree of pair-wise similarity between the GH1 promoter and its paralogs in macaque (mean, 92.0%) as compared to chimpanzee (93.5%) and human (94.0%), and hence provides further support for the idea of a threshold (perhaps around 92%) below which gene conversion is reduced or abolished.

Novel cytomegaloviruses in free-ranging and captive great apes: phylogenetic evidence for bidirectional horizontal transmission.

Wild great apes often suffer from diseases of unknown aetiology. This is among the causes of population declines. Because human cytomegalovirus (HCMV) is an important pathogen, especially in immunocompromised individuals, a search for cytomegaloviruses (CMVs) in deceased wild and captive chimpanzees, gorillas and orang-utans was performed. By using a degenerate PCR targeting four conserved genes (UL54-UL57), several distinct, previously unrecognized CMVs were found for each species. Sequences of up to 9 kb were determined for ten novel CMVs, located in the UL54-UL57 block. A phylogenetic tree was inferred for the ten novel CMVs, the previously characterized chimpanzee CMV, HCMV strains and Old World and New World monkey CMVs. The primate CMVs fell into four clades, containing New World monkey, Old World monkey, orang-utan and human CMVs, respectively, plus two clades that each contained both chimpanzee and gorilla isolates (termed CG1 and CG2). The tree loci of the first four clades mirrored those for their respective hosts in the primate tree, suggesting that these CMV lineages arose through cospeciation with host lineages. The CG1 and CG2 loci corresponded to those of the gorilla and chimpanzee hosts, respectively. This was interpreted as indicating that CG1 and CG2 represented CMV lineages that had arisen cospeciationally with the gorilla and chimpanzee lineages, respectively, with subsequent transfer within each clade between the host genera. Divergence dates were estimated and found to be consistent with overall cospeciational development of major primate CMV lineages. However, CMV transmission between chimpanzees and gorillas in both directions has also occurred.

Molecular diagnosis in dermatopathology: what makes sense, and what doesn't.

Molecular techniques have provided us with a wealth of information about biological events in healthy individual, and improved tremendously our understanding about the pathogenesis of a huge variety of cutaneous diseases. Those methods have originally been invented to support basic scientific investigations on a molecular level and are translated increasingly into sophisticated diagnostic tools changing the classic paradigm of diagnostic pathology; among them are immunohistochemistry (IHC), polymerase chain reaction (PCR), G-banding, loss of heterozygosity, fluorescence in situ hybridization (FISH), chromogen in situ hybridization (CISH), comparative genomic hybridization on chromosomes and microarray technology. Some of them such as IHC and PCR have already been standardized to a level that allows its utility in daily routine diagnostics for several dermatological diseases. For others like array-based technologies, their optimal indications await to be fully determined. These ancillary methods have the great potential to contribute important new information to challenging cases, and will help to improve diagnostic accuracy particularly in cases in which conventional histopathology is ambiguous. Thus, they will broaden our armamentarium for diagnostic pathology. Herein, some key techniques will be reviewed and their applicability towards the diagnosis of dermatological diseases critically discussed.

Evolutionary analysis of the highly dynamic CHEK2 duplicon in anthropoids.

BACKGROUND: Segmental duplications (SDs) are euchromatic portions of genomic DNA (> or = 1 kb) that occur at more than one site within the genome, and typically share a high level of sequence identity (>90%). Approximately 5% of the human genome is composed of such duplicated sequences. Here we report the detailed investigation of CHEK2 duplications. CHEK2 is a multiorgan cancer susceptibility gene encoding a cell cycle checkpoint kinase acting in the DNA-damage response signalling pathway. The continuous presence of the CHEK2 gene in all eukaryotes and its important role in maintaining genome stability prompted us to investigate the duplicative evolution and phylogeny of CHEK2 and its paralogs during anthropoid evolution. RESULTS: To study CHEK2 duplicon evolution in anthropoids we applied a combination of comparative FISH and in silico analyses. Our comparative FISH results with a CHEK2 fosmid probe revealed the single-copy status of CHEK2 in New World monkeys, Old World monkeys and gibbons. Whereas a single CHEK2 duplication was detected in orangutan, a multi-site signal pattern indicated a burst of duplication in African great apes and human. Phylogenetic analysis of paralogous and ancestral CHEK2 sequences in human, chimpanzee and rhesus macaque confirmed this burst of duplication, which occurred after the radiation of orangutan and African great apes. In addition, we used inter-species quantitative PCR to determine CHEK2 copy numbers. An amplification of CHEK2 was detected in African great apes and the highest CHEK2 copy number of all analysed species was observed in the human genome. Furthermore, we detected variation in CHEK2 copy numbers within the analysed set of human samples. CONCLUSION: Our detailed analysis revealed the highly dynamic nature of CHEK2 duplication during anthropoid evolution. We determined a burst of CHEK2 duplication after the radiation of orangutan and African great apes and identified the highest CHEK2 copy number in human. In conclusion, our analysis of CHEK2 duplicon evolution revealed that SDs contribute to inter-species variation. Furthermore, our qPCR analysis led us to presume CHEK2 copy number variation in human, and molecular diagnostics of the cancer susceptibility gene CHEK2 inside the duplicated region might be hampered by the individual-specific set of duplicons.

Evolution of the DAZ gene and the AZFc region on primate Y chromosomes.

BACKGROUND: The Azoospermia Factor c (AZFc) region of the human Y chromosome is a unique product of segmental duplication. It consists almost entirely of very long amplicons, represented by different colors, and is frequently deleted in subfertile men. Most of the AZFc amplicons have high sequence similarity with autosomal segments, indicating recent duplication and transposition to the Y chromosome. The Deleted in Azoospermia (DAZ) gene within the red-amplicon arose from an ancestral autosomal DAZ-like (DAZL) gene. It varies significantly between different men regarding to its copy number and the numbers of RNA recognition motif and DAZ repeat it encodes. We used Southern analyses to study the evolution of DAZ and AZFc amplicons on the Y chromosomes of primates. RESULTS: The Old World monkey rhesus macaque has only one DAZ gene. In contrast, the great apes have multiple copies of DAZ, ranging from 2 copies in bonobos and gorillas to at least 6 copies in orangutans, and these DAZ genes have polymorphic structures similar to those of their human counterparts. Sequences homologous to the various AZFc amplicons are present on the Y chromosomes of some but not all primates, indicating that they arrived on the Y chromosome at different times during primate evolution. CONCLUSION: The duplication and transposition of AZFc amplicons to the human Y chromosome occurred in three waves, i.e., after the branching of the New World monkey, the gorilla, and the chimpanzee/bonobo lineages, respectively. The red-amplicon, one of the first to arrive on the Y chromosome, amplified by inverted duplication followed by direct duplication after the separation of the Old World monkey and the great ape lineages. Subsequent duplication/deletion in the various lineages gave rise to a spectrum of DAZ gene structure and copy number found in today's great apes.

Chromosomal evolution of the PKD1 gene family in primates.

BACKGROUND: The autosomal dominant polycystic kidney disease (ADPKD) is mostly caused by mutations in the PKD1 (polycystic kidney disease 1) gene located in 16p13.3. Moreover, there are six pseudogenes of PKD1 that are located proximal to the master gene in 16p13.1. In contrast, no pseudogene could be detected in the mouse genome, only a single copy gene on chromosome 17. The question arises how the human situation originated phylogenetically. To address this question we applied comparative FISH-mapping of a human PKD1-containing genomic BAC clone and a PKD1-cDNA clone to chromosomes of a variety of primate species and the dog as a non-primate outgroup species. RESULTS: Comparative FISH with the PKD1-cDNA clone clearly shows that in all primate species studied distinct single signals map in subtelomeric chromosomal positions orthologous to the short arm of human chromosome 16 harbouring the master PKD1 gene. Only in human and African great apes, but not in orangutan, FISH with both BAC and cDNA clones reveals additional signal clusters located proximal of and clearly separated from the PKD1 master genes indicating the chromosomal position of PKD1 pseudogenes in 16p of these species, respectively. Indeed, this is in accordance with sequencing data in human, chimpanzee and orangutan. Apart from the master PKD1 gene, six pseudogenes are identified in both, human and chimpanzee, while only a single-copy gene is present in the whole-genome sequence of orangutan. The phylogenetic reconstruction of the PKD1-tree reveals that all human pseudogenes are closely related to the human PKD1 gene, and all chimpanzee pseudogenes are closely related to the chimpanzee PKD1 gene. However, our statistical analyses provide strong indication that gene conversion events may have occurred within the PKD1 family members of human and chimpanzee, respectively. CONCLUSION: PKD1 must have undergone amplification very recently in hominid evolution. Duplicative transposition of the PKD1 gene and further amplification and evolution of the PKD1 pseudogenes may have arisen in a common ancestor of Homo, Pan and Gorilla approximately 8 MYA. Reticulate evolutionary processes such as gene conversion and non-allelic homologous recombination (NAHR) may have resulted in concerted evolution of PKD1 family members in human and chimpanzee and, thus, simulate an independent evolution of the PKD1 pseudogenes from their master PKD1 genes in human and chimpanzee.

Epigenetic origin of evolutionary novel centromeres.

Most evolutionary new centromeres (ENC) are composed of large arrays of satellite DNA and surrounded by segmental duplications. However, the hypothesis is that ENCs are seeded in an anonymous sequence and only over time have acquired the complexity of "normal" centromeres. Up to now evidence to test this hypothesis was lacking. We recently discovered that the well-known polymorphism of orangutan chromosome 12 was due to the presence of an ENC. We sequenced the genome of an orangutan homozygous for the ENC, and we focused our analysis on the comparison of the ENC domain with respect to its wild type counterpart. No significant variations were found. This finding is the first clear evidence that ENC seedings are epigenetic in nature. The compaction of the ENC domain was found significantly higher than the corresponding WT region and, interestingly, the expression of the only gene embedded in the region was significantly repressed.

An excess of chromosome 1 breakpoints in male infertility.

In a search for potential infertility loci, which might be revealed by clustering of chromosomal breakpoints, we compiled 464 infertile males with a balanced rearrangement from Mendelian Cytogenetics Network database (MCNdb) and compared their karyotypes with those of a Danish nation-wide cohort. We excluded Robertsonian translocations, rearrangements involving sex chromosomes and common variants. We identified 10 autosomal bands, five of which were on chromosome 1, with a large excess of breakpoints in the infertility group. Some of these could potentially harbour a male-specific infertility locus. However, a general excess of breakpoints almost everywhere on chromosome 1 was observed among the infertile males: 26.5 versus 14.5% in the cohort. This excess was observed both for translocation and inversion carriers, especially pericentric inversions, both for published and unpublished cases, and was significantly associated with azoospermia. The largest number of breakpoints was reported in 1q21; FISH mapping of four of these breakpoints revealed that they did not involve the same region at the molecular level. We suggest that chromosome 1 harbours a critical domain whose integrity is essential for male fertility.

Widespread differences in cortex DNA methylation of the "language gene" CNTNAP2 between humans and chimpanzees.

CNTNAP2, one of the largest genes in the human genome, has been linked to human-specific language abilities and neurodevelopmental disorders. Our hypothesis is that epigenetic rather than genetic changes have accelerated the evolution of the human brain. To compare the cortex DNA methylation patterns of human and chimpanzee CNTNAP2 at ultra-high resolution, we combined methylated DNA immunoprecipitation (MeDIP) with NimbleGen tiling arrays for the orthologous gene and flanking sequences. Approximately 1.59 Mb of the 2.51 Mb target region could be aligned and analyzed with a customized algorithm in both species. More than one fifth (0.34 Mb) of the analyzed sequence throughout the entire gene displayed significant methylation differences between six human and five chimpanzee cortices. One of the most striking interspecies differences with 28% methylation in human and 59% in chimpanzee cortex (by bisulfite pyrosequencing) lies in a region 300 bp upstream of human SNP rs7794745 which has been associated with autism and parent-of-origin effects. Quantitative real-time RT PCR revealed that the protein-coding splice variant CNTNAP2-201 is 1.6-fold upregulated in human cortex, compared with the chimpanzee. Transcripts CNTNAP2-001, -002, and -003 did not show skewed allelic expression, which argues against CNTNAP2 imprinting, at least in adult human brain. Collectively, our results suggest widespread cortex DNA methylation changes in CNTNAP2 since the human-chimpanzee split, supporting a role for CNTNAP2 fine-regulation in human-specific language and communication traits.

Y-Chromosome variation in hominids: intraspecific variation is limited to the polygamous chimpanzee.

BACKGROUND: We have previously demonstrated that the Y-specific ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y) varied with respect to copy number and position among chimpanzees (Pan troglodytes). In comparison, seven Y-chromosomal lineages of the bonobo (Pan paniscus), the chimpanzee's closest living relative, showed no variation. We extend our earlier comparative investigation to include an analysis of the intraspecific variation of these genes in gorillas (Gorilla gorilla) and orangutans (Pongo pygmaeus), and examine the resulting patterns in the light of the species' markedly different social and mating behaviors. METHODOLOGY/PRINCIPAL FINDINGS: Fluorescence in situ hybridization analysis (FISH) of DAZ and CDY in 12 Y-chromosomal lineages of western lowland gorilla (G. gorilla gorilla) and a single lineage of the eastern lowland gorilla (G. beringei graueri) showed no variation among lineages. Similar findings were noted for the 10 Y-chromosomal lineages examined in the Bornean orangutan (Pongo pygmaeus), and 11 Y-chromosomal lineages of the Sumatran orangutan (P. abelii). We validated the contrasting DAZ and CDY patterns using quantitative real-time polymerase chain reaction (qPCR) in chimpanzee and bonobo. CONCLUSION/SIGNIFICANCE: High intraspecific variation in copy number and position of the DAZ and CDY genes is seen only in the chimpanzee. We hypothesize that this is best explained by sperm competition that results in the variant DAZ and CDY haplotypes detected in this species. In contrast, bonobos, gorillas and orangutans-species that are not subject to sperm competition-showed no intraspecific variation in DAZ and CDY suggesting that monoandry in gorillas, and preferential female mate choice in bonobos and orangutans, probably permitted the fixation of a single Y variant in each taxon. These data support the notion that the evolutionary history of a primate Y chromosome is not simply encrypted in its DNA sequences, but is also shaped by the social and behavioral circumstances under which the specific species has evolved.

Diversification and molecular evolution of ATOH8, a gene encoding a bHLH transcription factor.

ATOH8 is a bHLH domain transcription factor implicated in the development of the nervous system, kidney, pancreas, retina and muscle. In the present study, we collected sequence of ATOH8 orthologues from 18 vertebrate species and 24 invertebrate species. The reconstruction of ATOH8 phylogeny and sequence analysis showed that this gene underwent notable divergences during evolution. For those vertebrate species investigated, we analyzed the gene structure and regulatory elements of ATOH8. We found that the bHLH domain of vertebrate ATOH8 was highly conserved. Mammals retained some specific amino acids in contrast to the non-mammalian orthologues. Mammals also developed another potential isoform, verified by a human expressed sequence tag (EST). Comparative genomic analyses of the regulatory elements revealed a replacement of the ancestral TATA box by CpG-islands in the eutherian mammals and an evolutionary tendency for TATA box reduction in vertebrates in general. We furthermore identified the region of the effective promoter of human ATOH8 which could drive the expression of EGFP reporter in the chicken embryo. In the opossum, both the coding region and regulatory elements of ATOH8 have some special features, such as the unique extended C-terminus encoded by the third exon and absence of both CpG islands and TATA elements in the regulatory region. Our gene mapping data showed that in human, ATOH8 was hosted in one chromosome which is a fusion product of two orthologous chromosomes in non-human primates. This unique chromosomal environment of human ATOH8 probably subjects its expression to the regulation at chromosomal level. We deduce that the great interspecific differences found in both ATOH8 gene sequence and its regulatory elements might be significant for the fine regulation of its spatiotemporal expression and roles of ATOH8, thus orchestrating its function in different tissues and organisms.

Y chromosomal variation tracks the evolution of mating systems in chimpanzee and bonobo.

The male-specific regions of the Y chromosome (MSY) of the human and the chimpanzee (Pan troglodytes) are fully sequenced. The most striking difference is the dramatic rearrangement of large parts of their respective MSYs. These non-recombining regions include ampliconic gene families that are known to be important for male reproduction,and are consequently under significant selective pressure. However, whether the published Y-chromosomal pattern of ampliconic fertility genes is invariable within P. troglodytes is an open but fundamental question pertinent to discussions of the evolutionary fate of the Y chromosome in different primate mating systems. To solve this question we applied fluorescence in situ hybridisation (FISH) of testis-specific expressed ampliconic fertility genes to metaphase Y chromosomes of 17 chimpanzees derived from 11 wild-born males and 16 bonobos representing seven wild-born males. We show that of eleven P. troglodytes Y-chromosomal lines, ten Y-chromosomal variants were detected based on the number and arrangement of the ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y)-a so-far never-described variation of a species' Y chromosome. In marked contrast, no variation was evident among seven Y-chromosomal lines of the bonobo, P. paniscus, the chimpanzee's closest living relative. Although, loss of variation of the Y chromosome in the bonobo by a founder effect or genetic drift cannot be excluded, these contrasting patterns might be explained in the context of the species' markedly different social and mating behaviour. In chimpanzees, multiple males copulate with a receptive female during a short period of visible anogenital swelling, and this may place significant selection on fertility genes. In bonobos, however, female mate choice may make sperm competition redundant (leading to monomorphism of fertility genes), since ovulation in this species is concealed by the prolonged anogenital swelling, and because female bonobos can occupy high-ranking positions in the group and are thus able to determine mate choice more freely.

Complex evolution of a Y-chromosomal double homeobox 4 (DUX4)-related gene family in hominoids.

The human Y chromosome carries four human Y-chromosomal euchromatin/heterochromatin transition regions, all of which are characterized by the presence of interchromosomal segmental duplications. The Yq11.1/Yq11.21 transition region harbours a peculiar segment composed of an imperfectly organized tandem-repeat structure encoding four members of the double homeobox (DUX) gene family. By comparative fluorescence in situ hybridization (FISH) analysis we have documented the primary appearance of Y-chromosomal DUX genes (DUXY) on the gibbon Y chromosome. The major amplification and dispersal of DUXY paralogs occurred after the gibbon and hominid lineages had diverged. Orthologous DUXY loci of human and chimpanzee show a highly similar structural organization. Sequence alignment survey, phylogenetic reconstruction and recombination detection analyses of human and chimpanzee DUXY genes revealed the existence of all copies in a common ancestor. Comparative analysis of the circumjacent beta-satellites indicated that DUXY genes and beta-satellites evolved in concert. However, evolutionary forces acting on DUXY genes may have induced amino acid sequence differences in the orthologous chimpanzee and human DUXY open reading frames (ORFs). The acquisition of complete ORFs in human copies might relate to evolutionary advantageous functions indicating neo-functionalization. We propose an evolutionary scenario in which an ancestral tandem array DUX gene cassette transposed to the hominoid Y chromosome followed by lineage-specific chromosomal rearrangements paved the way for a species-specific evolution of the Y-chromosomal members of a large highly diverged homeobox gene family.

Evolutionary dynamics of segmental duplications from human Y-chromosomal euchromatin/heterochromatin transition regions.

Human chromosomal regions enriched in segmental duplications are subject to extensive genomic reorganization. Such regions are particularly informative for illuminating the evolutionary history of a given chromosome. We have analyzed 866 kb of Y-chromosomal non-palindromic segmental duplications delineating four euchromatin/heterochromatin transition regions (Yp11.2/Yp11.1, Yq11.1/Yq11.21, Yq11.23/Yq12, and Yq12/PAR2). Several computational methods were applied to decipher the segmental duplication architecture and identify the ancestral origin of the 41 different duplicons. Combining computational and comparative FISH analysis, we reconstruct the evolutionary history of these regions. Our analysis indicates a continuous process of transposition of duplicated sequences onto the evolving higher primate Y chromosome, providing unique insights into the development of species-specific Y-chromosomal and autosomal duplicons. Phylogenetic sequence comparisons show that duplicons of the human Yp11.2/Yp11.1 region were already present in the macaque-human ancestor as multiple paralogs located predominantly in subtelomeric regions. In contrast, duplicons from the Yq11.1/Yq11.21, Yq11.23/Yq12, and Yq12/PAR2 regions show no evidence of duplication in rhesus macaque, but map to the pericentromeric regions in chimpanzee and human. This suggests an evolutionary shift in the direction of duplicative transposition events from subtelomeric in Old World monkeys to pericentromeric in the human/ape lineage. Extensive chromosomal relocation of autosomal-duplicated sequences from euchromatin/heterochromatin transition regions to interstitial regions as demonstrated on the pygmy chimpanzee Y chromosome support a model in which substantial reorganization and amplification of duplicated sequences may contribute to speciation.

Complex patterns of copy number variation at sites of segmental duplications: an important category of structural variation in the human genome.

The structural diversity of the human genome is much higher than previously assumed although its full extent remains unknown. To investigate the association between segmental duplications that display constitutive copy number differences (CNDs) between humans and the great apes and those which exhibit polymorphic copy number variations (CNVs) between humans, we analysed a BAC array enriched with segmental duplications displaying such CNDs. This study documents for the first time that in addition to human-specific gains common to all humans, these duplication clusters (DCs) also exhibit polymorphic CNVs > 40 kb. Segmental duplication is known to have been a frequent event during human genome evolution. Importantly, among the CNV-associated genes identified here, those involved in transcriptional regulation were found to be significantly overrepresented. Complex patterns of variation were evident at sites of DCs, manifesting as inter-individual differentially sized copy number alterations at the same genomic loci. Thus, CNVs associated with segmental duplications do not simply represent insertion/deletion polymorphisms, but rather constitute a wide variety of rearrangements involving differential amplification and partial gains and losses with high inter-individual variability. Although the number of CNVs was not found to differ between Africans and Caucasians/Asians, the average number of variant patterns per locus was significantly lower in Africans. Thus, complex variation patterns characterizing segmental duplications result from relatively recent genomic rearrangements. The high number of these rearrangements, some of which are potentially recurrent, together with differences in population size and expansion dynamics, may account for the greater diversity of CNV in Caucasians/Asians as compared with Africans.