Accelerated evolution in the human lineage led to gain and loss of transcriptional enhancers in the RBFOX1 locus.

A long-standing goal of evolutionary biology is to decode how changes in gene regulatory networks contribute to human-specific traits. Human accelerated regions (HARs) are prime candidates for driving gene regulatory modifications in human development. The RBFOX1 locus is densely populated with HARs, providing a set of potential regulatory elements that could have changed its expression in the human lineage. Here, we examined the role of RBFOX1-HARs using transgenic zebrafish reporter assays and identified 15 transcriptional enhancers that are active in the developing nervous system, 9 of which displayed differential activity between the human and chimpanzee sequences. The engineered loss of two selected RBFOX1-HARs in knockout mouse models modified Rbfox1 expression at specific developmental stages and tissues in the brain, influencing the expression and splicing of a high number of Rbfox1 target genes. Our results provided insight into the spatial and temporal changes in gene expression driven by RBFOX1-HARs.

Novel Evolution of Mineralocorticoid Receptor in Humans Compared to Chimpanzees, Gorillas, and Orangutans.

We identified five distinct full-length human mineralocorticoid receptor (MR) genes containing either 984 amino acids (MR-984) or 988 amino acids (MR-988), which can be distinguished by the presence or absence of Lys, Cys, Ser, and Trp (KCSW) in their DNA-binding domain (DBD) and mutations at codons 180 and 241 in their amino-terminal domain (NTD). Two human MR-KCSW genes contain either (Val-180, Val-241) or (Ile-180, Val-241) in their NTD, and three human MR-984 genes contain either (Ile-180, Ala-241), (Val-180, Val-241), or (Ile-180, Val-241). Human MR-KCSW with (Ile-180, Ala-241) has not been cloned. In contrast, chimpanzees contain four MRs: two MR-988s with KCSW in their DBD, or two MR-984s without KCSW in their DBD. Chimpanzee MRs only contain (Ile180, Val-241) in their NTD. A chimpanzee MR with either (Val-180, Val-241) or (Ile-180, Ala-241) in the NTD has not been cloned. Gorillas and orangutans each contain one MR-988 with KCSW in the DBD and one MR-984 without KCSW, and these MRs only contain (Ile-180, Val-241) in their NTD. A gorilla MR or orangutan MR with either (Val-180, Val-241) or (Ile-180, Ala-241) in the NTD has not been cloned. Together, these data suggest that human MRs with (Val-180, Val-241) or (Ile-180, Ala-241) in the NTD evolved after humans and chimpanzees diverged from their common ancestor. Considering the multiple functions in human development of the MR in kidney, brain, heart, skin, and lungs, as well as MR activity in interaction with the glucocorticoid receptor, we suggest that the evolution of human MRs that are absent in chimpanzees may have been important in the evolution of humans from chimpanzees. Investigation of the physiological responses to corticosteroids mediated by the MR in humans, chimpanzees, gorillas, and orangutans may provide insights into the evolution of humans and their closest relatives.

Long non-coding RNAs expression and regulation across different brain regions in primates.

Human and non-human primates have strikingly similar genomes, but they strongly differ in many brain-based processes (e.g., behaviour and cognition). While the functions of protein-coding genes have been extensively studied, rather little is known about the role of non-coding RNAs such as long non-coding RNAs (lncRNAs). Here, we predicted lncRNAs and analysed their expression pattern across different brain regions of human and non-human primates (chimpanzee, gorilla, and gibbon). Our analysis identified shared orthologous and non-orthologous lncRNAs, showing striking differences in the genomic features. Differential expression analysis of the shared orthologous lncRNAs from humans and chimpanzees revealed distinct expression patterns in subcortical regions (striatum, hippocampus) and neocortical areas while retaining a homogeneous expression in the cerebellum. Co-expression analysis of lncRNAs and protein-coding genes revealed massive proportions of co-expressed pairs in neocortical regions of humans compared to chimpanzees. Network analysis of co-expressed pairs revealed the distinctive role of the hub-acting orthologous lncRNAs in a region- and species-specific manner. Overall, our study provides novel insight into lncRNA driven gene regulatory landscape, neural regulation, brain evolution, and constitutes a resource for primate's brain lncRNAs.

The modern human aryl hydrocarbon receptor is more active when ancestralized by genome editing.

The aryl hydrocarbon receptor (AHR) is a transcription factor that has many functions in mammals. Its best known function is that it binds aromatic hydrocarbons and induces the expression of cytochrome P450 genes, which encode enzymes that metabolize aromatic hydrocarbons and other substrates. All present-day humans carry an amino acid substitution at position 381 in the AHR that occurred after the divergence of modern humans from Neandertals and Denisovans. Previous studies that have expressed the ancestral and modern versions of AHR from expression vectors have yielded conflicting results with regard to their activities. Here, we use genome editing to modify the endogenous AHR gene so that it encodes to the ancestral, Neandertal-like AHR protein in human cells. In the absence of exogenous ligands, the expression of AHR target genes is higher in cells expressing the ancestral AHR than in cells expressing the modern AHR, and similar to the expression in chimpanzee cells. Furthermore, the modern human AHR needs higher doses of three ligands than the ancestral AHR to induce the expression of target genes. Thus, the ability of AHR to induce the expression of many of its target genes is reduced in modern humans.

The complete sequence and comparative analysis of ape sex chromosomes.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

The variation and evolution of complete human centromeres.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Comprehensive characterization of ERV-K (HML-8) in the chimpanzee genome revealed less genomic activity than humans.

Endogenous retroviruses (ERVs) originate from ancestral germline infections caused by exogenous retroviruses. Throughout evolution, they have become fixed within the genome of the animals into which they were integrated. As ERV elements coevolve with the host, they are normally epigenetically silenced and can become upregulated in a series of physiological and pathological processes. Generally, a detailed ERV profile in the host genome is critical for understanding the evolutionary history and functional performance of the host genome. We previously characterized and cataloged all the ERV-K subtype HML-8 loci in the human genome; however, this has not been done for the chimpanzee, the nearest living relative of humans. In this study, we aimed to catalog and characterize the integration of HML-8 in the chimpanzee genome and compare it with the integration of HML-8 in the human genome. We analyzed the integration of HML-8 and found that HML-8 pervasively invaded the chimpanzee genome. A total of 76 proviral elements were characterized on 23/24 chromosomes, including detailed elements distribution, structure, phylogeny, integration time, and their potential to regulate adjacent genes. The incomplete structure of HML-8 proviral LTRs will undoubtedly affect their activity. Moreover, the results indicated that HML-8 integration occurred before the divergence between humans and chimpanzees. Furthermore, chimpanzees include more HML-8 proviral elements (76 vs. 40) and fewer solo long terminal repeats (LTR) (0 vs. 5) than humans. These results suggested that chimpanzee genome activity is less than the human genome and that humans may have a better ability to shape and screen integrated proviral elements. Our work is informative in both an evolutionary and a functional context for ERVs.

TRAILS: Tree reconstruction of ancestry using incomplete lineage sorting.

Genome-wide genealogies of multiple species carry detailed information about demographic and selection processes on individual branches of the phylogeny. Here, we introduce TRAILS, a hidden Markov model that accurately infers time-resolved population genetics parameters, such as ancestral effective population sizes and speciation times, for ancestral branches using a multi-species alignment of three species and an outgroup. TRAILS leverages the information contained in incomplete lineage sorting fragments by modelling genealogies along the genome as rooted three-leaved trees, each with a topology and two coalescent events happening in discretized time intervals within the phylogeny. Posterior decoding of the hidden Markov model can be used to infer the ancestral recombination graph for the alignment and details on demographic changes within a branch. Since TRAILS performs posterior decoding at the base-pair level, genome-wide scans based on the posterior probabilities can be devised to detect deviations from neutrality. Using TRAILS on a human-chimp-gorilla-orangutan alignment, we recover speciation parameters and extract information about the topology and coalescent times at high resolution.

Tempo and mode of gene expression evolution in the brain across primates.

Primate evolution has led to a remarkable diversity of behavioral specializations and pronounced brain size variation among species (Barton, 2012; DeCasien and Higham, 2019; Powell et al., 2017). Gene expression provides a promising opportunity for studying the molecular basis of brain evolution, but it has been explored in very few primate species to date (e.g. Khaitovich et al., 2005; Khrameeva et al., 2020; Ma et al., 2022; Somel et al., 2009). To understand the landscape of gene expression evolution across the primate lineage, we generated and analyzed RNA-seq data from four brain regions in an unprecedented eighteen species. Here, we show a remarkable level of variation in gene expression among hominid species, including humans and chimpanzees, despite their relatively recent divergence time from other primates. We found that individual genes display a wide range of expression dynamics across evolutionary time reflective of the diverse selection pressures acting on genes within primate brain tissue. Using our samples that represent a 190-fold difference in primate brain size, we identified genes with variation in expression most correlated with brain size. Our study extensively broadens the phylogenetic context of what is known about the molecular evolution of the brain across primates and identifies novel candidate genes for the study of genetic regulation of brain evolution.

Dating ancient splits in phylogenetic trees, with application to the human-Neanderthal split.

BACKGROUND: We tackle the problem of estimating species TMRCAs (Time to Most Recent Common Ancestor), given a genome sequence from each species and a large known phylogenetic tree with a known structure (typically from one of the species). The number of transitions at each site from the first sequence to the other is assumed to be Poisson distributed, and only the parity of the number of transitions is observed. The detailed phylogenetic tree contains information about the transition rates in each site. We use this formulation to develop and analyze multiple estimators of the species' TMRCA. To test our methods, we use mtDNA substitution statistics from the well-established Phylotree as a baseline for data simulation such that the substitution rate per site mimics the real-world observed rates. RESULTS: We evaluate our methods using simulated data and compare them to the Bayesian optimizing software BEAST2, showing that our proposed estimators are accurate for a wide range of TMRCAs and significantly outperform BEAST2. We then apply the proposed estimators on Neanderthal, Denisovan, and Chimpanzee mtDNA genomes to better estimate their TMRCA with modern humans and find that their TMRCA is substantially later, compared to values cited recently in the literature. CONCLUSIONS: Our methods utilize the transition statistics from the entire known human mtDNA phylogenetic tree (Phylotree), eliminating the requirement to reconstruct a tree encompassing the specific sequences of interest. Moreover, they demonstrate notable improvement in both running speed and accuracy compared to BEAST2, particularly for earlier TMRCAs like the human-Chimpanzee split. Our results date the human - Neanderthal TMRCA to be [Formula: see text] years ago, considerably later than values cited in other recent studies.

Astrocytes Drive Divergent Metabolic Gene Expression in Humans and Chimpanzees.

The human brain utilizes ∼20% of all of the body's metabolic resources, while chimpanzee brains use <10%. Although previous work shows significant differences in metabolic gene expression between the brains of primates, we have yet to fully resolve the contribution of distinct brain cell types. To investigate cell type-specific interspecies differences in brain gene expression, we conducted RNA-seq on neural progenitor cells, neurons, and astrocytes generated from induced pluripotent stem cells from humans and chimpanzees. Interspecies differential expression analyses revealed that twice as many genes exhibit differential expression in astrocytes (12.2% of all genes expressed) than neurons (5.8%). Pathway enrichment analyses determined that astrocytes, rather than neurons, diverged in expression of glucose and lactate transmembrane transport, as well as pyruvate processing and oxidative phosphorylation. These findings suggest that astrocytes may have contributed significantly to the evolution of greater brain glucose metabolism with proximity to humans.

The KIR repertoire of a West African chimpanzee population is characterized by limited gene, allele, and haplotype variation.

Introduction: The killer cell immunoglobulin-like receptors (KIR) play a pivotal role in modulating the NK cell responses, for instance, through interaction with major histocompatibility complex (MHC) class I molecules. Both gene systems map to different chromosomes but co-evolved during evolution. The human KIR gene family is characterized by abundant allelic polymorphism and copy number variation. In contrast, our knowledge of the KIR repertoire in chimpanzees is limited to 39 reported alleles, with no available population data. Only three genomic KIR region configurations have been mapped, and seventeen additional ones were deduced by genotyping. Methods: Previously, we documented that the chimpanzee MHC class I repertoire has been skewed due to an ancient selective sweep. To understand the depth of the sweep, we set out to determine the full-length KIR transcriptome - in our MHC characterized pedigreed West African chimpanzee cohort - using SMRT sequencing (PacBio). In addition, the genomic organization of 14 KIR haplotypes was characterized by applying a Cas9-mediated enrichment approach in concert with long-read sequencing by Oxford Nanopore Technologies. Results: In the cohort, we discovered 35 undescribed and 15 already recorded Patr-KIR alleles, and a novel hybrid KIR gene. Some KIR transcripts are subject to evolutionary conserved alternative splicing events. A detailed insight on the KIR region dynamics (location and order of genes) was obtained, however, only five new KIR region configurations were detected. The population data allowed to investigate the distribution of the MHC-C1 and C2-epitope specificity of the inhibitory lineage III KIR repertoire, and appears to be skewed towards C2. Discussion: Although the KIR region is known to evolve fast, as observed in other primate species, our overall conclusion is that the genomic architecture and repertoire in West African chimpanzees exhibit only limited to moderate levels of variation. Hence, the ancient selective sweep that affected the chimpanzee MHC class I region may also have impacted the KIR system.

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human.

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

Disease association and inter-connectivity analysis of human brain specific co-expressed functional modules

BACKGROUND: In the recent studies, it is suggested that the analysis of transcriptomic change of functional modules instead of individual genes would be more effective for system-wide identification of cellular functions. This could also provide a new possibility for the better understanding of difference between human and chimpanzee. RESULTS: In this study, we analyzed to find molecular characteristics of human brain functions from the difference of transcriptome between human and chimpanzee's brain using the functional module-centric co-expression analysis. We performed analysis of brain disease association and systems-level connectivity of species-specific co-expressed functional modules. CONCLUSIONS: Throughout the analyses, we found human-specific functional modules and significant overlap between their genes in known brain disease genes, suggesting that human brain disorder could be mediated by the perturbation of modular activities emerged in human brain specialization. In addition, the human-specific modules having neurobiological functions exhibited higher networking than other functional modules. This finding suggests that the expression of neural functions are more connected than other functions, and the resulting high-order brain functions could be identified as a result of consolidated inter-modular gene activities. Our result also showed that the functional module based transcriptome analysis has a potential to expand molecular understanding of high-order complex functions like cognitive abilities and brain disorders.

Molecular evolution and expression profile of the chemerine encoding gene RARRES2 in baboon and chimpanzee

BACKGROUND: Chemerin, encoded by the retinoic acid receptor responder 2 (RARRES2) gene is an adipocytesecreted protein with autocrine/paracrine functions in adipose tissue, metabolism and inflammation with a recently described function in vascular tone regulation, liver, steatosis, etc. This molecule is believed to represent a critical endocrine signal linking obesity to diabetes. There are no data available regarding evolution of RARRES2 in non-human primates and great apes. Expression profile and orthology in RARRES2 genes are unknown aspects in the biology of this multigene family in primates. Thus; we attempt to describe expression profile and phylogenetic relationship as complementary knowledge in the function of this gene in primates. To do that, we performed A RT-PCR from different tissues obtained during necropsies. Also we tested the hypotheses of positive evolution, purifying selection, and neutrality. And finally a phylogenetic analysis was made between primates RARRES2 protein. RESULTS: RARRES2 transcripts were present in liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus and pituitary tissues. Expression in kidney and leukocytes were not detectable in either species. It was determined that the studied genes are orthologous. CONCLUSIONS: RARRES2 evolution fits the hypothesis of purifying selection. Expression profiles of the RARRES2 gene are similar in baboons and chimpanzees and are also phylogenetically related.