DNA sequence and comparative analysis of chimpanzee chromosome 22

Transposable elements that have recently been mobile in the human genome

Background

Transposable elements (TE) comprise nearly half of the human genome and their insertions have profound effects to human genetic diversification and as well as disease. Despite their abovementioned significance, there is no consensus on the TE subfamilies that remain active in the human genome. In this study, we therefore developed a novel statistical test for recently mobile subfamilies (RMSs), based on patterns of overlap with > 100,000 polymorphic indels.

Results

Our analysis produced a catalogue of 20 high-confidence RMSs, which excludes many false positives in public databases. Intriguingly though, it includes HERV-K, an LTR subfamily previously thought to be extinct. The RMS catalogue is strongly enriched for contributions to germline genetic disorders (P = 1.1e-10), and thus constitutes a valuable resource for diagnosing disorders of unknown aetiology using targeted TE-insertion screens. Remarkably, RMSs are also highly enriched for somatic insertions in diverse cancers (P = 2.8e-17), thus indicating strong correlations between germline and somatic TE mobility. Using CRISPR/Cas9 deletion, we show that an RMS-derived polymorphic TE insertion increased the expression of RPL17, a gene associated with lower survival in liver cancer. More broadly, polymorphic TE insertions from RMSs were enriched near genes with allele-specific expression, suggesting widespread effects on gene regulation.

Conclusions

By using a novel statistical test we have defined a catalogue of 20 recently mobile transposable element subfamilies. We illustrate the gene regulatory potential of RMS-derived polymorphic TE insertions, using CRISPR/Cas9 deletion in vitro on a specific candidate, as well as by genome wide analysis of allele-specific expression. Our study presents novel insights into TE mobility and regulatory potential and provides a key resource for human disease genetics and population history studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08085-0.

Multiple Genomic Events Altering Hominin SIGLEC Biology and Innate Immunity Predated the Common Ancestor of Humans and Archaic Hominins

Human-specific pseudogenization of the CMAH gene eliminated the mammalian sialic acid (Sia) Neu5Gc (generating an excess of its precursor Neu5Ac), thus changing ubiquitous cell surface “self-associated molecular patterns” that modulate innate immunity via engagement of CD33-related-Siglec receptors. The Alu-fusion-mediated loss-of-function of CMAH fixed ∼2–3 Ma, possibly contributing to the origins of the genus Homo. The mutation likely altered human self-associated molecular patterns, triggering multiple events, including emergence of human-adapted pathogens with strong preference for Neu5Ac recognition and/or presenting Neu5Ac-containing molecular mimics of human glycans, which can suppress immune responses via CD33-related-Siglec engagement. Human-specific alterations reported in some gene-encoding Sia-sensing proteins suggested a “hotspot” in hominin evolution. The availability of more hominid genomes including those of two extinct hominins now allows full reanalysis and evolutionary timing. Functional changes occur in 8/13 members of the human genomic cluster encoding CD33-related Siglecs, all predating the human common ancestor. Comparisons with great ape genomes indicate that these changes are unique to hominins. We found no evidence for strong selection after the Human–Neanderthal/Denisovan common ancestor, and these extinct hominin genomes include almost all major changes found in humans, indicating that these changes in hominin sialobiology predate the Neanderthal–human divergence ∼0.6 Ma. Multiple changes in this genomic cluster may also explain human-specific expression of CD33rSiglecs in unexpected locations such as amnion, placental trophoblast, pancreatic islets, ovarian fibroblasts, microglia, Natural Killer(NK) cells, and epithelia. Taken together, our data suggest that innate immune interactions with pathogens markedly altered hominin Siglec biology between 0.6 and 2 Ma, potentially affecting human evolution.

Social Networking of Quasi-Species Consortia drive Virolution via Persistence

The emergence of cooperative quasi-species consortia (QS-C) thinking from the more accepted quasispecies equations of Manfred Eigen, provides a conceptual foundation from which concerted action of RNA agents can now be understood. As group membership becomes a basic criteria for the emergence of living systems, we also start to understand why the history and context of social RNA networks become crucial for survival and function. History and context of social RNA networks also lead to the emergence of a natural genetic code. Indeed, this QS-C thinking can also provide us with a transition point between the chemical world of RNA replicators and the living world of RNA agents that actively differentiate self from non-self and generate group identity with membership roles. Importantly the social force of a consortia to solve complex, multilevel problems also depend on using opposing and minority functions. The consortial action of social networks of RNA stem-loops subsequently lead to the evolution of cellular organisms representing a tree of life.

Epigenetic regulation during human cortical development: Seq-ing answers from the brain to the organoid

Epigenetic regulation plays an important role in controlling gene expression during complex processes, such as development of the human brain. Mutations in genes encoding chromatin modifying proteins and in the non-protein coding sequences of the genome can potentially alter transcription factor binding or chromatin accessibility. Such mutations can frequently cause neurodevelopmental disorders, therefore understanding how epigenetic regulation shapes brain development is of particular interest. While epigenetic regulation of neural development has been extensively studied in murine models, significant species-specific differences in both the genome sequence and in brain development necessitate human models. However, access to human fetal material is limited and these tissues cannot be grown or experimentally manipulated ex vivo. Therefore, models that recapitulate particular aspects of human fetal brain development, such as the in vitro differentiation of human pluripotent stem cells (hPSCs), are instrumental for studying the epigenetic regulation of human neural development. Here, we examine recent studies that have defined changes in the epigenomic landscape during fetal brain development. We compare these studies with analogous data derived by in vitro differentiation of hPSCs into specific neuronal cell types or as three-dimensional cerebral organoids. Such comparisons can be informative regarding which aspects of fetal brain development are faithfully recapitulated by in vitro differentiation models and provide a foundation for using experimentally tractable in vitro models of human brain development to study neural gene regulation and the basis of its disruption to cause neurodevelopmental disorders.

Nucleotide sequencing of the HoxA gene cluster using Gorilla fosmid clones

We sequenced the western gorilla (Gorilla gorilla) HoxA cluster region using seven fosmid clones, and found that the total tiling path sequence was 214,185 bp from the 5' non-genic region of HoxA1 to the 3' non-genic region of Evx1. We compared the nucleotide sequence with the gorilla genome sequence in the NCBI database, and the overall proportion of nucleotide difference was estimated to be 0.0005-0.0007. These estimates are lower than overall genomic polymorphism in gorillas.

Two base pair deletion in IL2 receptor γ gene in NOD/SCID mice induces a highly severe immunodeficiency

Genome editing has recently emerged as a powerful tool for generating mutant mice. Small deletions of nucleotides in the target genes are frequently found in CRISPR/Cas9 mediated mutant mice. However, there are very few reports analyzing the phenotypes in small deleted mutant mice generated by CRISPR/Cas9. In this study, we generated a mutant by microinjecting sgRNAs targeting the IL2 receptor γ gene and Cas9 protein, into the cytoplasm of IVF-derived NOD.CB17/Prkdcscid/JKrb (NOD/SCID) mice embryos, and further investigated whether a 2 bp deletion of the IL2 receptor γ gene affects severe deficiency of immune cells as seen in NOD/LtSz-scid IL2 receptor γ^−/− (NSG) mice. Our results show that the thymus weight of mutant mice is significantly less than that of NOD/SCID mice, whereas the spleen weight was marginally less. T and B cells in the mutant mice were severely deficient, and NK cells were almost absent. In addition, tumor growth was exceedingly increased in the mutant mice transplanted with HepG2, Raji and A549 cells, but not in nude and NOD/SCID mice. These results suggest that the NOD/SCID mice with deletion of 2 bp in the IL2 receptor γ gene shows same phenotype as NSG mice. Taken together, our data indicates that small deletions by genome editing is sufficient to generate null mutant mice.

The Use of Non-Human Primates in Biological and Medical Research: Evidence Submitted by FRAME to the Academy of Medical Sciences/Medical Research Council/Royal Society/Wellcome Trust Working Group

The Academy of Medical Sciences, the Medical Research Council, the Royal Society and the Wellcome Trust are undertaking a study into the use of non-human primates in biological and medical research. An independent working group of scientific experts, led by Sir David Weatherall, aims to produce a report summarising the findings of this study, early in 2006. The trends in primate research, and the nature and effects of recent and proposed changes in the global use of non-human primates in research, will be investigated. The associated ethical, welfare and regulatory issues, and the role and impact of the Three Rs principles of refinement, reduction and replacement will also be reviewed. As part of this study, a call for evidence was made. The evidence submitted by FRAME emphasised that the use of non-human primates for fundamental research or for regulatory testing still fails to take into account the fact that, although non-human primates are anatomically and physiologically similar to humans, they are not necessarily relevant models for studies on human disease or human physiology. FRAME continues to believe that we have a duty to ensure that these animals are not used without overwhelming evidence that they are the only suitable and relevant models for use in work of undeniable significance.

A Japanese history of the Human Genome Project

The Human Genome Project (HGP) is one of the most important international achievements in life sciences, to which Japanese scientists made remarkable contributions. In the early 1980s, Akiyoshi Wada pioneered the first project for the automation of DNA sequencing technology. Ken-ichi Matsubara exhibited exceptional leadership to launch the comprehensive human genome program in Japan. Hideki Kambara made a major contribution by developing a key device for high-speed DNA sequencers, which enabled scientists to construct human genome draft sequences. The RIKEN team led by Yoshiyuki Sakaki (the author) played remarkable roles in the draft sequencing and completion of chromosomes 21, 18, and 11. Additionally, the Keio University team led by Nobuyoshi Shimizu made noteworthy contributions to the completion of chromosomes 22, 21, and 8. In April 2003, the Japanese team joined the international consortium in declaring the completion of the human genome sequence. Consistent with the HGP mandate, Japan has successfully developed a wide range of ambitious genomic sciences.

Differential Brain, Cognitive and Motor Profiles Associated with Partial Trisomy. Modeling Down Syndrome in Mice

We hypothesize that the trisomy 21 (Down syndrome) is the additive and interactive outcome of the triple copy of different regions of HSA21. Because of the small number of patients with partial trisomy 21, we addressed the question in the Mouse in which three chromosomal regions located on MMU10, MMU17 and MMU16 carries almost all the HSA21 homologs. Male mice from four segmental trisomic strains covering the D21S17-ETS2 (syntenic to MMU16) were examined with an exhaustive battery of cognitive tests, motor tasks and MRI and compared with TS65Dn that encompasses D21S17-ETS2. None of the four strains gather all the impairments (measured by the effect size) of TS65Dn strain. The 152F7 strain was close to TS65Dn for motor behavior and reference memory and the three other strains 230E8, 141G6 and 285E6 for working memory. Episodic memory was impaired only in strain 285E6. The hippocampus and cerebellum reduced sizes that were seen in all the strains indicate that trisomy 21 is not only a hippocampus syndrome but that it results from abnormal interactions between the two structures.

Application of Monte Carlo simulation in addressing key issues of complex coacervation formed by polyelectrolytes and oppositely charged colloids

This paper reviews the recent advance of Monte Carlo (MC) simulation in addressing key issues of complex coacervation between polyelectrolytes and oppositely charged colloids. Readers were first supplied with a brief overview of current knowledge and experimental strategies in the study of complex coacervation. In the next section, the general MC simulation procedures as well as representative strategies applied in complex coacervation were summarized. The unique contributions of MC simulation in either capturing delicate features, easing the experimental trials or proving the concept were then elucidated through the following aspects: i) identify phase boundary and decouple interaction contributions; ii) clarify composition distribution and internal structure; iii) predict the influences of physicochemical conditions on complex coacervation; iv) delineate the mechanisms for "binding on the wrong side of the isoelectric point". Finally, current challenges as well as prospects of MC simulation in complex coacervation are also discussed. The ultimate goal of this review is to provide readers with basic guideline for synergistic design of experiments in combination with MC simulation, and deliver convincing interpretation and reliable prediction for the structure and behavior in polyelectrolyte-macroion complex coacervation.

General continuous-time Markov model of sequence evolution via insertions/deletions: local alignment probability computation

Background

Insertions and deletions (indels) account for more nucleotide differences between two related DNA sequences than substitutions do, and thus it is imperative to develop a method to reliably calculate the occurrence probabilities of sequence alignments via evolutionary processes on an entire sequence. Previously, we presented a perturbative formulation that facilitates the ab initio calculation of alignment probabilities under a continuous-time Markov model, which describes the stochastic evolution of an entire sequence via indels with quite general rate parameters. And we demonstrated that, under some conditions, the ab initio probability of an alignment can be factorized into the product of an overall factor and contributions from regions (or local alignments) delimited by gapless columns.

Results

Here, using our formulation, we attempt to approximately calculate the probabilities of local alignments under space-homogeneous cases. First, for each of all types of local pairwise alignments (PWAs) and some typical types of local multiple sequence alignments (MSAs), we numerically computed the total contribution from all parsimonious indel histories and that from all next-parsimonious histories, and compared them. Second, for some common types of local PWAs, we derived two integral equation systems that can be numerically solved to give practically exact solutions. We compared the total parsimonious contribution with the practically exact solution for each such local PWA. Third, we developed an algorithm that calculates the first-approximate MSA probability by multiplying total parsimonious contributions from all local MSAs. Then we compared the first-approximate probability of each local MSA with its absolute frequency in the MSAs created via a genuine sequence evolution simulator, Dawg. In all these analyses, the total parsimonious contributions approximated the multiplication factors fairly well, as long as gap sizes and branch lengths are at most moderate. Examination of the accuracy of another indel probabilistic model in the light of our formulation indicated some modifications necessary for the model’s accuracy improvement.

Conclusions

At least under moderate conditions, the approximate methods can quite accurately calculate ab initio alignment probabilities under biologically more realistic models than before. Thus, our formulation will provide other indel probabilistic models with a sound reference point.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1167-6) contains supplementary material, which is available to authorized users.

General continuous-time Markov model of sequence evolution via insertions/deletions: are alignment probabilities factorable?

Background

Insertions and deletions (indels) account for more nucleotide differences between two related DNA sequences than substitutions do, and thus it is imperative to develop a stochastic evolutionary model that enables us to reliably calculate the probability of the sequence evolution through indel processes. Recently, indel probabilistic models are mostly based on either hidden Markov models (HMMs) or transducer theories, both of which give the indel component of the probability of a given sequence alignment as a product of either probabilities of column-to-column transitions or block-wise contributions along the alignment. However, it is not a priori clear how these models are related with any genuine stochastic evolutionary model, which describes the stochastic evolution of an entire sequence along the time-axis. Moreover, currently none of these models can fully accommodate biologically realistic features, such as overlapping indels, power-law indel-length distributions, and indel rate variation across regions.

Results

Here, we theoretically dissect the ab initio calculation of the probability of a given sequence alignment under a genuine stochastic evolutionary model, more specifically, a general continuous-time Markov model of the evolution of an entire sequence via insertions and deletions. Our model is a simple extension of the general “substitution/insertion/deletion (SID) model”. Using the operator representation of indels and the technique of time-dependent perturbation theory, we express the ab initio probability as a summation over all alignment-consistent indel histories. Exploiting the equivalence relations between different indel histories, we find a “sufficient and nearly necessary” set of conditions under which the probability can be factorized into the product of an overall factor and the contributions from regions separated by gapless columns of the alignment, thus providing a sort of generalized HMM. The conditions distinguish evolutionary models with factorable alignment probabilities from those without ones. The former category includes the “long indel” model (a space-homogeneous SID model) and the model used by Dawg, a genuine sequence evolution simulator.

Conclusions

With intuitive clarity and mathematical preciseness, our theoretical formulation will help further advance the ab initio calculation of alignment probabilities under biologically realistic models of sequence evolution via indels.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1105-7) contains supplementary material, which is available to authorized users.

No Distinction of Orthology/Paralogy between Human and Chimpanzee Rh Blood Group Genes

On human (Homo sapiens) chromosome 1, there is a tandem duplication encompassing Rh blood group genes (Hosa_RHD and Hosa_RHCE). This duplication occurred in the common ancestor of humans, chimpanzees (Pan troglodytes), and gorillas, after splitting from their common ancestor with orangutans. Although several studies have been conducted on ape Rh blood group genes, the clear genome structures of the gene clusters remain unknown. Here, we determined the genome structure of the gene cluster of chimpanzee Rh genes by sequencing five BAC (Bacterial Artificial Chromosome) clones derived from chimpanzees. We characterized three complete loci (Patr_RHα, Patr_RHβ, and Patr_RHγ). In the Patr_RHβ locus, a short version of the gene, which lacked the middle part containing exons 4-8, was observed. The Patr_RHα and Patr_RHβ genes were located on the locations corresponding to Hosa_RHD and Hosa_RHCE, respectively, and Patr_RHγ was in the immediate vicinity of Patr_RHβ. Sequence comparisons revealed high sequence similarity between Patr_RHβ and Hosa_RHCE, while the chimpanzee Rh gene closest to Hosa_RHD was not Patr_RHα but rather Patr_RHγ. The results suggest that rearrangements and gene conversions frequently occurred between these genes and that the classic orthology/paralogy dichotomy no longer holds between human and chimpanzee Rh blood group genes.

Regulatory Divergence of Transcript Isoforms in a Mammalian Model System

Phenotypic differences between species are driven by changes in gene expression and, by extension, by modifications in the regulation of the transcriptome. Investigation of mammalian transcriptome divergence has been restricted to analysis of bulk gene expression levels and gene-internal splicing. Using allele-specific expression analysis in inter-strain hybrids of Mus musculus, we determined the contribution of multiple cellular regulatory systems to transcriptome divergence, including: alternative promoter usage, transcription start site selection, cassette exon usage, alternative last exon usage, and alternative polyadenylation site choice. Between mouse strains, a fifth of genes have variations in isoform usage that contribute to transcriptomic changes, half of which alter encoded amino acid sequence. Virtually all divergence in isoform usage altered the post-transcriptional regulatory instructions in gene UTRs. Furthermore, most genes with isoform differences between strains contain changes originating from multiple regulatory systems. This result indicates widespread cross-talk and coordination exists among different regulatory systems. Overall, isoform usage diverges in parallel with and independently to gene expression evolution, and the cis and trans regulatory contribution to each differs significantly.

High-throughput RNA sequencing reveals structural differences of orthologous brain-expressed genes between western lowland gorillas and humans

The human brain and human cognitive abilities are strikingly different from those of other great apes despite relatively modest genome sequence divergence. However, little is presently known about the interspecies divergence in gene structure and transcription that might contribute to these phenotypic differences. To date, most comparative studies of gene structure in the brain have examined humans, chimpanzees, and macaque monkeys. To add to this body of knowledge, we analyze here the brain transcriptome of the western lowland gorilla (Gorilla gorilla gorilla), an African great ape species that is phylogenetically closely related to humans, but with a brain that is approximately one-third the size. Manual transcriptome curation from a sample of the planum temporale region of the neocortex revealed 12 protein-coding genes and one noncoding-RNA gene with exons in the gorilla unmatched by public transcriptome data from the orthologous human loci. These interspecies gene structure differences accounted for a total of 134 amino acids in proteins found in the gorilla that were absent from protein products of the orthologous human genes. Proteins varying in structure between human and gorilla were involved in immunity and energy metabolism, suggesting their relevance to phenotypic differences. This gorilla neocortical transcriptome comprises an empirical, not homology- or prediction-driven, resource for orthologous gene comparisons between human and gorilla. These findings provide a unique repository of the sequences and structures of thousands of genes transcribed in the gorilla brain, pointing to candidate genes that may contribute to the traits distinguishing humans from other closely related great apes.

Germline mutation rates and the long-term phenotypic effects of mutation accumulation in wild-type laboratory mice and mutator mice

The germline mutation rate is an important parameter that affects the amount of genetic variation and the rate of evolution. However, neither the rate of germline mutations in laboratory mice nor the biological significance of the mutation rate in mammalian populations is clear. Here we studied genome-wide mutation rates and the long-term effects of mutation accumulation on phenotype in more than 20 generations of wild-type C57BL/6 mice and mutator mice, which have high DNA replication error rates. We estimated the base-substitution mutation rate to be 5.4 × 10⁻⁹ (95% confidence interval = 4.6 × 10⁻⁹–6.5 × 10⁻⁹) per nucleotide per generation in C57BL/6 laboratory mice, about half the rate reported in humans. The mutation rate in mutator mice was 17 times that in wild-type mice. Abnormal phenotypes were 4.1-fold more frequent in the mutator lines than in the wild-type lines. After several generations, the mutator mice reproduced at substantially lower rates than the controls, exhibiting low pregnancy rates, lower survival rates, and smaller litter sizes, and many of the breeding lines died out. These results provide fundamental information about mouse genetics and reveal the impact of germline mutation rates on phenotypes in a mammalian population.

The pattern of DNA cleavage intensity around indels

Indels (insertions and deletions) are the second most common form of genetic variations in the eukaryotic genomes and are responsible for a multitude of genetic diseases. Despite its significance, detailed molecular mechanisms for indel generation are still unclear. Here we examined 2,656,597 small human and mouse germline indels, 16,742 human somatic indels, 10,599 large human insertions, and 5,822 large chimpanzee insertions and systematically analyzed the patterns of DNA cleavage intensities in the 200 base pair regions surrounding these indels. Our results show that DNA cleavage intensities close to the start and end points of indels are significantly lower than other regions, for both small human germline and somatic indels and also for mouse small indels. Compared to small indels, the patterns of DNA cleavage intensity around large indels are more complex, and there are two low intensity regions near each end of the indels that are approximately 13 bp apart from each other. Detailed analyses of a subset of indels show that there is slight difference in cleavage intensity distribution between insertion indels and deletion indels that could be contributed by their respective enrichment of different repetitive elements. These results will provide new insight into indel generation mechanisms.

The role of DNA insertions in phenotypic differentiation between humans and other primates

What makes us human is one of the most interesting and enduring questions in evolutionary biology. To assist in answering this question, we have identified insertions in the human genome which cannot be found in five comparison primate species: Chimpanzee, gorilla, orangutan, gibbon, and macaque. A total of 21,269 nonpolymorphic human-specific insertions were identified, of which only 372 were found in exons. Any function conferred by the remaining 20,897 is likely to be regulatory. Many of these insertions are likely to have been fitness neutral; however, a small number has been identified in genes showing signs of positive selection. Insertions found within positively selected genes show associations to neural phenotypes, which were also enriched in the whole data set. Other phenotypes that are found to be enriched in the data set include dental and sensory perception-related phenotypes, features which are known to differ between humans and other apes. The analysis provides several likely candidates, either genes or regulatory regions, which may be involved in the processes that differentiate humans from other apes.

Endogenous retrovirus-mediated genomic variations in chimpanzees

Transposable elements (TEs) have played a significant role in the evolution of host genome by triggering genomic rearrangements. TEs have been studied in various research fields, ranging from population genomics to personalized medicines. Human-specific TEs and TEs existing in the human genome have been well studied. Unlike them, non-human primate-specific TEs remain shrouded in mystery. However, the study of TE-mediated genomic or genetic variations through comparative genomics is essential to understand mechanisms which TEs utilize to modify species-specific genome architecture and to cause species-specific diseases, Therefore, we have studied chimpanzee-specific TEs as well as human-specific TEs. At first, we identified human-specific HERV-K integrated into the human genome after the divergence of human and chimpanzee. Then, for a comparative study of HERV-Ks and non-human ERVs, we extracted chimpanzee-specific endogenous retroviruses (PtERVs) from the chimpanzee genome. We identified 256 chimpanzee-specific PtERVs and characterized them, focusing on their estimated evolutionary age, polymorphism level in chimpanzee populations, and potential impact on the difference between the human and chimpanzee genomes.

Chimpanzee-specific endogenous retrovirus generates genomic variations in the chimpanzee genome

Endogenous retroviruses (ERVs), eukaryotic transposable elements, exist as proviruses in vertebrates including primates and contribute to genomic changes during the evolution of their host genomes. Many studies about ERVs have focused on the elements residing in the human genome but only a few studies have focused on the elements which exist in non-human primate genomes. In this study, we identified 256 chimpanzee-specific endogenous retrovirus copies (PtERVs: Pan troglodyte endogenous retroviruses) from the chimpanzee reference genome sequence through comparative genomics. Among the chimpanzee-specific ERV copies, 121 were full-length chimpanzee-specific ERV elements while 110 were chimpanzee-specific solitary LTR copies. In addition, we found eight potential retrotransposition-competent full-length chimpanzee-specific ERV copies containing an intact env gene, and two of them were polymorphic in chimpanzee individuals. Through computational analysis and manual inspection, we found that some of the chimpanzee-specific ERVs have propagated via non-classical PtERV insertion (NCPI), and at least one of the PtERVs may have played a role in creating an alternative transcript of a chimpanzee gene. Based on our findings in this study, we state that the chimpanzee-specific ERV element is one of the sources of chimpanzee genomic variations, some of which might be related to the alternative transcripts in the chimpanzee population.

Human-specific epigenetic variation in the immunological Leukotriene B4 Receptor (LTB4R/BLT1) implicated in common inflammatory diseases

BACKGROUND

Common human diseases are caused by the complex interplay of genetic susceptibility as well as environmental factors. Due to the environment's influence on the epigenome, and therefore genome function, as well as conversely the genome's facilitative effect on the epigenome, analysis of this level of regulation may increase our knowledge of disease pathogenesis.

METHODS

In order to identify human-specific epigenetic influences, we have performed a novel genome-wide DNA methylation analysis comparing human, chimpanzee and rhesus macaque.

RESULTS

We have identified that the immunological Leukotriene B4 receptor (LTB4R, BLT1 receptor) is the most epigenetically divergent human gene in peripheral blood in comparison with other primates. This difference is due to the co-ordinated active state of human-specific hypomethylation in the promoter and human-specific increased gene body methylation. This gene is significant in innate immunity and the LTB4/LTB4R pathway is involved in the pathogenesis of the spectrum of human inflammatory diseases. This finding was confirmed by additional neutrophil-only DNA methylome and lymphoblastoid H3K4me3 chromatin comparative data. Additionally we show through functional analysis that this receptor has increased expression and a higher response to the LTB4 ligand in human versus rhesus macaque peripheral blood mononuclear cells. Genome-wide we also find human species-specific differentially methylated regions (human s-DMRs) are more prevalent in CpG island shores than within the islands themselves, and within the latter are associated with the CTCF motif.

CONCLUSIONS

This result further emphasises the exclusive nature of the human immunological system, its divergent adaptation even from very closely related primates, and the power of comparative epigenomics to identify and understand human uniqueness.

Coelacanth genomes reveal signatures for evolutionary transition from water to land

Coelacanths are known as “living fossils,” as they show remarkable morphological resemblance to the fossil record and belong to the most primitive lineage of living Sarcopterygii (lobe-finned fishes and tetrapods). Coelacanths may be key to elucidating the tempo and mode of evolution from fish to tetrapods. Here, we report the genome sequences of five coelacanths, including four Latimeria chalumnae individuals (three specimens from Tanzania and one from Comoros) and one L. menadoensis individual from Indonesia. These sequences cover two African breeding populations and two known extant coelacanth species. The genome is ∼2.74 Gbp and contains a high proportion (∼60%) of repetitive elements. The genetic diversity among the individuals was extremely low, suggesting a small population size and/or a slow rate of evolution. We found a substantial number of genes that encode olfactory and pheromone receptors with features characteristic of tetrapod receptors for the detection of airborne ligands. We also found that limb enhancers of bmp7 and gli3, both of which are essential for limb formation, are conserved between coelacanth and tetrapods, but not ray-finned fishes. We expect that some tetrapod-like genes may have existed early in the evolution of primitive Sarcopterygii and were later co-opted to adapt to terrestrial environments. These coelacanth genomes will provide a cornerstone for studies to elucidate how ancestral aquatic vertebrates evolved into terrestrial animals.

Complexity in cancer biology: is systems biology the answer?

Complex phenotypes emerge from the interactions of thousands of macromolecules that are organized in multimolecular complexes and interacting functional modules. In turn, modules form functional networks in health and disease. Omics approaches collect data on changes for all genes and proteins and statistical analysis attempts to uncover the functional modules that perform the functions that characterize higher levels of biological organization. Systems biology attempts to transcend the study of individual genes/proteins and to integrate them into higher order information. Cancer cells exhibit defective genetic and epigenetic networks formed by altered complexes and network modules arising in different parts of tumor tissues that sustain autonomous cell behavior which ultimately lead tumor growth. We suggest that an understanding of tumor behavior must address not only molecular but also, and more importantly, tumor cell heterogeneity, by considering cancer tissue genetic and epigenetic networks, by characterizing changes in the types, composition, and interactions of complexes and networks in the different parts of tumor tissues, and by identifying critical hubs that connect them in time and space.

Transposable elements are a significant contributor to tandem repeats in the human genome

Sequence repeats are an important phenomenon in the human genome, playing important roles in genomic alteration often with phenotypic consequences. The two major types of repeat elements in the human genome are tandem repeats (TRs) including microsatellites, minisatellites, and satellites and transposable elements (TEs). So far, very little has been known about the relationship between these two types of repeats. In this study, we identified TRs that are derived from TEs either based on sequence similarity or overlapping genomic positions. We then analyzed the distribution of these TRs among TE families/subfamilies. Our study shows that at least 7,276 TRs or 23% of all minisatellites/satellites is derived from TEs, contributing ∼0.32% of the human genome. TRs seem to be generated more likely from younger/more active TEs, and once initiated they are expanded with time via local duplication of the repeat units. The currently postulated mechanisms for origin of TRs can explain only 6% of all TE-derived TRs, indicating the presence of one or more yet to be identified mechanisms for the initiation of such repeats. Our result suggests that TEs are contributing to genome expansion and alteration not only by transposition but also by generating tandem repeats.

Is the Y chromosome disappearing?--both sides of the argument

On August 31, 2011 at the 18th International Chromosome Conference in Manchester, Jenny Graves took on Jenn Hughes to debate the demise (or otherwise) of the mammalian Y chromosome. Sex chromosome evolution is an example of convergence; there are numerous examples of XY and ZW systems with varying degrees of differentiation and isolated examples of the Y disappearing in some lineages. It is agreed that the Y was once genetically identical to its partner and that the present-day human sex chromosomes retain only traces of their shared ancestry. The euchromatic portion of the male-specific region of the Y is ~1/6 of the size of the X and has only ~1/12 the number of genes. The big question however is whether this degradation will continue or whether it has reached a point of equilibrium. Jenny Graves argued that the Y chromosome is subject to higher rates of variation and inefficient selection and that Ys (and Ws) degrade inexorably. She argued that there is evidence that the Y in other mammals has undergone lineage-specific degradation and already disappeared in some rodent lineages. She also pointed out that there is practically nothing left of the original human Y and the added part of the human Y is degrading rapidly. Jenn Hughes on the other hand argued that the Y has not disappeared yet and it has been around for hundreds of millions of years. She stated that it has shown that it can outsmart genetic decay in the absence of "normal" recombination and that most of its genes on the human Y exhibit signs of purifying selection. She noted that it has added at least eight different genes, many of which have subsequently expanded in copy number, and that it has not lost any genes since the human and chimpanzee diverged ~6 million years ago. The issue was put to the vote with an exact 50/50 split among the opinion of the audience; an interesting (though perhaps not entirely unexpected) skew however was noted in the sex ratio of those for and against the notion.

Lessons from chimpanzee-based research on human disease: the implications of genetic differences

Assertions that the use of chimpanzees to investigate human diseases is valid scientifically are frequently based on a reported 98-99% genetic similarity between the species. Critical analyses of the relevance of chimpanzee studies to human biology, however, indicate that this genetic similarity does not result in sufficient physiological similarity for the chimpanzee to constitute a good model for research, and furthermore, that chimpanzee data do not translate well to progress in clinical practice for humans. Leading examples include the minimal citations of chimpanzee research that is relevant to human medicine, the highly different pathology of HIV/AIDS and hepatitis C virus infection in the two species, the lack of correlation in the efficacy of vaccines and treatments between chimpanzees and humans, and the fact that chimpanzees are not useful for research on human cancer. The major molecular differences underlying these inter-species phenotypic disparities have been revealed by comparative genomics and molecular biology - there are key differences in all aspects of gene expression and protein function, from chromosome and chromatin structure to post-translational modification. The collective effects of these differences are striking, extensive and widespread, and they show that the superficial similarity between human and chimpanzee genetic sequences is of little consequence for biomedical research. The extrapolation of biomedical data from the chimpanzee to the human is therefore highly unreliable, and the use of the chimpanzee must be considered of little value, particularly given the breadth and potential of alternative methods of enquiry that are currently available to science.

Human blood identification using the genome profiling method

In criminal investigations, usually it is necessary to identify whether blood spots found at crime scenes are from humans or not. Nowadays, immunohistochemical methods and DNA analysis are usually used for this purpose. However, such methods and DNA analysis are labor intensive and expensive, and require highly trained skilled technicians. Recently, the genome profiling method (GP method) was developed. However, its use as a human DNA analysis method has not been reported. In this report, an attempt was made to differentiate human blood samples from animal blood samples using the GP method for forensic purposes. DNA extracted from a rat, squirrel, cat, dog, cow, and antelope along with human blood samples were analyzed. Following cluster analysis the human samples clustered into a single group separate from the animal samples. Therefore, although the number of samples was small the results suggest that the GP method might enable us to differentiate human samples from various animal samples. It may become a powerful tool in the field of forensic science.

Characterization and potential functional significance of human-chimpanzee large INDEL variation

BACKGROUND

Although humans and chimpanzees have accumulated significant differences in a number of phenotypic traits since diverging from a common ancestor about six million years ago, their genomes are more than 98.5% identical at protein-coding loci. This modest degree of nucleotide divergence is not sufficient to explain the extensive phenotypic differences between the two species. It has been hypothesized that the genetic basis of the phenotypic differences lies at the level of gene regulation and is associated with the extensive insertion and deletion (INDEL) variation between the two species. To test the hypothesis that large INDELs (80 to 12,000 bp) may have contributed significantly to differences in gene regulation between the two species, we categorized human-chimpanzee INDEL variation mapping in or around genes and determined whether this variation is significantly correlated with previously determined differences in gene expression.

RESULTS

Extensive, large INDEL variation exists between the human and chimpanzee genomes. This variation is primarily attributable to retrotransposon insertions within the human lineage. There is a significant correlation between differences in gene expression and large human-chimpanzee INDEL variation mapping in genes or in proximity to them.

CONCLUSIONS

The results presented herein are consistent with the hypothesis that large INDELs, particularly those associated with retrotransposons, have played a significant role in human-chimpanzee regulatory evolution.

Genome-based analysis of the nonhuman primate Macaca fascicularis as a model for drug safety assessment

The long-tailed macaque, also referred to as cynomolgus monkey (Macaca fascicularis), is one of the most important nonhuman primate animal models in basic and applied biomedical research. To improve the predictive power of primate experiments for humans, we determined the genome sequence of a Macaca fascicularis female of Mauritian origin using a whole-genome shotgun sequencing approach. We applied a template switch strategy that uses either the rhesus or the human genome to assemble sequence reads. The sixfold sequence coverage of the draft genome sequence enabled discovery of about 2.1 million potential single-nucleotide polymorphisms based on occurrence of a dimorphic nucleotide at a given position in the genome sequence. Homology-based annotation allowed us to identify 17,387 orthologs of human protein-coding genes in the M. fascicularis draft genome, and the predicted transcripts enabled the design of a M. fascicularis-specific gene expression microarray. Using liver samples from 36 individuals of different geographic origin we identified 718 genes with highly variable expression in liver, whereas the majority of the transcriptome shows relatively stable and comparable expression. Knowledge of the M. fascicularis draft genome is an important contribution to both the use of this animal in disease models and the safety assessment of drugs and their metabolites. In particular, this information allows high-resolution genotyping and microarray-based gene-expression profiling for animal stratification, thereby allowing the use of well-characterized animals for safety testing. Finally, the genome sequence presented here is a significant contribution to the global "3R" animal welfare initiative, which has the goal to reduce, refine, and replace animal experiments.

Paired opposing leukocyte receptors recognizing rapidly evolving ligands are subject to homogenization of their ligand binding domains

Some leukocyte receptors come in groups of two or more where the partners share ligand(s) but transmit opposite signals. Some of the ligands, such as MHC class I, are fast evolving, raising the problem of how paired opposing receptors manage to change in step with respect to ligand binding properties and at the same time conserve opposite signaling functions. An example is the KLRC (NKG2) family, where opposing variants have been conserved in both rodents and primates. Phylogenetic analyses of the KLRC receptors within and between the two orders show that the opposing partners have been subject to post-speciation gene homogenization restricted mainly to the parts of the genes that encode the ligand binding domains. Concerted evolution similarly restricted is demonstrated also for the KLRI, KLRB (NKR-P1), KLRA (Ly49), and PIR receptor families. We propose the term merohomogenization for this phenomenon and discuss its significance for the evolution of immune receptors.

Naming 'junk': human non-protein coding RNA (ncRNA) gene nomenclature

Previously, the majority of the human genome was thought to be 'junk' DNA with no functional purpose. Over the past decade, the field of RNA research has rapidly expanded, with a concomitant increase in the number of non-protein coding RNA (ncRNA) genes identified in this 'junk'. Many of the encoded ncRNAs have already been shown to be essential for a variety of vital functions, and this wealth of annotated human ncRNAs requires standardised naming in order to aid effective communication. The HUGO Gene Nomenclature Committee (HGNC) is the only organisation authorised to assign standardised nomenclature to human genes. Of the 30,000 approved gene symbols currently listed in the HGNC database (http://www.genenames.org/search), the majority represent protein-coding genes; however, they also include pseudogenes, phenotypic loci and some genomic features. In recent years the list has also increased to include almost 3,000 named human ncRNA genes. HGNC is actively engaging with the RNA research community in order to provide unique symbols and names for each sequence that encodes an ncRNA. Most of the classical small ncRNA genes have now been provided with a unique nomenclature, and work on naming the long (>200 nucleotides) non-coding RNAs (lncRNAs) is ongoing.

Transcript catalogs of human chromosome 21 and orthologous chimpanzee and mouse regions

A comprehensive representation of the gene content of the long arm of human chromosome 21 (Hsa21q) remains of interest for the study of Down syndrome, its associated phenotypic features, and mouse models. Here we compare transcript catalogs for Hsa21q, chimpanzee chromosome 21 (Ptr21q), and orthologous regions of mouse chromosomes 16, 17, and 10 for open reading frame (ORF) characteristics and conservation. The Hsa21q and mouse catalogs contain 552 and 444 gene models, respectively, of which only 162 are highly conserved. Hsa21q transcripts were used to identify orthologous exons in Ptr21q and assemble 533 putative transcripts. Transcript catalogs for all three organisms are searchable for nucleotide and amino acid sequence features of ORF length, repeat content, experimental support, gene structure, and conservation. For human and mouse comparisons, three additional summaries are provided: (1) the chromosomal distribution of novel ORF transcripts versus potential functional RNAs, (2) the distribution of species-specific transcripts within Hsa21q and mouse models of Down syndrome, and (3) the organization of sense-antisense and putative sense-antisense structures defining potential regulatory mechanisms. Catalogs, summaries, and nucleotide and amino acid sequences of all composite transcripts are available and searchable at http://gfuncpathdb.ucdenver.edu/iddrc/chr21/home.php. These data sets provide comprehensive information useful for evaluation of candidate genes and mouse models of Down syndrome and for identification of potential functional RNA genes and novel regulatory mechanisms involving Hsa21q genes. These catalogs and search tools complement and extend information available from other gene annotation projects.

Mobilizing diversity: transposable element insertions in genetic variation and disease

Transposable elements (TEs) comprise a large fraction of mammalian genomes. A number of these elements are actively jumping in our genomes today. As a consequence, these insertions provide a source of genetic variation and, in rare cases, these events cause mutations that lead to disease. Yet, the extent to which these elements impact their host genomes is not completely understood. This review will summarize our current understanding of the mechanisms underlying transposon regulation and the contribution of TE insertions to genetic diversity in the germline and in somatic cells. Finally, traditional methods and emerging technologies for identifying transposon insertions will be considered.

Human-chimpanzee promoter comparisons: property-conserved evolution?

Identification of different functional elements and their properties is a fundamental need in biomedical research and phylogenetic comparisons of a growing number of sequenced genomes form a solid basis for this task. Most available phylogenetic approaches are focused on searching for individual sequence alterations, responsible for the observed phenotype, or statistically evaluate observed mutations to infer general trends. However, being applied to close genomes such methods suffer from poor statistics of rare mutations and give only (at its best) coarse results concerning the potential functional importance of the nucleotide differences. However, quantifying the changes in physical properties of DNA allows to see the strength of introduced mutations and hence to classify them for further investigations. In this work we present the comparative sequence analysis of two evolutionarily close species-human and chimpanzee. In contrast to previous studies we evaluate changes in melting enthalpy of DNA rather than count nucleotide mismatches. We find that nucleotide mismatches in promoters were apparently introduced in a correlated manner during the course of evolution, so that, for example, the DNA property "melting enthalpy" was retained. Such property conservation of promoters is significantly different from nucleotide conservation, shows significant positional and functional biases, and seems to represent a novel feature of gene regulation.

Variation in transcription factor binding among humans

Differences in gene expression may play a major role in speciation and phenotypic diversity. We examined genome-wide differences in transcription factor (TF) binding in several humans and a single chimpanzee by using chromatin immunoprecipitation followed by sequencing. The binding sites of RNA polymerase II (PolII) and a key regulator of immune responses, nuclear factor kappaB (p65), were mapped in 10 lymphoblastoid cell lines, and 25 and 7.5% of the respective binding regions were found to differ between individuals. Binding differences were frequently associated with single-nucleotide polymorphisms and genomic structural variants, and these differences were often correlated with differences in gene expression, suggesting functional consequences of binding variation. Furthermore, comparing PolII binding between humans and chimpanzee suggests extensive divergence in TF binding. Our results indicate that many differences in individuals and species occur at the level of TF binding, and they provide insight into the genetic events responsible for these differences.

Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content

The human Y chromosome began to evolve from an autosome hundreds of millions of years ago, acquiring a sex-determining function and undergoing a series of inversions that suppressed crossing over with the X chromosome^1,2. Little is known about the Y chromosome’s recent evolution because only the human Y chromosome has been fully sequenced. Prevailing theories hold that Y chromosomes evolve by gene loss, the pace of which slows over time, eventually leading to a paucity of genes, and stasis^3,4. These theories have been buttressed by partial sequence data from newly emergent plant and animal Y chromosomes^5-8, but they have not been tested in older, highly evolved Y chromosomes like that of humans. We therefore finished sequencing the male-specific region of the Y chromosome (MSY) in our closest living relative, the chimpanzee, achieving levels of accuracy and completion previously reached for the human MSY. We then compared the MSYs of the two species and found that they differ radically in sequence structure and gene content, implying rapid evolution during the past 6 million years. The chimpanzee MSY harbors twice as many massive palindromes as the human MSY, yet it has lost large fractions of the MSY protein-coding genes and gene families present in the last common ancestor. We suggest that the extraordinary divergence of the chimpanzee and human MSYs was driven by four synergistic factors: the MSY’s prominent role in sperm production, genetic hitchhiking effects in the absence of meiotic crossing over, frequent ectopic recombination within the MSY, and species differences in mating behavior. While genetic decay may be the principal dynamic in the evolution of newly emergent Y chromosomes, wholesale renovation is the paramount theme in the ongoing evolution of chimpanzee, human, and perhaps other older MSYs.

Comparative analysis of Alu repeats in primate genomes

Using bacteria artificial chromosome (BAC) end sequences (16.9 Mb) and high-quality alignments of genomic sequences (17.4 Mb), we performed a global assessment of the divergence distributions, phylogenies, and consensus sequences for Alu elements in primates including lemur, marmoset, macaque, baboon, and chimpanzee as compared to human. We found that in lemurs, Alu elements show a broader and more symmetric sequence divergence distribution, suggesting a steady rate of Alu retrotransposition activity among prosimians. In contrast, Alu elements in anthropoids show a skewed distribution shifted toward more ancient elements with continual declining rates in recent Alu activity along the hominoid lineage of evolution. Using an integrated approach combining mutation profile and insertion/deletion analyses, we identified nine novel lineage-specific Alu subfamilies in lemur (seven), marmoset (one), and baboon/macaque (one) containing multiple diagnostic mutations distinct from their human counterparts-Alu J, S, and Y subfamilies, respectively. Among these primates, we show that that the lemur has the lowest density of Alu repeats (55 repeats/Mb), while marmoset has the greatest abundance (188 repeats/Mb). We estimate that approximately 70% of lemur and 16% of marmoset Alu elements belong to lineage-specific subfamilies. Our analysis has provided an evolutionary framework for further classification and refinement of the Alu repeat phylogeny. The differences in the distribution and rates of Alu activity have played an important role in subtly reshaping the structure of primate genomes. The functional consequences of these changes among the diverse primate lineages over such short periods of evolutionary time are an important area of future investigation.

Gain of new exons and promoters by lineage-specific transposable elements-integration and conservation event on CHRM3 gene

The CHRM3 gene is a member of the muscarinic acetylcholine receptor family that plays important roles in the regulation of fundamental physiological functions. The evolutionary mechanism of exon-acquisition and alternative splicing of the CHRM3 gene in relation to transposable elements (TEs) were analyzed using experimental approaches and in silico analysis. Five different transcript variants (T1, T2, T3, T3-1, and T4) derived from three distinct promoter regions (T1: L1HS, T2, T4: original, T3, T3-1: THE1C) were identified. A placenta (T1) and testis (T3 and T3-1)-dominated expression pattern appeared to be controlled by different TEs (L1HS and THE1C) that were integrated into the common ancestor genome during primate evolution. Remarkably, the T1 transcript was formed by the integration event of the human specific L1HS element. Among the 12 different brain regions, the brain stem, olfactory region, and cerebellum showed decreased expression patterns. Evolutionary analysis of splicing sites and alternative splicing suggested that the exon-acquisition event was determined by a selection and conservation mechanism. Furthermore, continuous integration events of transposable elements could produce lineage specific alternative transcripts by providing novel promoters and splicing sites. Taken together, exon-acquisition and alternative splicing events of CHRM3 genes were shown to have occurred through the continuous integration of transposable elements following conservation.

Cytogenetic anchoring of radiation hybrid and virtual maps of sheep chromosome X and comparison of X chromosomes in sheep, cattle, and human

A comprehensive physical map was generated for Ovis aries chromosome X (OARX) based on a cytogenomics approach. DNA probes were prepared from bacterial artificial chromosome (BAC) clones from the CHORI-243 sheep library and were assigned to G-banded metaphase spreads via fluorescence in-situ hybridization (FISH). A total of 22 BACs gave a single hybridization signal to the X chromosome and were assigned out of 32 tested. The positioned BACs contained 16 genes and a microsatellite marker which represent new cytogenetically mapped loci in the sheep genome. The gene and microsatellite loci serve to anchor between the existing radiation hybrid (RH) and virtual sheep genome (VSG) maps to the cytogenetic OARX map, whilst the BACs themselves also serve as anchors between the VSG and the cytogenetic maps. An additional 17 links between the RH and cytogenetic maps are provided by BAC end sequence (BES) derived markers that have also been positioned on the RH map. Comparison of the map orders for the cytogenetic, RH, and virtual maps reveals that the orders for the cytogenetic and RH maps are most similar, with only one locus, represented by BAC CH243-330E18, mapping to relatively different positions. Several discrepancies, including an inverted segment are found when comparing both the cytogenetic and RH maps with the virtual map. These discrepancies highlight the value of using physical mapping methods to inform the process of future in silico map construction. A detailed comparative analysis of sheep, human, and cattle mapping data allowed the construction of a comparative map that confirms and expands the knowledge about evolutionary conservation and break points between the X chromosomes of the three mammalian species.

Analysis of newly identified low copy AluYj subfamily

Human specific AluY elements were investigated by comparative analysis between human chromosome 21 and chimpanzee chromosome 22. Human specific AluY element was identified on human chromosome 21q22 (accession no. AL163282), and then that was a new member of AluYj subfamily. From the bioinformatic analysis, AluYj subfamily was investigated in human whole genome using AluYj4 consensus sequence (accession no. AL163282). Thirteen members of the AluYj4 elements (4 diagnostic mutations) and eight members of the AluYj3 elements (3 diagnostic mutations) were identified with distinct diagnostic mutation from AluY consensus sequence. The results of the molecular clock calculation of non-CpG region substitution indicated that, AluYj4 elements (2.1 million years old) may be proliferated more recent time than AluYj3 elements (14.1 million years old). For the verification of recent insertion time, four of AluYj4 elements (ch2-AC017101, ch10-AC044786, ch12-AC007656 and ch21-AL163282) from human chromosomes 2, 10, 12, 21 were analyzed by PCR amplification using various human and primate DNA samples. Though, no polymorphism was detected in human population, we identified the new AluYj4 subfamily as the human specific elements.

Genome-wide analysis of chimpanzee genes with premature termination codons

Background

Premature termination codons (PTCs) cause mRNA degradation or a truncated protein and thereby contribute to the transcriptome and proteome divergence between species. Here we present the first genome-wide study of PTCs in the chimpanzee. By comparing the human and chimpanzee genome sequences we identify and characterize genes with PTCs, in order to understand the contribution of these mutations to the transcriptome diversity between the species.

Results

We have studied a total of 13,487 human-chimpanzee gene pairs and found that ~8% were affected by PTCs in the chimpanzee. A majority (764/1,109) of PTCs were caused by insertions or deletions and the remaining part was caused by substitutions. The distribution of PTC genes varied between chromosomes, with Y having the highest proportion. Furthermore, the density of PTC genes varied on a megabasepair scale within chromosomes and we found the density to be correlated both with indel divergence and proximity to the telomere. Within genes, PTCs were more common close to the 5' and 3' ends of the amino acid sequence. Gene Ontology classification revealed that olfactory receptor genes were over represented among the PTC genes.

Conclusion

Our results showed that the density of PTC genes fluctuated across the genome depending on the local genomic context. PTCs were preferentially located in the terminal parts of the transcript, which generally have a lower frequency of functional domains, indicating that selection was operating against PTCs at sites central to protein function. The enrichment of GO terms associated with olfaction suggests that PTCs may have influenced the difference in the repertoire of olfactory genes between humans and chimpanzees. In summary, 8% of the chimpanzee genes were affected by PTCs and this type of variation is likely to have an important effect on the transcript and proteomic divergence between humans and chimpanzees.

Genome and gene alterations by insertions and deletions in the evolution of human and chimpanzee chromosome 22

BACKGROUND

Understanding structure and function of human genome requires knowledge of genomes of our closest living relatives, the primates. Nucleotide insertions and deletions (indels) play a significant role in differentiation that underlies phenotypic differences between humans and chimpanzees. In this study, we evaluated distribution, evolutionary history, and function of indels found by comparing syntenic regions of the human and chimpanzee genomes.

RESULTS

Specifically, we identified 6,279 indels of 10 bp or greater in a ~33 Mb alignment between human and chimpanzee chromosome 22. After the exclusion of those in repetitive DNA, 1,429 or 23% of indels still remained. This group was characterized according to the local or genome-wide repetitive nature, size, location relative to genes, and other genomic features. We defined three major classes of these indels, using local structure analysis: (i) those indels found uniquely without additional copies of indel sequence in the surrounding (10 Kb) region, (ii) those with at least one exact copy found nearby, and (iii) those with similar but not identical copies found locally. Among these classes, we encountered a high number of exactly repeated indel sequences, most likely due to recent duplications. Many of these indels (683 of 1,429) were in proximity of known human genes. Coding sequences and splice sites contained significantly fewer of these indels than expected from random expectations, suggesting that selection is a factor in limiting their persistence. A subset of indels from coding regions was experimentally validated and their impacts were predicted based on direct sequencing in several human populations as well as chimpanzees, bonobos, gorillas, and two subspecies of orangutans.

CONCLUSION

Our analysis demonstrates that while indels are distributed essentially randomly in intergenic and intronic genomic regions, they are significantly under-represented in coding sequences. There are substantial differences in representation of indel classes among genomic elements, most likely caused by differences in their evolutionary histories. Using local sequence context, we predicted origins and phylogenetic relationships of gene-impacting indels in primate species. These results suggest that genome plasticity is a major force behind speciation events separating the great ape lineages.

Cooperative exonization of MaLR and AluJo elements contributed an alternative promoter and novel splice variants of RNF19

The RNF19 protein, which contains RING-finger and IBR motifs, acts as an E3 ubiquitin ligase localized to Lewy bodies. RNF19 is located on human chromosome 8q22.2, has a 4.4-kb transcript, and is ubiquitously expressed in various tissues. Here, we identified an alternative RNF19 promoter region and alternative RNF19 transcripts derived from MaLR (mammalian apparent LTR-retrotransposon) and AluJo elements. Comparative analyses indicated human-specific expression of the MaLR- and AluJo-related transcripts. From the expression analysis of 72 tissue samples including human normal, tumor, and primate tissues, three different isoforms (V1, V2, and V3) of MaLR-derived transcripts were identified. Quantitative RT-PCR analysis showed a dominant expression pattern of the V2 MaLR-derived transcript. A reporter gene assay for MaLR element promoter activity indicated that pGL2-RNF19/MaLR in the forward orientation is capable of driving luciferase gene expression in Cos7 and HCT116 cells. These findings suggest that RNF19 has acquired a new promoter and alternative exons via continuous retrotransposition events of MaLR and AluJo elements during mammalian and primate evolution, respectively.

High divergence in primate-specific duplicated regions: human and chimpanzee chorionic gonadotropin beta genes

Background

Low nucleotide divergence between human and chimpanzee does not sufficiently explain the species-specific morphological, physiological and behavioral traits. As gene duplication is a major prerequisite for the emergence of new genes and novel biological processes, comparative studies of human and chimpanzee duplicated genes may assist in understanding the mechanisms behind primate evolution. We addressed the divergence between human and chimpanzee duplicated genomic regions by using Luteinizing Hormone Beta (LHB)/Chorionic Gonadotropin Beta (CGB) gene cluster as a model. The placental CGB genes that are essential for implantation have evolved from an ancestral pituitary LHB gene by duplications in the primate lineage.

Results

We shotgun sequenced and compared the human (45,165 bp) and chimpanzee (39,876 bp) LHB/CGB regions and hereby present evidence for structural variation resulting in discordant number of CGB genes (6 in human, 5 in chimpanzee). The scenario of species-specific parallel duplications was supported (i) as the most parsimonious solution requiring the least rearrangement events to explain the interspecies structural differences; (ii) by the phylogenetic trees constructed with fragments of intergenic regions; (iii) by the sequence similarity calculations. Across the orthologous regions of LHB/CGB cluster, substitutions and indels contributed approximately equally to the interspecies divergence and the distribution of nucleotide identity was correlated with the regional repeat content. Intraspecies gene conversion may have shaped the LHB/CGB gene cluster. The substitution divergence (1.8–2.59%) exceeded two-three fold the estimates for single-copy loci and the fraction of transversional mutations was increased compared to the unique sequences (43% versus ~30%). Despite the high sequence identity among LHB/CGB genes, there are signs of functional differentiation among the gene copies. Estimates for d_n/d_s rate ratio suggested a purifying selection on LHB and CGB8, and a positive evolution of CGB1.

Conclusion

If generalized, our data suggests that in addition to species-specific deletions and duplications, parallel duplication events may have contributed to genetic differences separating humans from their closest relatives. Compared to unique genomic segments, duplicated regions are characterized by high divergence promoted by intraspecies gene conversion and species-specific chromosomal rearrangements, including the alterations in gene copy number.