Chimpanzee Orthologs Research Articles

Resolving the three-dimensional interactome of Human Accelerated Regions during human and chimpanzee neurodevelopment.

Human Accelerated Regions (HARs) are highly conserved across species but exhibit a significant excess of human-specific sequence changes, suggesting they may have gained novel functions in human evolution. HARs include transcriptional enhancers with human-specific activity and have been implicated in the evolution of the human brain. However, our understanding of how HARs contributed to uniquely human features of the brain is hindered by a lack of insight into the genes and pathways that HARs regulate. It is unclear whether HARs acted by altering the expression of gene targets conserved between HARs and their chimpanzee orthologs or by gaining new gene targets in human, a mechanism termed enhancer hijacking. We generated a high-resolution map of chromatin interactions for 1,590 HARs and their orthologs in human and chimpanzee neural stem cells (NSCs) to comprehensively identify gene targets in both species. HARs and their chimpanzee orthologs targeted a conserved set of 2,963 genes enriched for neurodevelopmental processes including neurogenesis and synaptic transmission. Changes in HAR enhancer activity were correlated with changes in conserved gene target expression. Conserved targets were enriched among genes differentially expressed between human and chimpanzee NSCs or between human and non-human primate developing and adult brain. Species-specific HAR gene targets did not converge on known biological functions and were not significantly enriched among differentially expressed genes, suggesting that HARs did not alter gene expression via enhancer hijacking. HAR gene targets, including differentially expressed targets, also showed cell type-specific expression patterns in the developing human brain, including outer radial glia, which are hypothesized to contribute to human cortical expansion. Our findings support that HARs influenced human brain evolution by altering the expression of conserved gene targets and provide the means to functionally link HARs with novel human brain features.

Enhancer Function and Evolutionary Roles of Human Accelerated Regions.

Human accelerated regions (HARs) are the fastest-evolving sequences in the human genome. When HARs were discovered in 2006, their function was mysterious due to scant annotation of the noncoding genome. Diverse technologies, from transgenic animals to machine learning, have consistently shown that HARs function as gene regulatory enhancers with significant enrichment in neurodevelopment. It is now possible to quantitatively measure the enhancer activity of thousands of HARs in parallel and model how each nucleotide contributes to gene expression. These strategies have revealed that many human HAR sequences function differently than their chimpanzee orthologs, though individual nucleotide changes in the same HAR may have opposite effects, consistent with compensatory substitutions. To fully evaluate the role of HARs in human evolution, it will be necessary to experimentally and computationally dissect them across more cell types and developmental stages.

Discovery of Novel Sequences in 1,000 Swedish Genomes

Novel sequences (NSs), not present in the human reference genome, are abundant and remain largely unexplored. Here, we utilize de novo assembly to study NS in 1,000 Swedish individuals first sequenced as part of the SweGen project revealing a total of 46 Mb in 61,044 distinct contigs of sequences not present in GRCh38. The contigs were aligned to recently published catalogs of Icelandic and Pan-African NSs, as well as the chimpanzee genome, revealing a great diversity of shared sequences. Analyzing the positioning of NS across the chimpanzee genome, we find that 2,807 NS align confidently within 143 chimpanzee orthologs of human genes. Aligning the whole genome sequencing data to the chimpanzee genome, we discover ancestral NS common throughout the Swedish population. The NSs were searched for repeats and repeat elements: revealing a majority of repetitive sequence (56%), and enrichment of simple repeats (28%) and satellites (15%). Lastly, we align the unmappable reads of a subset of the thousand genomes data to our collection of NS, as well as the previously published Pan-African NS: revealing that both the Swedish and Pan-African NS are widespread, and that the Swedish NSs are largely a subset of the Pan-African NS. Overall, these results highlight the importance of creating a more diverse reference genome and illustrate that significant amounts of the NS may be of ancestral origin.

Evolution of the exclusively human pathogen Neisseria gonorrhoeae: Human-specific engagement of immunoregulatory Siglecs.

Neisseria gonorrhoeae causes the sexually transmitted disease gonorrhea exclusively in humans and uses multiple strategies to infect, including acquisition of host sialic acids that cap and mask lipooligosaccharide termini, while restricting complement activation. We hypothesized that gonococci selectively target human anti‐inflammatory sialic acid‐recognizing Siglec receptors on innate immune cells to blunt host responses and that pro‐inflammatory Siglecs and SIGLEC pseudogene polymorphisms represent host evolutionary adaptations to counteract this interaction. N. gonorrhoeae can indeed engage multiple human but not chimpanzee CD33rSiglecs expressed on innate immune cells and in the genitourinary tract––including Siglec‐11 (inhibitory) and Siglec‐16 (activating), which we detected for the first time on human cervical epithelium. Surprisingly, in addition to LOS sialic acid, we found that gonococcal porin (PorB) mediated binding to multiple Siglecs. PorB also bound preferentially to human Siglecs and not chimpanzee orthologs, modulating host immune reactions in a human‐specific manner. Lastly, we studied the distribution of null SIGLEC polymorphisms in a Namibian cohort with a high prevalence of gonorrhea and found that uninfected women preferentially harbor functional SIGLEC16 alleles encoding an activating immune receptor. These results contribute to the understanding of the human specificity of N. gonorrhoeae and how it evolved to evade the human immune defense.

Functional test of PCDHB11, the most human-specific neuronal surface protein.

BackgroundBrain-expressed proteins that have undergone functional change during human evolution may contribute to human cognitive capacities, and may also leave us vulnerable to specifically human diseases, such as schizophrenia, autism or Alzheimer’s disease. In order to search systematically for those proteins that have changed the most during human evolution and that might contribute to brain function and pathology, all proteins with orthologs in chimpanzee, orangutan and rhesus macaque and annotated as being expressed on the surface of cells in the human central nervous system were ordered by the number of human-specific amino acid differences that are fixed in modern populations.ResultsPCDHB11, a beta-protocadherin homologous to murine cell adhesion proteins, stood out with 12 substitutions and maintained its lead after normalizing for protein size and applying weights for amino acid exchange probabilities. Human PCDHB11 was found to cause homophilic cell adhesion, but at lower levels than shown for other clustered protocadherins. Homophilic adhesion caused by a PCDHB11 with reversion of human-specific changes was as low as for modern human PCDHB11; while neither human nor reverted PCDHB11 adhered to controls, they did adhere to each other. A loss of function in PCDHB11 is unlikely because intra-human variability did not increase relative to the other human beta-protocadherins.ConclusionsThe brain-expressed protein with the highest number of human-specific substitutions is PCDHB11. In spite of its fast evolution and low intra-human variability, cell-based tests on the only proposed function for PCDHB11 did not indicate a functional change.

No Evidence that MicroRNAs Coevolve with Genes Located in Copy Number Regions.

MicroRNAs (miRNAs) are a widespread class of regulatory noncoding RNAs with key roles in physiology and development, conferring robustness to noise in regulatory networks. Consistent with this buffering function, it was recently suggested that human miRNAs coevolve with genes in copy number regions (copy number variation [CNV] genes) to reduce dosage imbalance. Here, I compare miRNA regulation between CNV and non-CNV genes in four model organisms. miRNA regulation of CNV genes is elevated in human and fly but reduced in nematode and zebrafish. By analyzing 31 human CNV data sets, careful analysis of human and chimpanzee orthologs, resampling genes within species and comparing structural variant types, I show that the apparent coevolution between CNV genes and miRNAs is due to the strong dependency between 3'-untranslated region length and miRNA target prediction. Deciphering the interplay between CNVs and miRNAs will likely require a deeper understanding of how miRNAs are embedded in regulatory circuits.

Allelic Frequency in Human SNPs Predicts the Rate of Non-Synonymous Nucleotide Substitutions between Human and Chimpanzee Genes

The combination of comparative genomics and population genetics may provide important clues regarding human evolution. We have hypothesized that the divergence between various human and chimpanzee orthologs will be reflected in the variability of single nucleotide polymorphisms (SNPs) that are localized in the vicinity of the corresponding loci in different human populations. Consistently with this notion, more diverged genes between humans and chimpanzees are more likely to be associated with human speciation and are anticipated to contain SNPs with reduced variability between different human populations. In order to test this hypothesis, we have compared the rate of non synonymous nucleotide substitutions (Ka) between 255 chimpanzee and human orthologs with the average deviation in the allelic frequencies of corresponding closely linked SNPs in two distinct human populations: The Yoruba people in Ibadan, Nigeria (YRI) and US residents with ancestry from Northern and Western Europe, collected in 1980 by the Centre d’Etude du Polymorphisme Humain (CEU). We found a significant (p

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

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

Functional impact of transposable elements using bioinformatic analysis and a comparative genomic approach.

A dual coding event, which is the translation of different isoforms from a single gene, is one of the special patterns among the alternative splicing events. This is an important mechanism for the regulation of protein diversity in human and mouse genomes. Although the regulation for dual coding events has been characterized in a few genes, the individual mechanism remains unclear. Numerous studies have described the exonization of transposable elements, which is the splicing mediated insertion of transposable element sequence fragments into mature mRNAs. Therefore, in this study, we investigated the number of transposable element (TE)-derived dual coding genes in human, chimpanzee and mouse genomes. TE fusion exons appeared in the dual coding regions of 309 human genes. Functional protein domain alterations by TE-derived dual coding events were observed in 129 human genes. Comparative TE-derived dual coding events were also analyzed in chimpanzee and mouse orthologs. Seventy chimpanzee orthologs had TE-derived dual coding events, but mouse orthologs did not have any TE-derived dual coding events. Taken together, our analyses listed the number of TE-derived dual coding genes which could be investigated by experimental analysis and suggested that TE-derived dual coding events were major sources for the functional diversity of human genes, but not mouse genes.

Transposable Elements: Insertion Pattern and Impact on Gene Expression Evolution in Hominids

Transposable elements (TEs) can affect the regulation of nearby genes through several mechanisms. Here, we examine to what extent recent TE insertions have contributed to the evolution of gene expression in hominids. We compare expression levels of human and chimpanzee orthologs and detect a weak increase in expression divergence (ED) for genes with species-specific TE insertions compared with unaffected genes. However, we show that genes with TE insertions predating the human-chimpanzee split also exhibit a similar increase in ED and therefore conclude that the increase is not due to the transcriptional influence of the TEs. These results are further confirmed by lineage-specific analysis of ED, using rhesus macaque as an outgroup: Human-chimpanzee ortholog pairs, where one ortholog has suffered TE insertion but not the other, do not show increased ED along the lineage where the insertion occurred, relative to the other lineage. We also show that genes with recent TE insertions tend to produce more alternative transcripts but find no evidence that the TEs themselves promote transcript diversity. Finally, we observe that TEs are enriched upstream relative to downstream of genes and show that this is due to insertional bias, rather than selection, because this bias is only observed in genes expressed in the germ line. This provides an alternative neutral explanation for the accumulation of TEs in upstream sequences.

Interrogating 11 Fast-Evolving Genes for Signatures of Recent Positive Selection in Worldwide Human Populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs in chimpanzee, mouse, rat, and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes was then further explored for signatures of recent positive selection using single nucleotide polymorphism (SNP) data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP-CEPH diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintenance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal of recent selection in this region may come from a different, neighboring gene (CD5) rather than from the candidate gene itself (VPS37C). For this set of positively selected genes in the human lineage, we find no indication that these genes maintained their rapid evolutionary pace among human populations. Based on these data, it therefore appears that adaptation for human-specific and for population-specific traits may have involved different genes.

Evolution of C2H2-zinc finger genes revisited

BackgroundA recent study by Tadepally et al. describes the clustering of zinc finger (ZF) genes in the human genome and traces their evolutionary history among several placental mammals with complete or draft genome sequences. One of the main conclusions from the paper is that there is a dramatic rate of gene duplication and gene loss, including the surprising result that 118 human ZF genes are absent in chimpanzee. The authors also present evidence concerning the ancestral order in which the ZF-associated KRAB and SCAN domains were recruited to ZF proteins.ResultsBased on our analysis of two of the largest human ZF gene clusters, we find that nearly all of the human genes have plausible orthologs in chimpanzee. The one exception may be a result of the incomplete sequence coverage in the draft chimpanzee genome. The discrepancy in gene content analysis may result from the authors' dependence on the preliminary NCBI gene prediction set for chimpanzee, which appears to either fail to predict or to mispredict many chimpanzee ZF genes. Similar problems may affect the authors' interpretation of the more divergent dog, mouse, and rat ZF gene complements. In addition, we present evidence that the KRAB domain was recruited to ZF genes before the SCAN domain, rather than the reverse as the authors suggest. This discrepancy appears to result from the fact that the SCAN domain did indeed arise before the KRAB domain but is present only in non-ZF genes until a much later date.ConclusionWhen comparing gene content among species, especially when using draft genome assemblies, dependence on preliminary gene prediction sets can be seriously misleading. In such studies, genic sequences must be identified in a manner that is as independent as possible of prediction sets. In addition, we present evidence that provides a more parsimonious explanation for the large proportion of mammalian KRAB-ZF genes without a SCAN domain.

NF90 binds novel human non‐coding RNAs

Proteins in the nuclear factor 90 (NF90) family of dsRNA binding proteins are abundantly expressed in a broad spectrum of vertebrate tissues and participate in many aspects of RNA metabolism. NF90, also known as DRBP76 and NFAR1, forms heterodimeric complexes with NF45 protein. We detected several small non‐protein‐coding RNAs that associate with NF90 in vivo. They include snaRs (small NF90‐associated RNAs), a novel family of RNAs that are unique to the African Great Apes and display a restricted tissue distribution in humans.The genes for 3 of these RNA species are located within the proximal promoters of protein‐coding genes and are transcribed by RNA polymerase III (Pol III) in the opposite direction to their ‘host’ Pol II genes. snaRs are highly structured RNAs, ∼117 nucleotides in length. The snaR‐G1 and snaR‐G2 genes are in the proximal promoters of chorionic gonadotropin genes CG □ 1 and 2, respectively, in locations shared by chimpanzee orthologs. The gene encoding a novel tRNAimet precursor is 242 bp upstream of the 5′‐end of NF45 mRNA, in a highly conserved position in vertebrates. The binding of NF90 to a pre‐tRNAimet transcript whose gene resides in the NF45 promoter is suggestive of an intricate feedback control mechanism. We speculate that NF90 regulates the transcription of specific protein‐coding genes through interaction with small non‐coding RNAs transcribed by Pol III.Supported by NIH grant AI034552.

Rapid Detection of Positive Selection in Genes and Genomes Through Variation Clusters

Positive selection in genes and genomes can point to the evolutionary basis for differences among species and among races within a species. The detection of positive selection can also help identify functionally important protein regions and thus guide protein engineering. Many existing tests for positive selection are excessively conservative, vulnerable to artifacts caused by demographic population history, or computationally very intensive. I here propose a simple and rapid test that is complementary to existing tests and that overcomes some of these problems. It relies on the null hypothesis that neutrally evolving DNA regions should show a Poisson distribution of nucleotide substitutions. The test detects significant deviations from this expectation in the form of variation clusters, highly localized groups of amino acid changes in a coding region. In applying this test to several thousand human-chimpanzee gene orthologs, I show that such variation clusters are not generally caused by relaxed selection. They occur in well-defined domains of a protein's tertiary structure and show a large excess of amino acid replacement over silent substitutions. I also identify multiple new human-chimpanzee orthologs subject to positive selection, among them genes that are involved in reproductive functions, immune defense, and the nervous system.

Detecting positive darwinian selection in brain-expressed genes during human evolution

To understand the genetic basis that underlies the phenotypic divergence between human and non-human primates, we screened a total of 7176 protein-coding genes expressed in the human brain and compared them with the chimpanzee orthologs to identify genes that show evidence of rapid evolution in the human lineage. Our results showed that the nonsynonymous/synonymous substitution (Ka/Ks) ratio for genes expressed in the brain of human and chimpanzee is 0.3854, suggesting that the brain-expressed genes are under functional constraint. The X-linked human brain-expressed genes evolved more rapidly than autosomal ones. We further dissected the molecular evolutionary patterns of 34 candidate genes by sequencing representative primate species to identify lineage-specific adaptive evolution. Fifteen out of the 34 candidate genes showed evidence of positive Darwinian selection in human and/or chimpanzee lineages. These genes are predicted to play diverse functional roles in embryonic development, spermatogenesis and male fertility, signal transduction, sensory nociception, and neural function. This study together with others demonstrated the usefulness and power of phylogenetic comparison of multiple closely related species in detecting lineage-specific adaptive evolution, and the identification of the positively selected brain-expressed genes may add new knowledge to the understanding of molecular mechanism of human origin.

Large-scale sequencing of the CD33-related Siglec gene cluster in five mammalian species reveals rapid evolution by multiple mechanisms.

Siglecs are a recently discovered family of animal lectins that belong to the Ig superfamily and recognize sialic acids (Sias). CD33-related Siglecs (CD33rSiglecs) are a subgroup with as-yet-unknown functions, characterized by sequence homology, expression on innate immune cells, conserved cytosolic tyrosine-based signaling motifs, and a clustered localization of their genes. To better understand the biology and evolution of CD33rSiglecs, we sequenced and compared the CD33rSiglec gene cluster from multiple mammalian species. Within the sequenced region, the segments containing CD33rSiglec genes showed a lower degree of sequence conservation. In contrast to the adjacent conserved kallikrein-like genes, the CD33rSiglec genes showed extensive species differences, including expansions of gene subsets; gene deletions, including one human-specific loss of a novel functional primate Siglec (Siglec-13); exon shuffling, generating hybrid genes; accelerated accumulation of nonsynonymous substitutions in the Sia-recognition domain; and multiple instances of mutations of an arginine residue essential for Sia recognition in otherwise intact Siglecs. Nonsynonymous differences between human and chimpanzee orthologs showed uneven distribution between the two beta sheets of the Sia-recognition domain, suggesting biased mutation accumulation. These data indicate that CD33rSiglec genes are undergoing rapid evolution via multiple genetic mechanisms, possibly due to an evolutionary "arms race" between hosts and pathogens involving Sia recognition. These studies, which reflect one of the most complete comparative sequence analyses of a rapidly evolving gene cluster, provide a clearer picture of the ortholog status of CD33rSiglecs among primates and rodents and also facilitate rational recommendations regarding their nomenclature.

Selective pressures on the olfactory receptor repertoire since the human-chimpanzee divergence.

The availability of the sequence of the chimpanzee genome provides an opportunity to examine human genes and their chimpanzee orthologs and to analyze selective pressures that have been shaping the olfactory receptor repertoire since the human-chimpanzee divergence. We determined the ratio of nonsynonymous to synonymous changes for each of 186 orthologous pairs and then examined how the distribution of these ratios compares with the distribution expected under neutral drift. Consistent with the diminishing importance of olfaction for these species, we find no evidence for positive selection and we find evidence of weak purifying selection affecting over half of the repertoire.

Testing Multiregionality of Modern Human Origins

In order to examine the possibility of multiple founding populations of anatomically modern Homo sapiens, we collected DNA sequence data from 10 X-chromosomal regions, 5 autosomal regions, and 1 Y-chromosomal region, in addition to mitochondrial DNA. Except for five regions which are genealogically uninformative and two other regions for which chimpanzee orthologs are not available, the ancestral sequence and population for each of the remaining regions were successfully inferred. Of these 10 ancestral sequences, 9 occurred in Africa and only 1 occurred in Asia during the Pleistocene. Computer simulation was carried out to quantify the multiregional hypothesis based solely on the premise that there was more than one founding population in the Pleistocene. Allowing the breeding size to vary among the founding populations, the hypothesis may account for the observed African ancestry in 90% of the genomic regions. However, it is required that the founding population in Africa was much larger than that outside Africa. Likelihood estimates of the breeding sizes in the founding populations were more than 9,000 in Africa and less than 1,000 in outside of Africa, although these estimates can be much less biased at the 1% significance level. If the number of African ancestral sequences further increases as more data accumulate in other genomic regions, the conclusion of a single founding population of modern H. sapiens is inevitable.