Comparative and demographic analysis of orang-utan genomes

junctionCounts: comprehensive alternative splicing analysis and prediction of isoform-level impacts to the coding sequence

Alternative splicing (AS) is emerging as an important regulatory process for complex biological processes. Transcriptomic studies therefore commonly involve the identification and quantification of alternative processing events, but the need for predicting the functional consequences of changes to the relative inclusion of alternative events remains largely unaddressed. Many tools exist for the former task, albeit each constrained to its own event type definitions. Few tools exist for the latter task; each with significant limitations. To address these issues we developed junctionCounts, which captures both simple and complex pairwise AS events and quantifies them with straightforward exon-exon and exon-intron junction reads in RNA-seq data, performing competitively among similar tools in terms of sensitivity, false discovery rate and quantification accuracy. Its partner utility, cdsInsertion, identifies transcript coding sequence (CDS) information via in silico translation from annotated start codons, including the presence of premature termination codons. Finally, findSwitchEvents connects AS events with CDS information to predict the impact of individual events to the isoform-level CDS. We used junctionCounts to characterize splicing dynamics and NMD regulation during neuronal differentiation across four primates, demonstrating junctionCounts' capacity to robustly characterize AS in a variety of organisms and to predict its effect on mRNA isoform fate.

Structural variation in humans and our primate kin in the era of telomere-to-telomere genomes and pangenomics

Structural variants (SVs) account for the majority of base pair differences both within and between primate species. However, our understanding of inter- and intra-species SV has been historically hampered by the quality of draft primate genomes and the absence of genome resources for key taxa. Recently, advances in long-read sequencing and genome assembly have begun to radically reshape our understanding of SVs. Two landmark achievements include the publication of a human telomere-to-telomere (T2T) genome as well as the development of the first human pangenome reference. In this review, we first look back to the major works laying the foundation for these projects. We then examine the ways in which T2T genome assemblies and pangenomes are transforming our understanding of and approach to primate SV. Finally, we discuss what the future of primate SV research may look like in the era of T2T genomes and pangenomics.

Complete sequencing of ape genomes

We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.

Novel Cascade Alpha Satellite HORs in Orangutan Chromosome 13 Assembly: Discovery of the 59mer HOR-The largest Unit in Primates-And the Missing Triplet 45/27/18 HOR in Human T2T-CHM13v2.0 Assembly

From the recent genome assembly NHGRI_mPonAbe1-v2.0_NCBI (GCF_028885655.2) of orangutan chromosome 13, we computed the precise alpha satellite higher-order repeat (HOR) structure using the novel high-precision GRM2023 algorithm with Global Repeat Map (GRM) and Monomer Distance (MD) diagrams. This study rigorously identified alpha satellite HORs in the centromere of orangutan chromosome 13, discovering a novel 59mer HOR-the longest HOR unit identified in any primate to date. Additionally, it revealed the first intertwined sequence of three HORs, 18mer/27mer/45mer HORs, with a common aligned "backbone" across all HOR copies. The major 7mer HOR exhibits a Willard's-type canonical copy, although some segments of the array display significant irregularities. In contrast, the 14mer HOR forms a regular Willard's-type HOR array. Surprisingly, the GRM2023 high-precision analysis of chromosome 13 of human genome assembly T2T-CHM13v2.0 reveals the presence of only a 7mer HOR, despite both the orangutan and human genome assemblies being derived from whole genome shotgun sequences.

Evidence for the smallest fossil Pongo in southern China

The rarity of Pongo fossils with precise absolute dating from the Middle Pleistocene hampers our understanding of the taxonomy and spatiotemporal distribution of Quaternary orangutans in southern China. Here, we report a newly discovered sample of 113 isolated teeth of fossil Pongo from Zhongshan Cave in the Bubing Basin, Guangxi, southern China. We describe the Pongo specimens from Zhongshan Cave and compare them metrically to other samples of fossil Pongo species (i.e., Pongo weidenreichi, Pongo devosi, Pongo duboisi, Pongo palaeosumatrensis, Pongo javensis, and Pongo sp.) and to extant orangutans (i.e., Pongo pygmaeus and Pongo abelii). The Zhongshan Pongo assemblage is dated using U-series and coupled electron spin resonance/U-series methods. Our results reasonably constrain the Zhongshan Pongo assemblage to 184 ± 16 ka, which is consistent with the biostratigraphic evidence. The Zhongshan Pongo teeth are only 6.5% larger on average than those of extant Pongo. The Zhongshan teeth are smaller overall than those of Pongo from all other cave sites in southern China, and they currently represent the smallest fossil orangutans in southern China. Based on their dental size, and the presence of a well-developed lingual pillar and lingual cingulum on the upper and lower incisors, an intermediate frequency of lingual cingulum remnants on the upper molars, and a higher frequency of moderate to heavy wrinkling on the upper and lower molars, we provisionally assign the Zhongshan fossils to P. devosi. Our results confirm earlier claims that P. weidenreichi is replaced by a smaller species in southern China, P. devosi, by the late Middle Pleistocene. The occurrence of P. devosi in Zhongshan Cave further extends its spatial and temporal distribution. The Pongo specimens from Zhongshan provide important new evidence to demonstrate that the dental morphological features of Pongo in southern China changed substantially during the late Middle Pleistocene.

Structurally divergent and recurrently mutated regions of primate genomes

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Relative leg-to-arm skeletal strength proportions in orangutans by species and sex

Among extant great apes, orangutans climb most frequently. However, Bornean orangutans (Pongo pygmaeus) exhibit higher frequencies of terrestrial locomotion than do Sumatran orangutans (Pongo abelii). Variation in long bone cross-sectional geometry is known to reflect differential loading of the limbs. Thus, Bornean orangutans should show greater relative leg-to-arm strength than their Sumatran counterparts. Using skeletal specimens from museum collections, we measured two cross-sectional geometric measures of bone strength: the polar section modulus (Zpol) and the ratio of maximum to minimum area moments of inertia (Imax/Imin), at the midshaft of long bones in Bornean (n = 19) and Sumatran adult orangutans (n = 12) using medical CT and peripheral quantitative CT scans, and compared results to published data of other great apes. Relative leg-to-arm strength was quantified using ratios of femur and tibia over humerus, radius, and ulna, respectively. Differences between orangutan species and between sexes in median ratios were assessed using Wilcoxon rank sum tests. The tibia of Bornean orangutans was stronger relative to the humerus and the ulna than in Sumatran orangutans (p = 0.008 and 0.025, respectively), consistent with behavioral studies that indicate higher frequencies of terrestrial locomotion in the former. In three Zpol ratios, adult female orangutans showed greater leg-to-arm bone strength compared to flanged males, which may relate to females using their legs more during arboreal locomotion than in adult flanged males. A greater amount of habitat discontinuity on Borneo compared to Sumatra has been posited as a possible explanation for observed interspecific differences in locomotor behaviors, but recent camera trap studies has called this into question. Alternatively, greater frequencies of terrestriality in Pongo pygmaeus may be due to the absence of tigers on Borneo. The results of this study are consistent with the latter explanation given that habitat continuity was greater a century ago when our study sample was collected.

Orangutan males make increased use of social learning opportunities, when resource availability is high

Humans' colonization of diverse habitats relied on our ancestors' abilities to innovate and share innovations with others. While ecological impacts on innovations are well studied, their effect on social learning remains poorly understood. We examined how food availability affects social learning in migrant orangutan unflanged males, who may learn from local orangutans through peering (i.e., observational social learning). We analyzed 1,384 dyadic associations, including 360 peering events, among 46 wild Sumatran orangutan and 25 Bornean orangutan males, collected over 18 years. Migrants' peering rates significantly increased with higher food availability and time spent in proximity to others. Furthermore, migrants in the more sociable Sumatran population exhibited significantly higher peering rates compared to the Borneans, suggesting intrinsic and/or developmental effects of food availability on social learning. These findings emphasize the importance of investigating ecological effects on social learning on the immediate, developmental, and intrinsic levels for our understanding of cultural evolution.

Generation of induced pluripotent stem cells from Bornean orangutans

Introduction: Orangutans, classified under the Pongo genus, are an endangered non-human primate (NHP) species. Derivation of induced pluripotent stem cells (iPSCs) represents a promising avenue for conserving the genetic resources of these animals. Earlier studies focused on deriving orangutan iPSCs (o-iPSCs) from Sumatran orangutans (Pongo abelii). To date, no reports specifically target the other Critically Endangered species in the Pongo genus, the Bornean orangutans (Pongo pygmaeus). Methods: Using Sendai virus-mediated Yamanaka factor-based reprogramming of peripheral blood mononuclear cells to generate iPSCs (bo-iPSCs) from a female captive Bornean orangutan. In this study, we evaluate the colony morphology, pluripotent markers, X chromosome activation status, and transcriptomic profile of the bo-iPSCs to demonstrate the pluripotency of iPSCs from Bornean orangutans. Results: The bo-iPSCs were successfully derived from Bornean orangutans, using Sendai virus-mediated Yamanaka factor-based reprogramming of peripheral blood mononuclear cells. When a modified 4i/L/A (m4i/L/A) culture system was applied to activate the WNT signaling pathway in these bo-iPSCs, the derived cells (m-bo-iPSCs) manifested characteristics akin to human naive pluripotent stem cells, including high expression levels of KLF17, DNMT3L, and DPPA3/5, as well as the X chromosome reactivation. Comparative RNA-seq analysis positioned the m-bo-iPSCs between human naive and formative pluripotent states. Furthermore, the m-bo-iPSCs express differentiation capacity into all three germlines, evidenced by controlled in vitro embryoid body formation assay. Discussion: Our work establishes a novel approach to preserve the genetic diversity of endangered Bornean orangutans while offering insights into primate stem cell pluripotency. In the future, derivation of the primordial germ cell-like cells (PGCLCs) from m-bo-iPSCs is needed to demonstrate the further specific application in species preservation and broaden the knowledge of primordial germ cell specification across species.

The haplotype-resolved telomere-to-telomere carnation (Dianthus caryophyllus) genome reveals the correlation between genome architecture and gene expression

Carnation (Dianthus caryophyllus) is one of the most valuable commercial flowers, due to its richness of color and form, and its excellent storage and vase life. The diverse demands of the market require faster breeding in carnations. A full understanding of carnations is therefore required to guide the direction of breeding. Hence, we assembled the haplotype-resolved gap-free carnation genome of the variety 'Baltico', which is the most common white standard variety worldwide. Based on high-depth HiFi, ultra-long nanopore, and Hi-C sequencing data, we assembled the telomere-to-telomere (T2T) genome to be 564 479 117 and 568 266 215 bp for the two haplotypes Hap1 and Hap2, respectively. This T2T genome exhibited great improvement in genome assembly and annotation results compared with the former version. The improvements were seen when different approaches to evaluation were used. Our T2T genome first informs the analysis of the telomere and centromere region, enabling us to speculate about specific centromere characteristics that cannot be identified by high-order repeats in carnations. We analyzed allele-specific expression in three tissues and the relationship between genome architecture and gene expression in the haplotypes. This demonstrated that the length of the genes, coding sequences, and introns, the exon numbers and the transposable element insertions correlate with gene expression ratios and levels. The insertions of transposable elements repress expression in gene regulatory networks in carnation. This gap-free finished T2T carnation genome provides a valuable resource to illustrate the genome characteristics and for functional genomics analysis in further studies and molecular breeding.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Centromere innovations within a mouse species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying centromere protein-A (CENP-A) nucleosomes at the nexus of a satellite repeat that we identified and termed π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 mega-base pairs of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance promotes accumulation of microtubule-binding components of the kinetochore and a microtubule-destabilizing kinesin of the inner centromere. We propose that the balance of pro- and anti-microtubule binding by the new centromere is what permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

An expedition in the jungle of pluripotent stem cells of non-human primates

For nearly three decades, more than 80 embryonic stem cell lines and more than 100 induced pluripotent stem cell lines have been derived from New World monkeys, Old World monkeys, and great apes. In this comprehensive review, we examine these cell lines originating from marmoset, cynomolgus macaque, rhesus macaque, pig-tailed macaque, Japanese macaque, African green monkey, baboon, chimpanzee, bonobo, gorilla, and orangutan. We outline the methodologies implemented for their establishment, the culture protocols for their long-term maintenance, and their basic molecular characterization. Further, we spotlight any cell lines that express fluorescent reporters. Additionally, we compare these cell lines with human pluripotent stem cell lines, and we discuss cell lines reprogrammed into a pluripotent naive state, detailing the processes used to attain this. Last, we present the findings from the application of these cell lines in two emerging fields: intra- and interspecies embryonic chimeras and blastoids.

Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes

Y chromosomal ampliconic genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been studied in great apes; however, the diversity of splicing variants remains unexplored. Here, we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this data set resulted in several findings. First, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Second, our results suggest that BPY2 transcripts and proteins originate from separate genomic regions in bonobo versus human, which is possibly facilitated by acquiring new promoters. Third, our analysis indicates that the PRY gene family, having the highest representation of noncoding transcripts, has been undergoing pseudogenization. Fourth, we have not detected signatures of selection in the five YAG families shared among great apes, even though we identified many species-specific protein-coding transcripts. Fifth, we predicted consensus disorder regions across most gene families and species, which could be used for future investigations of male infertility. Overall, our work illuminates the YAG isoform landscape and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Decoding DNA sequence-driven evolution of the human brain epigenome at cellular resolution

DNA-based evolutionary comparisons of regulatory genomic elements enable insight into functional changes, overcoming tissue inaccessibility. Here, we harnessed adult and fetal cortex single-cell ATAC-seq datasets to uncover DNA substitutions specific to the human and human-ancestral lineages within apes. We found that fetal microglia identity is evolutionarily divergent in all lineages, whereas other cell types are conserved. Using multiomic datasets, we further identified genes linked to multiple lineage-divergent gene regulatory elements and implicated biological pathways associated with these divergent features. We also uncovered patterns of transcription factor binding site evolution across lineages and identified expansion of bHLH-PAS factor targets in human-hominin lineages, and MEF2 factor targets in the ape lineage. Finally, conserved features were more enriched in brain disease variants, whereas there was no distinct enrichment on the human lineage compared to its ancestral lineages. Our study identifies major evolutionary patterns in the human brain epigenome at cellular resolution.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

A global catalog of whole-genome diversity from 233 primate species

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

Phylogenomic analyses provide insights into primate evolution

Comparative analysis of primate genomes within a phylogenetic context is essential for understanding the evolution of human genetic architecture and primate diversity. We present such a study of 50 primate species spanning 38 genera and 14 families, including 27 genomes first reported here, with many from previously less well represented groups, the New World monkeys and the Strepsirrhini. Our analyses reveal heterogeneous rates of genomic rearrangement and gene evolution across primate lineages. Thousands of genes under positive selection in different lineages play roles in the nervous, skeletal, and digestive systems and may have contributed to primate innovations and adaptations. Our study reveals that many key genomic innovations occurred in the Simiiformes ancestral node and may have had an impact on the adaptive radiation of the Simiiformes and human evolution.

Centromere Innovations Within a Mouse Species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying CENP-A nucleosomes at the nexus of a satellite repeat that we identified and term π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 Mbp of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance drives accumulation of microtubule-binding components of the kinetochore, as well as a microtubule-destabilizing kinesin of the inner centromere. The balance of pro- and anti-microtubule-binding by the new centromere permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

Harnessing the omics revolution to address the global biodiversity crisis

Human disturbances are altering global biodiversity in unprecedented ways. We identify three fundamental challenges underpinning our understanding of global biodiversity (namely discovery, loss, and preservation), and discuss how the omics revolution (e.g. genomics, transcriptomics, proteomics, metabolomics, and meta-omics) can help address these challenges. We also discuss how omics tools can illuminate the major drivers of biodiversity loss, including invasive species, pollution, urbanization, overexploitation, and climate change, with a special focus on highly diverse tropical environments. Although omics tools are transforming the traditional toolkit of biodiversity research, their application to addressing the current biodiversity crisis remains limited and may not suffice to offset current rates of biodiversity loss. Despite technical and logistical challenges, omics tools need to be fully integrated into global biodiversity research, and better strategies are needed to improve their translation into biodiversity policy and practice. It is also important to recognize that although the omics revolution can be considered the biologist's dream, socioeconomic disparity limits their application in biodiversity research.

Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes

Y-chromosomal Ampliconic Genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been recently studied in great apes, however, the diversity of splicing variants remains unexplored. Here we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture-probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this dataset resulted in several findings. First, we uncovered a high diversity of YAG transcripts across great apes. Second, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Our results suggest that BPY2 transcripts and predicted proteins in several great ape species (bonobo and the two orangutans) have independent evolutionary origins and are not homologous to human reference transcripts and proteins. In contrast, our results suggest that the PRY gene family, having the highest representation of transcripts without open reading frames, has been undergoing pseudogenization. Third, even though we have identified many species-specific protein-coding YAG transcripts, we have not detected any signatures of positive selection. Overall, our work illuminates the YAG isoform landscape and its evolutionary history, and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.

Structurally divergent and recurrently mutated regions of primate genomes

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

A composite menstrual cycle of captive orangutans, with associated hormonal and behavioral variability

Knowledge of species-typical reproductive endocrinology profiles is crucial for testing hypotheses pertaining to the evolutionary history, reproductive parameters, and life history of a species, and for managing the well-being of individual animals in human care. Large-scale empirical measurements of ovarian hormones, however, are rare for most primate species, including orangutans. In this study, we used enzyme immunoassays (EIA) to quantify estrogen (estrone conjugates; E₁ C) and progesterone (pregnanediol-3-glucuronide; PdG) levels for 98 cycles in 7 cycling zoo-housed female orangutans (10-43 years old). We use a subset of these cycles (N = 44) to create the first composite menstrual cycle for orangutans, which serves as a valuable baseline for future comparative analyses and veterinary considerations. Similar to previous studies, we determined the mean ovarian cycle length of orangutans to be 29.7 days (N = 98 cycles), although we illustrate evidence of both intra- and interindividual variation in ovarian steroid production. Given that this study took place in captivity, we consider how energetic and psychosocial aspects of the zoo environment, such as greater food availability and potential stress, may affect the reproductive physiology and sexual behavior of these females. Furthermore, we discuss the role that age and genetic background may play in producing variability. Finally, we test whether ovarian hormone levels correlate with the reproductive behaviors of these female orangutans using associated behavioral data. Our results suggest that matings are more common during the periovulatory period than outside of it, but do not support a consistent link between hormonal indices of fecundability and mating behaviors in these individuals.

GADMA2: more efficient and flexible demographic inference from genetic data

BACKGROUND

Inference of complex demographic histories is a source of information about events that happened in the past of studied populations. Existing methods for demographic inference typically require input from the researcher in the form of a parameterized model. With an increased variety of methods and tools, each with its own interface, the model specification becomes tedious and error-prone. Moreover, optimization algorithms used to find model parameters sometimes turn out to be inefficient, for instance, by being not properly tuned or highly dependent on a user-provided initialization. The open-source software GADMA addresses these problems, providing automatic demographic inference. It proposes a common interface for several likelihood engines and provides global parameters optimization based on a genetic algorithm.

RESULTS

Here, we introduce the new GADMA2 software and provide a detailed description of the added and expanded features. It has a renovated core code base, new likelihood engines, an updated optimization algorithm, and a flexible setup for automatic model construction. We provide a full overview of GADMA2 enhancements, compare the performance of supported likelihood engines on simulated data, and demonstrate an example of GADMA2 usage on 2 empirical datasets.

CONCLUSIONS

We demonstrate the better performance of a genetic algorithm in GADMA2 by comparing it to the initial version and other existing optimization approaches. Our experiments on simulated data indicate that GADMA2's likelihood engines are able to provide accurate estimations of demographic parameters even for misspecified models. We improve model parameters for 2 empirical datasets of inbred species.

Ape femoral-humeral rigidities and arboreal locomotion

OBJECTIVES

This study investigates patterns of bone functional adaptations in extant apes through comparing hindlimb to forelimb bone rigidity ratios in groups with varying levels of arboreality.

MATERIALS AND METHODS

Using CT scans, bone rigidity (J) was calculated at three regions of interest (ROI) along femoral and humeral diaphyses in Homo, Pongo, Pan, and Gorilla with further comparisons made between species and subspecies divisions within Pan and Gorilla.

RESULTS

Consistent with previous work on extant hominoids, species exhibited differences in midshaft femoral to humeral (F/H) rigidity ratios. Results of the present study confirm that these midshaft differences extend to 35% and 65% diaphyseal ROIs. Modern humans, exhibiting larger ratios, and orangutans, exhibiting smaller ratios, bracketed the intermediate African apes in comparisons. Within some African apes, limb rigidity ratios varied significantly between taxonomic groups. Eastern gorillas exhibited the highest mean ratios and chimpanzees the lowest at all three ROIs. In posthoc comparisons, chimpanzees and bonobos did not differ in relative limb rigidity ratios at any of the three ROIs. However, western gorillas were more similar to bonobos than eastern gorillas at 50% and 35% ROIs, but not at the 65% ROI.

CONCLUSION

Species, and to a lesser extent subspecies, can be distinguished by F/H limb rigidity ratios according to broad positional behavior patterns at multiple regions of interest along the diaphyses. Similarity of bonobos and western gorillas is in line with behavioral data of bonobos being the most terrestrial of Pan species, and western gorillas the most arboreal of the Gorilla groups.

The coevolution between APOBEC3 and retrotransposons in primates

Retrotransposons are genetic elements with the ability to replicate in the genome using reverse transcriptase: they have been associated with the development of different biological structures, such as the Central Nervous System (CNS), and their high mutagenic potential has been linked to various diseases, including cancer and neurological disorders. Throughout evolution and over time, Primates and Homo had to cope with infections from viruses and bacteria, and also with endogenous retroelements. Therefore, host genomes have evolved numerous methods to counteract the activity of endogenous and exogenous pathogens, and the APOBEC3 family of mutators is a prime example of a defensive mechanism in this context.In most Primates, there are seven members of the APOBEC3 family of deaminase proteins: among their functions, there is the ability to inhibit the mobilization of retrotransposons and the functionality of viruses. The evolution of the APOBEC3 proteins found in Primates is correlated with the expansion of two major families of retrotransposons, i.e. ERV and LINE-1.In this review, we will discuss how the rapid expansion of the APOBEC3 family is linked to the evolution of retrotransposons, highlighting the strong evolutionary arms race that characterized the history of APOBEC3s and endogenous retroelements in Primates. Moreover, the possible role of this relationship will be assessed in the context of embryonic development and brain-associated diseases.

Human centromere repositioning activates transcription and opens chromatin fibre structure

Human centromeres appear as constrictions on mitotic chromosomes and form a platform for kinetochore assembly in mitosis. Biophysical experiments led to a suggestion that repetitive DNA at centromeric regions form a compact scaffold necessary for function, but this was revised when neocentromeres were discovered on non-repetitive DNA. To test whether centromeres have a special chromatin structure we have analysed the architecture of a neocentromere. Centromere repositioning is accompanied by RNA polymerase II recruitment and active transcription to form a decompacted, negatively supercoiled domain enriched in ‘open’ chromatin fibres. In contrast, centromerisation causes a spreading of repressive epigenetic marks to surrounding regions, delimited by H3K27me3 polycomb boundaries and divergent genes. This flanking domain is transcriptionally silent and partially remodelled to form ‘compact’ chromatin, similar to satellite-containing DNA sequences, and exhibits genomic instability. We suggest transcription disrupts chromatin to provide a foundation for kinetochore formation whilst compact pericentromeric heterochromatin generates mechanical rigidity.

Evidence for the latest fossil Pongo in southern China

Pongo fossils with precise absolute age brackets are rare, limiting our understanding of their taxonomy and spatiotemporal distribution in southern China during the Late Pleistocene. Twenty-four isolated teeth of fossil orangutans were recently discovered during excavations at Yicun Cave in Guangxi Zhuang Autonomous Region, southern China. Here, we dated the fossil-bearing layer using Uranium-series dating of the associated flowstone and soda straw stalactites. Our results date the Yicun orangutan fossils to between 66 ± 0.32 ka and 57 ± 0.26 ka; thus, these fossils currently represent the last appearance datum of Pongo in southern China. We further conducted a detailed morphological comparison of the Yicun fossil teeth with large samples of fossil (n = 2454) and extant (n = 441) orangutans from mainland and island Southeast Asia to determine their taxonomic position. Compared to other fossil and extant orangutan samples, the Yicun Pongo assemblage has larger teeth and displays greater variation in occlusal structure. Based on the high frequency of cingular remnants and light to moderate enamel wrinkling of the molars, we assigned the Yicun fossils to Pongo weidenreichi, a species that was widespread in southern China throughout the Pleistocene. Lastly, we used published stable carbon isotope data of Early to Late Pleistocene mammalian fossil teeth from mainland Southeast Asia to reconstruct changes in the paleoenvironment and to interpret dental size variation of Pongo assemblages in a broader temporal and environmental context. The carbon isotope data show that dental size reduction in Pongo is associated with environmental changes. These morphological changes in Pongo appear to coincide with the expansion of savannah biomes and the contraction of forest habitats from the Middle Pleistocene onward. The variation in dental size of forest-dwelling Pongo in mainland Southeast Asia may have resulted from habitat differentiation during the Pleistocene.

Nine out of ten samples were mistakenly switched by The Orang-utan Genome Consortium

The Sumatran orang-utan (Pongo abelii) reference genome was first published in 2011, in conjunction with ten re-sequenced genomes from unrelated wild-caught individuals. Together, these published data have been utilized in almost all great ape genomic studies, plus in much broader comparative genomic research. Here, we report that the original sequencing Consortium inadvertently switched nine of the ten samples and/or resulting re-sequenced genomes, erroneously attributing eight of these to the wrong source individuals. Among them is a genome from the recently identified Tapanuli (P. tapanuliensis) species: thus, this genome was sequenced and published a full six years prior to the species’ description. Sex was wrongly assigned to five known individuals; the numbers in one sample identifier were swapped; and the identifier for another sample most closely resembles that of a sample from another individual entirely. These errors have been reproduced in countless subsequent manuscripts, with noted implications for studies reliant on data from known individuals.

Evolution of eukaryotic centromeres by drive and suppression of selfish genetic elements

Despite the universal requirement for faithful chromosome segregation, eukaryotic centromeres are rapidly evolving. It is hypothesized that rapid centromere evolution represents an evolutionary arms race between selfish genetic elements that drive, or propagate at the expense of organismal fitness, and mechanisms that suppress fitness costs. Selfish centromere DNA achieves preferential inheritance in female meiosis by recruiting more effector proteins that alter spindle microtubule interaction dynamics. Parallel pathways for effector recruitment are adaptively evolved to suppress functional differences between centromeres. Opportunities to drive are not limited to female meiosis, and selfish transposons, plasmids and B chromosomes also benefit by maximizing their inheritance. Rapid evolution of selfish genetic elements can diversify suppressor mechanisms in different species that may cause hybrid incompatibility.

The evolutionary history of human spindle genes includes back-and-forth gene flow with Neandertals

Proteins associated with the spindle apparatus, a cytoskeletal structure that ensures the proper segregation of chromosomes during cell division, experienced an unusual number of amino acid substitutions in modern humans after the split from the ancestors of Neandertals and Denisovans. Here, we analyze the history of these substitutions and show that some of the genes in which they occur may have been targets of positive selection. We also find that the two changes in the kinetochore scaffold 1 (KNL1) protein, previously believed to be specific to modern humans, were present in some Neandertals. We show that the KNL1 gene of these Neandertals shared a common ancestor with present-day Africans about 200,000 years ago due to gene flow from the ancestors (or relatives) of modern humans into Neandertals. Subsequently, some non-Africans inherited this modern human-like gene variant from Neandertals, but none inherited the ancestral gene variants. These results add to the growing evidence of early contacts between modern humans and archaic groups in Eurasia and illustrate the intricate relationships among these groups.

Evolution of holocentric chromosomes: Drivers, diversity, and deterrents

Centromeres are specialized chromosomal regions that recruit kinetochore proteins and mediate spindle microtubule attachment to ensure faithful chromosome segregation during mitosis and meiosis. Centromeres can be restricted to one region of the chromosome. Named "monocentromere", this type represents the most commonly found centromere organization across eukaryotes. Alternatively, centromeres can also be assembled at sites chromosome-wide. This second type is called "holocentromere". Despite their early description over 100 years ago, research on holocentromeres has lagged behind that of monocentromeres. Nevertheless, the application of next generation sequencing approaches and advanced microscopic technologies enabled recent advances understanding the molecular organization and regulation of holocentromeres in different organisms. Here we review the current state of research on holocentromeres focusing on evolutionary considerations. First, we provide a brief historical perspective on the discovery of holocentric chromosomes. We then discuss models/drivers that have been proposed over the years to explain the evolutionary transition from mono- to holocentric chromosomes. We continue to review the description of holocentric chromosomes in diverse eukaryotic groups and then focus our discussion on a specific and recently characterized type of holocentromere organization in insects that functions independently of the otherwise essential centromeric marker protein CenH3, thus providing novel insights into holocentromere evolution in insects. Finally, we propose reasons to explain why the holocentric trait is not more frequent across eukaryotes despite putative selective advantages.

Phyloepigenetics

Traditionally, phylogenetic interspecies relationships are estimated based on genetic diversity, since it is assumed that the more recently diverged a species, with comparable constancy of development, the more similar their genetic material and proteins should be. However, occasional controversies in the field may reflect limited resolution and accuracy of this approach. Epigenetics has, meanwhile, provided significant evidence that CpG dinucleotides (CpGs) within genetic material are of particular importance for the annotation and function of the genome and the formation of the phenotype, which is continuously shaped by evolutionary interaction with environmental factors. Based on this, it can be concluded that CpGs follow a distinct rate of evolution, compared to all other nucleotide positions. Evidence is provided that supports this conclusion. Therefore, using CpGs to fathom evolutionary relationships between species could turn out to be a valuable approach to achieve, in some cases, an improved understanding of evolutionary development.

Molecular Dynamics and Evolution of Centromeres in the Genus Equus

The centromere is the chromosomal locus essential for proper chromosome segregation. While the centromeric function is well conserved and epigenetically specified, centromeric DNA sequences are typically composed of satellite DNA and represent the most rapidly evolving sequences in eukaryotic genomes. The presence of satellite sequences at centromeres hampered the comprehensive molecular analysis of these enigmatic loci. The discovery of functional centromeres completely devoid of satellite repetitions and fixed in some animal and plant species represented a turning point in centromere biology, definitively proving the epigenetic nature of the centromere. The first satellite-free centromere, fixed in a vertebrate species, was discovered in the horse. Later, an extraordinary number of satellite-free neocentromeres had been discovered in other species of the genus Equus, which remains the only mammalian genus with numerous satellite-free centromeres described thus far. These neocentromeres arose recently during evolution and are caught in a stage of incomplete maturation. Their presence made the equids a unique model for investigating, at molecular level, the minimal requirements for centromere seeding and evolution. This model system provided new insights on how centromeres are established and transmitted to the progeny and on the role of satellite DNA in different aspects of centromere biology.

The contribution of evolutionarily volatile promoters to molecular phenotypes and human trait variation

BACKGROUND

Promoters are sites of transcription initiation that harbour a high concentration of phenotype-associated genetic variation. The evolutionary gain and loss of promoters between species (collectively, termed turnover) is pervasive across mammalian genomes and may play a prominent role in driving human phenotypic diversity.

RESULTS

We classified human promoters by their evolutionary history during the divergence of mouse and human lineages from a common ancestor. This defined conserved, human-inserted and mouse-deleted promoters, and a class of functional-turnover promoters that align between species but are only active in humans. We show that promoters of all evolutionary categories are hotspots for substitution and often, insertion mutations. Loci with a history of insertion and deletion continue that mode of evolution within contemporary humans. The presence of an evolutionary volatile promoter within a gene is associated with increased expression variance between individuals, but only in the case of human-inserted and mouse-deleted promoters does that correspond to an enrichment of promoter-proximal genetic effects. Despite the enrichment of these molecular quantitative trait loci (QTL) at evolutionarily volatile promoters, this does not translate into a corresponding enrichment of phenotypic traits mapping to these loci.

CONCLUSIONS

Promoter turnover is pervasive in the human genome, and these promoters are rich in molecularly quantifiable but phenotypically inconsequential variation in gene expression. However, since evolutionarily volatile promoters show evidence of selection, coupled with high mutation rates and enrichment of QTLs, this implicates them as a source of evolutionary innovation and phenotypic variation, albeit with a high background of selectively neutral expression variation.

Mosaic cis-regulatory evolution drives transcriptional partitioning of HERVH endogenous retrovirus in the human embryo

The human endogenous retrovirus type-H (HERVH) family is expressed in the preimplantation embryo. A subset of these elements are specifically transcribed in pluripotent stem cells where they appear to exert regulatory activities promoting self-renewal and pluripotency. How HERVH elements achieve such transcriptional specificity remains poorly understood. To uncover the sequence features underlying HERVH transcriptional activity, we performed a phyloregulatory analysis of the long terminal repeats (LTR7) of the HERVH family, which harbor its promoter, using a wealth of regulatory genomics data. We found that the family includes at least eight previously unrecognized subfamilies that have been active at different timepoints in primate evolution and display distinct expression patterns during human embryonic development. Notably, nearly all HERVH elements transcribed in ESCs belong to one of the youngest subfamilies we dubbed LTR7up. LTR7 sequence evolution was driven by a mixture of mutational processes, including point mutations, duplications, and multiple recombination events between subfamilies, that led to transcription factor binding motif modules characteristic of each subfamily. Using a reporter assay, we show that one such motif, a predicted SOX2/3 binding site unique to LTR7up, is essential for robust promoter activity in induced pluripotent stem cells. Together these findings illuminate the mechanisms by which HERVH diversified its expression pattern during evolution to colonize distinct cellular niches within the human embryo.

I Wanna Draw Like You: Inter- and Intra-Individual Differences in Orang-Utan Drawings

This study analyses 749 drawings by five female Bornean orang-utans (Pongo pygmaeus) at Tama Zoological Park in Japan. We searched for differences between individuals but also tried to identify possible temporal changes among the drawings of one individual, Molly, who drew almost 1300 drawings from 2006 to 2011. An analysis of the drawings was carried out after collecting quantitative and qualitative variables. Our findings reveal evidence of differences in the drawing style of the five individuals as well as creative changes in Molly's drawing style throughout her lifetime. Individuals differed in terms of the colours used, the space they filled, and the shapes (fan patterns, circles, or loops) they drew. Molly drew less and less as she grew older, and we found a significant difference between drawings produced in winter, when orang-utans were kept inside and had less activity, and those produced during other seasons. Our results suggest that the drawing behaviour of these five orang-utans is not random and that differences among individuals might reflect differences of styles, states of mind, and motivation to draw.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

Gene replacement strategies validate the use of functional tags on centromeric chromatin and invalidate an essential role for CENP-AK124ub

Functional tags are ubiquitous in cell biology, and for studies of one chromosomal locus, the centromere, tags have been remarkably useful. The centromere directs chromosome inheritance at cell division. The location of the centromere is defined by a histone H3 variant, CENP-A. The regulation of the chromatin assembly pathway essential for centromere inheritance and function includes posttranslational modification (PTM) of key components, including CENP-A itself. Others have recently called into question the use of functional tags, with the claim that at least two widely used tags obscured the essentiality of one particular PTM, CENP-A^K124 ubiquitination (ub). Here, we employ three independent gene replacement strategies that eliminate large, lysine-containing tags to interrogate these claims. Using these approaches, we find no evidence to support an essential function of CENP-A^K124ub. Our general methodology will be useful to validate discoveries permitted by powerful functional tagging schemes at the centromere and other cellular locations.

Using three gene replacement strategies, Salinas-Luypaert et al. demonstrate that CENP-A^K124ub is not essential for CENP-A function at centromeres. Thus, functional tags do not mask the role of K124 when it is mutated. These strategies can be employed to interrogate posttranslational modifications at the centromere and other cellular locations.

An etiology of human modernity

Following the refutation of the replacement hypothesis, which had proposed that a ‘superior’ hominin species arose in Africa and replaced all other humans existing at the time, the auto-domestication hypothesis remains the only viable explanation for the relatively abrupt change from robust to gracile humans in the Late Pleistocene. It invokes the incidental institution of the domestication syndrome in humans, most probably by newly introduced cultural practices. It also postulates that the induction of exograms compensated for the atrophy of the brain caused by domestication. This new explanation of the origins of modernity in humans elucidates practically all its many aspects, in stark contrast to the superseded replacement hypothesis, which explained virtually nothing. The first results of the domestication syndrome’s genetic exploration have become available in recent years, and they endorse the human self-domestication hypothesis.

Current Progress in Evolutionary Comparative Genomics of Great Apes

The availability of high-quality genome sequences of great ape species provides unprecedented opportunities for genomic analyses. Herein, we reviewed the recent progress in evolutionary comparative genomic studies of the existing great ape species, including human, chimpanzee, bonobo, gorilla, and orangutan. We elaborate discovery on evolutionary history, natural selection, structural variations, and new genes of these species, which is informative for understanding the origin of human-specific phenotypes.

Short Tandem Repeats as a High-Resolution Marker for Capturing Recent Orangutan Population Evolution

The genus Pongo is ideal to study population genetics adaptation, given its remarkable phenotypic divergence and the highly contrasting environmental conditions it’s been exposed to. Studying its genetic variation bears the promise to reveal a motion picture of these great apes’ evolutionary and adaptive history, and also helps us expand our knowledge of the patterns of adaptation and evolution. In this work, we advance the understanding of the genetic variation among wild orangutans through a genome-wide study of short tandem repeats (STRs). Their elevated mutation rate makes STRs ideal markers for the study of recent evolution within a given population. Current technological and algorithmic advances have rendered their sequencing and discovery more accurate, therefore their potential can be finally leveraged in population genetics studies. To study patterns of population variation within the wild orangutan population, we genotyped the short tandem repeats in a population of 21 individuals spanning four Sumatran and Bornean (sub-) species and eight Southeast Asian regions. We studied the impact of sequencing depth on our ability to genotype STRs and found that the STR copy number changes function as a powerful marker, correctly capturing the demographic history of these populations, even the divergences as recent as 10 Kya. Moreover, gene ontology enrichments for genes close to STR variants are aligned with local adaptations in the two islands. Coupled with more advanced STR-compatible population models, and selection tests, genomic studies based on STRs will be able to reduce the gap caused by the missing heritability for species with recent adaptations.

MSABrowser: dynamic and fast visualization of sequence alignments, variations and annotations

Summary

Sequence alignment is an excellent way to visualize the similarities and differences between DNA, RNA or protein sequences, yet it is currently difficult to jointly view sequence alignment data with genetic variations, modifications such as post-translational modifications and annotations (i.e. protein domains). Here, we present the MSABrowser tool that makes it easy to co-visualize genetic variations, modifications and annotations on the respective positions of amino acids or nucleotides in pairwise or multiple sequence alignments. MSABrowser is developed entirely in JavaScript and works on any modern web browser at any platform, including Linux, Mac OS X and Windows systems without any installation. MSABrowser is also freely available for the benefit of the scientific community.

Availability and implementation

MSABrowser is released as open-source and web-based software under MIT License. The visualizer, documentation, all source codes and examples are available at https://thekaplanlab.github.io/ and GitHub repository https://github.com/thekaplanlab/msabrowser.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Estimating the age of single nucleotide polymorphic sites in humans

BACKGROUND

Determining the ages of polymorphic sites in human genomes needs to be carried out in a careful balance between the degree of complexity of the approach and the desired accuracy.

OBJECTIVE

We provide evidence that a simpler approach where age determination is based upon the degree to which the alternative allele is spread among the population can be competitive with more complex methods.

METHODS

The information contained in the vcf files of Phase 1 of the 1000 Genomes Project combined with the genomic sequences of seven non-human primate species was analyzed. The analyses were supplemented by a computer simulation of the mutational changes in 10,000 model chromosomes with a length of 10,000 nucleotides over a period of 16 million years. The list of the birth dates of the derived alleles of homozygous and heterozygous components of the derived alleles served as a yardstick to estimate the ages of human alternative alleles.

RESULTS

The age of the derived alleles born in Africa before the "Out of Africa" event and not brought to other continents are estimated to be 0.17 Ma, the derived alleles present simultaneously on all continents are estimated to be 1.3 Ma old and the age of alleles arising in Europe or Asia is 0.06 Ma.

CONCLUSION

Our approach functions with polymorphisms that respect the "more frequent means older" principle. However, this shortcoming only leads to disagreement with the Atlas of Variant Age in about 20% of cases.

Ampliconic Genes on the Great Ape Y Chromosomes: Rapid Evolution of Copy Number but Conservation of Expression Levels

Multicopy ampliconic gene families on the Y chromosome play an important role in spermatogenesis. Thus, studying their genetic variation in endangered great ape species is critical. We estimated the sizes (copy number) of nine Y ampliconic gene families in population samples of chimpanzee, bonobo, and orangutan with droplet digital polymerase chain reaction, combined these estimates with published data for human and gorilla, and produced genome-wide testis gene expression data for great apes. Analyzing this comprehensive data set within an evolutionary framework, we, first, found high inter- and intraspecific variation in gene family size, with larger families exhibiting higher variation as compared with smaller families, a pattern consistent with random genetic drift. Second, for four gene families, we observed significant interspecific size differences, sometimes even between sister species—chimpanzee and bonobo. Third, despite substantial variation in copy number, Y ampliconic gene families’ expression levels did not differ significantly among species, suggesting dosage regulation. Fourth, for three gene families, size was positively correlated with gene expression levels across species, suggesting that, given sufficient evolutionary time, copy number influences gene expression. Our results indicate high variability in size but conservation in gene expression levels in Y ampliconic gene families, significantly advancing our understanding of Y-chromosome evolution in great apes.

Selection and thermostability suggest G-quadruplexes are novel functional elements of the human genome

Approximately 1% of the human genome has the ability to fold into G-quadruplexes (G4s)-noncanonical strand-specific DNA structures forming at G-rich motifs. G4s regulate several key cellular processes (e.g., transcription) and have been hypothesized to participate in others (e.g., firing of replication origins). Moreover, G4s differ in their thermostability, and this may affect their function. Yet, G4s may also hinder replication, transcription, and translation and may increase genome instability and mutation rates. Therefore, depending on their genomic location, thermostability, and functionality, G4 loci might evolve under different selective pressures, which has never been investigated. Here we conducted the first genome-wide analysis of G4 distribution, thermostability, and selection. We found an overrepresentation, high thermostability, and purifying selection for G4s within genic components in which they are expected to be functional-promoters, CpG islands, and 5' and 3' UTRs. A similar pattern was observed for G4s within replication origins, enhancers, eQTLs, and TAD boundary regions, strongly suggesting their functionality. In contrast, G4s on the nontranscribed strand of exons were underrepresented, were unstable, and evolved neutrally. In general, G4s on the nontranscribed strand of genic components had lower density and were less stable than those on the transcribed strand, suggesting that the former are avoided at the RNA level. Across the genome, purifying selection was stronger at stable G4s. Our results suggest that purifying selection preserves the sequences of functional G4s, whereas nonfunctional G4s are too costly to be tolerated in the genome. Thus, G4s are emerging as fundamental, functional genomic elements.

The genome of the Pyrenean desman and the effects of bottlenecks and inbreeding on the genomic landscape of an endangered species

The Pyrenean desman (Galemys pyrenaicus) is a small semiaquatic mammal endemic to the Iberian Peninsula. Despite its limited range, this species presents a strong genetic structure due to past isolation in glacial refugia and subsequent bottlenecks. Additionally, some populations are highly fragmented today as a consequence of river barriers, causing substantial levels of inbreeding. These features make the Pyrenean desman a unique model in which to study the genomic footprints of differentiation, bottlenecks and extreme isolation in an endangered species. To understand these processes, the complete genome of the Pyrenean desman was sequenced and assembled using a Bloom filter-based approach. An analysis of the 1.83 Gb reference genome and the sequencing of five additional individuals from different evolutionary units allowed us to detect its main genomic characteristics. The population differentiation of the species was reflected in highly distinctive demographic trajectories. In addition, a severe population bottleneck during the postglacial recolonization of the eastern Pyrenees created one of the lowest genomic heterozygosity values recorded in a mammal. Moreover, isolation and inbreeding gave rise to a high proportion of runs of homozygosity (ROH). Despite these extremely low levels of genetic diversity, two key multigene families from an eco-evolutionary perspective, the major histocompatibility complex and olfactory receptor genes, showed heterozygosity excess in the majority of individuals, revealing that functional diversity can be maintained up to a certain extent. Furthermore, these two classes of genes were significantly less abundant than expected within ROH. In conclusion, the genomic landscape of each analysed Pyrenean desman turned out to be strikingly distinctive and was a clear reflection of its recent ancestry and current conservation conditions. These results may help characterize the genomic health of each individual, and can be crucial for the conservation and management of the species.