Molecular refinement of gibbon genome rearrangements

Comparative analyses of gibbon centromeres reveal dynamic genus specific shifts in repeat composition

Centromeres are functionally conserved chromosomal loci essential for proper chromosome segregation during cell division, yet they show high sequence diversity across species. Despite their variation, a near universal feature of centromeres is the presence of repetitive sequences, such as DNA satellites and transposable elements (TEs). Because of their rapidly evolving karyotypes, gibbons represent a compelling model to investigate divergence of functional centromere sequences across short evolutionary timescales. In this study, we use ChIP-seq, RNA-seq, and fluorescence in situ hybridization to comprehensively investigate the centromeric repeat content of the four extant gibbon genera (Hoolock, Hylobates, Nomascus, and Siamang). In all gibbon genera, we find that CENP-A nucleosomes and the DNA-proteins that interface with the inner kinetochore preferentially bind retroelements of broad classes rather than satellite DNA. A previously identified gibbon-specific composite retrotransposon, LAVA, known to be expanded within the centromere regions of one gibbon genus (Hoolock), displays centromere- and species-specific sequence differences, potentially as a result of its co-option to a centromeric function. When dissecting centromere satellite composition, we discovered the presence of the retroelement-derived macrosatellite SST1 in multiple centromeres of Hoolock, whereas alpha-satellites represent the predominate satellite in the other genera, further suggesting an independent evolutionary trajectory for Hoolock centromeres. Finally, using de novo assembly of centromere sequences, we determined that transcripts originating from gibbon centromeres recapitulate the species-specific TE composition. Combined, our data reveal dynamic shifts in the repeat content that define gibbon centromeres and coincide with the extensive karyotypic diversity within this lineage.

Coalescent-Based Analyses of Genomic Sequence Data Provide a Robust Resolution of Phylogenetic Relationships among Major Groups of Gibbons

The phylogenetic relationships among extant gibbon species remain unresolved despite numerous efforts using morphological, behavorial, and genetic data and the sequencing of whole genomes. A major challenge in reconstructing the gibbon phylogeny is the radiative speciation process, which resulted in extremely short internal branches in the species phylogeny and extensive incomplete lineage sorting with extensive gene-tree heterogeneity across the genome. Here, we analyze two genomic-scale data sets, with ∼10,000 putative noncoding and exonic loci, respectively, to estimate the species tree for the major groups of gibbons. We used the Bayesian full-likelihood method bpp under the multispecies coalescent model, which naturally accommodates incomplete lineage sorting and uncertainties in the gene trees. For comparison, we included three heuristic coalescent-based methods (mp-est, SVDQuartets, and astral) as well as concatenation. From both data sets, we infer the phylogeny for the four extant gibbon genera to be (Hylobates, (Nomascus, (Hoolock, Symphalangus))). We used simulation guided by the real data to evaluate the accuracy of the methods used. Astral, while not as efficient as bpp, performed well in estimation of the species tree even in presence of excessive incomplete lineage sorting. Concatenation, mp-est and SVDQuartets were unreliable when the species tree contains very short internal branches. Likelihood ratio test of gene flow suggests a small amount of migration from Hylobates moloch to H. pileatus, while cross-genera migration is absent or rare. Our results highlight the utility of coalescent-based methods in addressing challenging species tree problems characterized by short internal branches and rampant gene tree-species tree discordance.

Accelerated pseudogenization of trace amine-associated receptor genes in primates

Trace amines (TAs) in the mammalian brain have been investigated for four decades. Trace amine‐associated receptors (TAARs) were discovered during the search for receptors activated by TAs. TAARs are considered a second class of vertebrate olfactory receptors and successfully proliferated in conjunction with adaptation to living on the ground to detect carnivore odors. Thus, therian mammals have a high number of TAAR genes due to rapid species‐specific gene duplications. In primate lineages, however, their genomes have significantly smaller numbers of TAAR genes than do other mammals. To elucidate the evolutionary force driving these patterns, exhaustive data mining of TAAR genes was performed for 13 primate genomes (covering all four infraorders) and two nonprimate euarchontan genomes. This study identified a large number of pseudogenes in many of these primate genomes and thus investigated the pseudogenization event process for the TAAR repertoires. The degeneration of TAARs is likely associated with arboreal inhabitants reducing their exposure to carnivores, and this was accelerated by the change in the nose shape of haplorhines after their divergence from strepsirrhines. Arboreal life may have decreased the reliance on the chemosensing of predators, suggestive of leading to the depauperation of TAAR subfamilies. The evolutionary deterioration of TAARs in primates has been reestablished in recently derived primates due to high selection pressure and probably functional diversity.

Epigenetic maintenance of topological domains in the highly rearranged gibbon genome

The relationship between evolutionary genome remodeling and the three-dimensional structure of the genome remain largely unexplored. Here, we use the heavily rearranged gibbon genome to examine how evolutionary chromosomal rearrangements impact genome-wide chromatin interactions, topologically associating domains (TADs), and their epigenetic landscape. We use high-resolution maps of gibbon–human breaks of synteny (BOS), apply Hi-C in gibbon, measure an array of epigenetic features, and perform cross-species comparisons. We find that gibbon rearrangements occur at TAD boundaries, independent of the parameters used to identify TADs. This overlap is supported by a remarkable genetic and epigenetic similarity between BOS and TAD boundaries, namely presence of CpG islands and SINE elements, and enrichment in CTCF and H3K4me3 binding. Cross-species comparisons reveal that regions orthologous to BOS also correspond with boundaries of large (400–600 kb) TADs in human and other mammalian species. The colocalization of rearrangement breakpoints and TAD boundaries may be due to higher chromatin fragility at these locations and/or increased selective pressure against rearrangements that disrupt TAD integrity. We also examine the small portion of BOS that did not overlap with TAD boundaries and gave rise to novel TADs in the gibbon genome. We postulate that these new TADs generally lack deleterious consequences. Last, we show that limited epigenetic homogenization occurs across breakpoints, irrespective of their time of occurrence in the gibbon lineage. Overall, our findings demonstrate remarkable conservation of chromatin interactions and epigenetic landscape in gibbons, in spite of extensive genomic shuffling.

Breaking Good: Accounting for Fragility of Genomic Regions in Rearrangement Distance Estimation

Models of evolution by genome rearrangements are prone to two types of flaws: One is to ignore the diversity of susceptibility to breakage across genomic regions, and the other is to suppose that susceptibility values are given. Without necessarily supposing their precise localization, we call "solid" the regions that are improbably broken by rearrangements and "fragile" the regions outside solid ones. We propose a model of evolution by inversions where breakage probabilities vary across fragile regions and over time. It contains as a particular case the uniform breakage model on the nucleotidic sequence, where breakage probabilities are proportional to fragile region lengths. This is very different from the frequently used pseudouniform model where all fragile regions have the same probability to break. Estimations of rearrangement distances based on the pseudouniform model completely fail on simulations with the truly uniform model. On pairs of amniote genomes, we show that identifying coding genes with solid regions yields incoherent distance estimations, especially with the pseudouniform model, and to a lesser extent with the truly uniform model. This incoherence is solved when we coestimate the number of fragile regions with the rearrangement distance. The estimated number of fragile regions is surprisingly small, suggesting that a minority of regions are recurrently used by rearrangements. Estimations for several pairs of genomes at different divergence times are in agreement with a slowly evolvable colocalization of active genomic regions in the cell.

Integration of molecular cytogenetics, dated molecular phylogeny, and model-based predictions to understand the extreme chromosome reorganization in the Neotropical genus Tonatia (Chiroptera: Phyllostomidae)

Background

Defining factors that contributed to the fixation of a high number of underdominant chromosomal rearrangements is a complex task because not only molecular mechanisms must be considered, but also the uniqueness of natural history attributes of each taxon. Ideally, detailed investigation of the chromosome architecture of an organism and related groups, placed within a phylogenetic context, is required. We used multiple approaches to investigate the dynamics of chromosomal evolution in lineages of bats with considerable karyotypic variation, focusing on the different facets contributing to fixation of the exceptional chromosomal changes in Tonatia saurophila. Integration of empirical data with proposed models of chromosome evolution was performed to understand the probable conditions for Tonatia’s karyotypic evolution.

Results

The trajectory of reorganization of chromosome blocks since the common ancestor of Glossophaginae and Phyllostominae subfamilies suggests that multiple tandem fusions, as well as disruption and fusions of conserved phyllostomid chromosomes were major drivers of karyotypic reshuffling in Tonatia. Considerable variation in the rates of chromosomal evolution between phyllostomid lineages was observed. Thirty–nine unique fusions and fission events reached fixation in Tonatia over a short period of time, followed by ~12 million years of chromosomal stasis. Physical mapping of repetitive DNA revealed an unusual accumulation of LINE-1 sequences on centromeric regions, probably associated with the chromosomal dynamics of this genus.

Conclusions

Multiple rearrangements have reached fixation in a wave-like fashion in phyllostomid bats. Different biological features of Tonatia support distinct models of rearrangement fixation, and it is unlikely that the fixations were a result of solely stochastic processes in small ancient populations. Increased recombination rates were probably facilitated by expansion of repetitive DNA, reinforced by aspects of taxon reproduction and ecology.

Electronic supplementary material

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

Chromosomal Rearrangements as Barriers to Genetic Homogenization between Archaic and Modern Humans

Chromosomal rearrangements, which shuffle DNA throughout the genome, are an important source of divergence across taxa. Using a paired-end read approach with Illumina sequence data for archaic humans, I identify changes in genome structure that occurred recently in human evolution. Hundreds of rearrangements indicate genomic trafficking between the sex chromosomes and autosomes, raising the possibility of sex-specific changes. Additionally, genes adjacent to genome structure changes in Neanderthals are associated with testis-specific expression, consistent with evolutionary theory that new genes commonly form with expression in the testes. I identify one case of new-gene creation through transposition from the Y chromosome to chromosome 10 that combines the 5'-end of the testis-specific gene Fank1 with previously untranscribed sequence. This new transcript experienced copy number expansion in archaic genomes, indicating rapid genomic change. Among rearrangements identified in Neanderthals, 13% are transposition of selfish genetic elements, whereas 32% appear to be ectopic exchange between repeats. In Denisovan, the pattern is similar but numbers are significantly higher with 18% of rearrangements reflecting transposition and 40% ectopic exchange between distantly related repeats. There is an excess of divergent rearrangements relative to polymorphism in Denisovan, which might result from nonuniform rates of mutation, possibly reflecting a burst of transposable element activity in the lineage that led to Denisovan. Finally, loci containing genome structure changes show diminished rates of introgression from Neanderthals into modern humans, consistent with the hypothesis that rearrangements serve as barriers to gene flow during hybridization. Together, these results suggest that this previously unidentified source of genomic variation has important biological consequences in human evolution.

Recent advances in primate phylogenomics

The world of primate genomics is expanding rapidly in new and exciting ways owing to lowered costs and new technologies in molecular methods and bioinformatics. The primate order is composed of 78 genera and 478 species, including human. Taxonomic inferences are complex and likely a consequence of ongoing hybridization, introgression, and reticulate evolution among closely related taxa. Recently, we applied large-scale sequencing methods and extensive taxon sampling to generate a highly resolved phylogeny that affirms, reforms, and extends previous depictions of primate speciation. The next stage of research uses this phylogeny as a foundation for investigating genome content, structure, and evolution across primates. Ongoing and future applications of a robust primate phylogeny are discussed, highlighting advancements in adaptive evolution of genes and genomes, taxonomy and conservation management of endangered species, next-generation genomic technologies, and biomedicine.

Comprehensive Analyses of White-Handed Gibbon Chromosomes Enables Access to 92 Evolutionary Conserved Breakpoints Compared to the Human Genome

Gibbon species (Hylobatidae) impress with an unusually high number of numerical and structural chromosomal changes within the family itself as well as compared to other Hominoidea including humans. In former studies applying molecular cytogenetic methods, 86 evolutionary conserved breakpoints (ECBs) were reported in the white-handed gibbon (Hylobates lar, HLA) with respect to the human genome. To analyze those ECBs in more detail and also to achieve a better understanding of the fast karyotype evolution in Hylobatidae, molecular data for these regions are indispensably necessary. In the present study, we obtained whole chromosome-specific probes by microdissection of all 21 HLA autosomes and prepared them for aCGH. Locus-specific DNA probes were also used for further molecular cytogenetic characterization of selected regions. Thus, we could map 6 yet unreported ECBs in HLA with respect to the human genome. Additionally, in 26 of the 86 previously reported ECBs, the present approach enabled a more precise breakpoint mapping. Interestingly, a preferred localization of ECBs within segmental duplications, copy number variant regions, and fragile sites was observed.

A comprehensive molecular cytogenetic analysis of chromosome rearrangements in gibbons

Chromosome rearrangements in small apes are up to 20 times more frequent than in most mammals. Because of their complexity, the full extent of chromosome evolution in these hominoids is not yet fully documented. However, previous work with array painting, BAC-FISH, and selective sequencing in two of the four karyomorphs has shown that high-resolution methods can precisely define chromosome breakpoints and map the complex flow of evolutionary chromosome rearrangements. Here we use these tools to precisely define the rearrangements that have occurred in the remaining two karyomorphs, genera Symphalangus (2n = 50) and Hoolock (2n = 38). This research provides the most comprehensive insight into the evolutionary origins of chromosome rearrangements involved in transforming small apes genome. Bioinformatics analyses of the human–gibbon synteny breakpoints revealed association with transposable elements and segmental duplications, providing some insight into the mechanisms that might have promoted rearrangements in small apes. In the near future, the comparison of gibbon genome sequences will provide novel insights to test hypotheses concerning the mechanisms of chromosome evolution. The precise definition of synteny block boundaries and orientation, chromosomal fusions, and centromere repositioning events presented here will facilitate genome sequence assembly for these close relatives of humans.

Evolutionary molecular cytogenetics of catarrhine primates: past, present and future

The catarrhine primates were the first group of species studied with comparative molecular cytogenetics. Many of the fundamental techniques and principles of analysis were initially applied to comparisons in these primates, including interspecific chromosome painting, reciprocal chromosome painting and the extensive use of cloned DNA probes for evolutionary analysis. The definition and importance of chromosome syntenies and associations for a correct cladistics analysis of phylogenomic relationships were first applied to catarrhines. These early chromosome painting studies vividly illustrated a striking conservation of the genome between humans and macaques. Contemporarily, it also revealed profound differences between humans and gibbons, a group of species more closely related to humans, making it clear that chromosome evolution did not follow a molecular clock. Chromosome painting has now been applied to more that 60 primate species and the translocation history has been mapped onto the major taxonomic divisions in the tree of primate evolution. In situ hybridization of cloned DNA probes, primarily BAC-FISH, also made it possible to more precisely map breakpoints with spanning and flanking BACs. These studies established marker order and disclosed intrachromosomal rearrangements. When applied comparatively to a range of primate species, they led to the discovery of evolutionary new centromeres as an important new category of chromosome evolution. BAC-FISH studies are intimately connected to genome sequencing, and probes can usually be assigned to a precise location in the genome assembly. This connection ties molecular cytogenetics securely to genome sequencing, assuring that molecular cytogenetics will continue to have a productive future in the multidisciplinary science of phylogenomics.

Centromere repositioning in mammals

The evolutionary history of chromosomes can be tracked by the comparative hybridization of large panels of bacterial artificial chromosome clones. This approach has disclosed an unprecedented phenomenon: 'centromere repositioning', that is, the movement of the centromere along the chromosome without marker order variation. The occurrence of evolutionary new centromeres (ENCs) is relatively frequent. In macaque, for instance, 9 out of 20 autosomal centromeres are evolutionarily new; in donkey at least 5 such neocentromeres originated after divergence from the zebra, in less than 1 million years. Recently, orangutan chromosome 9, considered to be heterozygous for a complex rearrangement, was discovered to be an ENC. In humans, in addition to neocentromeres that arise in acentric fragments and result in clinical phenotypes, 8 centromere-repositioning events have been reported. These 'real-time' repositioned centromere-seeding events provide clues to ENC birth and progression. In the present paper, we provide a review of the centromere repositioning. We add new data on the population genetics of the ENC of the orangutan, and describe for the first time an ENC on the X chromosome of squirrel monkeys. Next-generation sequencing technologies have started an unprecedented, flourishing period of rapid whole-genome sequencing. In this context, it is worth noting that these technologies, uncoupled from cytogenetics, would miss all the biological data on evolutionary centromere repositioning. Therefore, we can anticipate that classical and molecular cytogenetics will continue to have a crucial role in the identification of centromere movements. Indeed, all ENCs and human neocentromeres were found following classical and molecular cytogenetic investigations.

The genome diversity and karyotype evolution of mammals

The past decade has witnessed an explosion of genome sequencing and mapping in evolutionary diverse species. While full genome sequencing of mammals is rapidly progressing, the ability to assemble and align orthologous whole chromosome regions from more than a few species is still not possible. The intense focus on building of comparative maps for companion (dog and cat), laboratory (mice and rat) and agricultural (cattle, pig, and horse) animals has traditionally been used as a means to understand the underlying basis of disease-related or economically important phenotypes. However, these maps also provide an unprecedented opportunity to use multispecies analysis as a tool for inferring karyotype evolution. Comparative chromosome painting and related techniques are now considered to be the most powerful approaches in comparative genome studies. Homologies can be identified with high accuracy using molecularly defined DNA probes for fluorescence in situ hybridization (FISH) on chromosomes of different species. Chromosome painting data are now available for members of nearly all mammalian orders. In most orders, there are species with rates of chromosome evolution that can be considered as 'default' rates. The number of rearrangements that have become fixed in evolutionary history seems comparatively low, bearing in mind the 180 million years of the mammalian radiation. Comparative chromosome maps record the history of karyotype changes that have occurred during evolution. The aim of this review is to provide an overview of these recent advances in our endeavor to decipher the karyotype evolution of mammals by integrating the published results together with some of our latest unpublished results.

A comparative physical map reveals the pattern of chromosomal evolution between the turkey (Meleagris gallopavo) and chicken (Gallus gallus) genomes

Background

A robust bacterial artificial chromosome (BAC)-based physical map is essential for many aspects of genomics research, including an understanding of chromosome evolution, high-resolution genome mapping, marker-assisted breeding, positional cloning of genes, and quantitative trait analysis. To facilitate turkey genetics research and better understand avian genome evolution, a BAC-based integrated physical, genetic, and comparative map was developed for this important agricultural species.

Results

The turkey genome physical map was constructed based on 74,013 BAC fingerprints (11.9 × coverage) from two independent libraries, and it was integrated with the turkey genetic map and chicken genome sequence using over 41,400 BAC assignments identified by 3,499 overgo hybridization probes along with > 43,000 BAC end sequences. The physical-comparative map consists of 74 BAC contigs, with an average contig size of 13.6 Mb. All but four of the turkey chromosomes were spanned on this map by three or fewer contigs, with 14 chromosomes spanned by a single contig and nine chromosomes spanned by two contigs. This map predicts 20 to 27 major rearrangements distinguishing turkey and chicken chromosomes, despite up to 40 million years of separate evolution between the two species. These data elucidate the chromosomal evolutionary pattern within the Phasianidae that led to the modern turkey and chicken karyotypes. The predominant rearrangement mode involves intra-chromosomal inversions, and there is a clear bias for these to result in centromere locations at or near telomeres in turkey chromosomes, in comparison to interstitial centromeres in the orthologous chicken chromosomes.

Conclusion

The BAC-based turkey-chicken comparative map provides novel insights into the evolution of avian genomes, a framework for assembly of turkey whole genome shotgun sequencing data, and tools for enhanced genetic improvement of these important agricultural and model species.

Mitochondrial evidence for multiple radiations in the evolutionary history of small apes

BACKGROUND

Gibbons or small apes inhabit tropical and subtropical rain forests in Southeast Asia and adjacent regions, and are, next to great apes, our closest living relatives. With up to 16 species, gibbons form the most diverse group of living hominoids, but the number of taxa, their phylogenetic relationships and their phylogeography is controversial. To further the discussion of these issues we analyzed the complete mitochondrial cytochrome b gene from 85 individuals representing all gibbon species, including most subspecies.

RESULTS

Based on phylogenetic tree reconstructions, several monophyletic clades were detected, corresponding to genera, species and subspecies. A significantly supported branching pattern was obtained for members of the genus Nomascus but not for the genus Hylobates. The phylogenetic relationships among the four genera were also not well resolved. Nevertheless, the new data permitted the estimation of divergence ages for all taxa for the first time and showed that most lineages emerged during four short time periods. In the first, between approximately 6.7 and approximately 8.3 mya, the four gibbon genera diverged from each other. In the second (approximately 3.0 - approximately 3.9 mya) and in the third period (approximately 1.3 - approximately 1.8 mya), Hylobates and Hoolock differentiated. Finally, between approximately 0.5 and approximately 1.1 mya, Hylobates lar diverged into subspecies. In contrast, differentiation of Nomascus into species and subspecies was a continuous and prolonged process lasting from approximately 4.2 until approximately 0.4 mya.

CONCLUSIONS

Although relationships among gibbon taxa on various levels remain unresolved, the present study provides a more complete view of the evolutionary and biogeographic history of the hylobatid family, and a more solid genetic basis for the taxonomic classification of the surviving taxa. We also show that mtDNA constitutes a useful marker for the accurate identification of individual gibbons, a tool which is urgently required to locate hunting hotspots and select individuals for captive breeding programs. Further studies including nuclear sequence data are necessary to completely understand the phylogeny and phylogeography of gibbons.

A small, variable, and irregular killer cell Ig-like receptor locus accompanies the absence of MHC-C and MHC-G in gibbons

The killer cell Ig-like receptors (KIRs) of NK cells recognize MHC class I ligands and function in placental reproduction and immune defense against pathogens. During the evolution of monkeys, great apes, and humans, an ancestral KIR3DL gene expanded to become a diverse and rapidly evolving gene family of four KIR lineages. Characterizing the KIR locus are three framework regions, defining two intervals of variable gene content. By analysis of four KIR haplotypes from two species of gibbon, we find that the smaller apes do not conform to these rules. Although diverse and irregular in structure, the gibbon haplotypes are unusually small, containing only two to five functional genes. Comparison with the predicted ancestral hominoid KIR haplotype indicates that modern gibbon KIR haplotypes were formed by a series of deletion events, which created new hybrid genes as well as eliminating ancestral genes. Of the three framework regions, only KIR3DL3 (lineage V), defining the 5' end of the KIR locus, is present and intact on all gibbon KIR haplotypes. KIR2DL4 (lineage I) defining the central framework region has been a major target for elimination or inactivation, correlating with the absence of its putative ligand, MHC-G, in gibbons. Similarly, the MHC-C-driven expansion of lineage III KIR genes in great apes has not occurred in gibbons because they lack MHC-C. Our results indicate that the selective forces shaping the size and organization of the gibbon KIR locus differed from those acting upon the KIR of other hominoid species.

Extensive conserved synteny of genes between the karyotypes of Manduca sexta and Bombyx mori revealed by BAC-FISH mapping

Background

Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes.

Methodology/Principal Findings

We isolated 214 clones containing 159 orthologs of well conserved single-copy genes of a sequenced lepidopteran model, the silkworm, Bombyx mori, from a BAC library of a sphingid with an unexplored genome, the tobacco hornworm, Manduca sexta. We then constructed a BAC-FISH karyotype identifying all 28 chromosomes of M. sexta by mapping 124 loci using the corresponding BAC clones. BAC probes from three M. sexta chromosomes also generated clear signals on the corresponding chromosomes of the convolvulus hawk moth, Agrius convolvuli, which belongs to the same subfamily, Sphinginae, as M. sexta.

Conclusions/Significance

Comparison of the M. sexta BAC physical map with the linkage map and genome sequence of B. mori pointed to extensive conserved synteny including conserved gene order in most chromosomes. Only a few rearrangements, including three inversions, three translocations, and two fission/fusion events were estimated to have occurred after the divergence of Bombycidae and Sphingidae. These results add to accumulating evidence for the stability of lepidopteran genomes. Generating signals on A. convolvuli chromosomes using heterologous M. sexta probes demonstrated that BAC-FISH with orthologous sequences can be used for karyotyping a wide range of related and genetically uncharacterized species, significantly extending the ability to develop synteny maps for comparative and functional genomics.

Sequencing primate genomes: what have we learned?

We summarize the progress in whole-genome sequencing and analyses of primate genomes. These emerging genome datasets have broadened our understanding of primate genome evolution revealing unexpected and complex patterns of evolutionary change. This includes the characterization of genome structural variation, episodic changes in the repeat landscape, differences in gene expression, new models regarding speciation, and the ephemeral nature of the recombination landscape. The functional characterization of genomic differences important in primate speciation and adaptation remains a significant challenge. Limited access to biological materials, the lack of detailed phenotypic data and the endangered status of many critical primate species have significantly attenuated research into the genetic basis of primate evolution. Next-generation sequencing technologies promise to greatly expand the number of available primate genome sequences; however, such draft genome sequences will likely miss critical genetic differences within complex genomic regions unless dedicated efforts are put forward to understand the full spectrum of genetic variation.

Analysis of fine-scale mammalian evolutionary breakpoints provides new insight into their relation to genome organisation

Background

The Intergenic Breakage Model, which is the current model of structural genome evolution, considers that evolutionary rearrangement breakages happen with a uniform propensity along the genome but are selected against in genes, their regulatory regions and in-between. However, a growing body of evidence shows that there exists regions along mammalian genomes that present a high susceptibility to breakage. We reconsidered this question taking advantage of a recently published methodology for the precise detection of rearrangement breakpoints based on pairwise genome comparisons.

Results

We applied this methodology between the genome of human and those of five sequenced eutherian mammals which allowed us to delineate evolutionary breakpoint regions along the human genome with a finer resolution (median size 26.6 kb) than obtained before. We investigated the distribution of these breakpoints with respect to genome organisation into domains of different activity. In agreement with the Intergenic Breakage Model, we observed that breakpoints are under-represented in genes. Surprisingly however, the density of breakpoints in small intergenes (1 per Mb) appears significantly higher than in gene deserts (0.1 per Mb).

More generally, we found a heterogeneous distribution of breakpoints that follows the organisation of the genome into isochores (breakpoints are more frequent in GC-rich regions). We then discuss the hypothesis that regions with an enhanced susceptibility to breakage correspond to regions of high transcriptional activity and replication initiation.

Conclusion

We propose a model to describe the heterogeneous distribution of evolutionary breakpoints along human chromosomes that combines natural selection and a mutational bias linked to local open chromatin state.

New insights into centromere organization and evolution from the white-cheeked gibbon and marmoset

The evolutionary history of alpha-satellite DNA, the major component of primate centromeres, is hardly defined because of the difficulty in its sequence assembly and its rapid evolution when compared with most genomic sequences. By using several approaches, we have cloned, sequenced, and characterized alpha-satellite sequences from two species representing critical nodes in the primate phylogeny: the white-cheeked gibbon, a lesser ape, and marmoset, a New World monkey. Sequence analyses demonstrate that white-cheeked gibbon and marmoset alpha-satellite sequences are formed by units of approximately 171 and approximately 342 bp, respectively, and they both lack the high-order structure found in humans and great apes. Fluorescent in situ hybridization characterization shows a broad dispersal of alpha-satellite in the white-cheeked gibbon genome including centromeric, telomeric, and chromosomal interstitial localizations. On the other hand, centromeres in marmoset appear organized in highly divergent dimers roughly of 342 bp that show a similarity between monomers much lower than previously reported dimers, thus representing an ancient dimeric structure. All these data shed light on the evolution of the centromeric sequences in Primates. Our results suggest radical differences in the structure, organization, and evolution of alpha-satellite DNA among different primate species, supporting the notion that 1) all the centromeric sequence in Primates evolved by genomic amplification, unequal crossover, and sequence homogenization using a 171 bp monomer as the basic seeding unit and 2) centromeric function is linked to relatively short repeated elements, more than higher-order structure. Moreover, our data indicate that complex higher-order repeat structures are a peculiarity of the hominid lineage, showing the more complex organization in humans.

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.

Gene synteny comparisons between different vertebrates provide new insights into breakage and fusion events during mammalian karyotype evolution

BACKGROUND

Genome comparisons have made possible the reconstruction of the eutherian ancestral karyotype but also have the potential to provide new insights into the evolutionary inter-relationship of the different eutherian orders within the mammalian phylogenetic tree. Such comparisons can additionally reveal (i) the nature of the DNA sequences present within the evolutionary breakpoint regions and (ii) whether or not the evolutionary breakpoints occur randomly across the genome. Gene synteny analysis (E-painting) not only greatly reduces the complexity of comparative genome sequence analysis but also extends its evolutionary reach.

RESULTS

E-painting was used to compare the genome sequences of six different mammalian species and chicken. A total of 526 evolutionary breakpoint intervals were identified and these were mapped to a median resolution of 120 kb, the highest level of resolution so far obtained. A marked correlation was noted between evolutionary breakpoint frequency and gene density. This correlation was significant not only at the chromosomal level but also sub-chromosomally when comparing genome intervals of lengths as short as 40 kb. Contrary to previous findings, a comparison of evolutionary breakpoint locations with the chromosomal positions of well mapped common fragile sites and cancer-associated breakpoints failed to reveal any evidence for significant co-location. Primate-specific chromosomal rearrangements were however found to occur preferentially in regions containing segmental duplications and copy number variants.

CONCLUSION

Specific chromosomal regions appear to be prone to recurring rearrangement in different mammalian lineages ('breakpoint reuse') even if the breakpoints themselves are likely to be non-identical. The putative ancestral eutherian genome, reconstructed on the basis of the synteny analysis of 7 vertebrate genome sequences, not only confirmed the results of previous molecular cytogenetic studies but also increased the definition of the inferred structure of ancestral eutherian chromosomes. For the first time in such an analysis, the opossum was included as an outgroup species. This served to confirm our previous model of the ancestral eutherian genome since all ancestral syntenic segment associations were also noted in this marsupial.

A chromosomal inversion unique to the northern white-cheeked gibbon

The gibbon family belongs to the superfamily Hominoidea and includes 15 species divided into four genera. Each genus possesses a distinct karyotype with chromosome numbers varying from 38 to 52. This diversity is the result of numerous chromosomal changes that have accumulated during the evolution of the gibbon lineage, a quite unique feature in comparison with other hominoids and most of the other primates. Some gibbon species and subspecies rank among the most endangered primates in the world. Breeding programs can be extremely challenging and hybridization plays an important role within the factors responsible for the decline of captive gibbons. With less than 500 individuals left in the wild, the northern white-cheeked gibbon (Nomascus leucogenys leucogenys, NLE) is the most endangered primate in a successful captive breeding program. We present here the analysis of an inversion that we show being specific for the northern white-cheeked gibbon and can be used as one of the criteria to distinguish this subspecies from other gibbon taxa. The availability of the sequence spanning for one of the breakpoints of the inversion allows detecting it by a simple PCR test also on low quality DNA. Our results demonstrate the important role of genomics in providing tools for conservation efforts.

Evolutionary new centromeres in primates

The centromere has a pivotal role in structuring chromosomal architecture, but remains a poorly understood and seemingly paradoxical "black hole." Centromeres are a very rapidly evolving segment of the genome and it is now known that centromere shifts in evolution are not rare and must be considered on a par with other chromosome rearrangements. Recently, unprecedented findings on neocentromeres and evolutionary new centromeres (ENC) have helped clarify the relationship of the centromere within the genome and shown that these two phenomena are two faces of the same coin. No prominent sequence features are known that promote centromere formation and both types of new centromeres are formed epigenetically, both clinical neocentromeres and ENC cluster at chromosomal "hotspots." The clustering of neocentromeres in 8p is probably the result of the relatively high frequency of noncanonical pairing. Studies on the evolution of the chromosomes 3, 13, and 15 help explain why there are clusters of neocentromeres. These domains often correspond to ancestral inactivated centromeres and some regions can preserve features that trigger neocentromere emergence over tens of millions of years. Neocentromeres may be correlated with the distribution of segmental duplications (SDs) in regions of extreme plasticity that often can be characterized as gene deserts. Further, because centromeres and associated pericentric regions are dynamically complex, centromere shifts may turbocharge genome reorganization by influencing the distribution of heterochromatin. The "reuse" of regions as centromere seeding-points in evolution and in human clinical cases further extends the concept of "reuse" of specific domains for "chromosomal events."

Evolutionary-new centromeres preferentially emerge within gene deserts

BACKGROUND

Evolutionary-new centromeres (ENCs) result from the seeding of a centromere at an ectopic location along the chromosome during evolution. The novel centromere rapidly acquires the complex structure typical of eukaryote centromeres. This phenomenon has played an important role in shaping primate karyotypes. A recent study on the evolutionary-new centromere of macaque chromosome 4 (human 6) showed that the evolutionary-new centromere domain was deeply restructured, following the seeding, with respect to the corresponding human region assumed as ancestral. It was also demonstrated that the region was devoid of genes. We hypothesized that these two observations were not merely coincidental and that the absence of genes in the seeding area constituted a crucial condition for the evolutionary-new centromere fixation in the population.

RESULTS

To test our hypothesis, we characterized 14 evolutionary-new centromeres selected according to conservative criteria. Using different experimental approaches, we assessed the extent of genomic restructuring. We then determined the gene density in the ancestral domain where each evolutionary-new centromere was seeded.

CONCLUSIONS

Our study suggests that restructuring of the seeding regions is an intrinsic property of novel evolutionary centromeres that could be regarded as potentially detrimental to the normal functioning of genes embedded in the region. The absence of genes, which was found to be of high statistical significance, appeared as a unique favorable scenario permissive of evolutionary-new centromere fixation in the population.

Sequencing human-gibbon breakpoints of synteny reveals mosaic new insertions at rearrangement sites

The gibbon genome exhibits extensive karyotypic diversity with an increased rate of chromosomal rearrangements during evolution. In an effort to understand the mechanistic origin and implications of these rearrangement events, we sequenced 24 synteny breakpoint regions in the white-cheeked gibbon (Nomascus leucogenys, NLE) in the form of high-quality BAC insert sequences (4.2 Mbp). While there is a significant deficit of breakpoints in genes, we identified seven human gene structures involved in signaling pathways (DEPDC4, GNG10), phospholipid metabolism (ENPP5, PLSCR2), beta-oxidation (ECH1), cellular structure and transport (HEATR4), and transcription (ZNF461), that have been disrupted in the NLE gibbon lineage. Notably, only three of these genes show the expected evolutionary signatures of pseudogenization. Sequence analysis of the breakpoints suggested both nonclassical nonhomologous end-joining (NHEJ) and replication-based mechanisms of rearrangement. A substantial number (11/24) of human-NLE gibbon breakpoints showed new insertions of gibbon-specific repeats and mosaic structures formed from disparate sequences including segmental duplications, LINE, SINE, and LTR elements. Analysis of these sites provides a model for a replication-dependent repair mechanism for double-strand breaks (DSBs) at rearrangement sites and insights into the structure and formation of primate segmental duplications at sites of genomic rearrangements during evolution.

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.

How segmental duplications shape our genome: recent evolution of ABCC6 and PKD1 Mendelian disease genes

The completion of the Human Genome Project has brought the understanding that our genome contains an unexpectedly large proportion of segmental duplications. This poses the challenge of elucidating the consequences of recent duplications on physiology. We have conducted an in-depth study of a subset of segmental duplications on chromosome 16. We focused on PKD1 and ABCC6 duplications because mutations affecting these genes are responsible for the Mendelian disorders autosomal dominant polycystic kidney disease and pseudoxanthoma elasticum, respectively. We establish that duplications of PKD1 and ABCC6 are associated to low-copy repeat 16a and show that such duplications have occurred several times independently in different primate species. We demonstrate that partial duplication of PKD1 and ABCC6 has numerous consequences: the pseudogenes give rise to new transcripts and mediate gene conversion, which not only results in disease-causing mutations but also serves as a reservoir for sequence variation. The duplicated segments are also involved in submicroscopic and microscopic genomic rearrangements, contributing to structural variation in human and chromosomal break points in the gibbon. In conclusion, our data shed light on the recent and ongoing evolution of chromosome 16 mediated by segmental duplication and deepen our understanding of the history of two Mendelian disorder genes.

Evolution of a complex minisatellite DNA sequence

Minisatellites are tandem repeats of short DNA units widely distributed in genomes. However, the information on their dynamics in a phylogenetic context is very limited. Here we have studied the organization of the MsH43 locus in several species of primates and from these data we have reconstructed the evolutionary history of this complex minisatellite. Overall, with the exception of gibbon, MsH43 has an organization that is asymmetric, since the distribution of repeats is distinct between the 5' and 3' halves, and heterogeneous since there are many different repeats, some of them characteristic of each species. Inspection of the MsH43 arrays showed the existence of many duplications and deletions, suggesting the implication of slippage processes in the generation of polymorphism. Concerning the evolutionary history of this minisatellite, we propose that the birth of MsH43 may be situated before the divergence of Old World Monkeys since we found the existence of some MsH43 repeat motifs in prosimians and New World Monkeys. The analysis of MsH43 in apes revealed the existence of an evolutionary breakpoint in the pathway that originated African great apes and humans. Remarkably, human MsH43 is more homologous to orang-utan than to the corresponding sequence in gorilla and chimpanzee. This finding does not comply with the evolutionary paradigm that continuous alterations occur during the course of genome evolution. To adjust our results to the standard phylogeny of primates, we propose the existence of a wandering allele that was maintained almost unaltered during the period that extends between orang-utan and humans.

Construction, characterization, and chromosomal mapping of a fosmid library of the white-cheeked gibbon (Nomascus leucogenys)

Gibbons have experienced extensive karyotype rearrangements during evolution and represent an ideal model for studying the underlying molecular mechanism of evolutionary chromosomal rearrangements. It is anticipated that the cloning and sequence characterization of evolutionary chromosomal breakpoints will provide vital insights into the molecular force that has driven such a radical karyotype reshuffle in gibbons. We constructed and characterized a high-quality fosmid library of the white-cheeked gibbon (Nomascus leucogenys) containing 192,000 non-redundant clones with an average insert size of 38 kb and 2.5-fold genome coverage. By end sequencing of 100 randomly selected fosmid clones, we generated 196 sequence tags for the library. These end-sequenced fosmid clones were then mapped onto the chromosomes of the white-cheeked gibbon by fluorescence in situ hybridization, and no spurious chimeric clone was detected. BLAST search against the human genome showed a good correlation between the number of hit clones and the number of chromosomes, an indication of unbiased chromosomal distribution of the fosmid library. The chromosomal distribution of the mapped clones is also consistent with the BLAST search result against human and white-cheeked gibbon genomes. The fosmid library and the mapped clones will serve as a valuable resource for further studying gibbons’ chromosomal rearrangements and the underlying molecular mechanism as well as for comparative genomic study in the lesser apes.

Tracking the complex flow of chromosome rearrangements from the Hominoidea Ancestor to extant Hylobates and Nomascus Gibbons by high-resolution synteny mapping

In this study we characterized the extension, reciprocal arrangement, and orientation of syntenic chromosomal segments in the lar gibbon (Hylobates lar, HLA) by hybridization of a panel of approximately 1000 human BAC clones. Each lar gibbon rearrangement was defined by a splitting BAC clone or by two overlapping clones flanking the breakpoint. A reconstruction of the synteny arrangement of the last common ancestor of all living lesser apes was made by combining these data with previous results in Nomascus leucogenys, Hoolock hoolock, and Symphalangus syndactylus. The definition of the ancestral synteny organization facilitated tracking the cascade of chromosomal changes from the Hominoidea ancestor to the present day karyotype of Hylobates and Nomascus. Each chromosomal rearrangement could be placed within an approximate phylogenetic and temporal framework. We identified 12 lar-specific rearrangements and five previously undescribed rearrangements that occurred in the Hylobatidae ancestor. The majority of the chromosomal differences between lar gibbons and humans are due to rearrangements that occurred in the Hylobatidae ancestor (38 events), consistent with the hypothesis that the genus Hylobates is the most recently evolved lesser ape genus. The rates of rearrangements in gibbons are 10 to 20 times higher than the mammalian default rate. Segmental duplication may be a driving force in gibbon chromosome evolution, because a consistent number of rearrangements involves pericentromeric regions (10 events) and centromere inactivation (seven events). Both phenomena can be reasonably supposed to have strongly contributed to the euchromatic dispersal of segmental duplications typical of pericentromeric regions. This hypothesis can be more fully tested when the sequence of this gibbon species becomes available. The detailed synteny map provided here will, in turn, substantially facilitate sequence assembly efforts.

Molecular mechanisms of chromosomal rearrangement during primate evolution

Breakpoint analysis of the large chromosomal rearrangements which have occurred during primate evolution promises to yield new insights into the underlying mechanisms of mutagenesis. Comparison of these evolutionary breakpoints with those that are disease-associated in humans, and which occur during either meiotic or mitotic cell division, should help to identify basic mechanistic similarities as well as differences. It has recently become clear that segmental duplications (SDs) have had a very significant impact on genome plasticity during primate evolution. In comparisons of the human and chimpanzee genomes, SDs have been found in flanking regions of 70-80% of inversions and approximately 40% of deletions/duplications. A strong spatial association between primate-specific breakpoints and SDs has also become evident from comparisons of human with other mammalian genomes. The lineage-specific hyperexpansion of certain SDs observed in the genomes of human, chimpanzee, gorilla and gibbon is indicative of the intrinsic instability of some SDs in primates. However, since many primate-specific breakpoints map to regions lacking SDs, but containing interspersed high-copy repetitive sequence elements such as SINEs, LINEs, LTRs, alpha-satellites and (AT)( n ) repeats, we may infer that a range of different molecular mechanisms have probably been involved in promoting chromosomal breakage during the evolution of primate genomes.

Primate chromosome evolution: ancestral karyotypes, marker order and neocentromeres

In 1992 the Japanese macaque was the first species for which the homology of the entire karyotype was established by cross-species chromosome painting. Today, there are chromosome painting data on more than 50 species of primates. Although chromosome painting is a rapid and economical method for tracking translocations, it has limited utility for revealing intrachromosomal rearrangements. Fortunately, the use of BAC-FISH in the last few years has allowed remarkable progress in determining marker order along primate chromosomes and there are now marker order data on an array of primate species for a good number of chromosomes. These data reveal inversions, but also show that centromeres of many orthologous chromosomes are embedded in different genomic contexts. Even if the mechanisms of neocentromere formation and progression are just beginning to be understood, it is clear that these phenomena had a significant impact on shaping the primate genome and are fundamental to our understanding of genome evolution. In this report we complete and integrate the dataset of BAC-FISH marker order for human syntenies 1, 2, 4, 5, 8, 12, 17, 18, 19, 21, 22 and the X. These results allowed us to develop hypotheses about the content, marker order and centromere position in ancestral karyotypes at five major branching points on the primate evolutionary tree: ancestral primate, ancestral anthropoid, ancestral platyrrhine, ancestral catarrhine and ancestral hominoid. Current models suggest that between-species structural rearrangements are often intimately related to speciation. Comparative primate cytogenetics has become an important tool for elucidating the phylogeny and the taxonomy of primates. It has become increasingly apparent that molecular cytogenetic data in the future can be fruitfully combined with whole-genome assemblies to advance our understanding of primate genome evolution as well as the mechanisms and processes that have led to the origin of the human genome.

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.

Construction, alignment and analysis of twelve framework physical maps that represent the ten genome types of the genus Oryza

Bacterial artificial chromosome (BAC) fingerprint and end-sequenced physical maps representing the ten genome types of Oryza are presented

We describe the establishment and analysis of a genus-wide comparative framework composed of 12 bacterial artificial chromosome fingerprint and end-sequenced physical maps representing the 10 genome types of Oryza aligned to the O. sativa ssp. japonica reference genome sequence. Over 932 Mb of end sequence was analyzed for repeats, simple sequence repeats, miRNA and single nucleotide variations, providing the most extensive analysis of Oryza sequence to date.

A most distant intergeneric hybrid offspring (Larcon) of lesser apes, Nomascus leucogenys and Hylobates lar

Unlike humans, which are the sole remaining representatives of a once larger group of bipedal apes (hominins), the "lesser apes" (hylobatids) are a diverse radiation with numerous extant species. Consequently, the lesser apes can provide a valuable evolutionary window onto the possible interactions (e.g., interbreeding) of hominin lineages coexisting in the same time and place. In the present work, we employ chromosomal analyses to verify the hybrid ancestry of an individual (Larcon) produced by two of the most distant genera of lesser apes, Hylobates (lar-group gibbons) and Nomascus (concolor-group gibbons). In addition to a mixed pelage pattern, the hybrid animal carries a 48-chromosome karyotype that consists of the haploid complements of each parental species: Hylobates lar (n = 22) and Nomascus leucogenys leucogenys (n = 26). Studies of this animal's karyotype shed light onto the processes of speciation and genus-level divergence in the lesser apes and, by extension, across the Hominoidea.

Primate comparative genomics: lemur biology and evolution

Comparative genome sequencing projects are providing insight into aspects of genome biology that raise new questions and challenge existing paradigms. Placement in the phylogenetic tree can often be a major determinant of which organism to choose for study. Lemurs hold a key position at the base of the primate evolutionary tree and will be highly informative for the genomics community by offering comparisons of primate-specific characteristics and processes. Combining research in chromosome evolution, genome evolution and behavior with lemur comparative genomic sequencing will offer insights into many levels of primate evolution. We discuss the current state of lemur cytogenetic and phylogenetic analyses, and suggest how focusing more genomic efforts on lemurs will be beneficial to understanding human and primate evolution, as well as disease, and will contribute to conservation efforts.