Genome-wide characterization of centromeric satellites from multiple mammalian genomes

Despite its importance in cell biology and evolution, the centromere has remained the final frontier in genome assembly and annotation due to its complex repeat structure. However, isolation and characterization of the centromeric repeats from newly sequenced species are necessary for a complete understanding of genome evolution and function. In recent years, various genomes have been sequenced, but the characterization of the corresponding centromeric DNA has lagged behind. Here, we present a computational method (RepeatNet) to systematically identify higher-order repeat structures from unassembled whole-genome shotgun sequence and test whether these sequence elements correspond to functional centromeric sequences. We analyzed genome datasets from six species of mammals representing the diversity of the mammalian lineage, namely, horse, dog, elephant, armadillo, opossum, and platypus. We define candidate monomer satellite repeats and demonstrate centromeric localization for five of the six genomes. Our analysis revealed the greatest diversity of centromeric sequences in horse and dog in contrast to elephant and armadillo, which showed high-centromeric sequence homogeneity. We could not isolate centromeric sequences within the platypus genome, suggesting that centromeres in platypus are not enriched in satellite DNA. Our method can be applied to the characterization of thousands of other vertebrate genomes anticipated for sequencing in the near future, providing an important tool for annotation of centromeres.

Species identification and chromosome variation of captive two-toed sloths

Two-toed sloth species, Linnaeus's and Hoffmman's, are frequent residents of zoo collections in North America. However, species identification has always been problematic because of their large overlap in external morphology, which represents an obstacle to the captive breeding program. We describe here a PCR-based technique that allows species identification of two-toed sloths without requiring sequencing, by using a mitochondrial marker (COI gene) and restriction enzyme assay. We also report intra- and inter-specific patterns of chromosome variation in captive two-toed sloths. Molecularly, we identified 22 samples of Linnaeus's and Hoffmman's two-toed sloths corresponding to 14 and 8 individuals, respectively. One animal was identified as a hybrid using the nuclear gene Enam having alleles derived from both species. The chromosome number in Hoffman's two-toed sloths showed low variation ranging only between 50 and 51. In contrast, Linnaeus's two-toed sloths appeared to vary widely, with diploid numbers ranging from 53 to 67, suggesting distinct geographic groups. The species identification method presented here represents a low-cost easy-to-use tool that will help to improve management of the captive population of two-toed sloths.

Identification of differences in human and great ape phytanic acid metabolism that could influence gene expression profiles and physiological functions

Background

It has been proposed that anatomical differences in human and great ape guts arose in response to species-specific diets and energy demands. To investigate functional genomic consequences of these differences, we compared their physiological levels of phytanic acid, a branched chain fatty acid that can be derived from the microbial degradation of chlorophyll in ruminant guts. Humans who accumulate large stores of phytanic acid commonly develop cerebellar ataxia, peripheral polyneuropathy, and retinitis pigmentosa in addition to other medical conditions. Furthermore, phytanic acid is an activator of the PPAR-alpha transcription factor that influences the expression of genes relevant to lipid metabolism.

Results

Despite their trace dietary phytanic acid intake, all great ape species had elevated red blood cell (RBC) phytanic acid levels relative to humans on diverse diets. Unlike humans, chimpanzees showed sexual dimorphism in RBC phytanic acid levels, which were higher in males relative to females. Cultured skin fibroblasts from all species had a robust capacity to degrade phytanic acid. We provide indirect evidence that great apes, in contrast to humans, derive significant amounts of phytanic acid from the hindgut fermentation of plant materials. This would represent a novel reduction of metabolic activity in humans relative to the great apes.

Conclusion

We identified differences in the physiological levels of phytanic acid in humans and great apes and propose this is causally related to their gut anatomies and microbiomes. Phytanic acid levels could contribute to cross-species and sex-specific differences in human and great ape transcriptomes, especially those related to lipid metabolism. Based on the medical conditions caused by phytanic acid accumulation, we suggest that differences in phytanic acid metabolism could influence the functions of human and great ape nervous, cardiovascular, and skeletal systems.

Detection of two distinct forms of apoC-I in great apes

ApoC-I, the smallest of the soluble apolipoproteins, associates with both TG-rich lipoproteins and HDL. Mass spectral analyses of human apoC-I previously had demonstrated that in the circulation there are two forms, either a 57 amino acid protein or a 55 amino acid protein, due to the loss of two amino acids from the N-terminus. In our analyses of the apolipoproteins of the other great apes by mass spectrometry, four forms of apoC-I were detected. Two of these showed a high degree of identity to the mature and truncated forms of human apoC-I. The other two were homologous to the virtual protein and its truncated form that are encoded by a human pseudogene. In humans, the genes for apoC-I and its pseudogene are located on chromosome 19, the pseudogene being 2.5 kb downstream from the apoC-I gene. Based on the similarity between the apoC-I gene and the pseudogene, it has been concluded that the latter arose from the former as a result of gene duplication approximately 35 million years ago. Interestingly, the virtual protein encoded by the pseudogene is acidic, not basic like apoC-I. In the chimpanzee, there also are two genes for apoC-I, the one upstream encodes a basic protein and the downstream gene, rather than being a pseudogene, encodes an acidic protein (P86336). In addition to reporting on the molecular masses of great ape apoC-I, we were able to clearly demonstrate by "Top-down" sequencing that the acidic form arose from a separate gene. In our analyses, we have measured the molecular masses of apoC-I associated with the HDL of the following great apes: bonobo (Pan paniscus), chimpanzee (Pan troglodytes), and the Sumatran orangutan (Pongo abelii). Genomic variations in chromosome 19 among great apes, baboons and macaques as they relate to both genes for apoC-I and the pseudogene are compared and discussed.

Phylogeny of the macaques (Cercopithecidae: Macaca) based on Alu elements

Genus Macaca (Cercopithecidae: Papionini) is one of the most successful primate radiations. Despite previous studies on morphology and mitochondrial DNA analysis, a number of issues regarding the details of macaque evolution remain unsolved. Alu elements are a class of non-autonomous retroposons belonging to short interspersed elements that are specific to the primate lineage. Because retroposon insertions show very little homoplasy, and because the ancestral state (absence of the SINE) is known, Alu elements are useful genetic markers and have been utilized for analyzing primate phylogenentic relationships and human population genetic relationships. Using PCR display methodology, 298 new Alu insertions have been identified from ten species of macaques. Together with 60 loci reported previously, a total of 358 loci are used to infer the phylogenetic relationships of genus Macaca. With regard to earlier unresolved issues on the macaque evolution, the topology of our tree suggests that: 1) genus Macaca contains four monophyletic species groups; 2) within the Asian macaques, the silenus group diverged first, and members of the sinica and fascicularis groups share a common ancestor; 3) Macaca arctoides are classified in the sinica group. Our results provide a robust molecular phylogeny for genus Macaca with stronger statistical support than previous studies. The present study also illustrates that SINE-based approaches are a powerful tool in primate phylogenetic studies and can be used to successfully resolve evolutionary relationships between taxa at scales from the ordinal level to closely related species within one genus.

The sequence and de novo assembly of the giant panda genome

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

Phylogeny of horse chromosome 5q in the genus Equus and centromere repositioning

Horses, asses and zebras belong to the genus Equus and are the only extant species of the family Equidae in the order Perissodactyla. In a previous work we demonstrated that a key factor in the rapid karyotypic evolution of this genus was evolutionary centromere repositioning, that is, the shift of the centromeric function to a new position without alteration of the order of markers along the chromosome. In search of previously undiscovered evolutionarily new centromeres, we traced the phylogeny of horse chromosome 5, analyzing the order of BAC markers, derived from a horse genomic library, in 7 Equus species (E. caballus, E. hemionus onager, E. kiang, E. asinus, E. grevyi, E. burchelli and E. zebra hartmannae). This analysis showed that repositioned centromeres are present in E. asinus (domestic donkey, EAS) chromosome 16 and in E. burchelli (Burchell's zebra, EBU) chromosome 17, confirming that centromere repositioning is a strikingly frequent phenomenon in this genus. The observation that the neocentromeres in EAS16 and EBU17 are in the same chromosomal position suggests that they may derive from the same event and therefore, E. asinus and E. burchelli may be more closely related than previously proposed; alternatively, 2 centromere repositioning events, involving the same chromosomal region, may have occurred independently in different lineages, pointing to the possible existence of hot spots for neocentromere formation. Our comparative analysis also showed that, while E. caballus chromosome 5 seems to represent the ancestral configuration, centric fission followed by independent fusion events gave rise to 3 different submetacentric chromosomes in other Equus lineages.

Mass spectral analyses of the two major apolipoproteins of great ape high density lipoproteins

The two major apolipoproteins associated with human and chimpanzee (Pan troglodytes) high density lipoproteins (HDL) are apoA-I and dimeric apoA-II. Although humans are closely related to great apes, apolipoprotein data do not exist for bonobos (Pan paniscus), western lowland gorillas (Gorilla gorilla gorilla) and the Sumatran orangutans (Pongo abelii). In the absence of any data, other great apes simply have been assumed to have dimeric apoA-II while other primates and most other mammals have been shown to have monomeric apoA-II. Using mass spectrometry, we have measured the molecular masses of apoA-I and apoA-II associated with the HDL of these great apes. Each was observed to have dimeric apoA-II. Being phylogenetically related, one would anticipate these apolipoproteins having a high percentage of invariant sequences when compared with human apolipoproteins. However, the orangutan, which diverged from the human lineage between 16 and 21 million years ago, had an apoA-II with the lowest monomeric mass, 8031.3 Da and the highest apoA-I value, 28,311.7 Da, currently reported for various mammals. Interestingly, the gorilla that diverged from the lineage leading to the human–chimpanzee branch after the orangutan had almost identical mass values to those reported for human apoA-I and apoA-II. But chimpanzee and the bonobo that diverged more recently had identical apoA-II mass values that were slightly larger than reported for the human apolipoprotein. The chimpanzee A-I mass values were very close to those of humans; however, the bonobo had values intermediate to the molecular masses of orangutan and the other great apes. With the already existing genomic data for chimpanzee and the recent entries for the orangutan and gorilla, we were able to demonstrate a close agreement between our mass spectral data and the calculated molecular weights determined from the predicted primary sequences of the respective apolipoproteins. Post-translational modification of these apolipoproteins, involving truncation and oxidation of methionine, are also reported.

Genome sequence, comparative analysis, and population genetics of the domestic horse

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

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.

Molecular decay of the tooth gene Enamelin (ENAM) mirrors the loss of enamel in the fossil record of placental mammals

Vestigial structures occur at both the anatomical and molecular levels, but studies documenting the co-occurrence of morphological degeneration in the fossil record and molecular decay in the genome are rare. Here, we use morphology, the fossil record, and phylogenetics to predict the occurrence of "molecular fossils" of the enamelin (ENAM) gene in four different orders of placental mammals (Tubulidentata, Pholidota, Cetacea, Xenarthra) with toothless and/or enamelless taxa. Our results support the "molecular fossil" hypothesis and demonstrate the occurrence of frameshift mutations and/or stop codons in all toothless and enamelless taxa. We then use a novel method based on selection intensity estimates for codons (omega) to calculate the timing of iterated enamel loss in the fossil record of aardvarks and pangolins, and further show that the molecular evolutionary history of ENAM predicts the occurrence of enamel in basal representatives of Xenarthra (sloths, anteaters, armadillos) even though frameshift mutations are ubiquitous in ENAM sequences of living xenarthrans. The molecular decay of ENAM parallels the morphological degeneration of enamel in the fossil record of placental mammals and provides manifest evidence for the predictive power of Darwin's theory.

The value of avian genomics to the conservation of wildlife

Background

Genomic studies in non-domestic avian models, such as the California condor and white-throated sparrow, can lead to more comprehensive conservation plans and provide clues for understanding mechanisms affecting genetic variation, adaptation and evolution.

Developing genomic tools and resources including genomic libraries and a genetic map of the California condor is a prerequisite for identification of candidate loci for a heritable embryonic lethal condition. The white-throated sparrow exhibits a stable genetic polymorphism (i.e. chromosomal rearrangements) associated with variation in morphology, physiology, and behavior (e.g., aggression, social behavior, sexual behavior, parental care).

In this paper we outline the utility of these species as well as report on recent advances in the study of their genomes.

Results

Genotyping of the condor resource population at 17 microsatellite loci provided a better assessment of the current population's genetic variation. Specific New World vulture repeats were found in the condor genome. Using condor BAC library and clones, chicken-condor comparative maps were generated. A condor fibroblast cell line transcriptome was characterized using the 454 sequencing technology.

Our karyotypic analyses of the sparrow in combination with other studies indicate that the rearrangements in both chromosomes 2^m and 3^a are complex and likely involve multiple inversions, interchromosomal linkage, and pleiotropy. At least a portion of the rearrangement in chromosome 2^m existed in the common ancestor of the four North American species of Zonotrichia, but not in the one South American species, and that the 2^m form, originally thought to be the derived condition, might actually be the ancestral one.

Conclusion

Mining and characterization of candidate loci in the California condor using molecular genetic and genomic techniques as well as linkage and comparative genomic mapping will eventually enable the identification of carriers of the chondrodystrophy allele, resulting in improved genetic management of this disease.

In the white-throated sparrow, genomic studies, combined with ecological data, will help elucidate the basis of genic selection in a natural population. Morphs of the sparrow provide us with a unique opportunity to study intraspecific genomic differences, which have resulted from two separate yet linked evolutionary trajectories. Such results can transform our understanding of evolutionary and conservation biology.

Isolation and characterization of 23 polymorphic microsatellite loci for a West Indian iguana (Cyclura pinguis) from the British Virgin Islands

Twenty-three polymorphic microsatellite markers were identified and characterized for Cyclura pinguis, a critically endangered species of lizard (Sauria: Iguanidae) native to Anegada Island in the British Virgin Islands. We examined variation at these loci for 39 C. pinguis, finding up to five alleles per locus and an average expected heterozygosity of 0.55. Allele frequency estimates for these microsatellite loci will be used to characterize genetic diversity of captive and wild C. pinguis populations and to estimate relatedness among adult iguanas at the San Diego Zoo that form the nucleus of a captive breeding programme for this critically endangered species.

New insights into the evolution of intronic sequences of the beta-fibrinogen gene and their application in reconstructing mustelid phylogeny

Mustelidae is the largest and most diverse family in the order Carnivora. The phylogenetic relationships among the subfamilies have especially long been a focus of study. Herein we are among the first to employ two new introns (4 and 7) of the nuclear beta-fibrinogen gene to clarify these enigmatic problems. In addition, two previously available nuclear (IRBP exon 1 and TTR intron 1) and one mt (ND2) data sets were also combined and analyzed simultaneously with the newly obtained sequence data in this study. Detailed characterizations of the two intronic regions not only reveal the remarkable occurrences of short interspersed element (SINE) insertion events, providing a new example supporting the attractive hypothesis that attrition of an earlier retroposition may offer a proper environment for successive retropositions by forming a "dimer-like" structure, but also demonstrate their utility in the resolution of mustelid phylogeny. All of our analyses confirm the assemblage of Mustelinae, Lutrinae, and Melinae with confidence; moreover, two clades within Mustelinae were clearly recognized, i.e., genera Mustela and Martes. Notably, genus Martes of Mustelinae was found to branch off first, followed by Melinae and then a clade containing Lutrinae and genus Mustela of Mustelinae, indicating paraphyly of Mustelinae. In addition, Mephitinae diverges before the other mustelids and the monophyletic Procyonidae in all cases, supporting its elevation to a separate family. Additional independent genetic markers are still in need to resolve the trichotomy among Mephitinae and the other two carnivoran clades, Ailuridae and Procyonidae/non-mephitine Mustelidae.

Horse domestication and conservation genetics of Przewalski's horse inferred from sex chromosomal and autosomal sequences

Despite their ability to interbreed and produce fertile offspring, there is continued disagreement about the genetic relationship of the domestic horse (Equus caballus) to its endangered wild relative, Przewalski's horse (Equus przewalskii). Analyses have differed as to whether or not Przewalski's horse is placed phylogenetically as a separate sister group to domestic horses. Because Przewalski's horse and domestic horse are so closely related, genetic data can also be used to infer domestication-specific differences between the two. To investigate the genetic relationship of Przewalski's horse to the domestic horse and to address whether evolution of the domestic horse is driven by males or females, five homologous introns (a total of approximately 3 kb) were sequenced on the X and Y chromosomes in two Przewalski's horses and three breeds of domestic horses: Arabian horse, Mongolian domestic horse, and Dartmoor pony. Five autosomal introns (a total of approximately 6 kb) were sequenced for these horses as well. The sequences of sex chromosomal and autosomal introns were used to determine nucleotide diversity and the forces driving evolution in these species. As a result, X chromosomal and autosomal data do not place Przewalski's horses in a separate clade within phylogenetic trees for horses, suggesting a close relationship between domestic and Przewalski's horses. It was also found that there was a lack of nucleotide diversity on the Y chromosome and higher nucleotide diversity than expected on the X chromosome in domestic horses as compared with the Y chromosome and autosomes. This supports the hypothesis that very few male horses along with numerous female horses founded the various domestic horse breeds. Patterns of nucleotide diversity among different types of chromosomes were distinct for Przewalski's in contrast to domestic horses, supporting unique evolutionary histories of the two species.

Analysis of complete mitochondrial genome sequences increases phylogenetic resolution of bears (Ursidae), a mammalian family that experienced rapid speciation

Background

Despite the small number of ursid species, bear phylogeny has long been a focus of study due to their conservation value, as all bear genera have been classified as endangered at either the species or subspecies level. The Ursidae family represents a typical example of rapid evolutionary radiation. Previous analyses with a single mitochondrial (mt) gene or a small number of mt genes either provide weak support or a large unresolved polytomy for ursids. We revisit the contentious relationships within Ursidae by analyzing complete mt genome sequences and evaluating the performance of both entire mt genomes and constituent mtDNA genes in recovering a phylogeny of extremely recent speciation events.

Results

This mitochondrial genome-based phylogeny provides strong evidence that the spectacled bear diverged first, while within the genus Ursus, the sloth bear is the sister taxon of all the other five ursines. The latter group is divided into the brown bear/polar bear and the two black bears/sun bear assemblages. These findings resolve the previous conflicts between trees using partial mt genes. The ability of different categories of mt protein coding genes to recover the correct phylogeny is concordant with previous analyses for taxa with deep divergence times. This study provides a robust Ursidae phylogenetic framework for future validation by additional independent evidence, and also has significant implications for assisting in the resolution of other similarly difficult phylogenetic investigations.

Conclusion

Identification of base composition bias and utilization of the combined data of whole mitochondrial genome sequences has allowed recovery of a strongly supported phylogeny that is upheld when using multiple alternative outgroups for the Ursidae, a mammalian family that underwent a rapid radiation since the mid- to late Pliocene. It remains to be seen if the reliability of mt genome analysis will hold up in studies of other difficult phylogenetic issues. Although the whole mitochondrial DNA sequence based phylogeny is robust, it remains in conflict with phylogenetic relationships suggested by analysis of limited nuclear-encoded data, a situation that will require gathering more nuclear DNA sequence information.

Primate TNF promoters reveal markers of phylogeny and evolution of innate immunity

Background

Tumor necrosis factor (TNF) is a critical cytokine in the immune response whose transcriptional activation is controlled by a proximal promoter region that is highly conserved in mammals and, in particular, primates. Specific single nucleotide polymorphisms (SNPs) upstream of the proximal human TNF promoter have been identified, which are markers of human ancestry.

Methodology/Principal findings

Using a comparative genomics approach we show that certain fixed genetic differences in the TNF promoter serve as markers of primate speciation. We also demonstrate that distinct alleles of most human TNF promoter SNPs are identical to fixed nucleotides in primate TNF promoters. Furthermore, we identify fixed genetic differences within the proximal TNF promoters of Asian apes that do not occur in African ape or human TNF promoters. Strikingly, protein-DNA binding assays and gene reporter assays comparing these Asian ape TNF promoters to African ape and human TNF promoters demonstrate that, unlike the fixed differences that we define that are associated with primate phylogeny, these Asian ape-specific fixed differences impair transcription factor binding at an Sp1 site and decrease TNF transcription induced by bacterial stimulation of macrophages.

Conclusions/significance

Here, we have presented the broadest interspecies comparison of a regulatory region of an innate immune response gene to date. We have characterized nucleotide positions in Asian ape TNF promoters that underlie functional changes in cell type- and stimulus-specific activation of the TNF gene. We have also identified ancestral TNF promoter nucleotide states in the primate lineage that correspond to human SNP alleles. These findings may reflect evolution of Asian and African apes under a distinct set of infectious disease pressures involving the innate immune response and TNF.

Construction of a California condor BAC library and first-generation chicken–condor comparative physical map as an endangered species conservation genomics resource

To support genomic analysis of the endangered California condor (Gymnogyps californianus), a BAC library (CHORI-262) was generated using DNA from the blood of a female. The library consists of 89,665 recombinant BAC clones providing approximately 14-fold coverage of the presumed approximately 1.48-Gb genome. Taking advantage of recent progress in chicken genomics, we developed a first-generation comparative chicken-condor physical map using an overgo hybridization approach. The overgos were derived from chicken (164 probes) and New World vulture (8 probes) sequences. Screening a 2.8x subset of the total library resulted in 236 BAC-gene assignments with 2.5 positive BAC clones per successful probe. A preliminary comparative chicken-condor BAC-based map included 93 genes. Comparison of selected condor BAC sequences with orthologous chicken sequences suggested a high degree of conserved synteny between the two avian genomes. This work will aid in identification and characterization of candidate loci for the chondrodystrophy mutation to advance genetic management of this disease.

Sonic Hedgehog, a key development gene, experienced intensified molecular evolution in primates

Sonic Hedgehog (SHH) is one of the most intensively studied genes in developmental biology. It is a highly conserved gene, found in species as diverse as arthropods and mammals. The mammalian SHH encodes a signaling molecule that plays a central role in developmental patterning, especially of the nervous system and the skeletal system. Here, we show that the molecular evolution of SHH is markedly accelerated in primates relative to other mammals. We further show that within primates, the acceleration is most prominent along the lineage leading to humans. Finally, we show that the acceleration in the lineage leading to humans is coupled with signatures of adaptive evolution. In particular, the lineage leading to humans is characterized by a rampant and statistically highly non-random gain of serines and threonines, residues that are potential substrates of post-translational modifications. This suggests that SHH might have evolved more complex post-translational regulation in the lineage leading to humans. Collectively, these findings implicate SHH as a potential contributor to the evolution of primate- or human-specific morphological traits in the nervous and/or skeletal systems and provide the impetus for additional studies aimed at identifying the primate- or human-specific functions of this key development gene.