Insights into hominid evolution from the gorilla genome sequence

Estimating ancestral ranges and biogeographical processes in early hominins

Historical biogeography provides crucial insights into understanding the evolutionary history of hominins. We applied maximum-likelihood and biogeographical stochastic mapping to infer the ancestral ranges of hominins and estimate the frequency of biogeographical events. These events were inferred using two time-calibrated phylogenetic trees that differ in the position of Australopithecus sediba. Results suggest that regardless of which phylogeny was selected, Northcentral Africa was the preferred ancestral region for the ancestor of the Homo-Pan clade, as well as the ancestor of Sahelanthropus and later hominins. The northern and middle part of eastern Africa was the preferred ancestral region for several clades originating at subsequent deep nodes of the trees (∼5-4 Ma). The choice of tree topology had one important effect on results: whether hominin ancestors appearing after ∼4 Ma were widespread or endemic. These different patterns highlight the biogeographic significance of the phylogenetic relationships of A. sediba. Overall, the results showed that dispersal, local extinction, and sympatry played vital roles in creating the hominin distribution, whereas vicariance and jump dispersal were not as common. The results suggested symmetry in the directionality of dispersals. Distance probably influenced how rapidly taxa colonized a new region, and dispersals often followed the closest path. These findings are potentially impacted by the imperfection of the fossil record, suggesting that the results should be interpreted cautiously.

Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

Deep-time phylogenetic inference by paleoproteomic analysis of dental enamel

In temperate and subtropical regions, ancient proteins are reported to survive up to about 2 million years, far beyond the known limits of ancient DNA preservation in the same areas. Accordingly, their amino acid sequences currently represent the only source of genetic information available to pursue phylogenetic inference involving species that went extinct too long ago to be amenable for ancient DNA analysis. Here we present a complete workflow, including sample preparation, mass spectrometric data acquisition and computational analysis, to recover and interpret million-year-old dental enamel protein sequences. During sample preparation, the proteolytic digestion step, usually an integral part of conventional bottom-up proteomics, is omitted to increase the recovery of the randomly degraded peptides spontaneously generated by extensive diagenetic hydrolysis of ancient proteins over geological time. Similarly, we describe other solutions we have adopted to (1) authenticate the endogenous origin of the protein traces we identify, (2) detect and validate amino acid variation in the ancient protein sequences and (3) attempt phylogenetic inference. Sample preparation and data acquisition can be completed in 3-4 working days, while subsequent data analysis usually takes 2-5 days. The workflow described requires basic expertise in ancient biomolecules analysis, mass spectrometry-based proteomics and molecular phylogeny. Finally, we describe the limits of this approach and its potential for the reconstruction of evolutionary relationships in paleontology and paleoanthropology.

A human-specific insertion promotes cell proliferation and migration by enhancing TBC1D8B expression

Human-specific insertions play important roles in human phenotypes and diseases. Here we reported a 446-bp insertion (Insert-446) in intron 11 of the TBC1D8B gene, located on chromosome X, and traced its origin to a portion of intron 6 of the EBF1 gene on chromosome 5. Interestingly, Insert-446 was present in the human Neanderthal and Denisovans genomes, and was fixed in humans after human-chimpanzee divergence. We have demonstrated that Insert-446 acts as an enhancer through binding transcript factors that promotes a higher expression of human TBC1D8B gene as compared with orthologs in macaques. In addition, over-expression TBC1D8B promoted cell proliferation and migration through "a dual finger" catalytic mechanism (Arg538 and Gln573) in the TBC domain in vitro and knockdown of TBC1D8B attenuated tumorigenesis in vivo. Knockout of Insert-446 prevented cell proliferation and migration in cancer and normal cells. Our results reveal that the human-specific Insert-446 promotes cell proliferation and migration by upregulating the expression of TBC1D8B gene. These findings provide a significant insight into the effects of human-specific insertions on evolution.

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.

Primate archaeology 3.0

The new field of primate archaeology investigates the technological behavior and material record of nonhuman primates, providing valuable comparative data on our understanding of human technological evolution. Yet, paralleling hominin archaeology, the field is largely biased toward the analysis of lithic artifacts. While valuable comparative data have been gained through an examination of extant nonhuman primate tool use and its archaeological record, focusing on this one single aspect provides limited insights. It is therefore necessary to explore to what extent other non-technological activities, such as non-tool aided feeding, traveling, social behaviors or ritual displays, leave traces that could be detected in the archaeological record. Here we propose four new areas of investigation which we believe have been largely overlooked by primate archaeology and that are crucial to uncovering the full archaeological potential of the primate behavioral repertoire, including that of our own: (1) Plant technology; (2) Archaeology beyond technology; (3) Landscape archaeology; and (4) Primate cultural heritage. We discuss each theme in the context of the latest developments and challenges, as well as propose future directions. Developing a more "inclusive" primate archaeology will not only benefit the study of primate evolution in its own right but will aid conservation efforts by increasing our understanding of changes in primate-environment interactions over time.

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.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000-180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer's gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer's gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

Functional host-specific adaptation of the intestinal microbiome in hominids

Fine-scale knowledge of the changes in composition and function of the human gut microbiome compared that of our closest relatives is critical for understanding the evolutionary processes underlying its developmental trajectory. To infer taxonomic and functional changes in the gut microbiome across hominids at different timescales, we perform high-resolution metagenomic-based analyzes of the fecal microbiome from over two hundred samples including diverse human populations, as well as wild-living chimpanzees, bonobos, and gorillas. We find human-associated taxa depleted within non-human apes and patterns of host-specific gut microbiota, suggesting the widespread acquisition of novel microbial clades along the evolutionary divergence of hosts. In contrast, we reveal multiple lines of evidence for a pervasive loss of diversity in human populations in correlation with a high Human Development Index, including evolutionarily conserved clades. Similarly, patterns of co-phylogeny between microbes and hosts are found to be disrupted in humans. Together with identifying individual microbial taxa and functional adaptations that correlate to host phylogeny, these findings offer insights into specific candidates playing a role in the diverging trajectories of the gut microbiome of hominids. We find that repeated horizontal gene transfer and gene loss, as well as the adaptation to transient microaerobic conditions appear to have played a role in the evolution of the human gut microbiome.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

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.

The promise of great apes as model organisms for understanding the downstream consequences of early life experiences

Early life experiences have a significant influence on adult health and aging processes in humans. Despite widespread interest in the evolutionary roots of this phenomenon, very little research on this topic has been conducted in humans' closest living relatives, the great apes. The longitudinal data sets that are now available on wild and captive great ape populations hold great promise to clarify the nature, evolutionary function, and mechanisms underlying these connections in species which share key human life history characteristics. Here, we explain features of great ape life history and socioecologies that make them of particular interest for this topic, as well as those that may limit their utility as comparative models; outline the ways in which available data are complementary to and extend the kinds of data that are available for humans; and review what is currently known about the connections among early life experiences, social behavior, and adult physiology and biological fitness in our closest living relatives. We conclude by highlighting key next steps for this emerging area of research.

Ghost admixture in eastern gorillas

Archaic admixture has had a substantial impact on human evolution with multiple events across different clades, including from extinct hominins such as Neanderthals and Denisovans into modern humans. In great apes, archaic admixture has been identified in chimpanzees and bonobos but the possibility of such events has not been explored in other species. Here, we address this question using high-coverage whole-genome sequences from all four extant gorilla subspecies, including six newly sequenced eastern gorillas from previously unsampled geographic regions. Using approximate Bayesian computation with neural networks to model the demographic history of gorillas, we find a signature of admixture from an archaic 'ghost' lineage into the common ancestor of eastern gorillas but not western gorillas. We infer that up to 3% of the genome of these individuals is introgressed from an archaic lineage that diverged more than 3 million years ago from the common ancestor of all extant gorillas. This introgression event took place before the split of mountain and eastern lowland gorillas, probably more than 40 thousand years ago and may have influenced perception of bitter taste in eastern gorillas. When comparing the introgression landscapes of gorillas, humans and bonobos, we find a consistent depletion of introgressed fragments on the X chromosome across these species. However, depletion in protein-coding content is not detectable in eastern gorillas, possibly as a consequence of stronger genetic drift in this species.

The Bovine Pangenome Consortium: democratizing production and accessibility of genome assemblies for global cattle breeds and other bovine species

The Bovine Pangenome Consortium (BPC) is an international collaboration dedicated to the assembly of cattle genomes to develop a more complete representation of cattle genomic diversity. The goal of the BPC is to provide genome assemblies and a community-agreed pangenome representation to replace breed-specific reference assemblies for cattle genomics. The BPC invites partners sharing our vision to participate in the production of these assemblies and the development of a common, community-approved, pangenome reference as a public resource for the research community ( https://bovinepangenome.github.io/ ). This community-driven resource will provide the context for comparison between studies and the future foundation for cattle genomic selection.

Pervasive incomplete lineage sorting illuminates speciation and selection in primates

Incomplete lineage sorting (ILS) causes the phylogeny of some parts of the genome to differ from the species tree. In this work, we investigate the frequencies and determinants of ILS in 29 major ancestral nodes across the entire primate phylogeny. We find up to 64% of the genome affected by ILS at individual nodes. We exploit ILS to reconstruct speciation times and ancestral population sizes. Estimated speciation times are much more recent than genomic divergence times and are in good agreement with the fossil record. We show extensive variation of ILS along the genome, mainly driven by recombination but also by the distance to genes, highlighting a major impact of selection on variation along the genome. In many nodes, ILS is reduced more on the X chromosome compared with autosomes than expected under neutrality, which suggests higher impacts of natural selection on the X chromosome. Finally, we show an excess of ILS in genes with immune functions and a deficit of ILS in housekeeping genes. The extensive ILS in primates discovered in this study provides insights into the speciation times, ancestral population sizes, and patterns of natural selection that shape primate evolution.

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.

Distinct Traits of Structural and Regulatory Evolutional Conservation of Human Genes with Specific Focus on Major Cancer Molecular Pathways

The evolution of protein-coding genes has both structural and regulatory components. The first can be assessed by measuring the ratio of non-synonymous to synonymous nucleotide substitutions. The second component can be measured as the normalized proportion of transposable elements that are used as regulatory elements. For the first time, we characterized in parallel the regulatory and structural evolutionary profiles for 10,890 human genes and 2972 molecular pathways. We observed a ~0.1 correlation between the structural and regulatory metrics at the gene level, which appeared much higher (~0.4) at the pathway level. We deposited the data in the publicly available database RetroSpect. We also analyzed the evolutionary dynamics of six cancer pathways of two major axes: Notch/WNT/Hedgehog and AKT/mTOR/EGFR. The Hedgehog pathway had both components slower, whereas the Akt pathway had clearly accelerated structural evolution. In particular, the major hub nodes Akt and beta-catenin showed both components strongly decreased, whereas two major regulators of Akt TCL1 and CTMP had outstandingly high evolutionary rates. We also noticed structural conservation of serine/threonine kinases and the genes related to guanosine metabolism in cancer signaling: GPCRs, G proteins, and small regulatory GTPases (Src, Rac, Ras); however, this was compensated by the accelerated regulatory evolution.

Extraordinary selection on the human X chromosome associated with archaic admixture

The X chromosome in non-African humans shows less diversity and less Neanderthal introgression than expected from neutral evolution. Analyzing 162 human male X chromosomes worldwide, we identified fourteen chromosomal regions where nearly identical haplotypes spanning several hundred kilobases are found at high frequencies in non-Africans. Genetic drift alone cannot explain the existence of these haplotypes, which must have been associated with strong positive selection in partial selective sweeps. Moreover, the swept haplotypes are entirely devoid of archaic ancestry as opposed to the non-swept haplotypes in the same genomic regions. The ancient Ust'-Ishim male dated at 45,000 before the present (BP) also carries the swept haplotypes, implying that selection on the haplotypes must have occurred between 45,000 and 55,000 years ago. Finally, we find that the chromosomal positions of sweeps overlap previously reported hotspots of selective sweeps in great ape evolution, suggesting a mechanism of selection unique to X chromosomes.

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.

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.

On the effects of selection and mutation on species tree inference

The effect of selection acting on regions of the genome on the accuracy of species-level phylogenetic inference using methods that do not explicitly model selection is an open question that is relevant to most, if not all, phylogenomic studies. To address this, we derive a mathematical approximation to the Wright-Fisher model with mutation and selection in the limit as the population size becomes large. In contrast to previous approximations based on diffusion processes, our approximation can be used to study the distribution of coalescent times for an arbitrary number of lineages, allowing calculation of the probability distribution of gene genealogies under the coalescent model. We use these calculations to show that direct selection at strengths typically encountered in practice has only a small effect on the distribution of coalescent times, and hence on the distribution of gene trees. This implies that many coalescent-based methods for estimating the species tree topology will be robust to the presence of selection in a subset of the underlying genes. Selection will, however, bias the estimation of speciation times, causing them to underestimate the true speciation times. Our model captures the effects of selection on the genealogies that generate the observed sequence data, but does not model selective pressures that act only on the subsequent sequences or that negatively impact gene tree estimation.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

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.

Evolutionary dynamics of pseudoautosomal region 1 in humans and great apes

Background

The pseudoautosomal region 1 (PAR1) is a 2.7 Mb telomeric region of human sex chromosomes. PAR1 has a crucial role in ensuring proper segregation of sex chromosomes during male meiosis, exposing it to extreme recombination and mutation processes. We investigate PAR1 evolution using population genomic datasets of extant humans, eight populations of great apes, and two archaic human genome sequences.

Results

We find that PAR1 is fast evolving and closer to evolutionary nucleotide equilibrium than autosomal telomeres. We detect a difference between substitution patterns and extant diversity in PAR1, mainly driven by the conflict between strong mutation and recombination-associated fixation bias at CpG sites. We detect excess C-to-G mutations in PAR1 of all great apes, specific to the mutagenic effect of male recombination. Despite recent evidence for Y chromosome introgression from humans into Neanderthals, we find that the Neanderthal PAR1 retained similarity to the Denisovan sequence. We find differences between substitution spectra of these archaics suggesting rapid evolution of PAR1 in recent hominin history. Frequency analysis of alleles segregating in females and males provided no evidence for recent sexual antagonism in this region. We study repeat content and double-strand break hotspot regions in PAR1 and find that they may play roles in ensuring the obligate X-Y recombination event during male meiosis.

Conclusions

Our study provides an unprecedented quantification of population genetic forces governing PAR1 biology across extant and extinct hominids. PAR1 evolutionary dynamics are predominantly governed by recombination processes with a strong impact on mutation patterns across all species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02784-x.

Revised time estimation of the ancestral human chromosome 2 fusion

Background

The reduction of the chromosome number from 48 in the Great Apes to 46 in modern humans is thought to result from the end-to-end fusion of two ancestral non-human primate chromosomes forming the human chromosome 2 (HSA2). Genomic signatures of this event are the presence of inverted telomeric repeats at the HSA2 fusion site and a block of degenerate satellite sequences that mark the remnants of the ancestral centromere. It has been estimated that this fusion arose up to 4.5 million years ago (Mya).

Results

We have developed an enhanced algorithm for the detection and efficient counting of the locally over-represented weak-to-strong (AT to GC) substitutions. By analyzing the enrichment of these substitutions around the fusion site of HSA2 we estimated its formation time at 0.9 Mya with a 95% confidence interval of 0.4-1.5 Mya. Additionally, based on the statistics derived from our algorithm, we have reconstructed the evolutionary distances among the Great Apes (Hominoidea).

Conclusions

Our results shed light on the HSA2 fusion formation and provide a novel computational alternative for the estimation of the speciation chronology.

Enabling Inference for Context-Dependent Models of Mutation by Bounding the Propagation of Dependency

Although the rates at which positions in the genome mutate are known to depend not only on the nucleotide to be mutated, but also on neighboring nucleotides, it remains challenging to do phylogenetic inference using models of context-dependent mutation. In these models, the effects of one mutation may in principle propagate to faraway locations, making it difficult to compute exact likelihoods. This article shows how to use bounds on the propagation of dependency to compute likelihoods of mutation of a given segment of genome by marginalizing over sufficiently long flanking sequence. This can be used for maximum likelihood or Bayesian inference. Protocols examining residuals and iterative model refinement are also discussed. Tools for efficiently working with these models are provided in an R package, which could be used in other applications. The method is used to examine context dependence of mutations since the common ancestor of humans and chimpanzee.

The contribution of incomplete lineage sorting and introgression to the evolutionary history of the fast-evolving genus Ctenomys (Rodentia, Ctenomyidae)

Incomplete lineage sorting and introgression have been increasingly recognized as important processes involved in biological differentiation. Both incomplete lineage sorting and introgression result in incongruences between gene trees and species trees, consequently causing difficulties in phylogenetic reconstruction. This is particularly the case for rapid radiations, as short internodal distances and incomplete reproductive isolation increase the likelihood of both ILS and introgression. Estimation of the relative frequency of these processes requires assessments across many genomic regions. We use transcriptomics to test for introgression and estimate the frequency of incomplete lineage sorting in a set of three closely related and geographically adjacent South American tuco-tucos species (Ctenomys), a genus comprising 64 species resulting from recent, rapid radiation. After cleaning and filtering, 5764 orthologous genes strongly support paraphyly of C. pearsoni relative to C. brasiliensis (putatively represented by the population of Villa Serrana). In line with earlier phylogenetic work, the C. pearsoni - C. brasiliensis pair is closely related to C. torquatus, whereas C. rionegrensis is more distantly related to these three nominal species. Classical Patterson's D-statistic shows significant signals of introgression from C. torquatus into C. brasiliensis. However, a 5-taxon test shows no significant results. Incomplete lineage sorting was estimated to have involved about 9% of the loci, suggesting it represents an important process in the incipient diversification of tuco-tucos.

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.

Genome evolution and diversity of wild and cultivated potatoes

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production^1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

Population dynamics and genetic connectivity in recent chimpanzee history

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

Evolution of Placental Hormones: Implications for Animal Models

Human placenta secretes a variety of hormones, some of them in large amounts. Their effects on maternal physiology, including the immune system, are poorly understood. Not one of the protein hormones specific to human placenta occurs outside primates. Instead, laboratory and domesticated species have their own sets of placental hormones. There are nonetheless several examples of convergent evolution. Thus, horse and human have chorionic gonadotrophins with similar functions whilst pregnancy-specific glycoproteins have evolved in primates, rodents, horses, and some bats, perhaps to support invasive placentation. Placental lactogens occur in rodents and ruminants as well as primates though evolved through duplication of different genes and with functions that only partially overlap. There are also placental hormones, such as the pregnancy-associated glycoproteins of ruminants, that have no equivalent in human gestation. This review focusses on the evolution of placental hormones involved in recognition and maintenance of pregnancy, in maternal adaptations to pregnancy and lactation, and in facilitating immune tolerance of the fetal semiallograft. The contention is that knowledge gained from laboratory and domesticated mammals can translate to a better understanding of human placental endocrinology, but only if viewed in an evolutionary context.

Morphological variation of the maxilla in modern humans and African apes

Differences in morphology among modern humans and African apes are frequently used when assessing whether hominin fossils should be attributed to a single species or represent evidence for taxic diversity. A good understanding of the degree and structure of the intergeneric, interspecific, and intraspecific variation, including aspects such as sexual dimorphism and age, are key in this context. Here we explore the variation and differences shown by the maxilla of extant hominines, as maxillary morphology is central in the diagnosis of several hominin taxa. Our sample includes adults of all currently recognized hominine species and subspecies, with a balanced species sex ratio. In addition, we compared the adults with a small sample of late juveniles. The morphology of the maxillae was captured using three-dimensional landmarks, and the size and shape were analyzed using geometric morphometric methods. Key observations are that 1) the maxillae of all extant hominine species and subspecies show statistically significant differences, but complete separation in shape is only seen at the genus level; 2) the degree of variation is not consistent between genera, with subspecies of Gorilla being more different from each other than are species of Pan; 3) the pattern of sexual shape dimorphism is different in Pan, Gorilla, and Homo, often showing opposite trends; and 4) differentiation between maxillary shapes is increased after adjustment for static intraspecific allometry. These results provide a taxonomically up-to-date comparative morphological framework to help interpret the hominin fossil record, and we discuss the practical implications in that context.

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.

Incomplete lineage sorting and phenotypic evolution in marsupials

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Predicting Archaic Hominin Phenotypes from Genomic Data

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Population Genomic Sequencing Delineates Global Landscape of Copy Number Variations that Drive Domestication and Breed Formation of in Chicken

Copy number variation (CNV) is an important genetic mechanism that drives evolution and generates new phenotypic variations. To explore the impact of CNV on chicken domestication and breed shaping, the whole-genome CNVs were detected via multiple methods. Using the whole-genome sequencing data from 51 individuals, corresponding to six domestic breeds and wild red jungle fowl (RJF), we determined 19,329 duplications and 98,736 deletions, which covered 11,123 copy number variation regions (CNVRs) and 2,636 protein-coding genes. The principal component analysis (PCA) showed that these individuals could be divided into four populations according to their domestication and selection purpose. Seventy-two highly duplicated CNVRs were detected across all individuals, revealing pivotal roles of nervous system (NRG3, NCAM2), sensory (OR), and follicle development (VTG2) in chicken genome. When contrasting the CNVs of domestic breeds to those of RJFs, 235 CNVRs harboring 255 protein-coding genes, which were predominantly involved in pathways of nervous, immunity, and reproductive system development, were discovered. In breed-specific CNVRs, some valuable genes were identified, including HOXB7 for beard trait in Beijing You chicken; EDN3, SLMO2, TUBB1, and GFPT1 for melanin deposition in Silkie chicken; and SORCS2 for aggressiveness in Luxi Game fowl. Moreover, CSMD1 and NTRK3 with high duplications found exclusively in White Leghorn chicken, and POLR3H, MCM9, DOCK3, and AKR1B1L found in Recessive White Rock chicken may contribute to high egg production and fast-growing traits, respectively. The candidate genes of breed characteristics are valuable resources for further studies on phenotypic variation and the artificial breeding of chickens.

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.

Genetics Matters: Voyaging from the Past into the Future of Humanity and Sustainability

The understanding of how genetic information may be inherited through generations was established by Gregor Mendel in the 1860s when he developed the fundamental principles of inheritance. The science of genetics, however, began to flourish only during the mid-1940s when DNA was identified as the carrier of genetic information. The world has since then witnessed rapid development of genetic technologies, with the latest being genome-editing tools, which have revolutionized fields from medicine to agriculture. This review walks through the historical timeline of genetics research and deliberates how this discipline might furnish a sustainable future for humanity.

Body proportions and environmental adaptation in gorillas

OBJECTIVES

Limb length and trunk proportions are determined in a large, taxonomically and environmentally diverse sample of gorillas and related to variation in locomotion, climate, altitude, and diet.

MATERIALS AND METHODS

The sample includes 299 gorilla skeletons, 115 of which are infants and juveniles, distributed between western lowland (G. gorilla gorilla), low and high elevation grauer (G. beringei graueri), and Virunga mountain gorillas (G. b. beringei). Limb bone and vertebral column lengths scaled to body mass are compared between subgroups by age group.

RESULTS

All G. beringei have relatively short 3rd metapodials and manual proximal phalanges compared to G. gorilla, and this difference is apparent in infancy. All G. beringei also have shortened total limb lengths relative to either body mass or vertebral column length, although patterns of variation in individual skeletal elements are more complex, and infants do not display the same patterns as adults. Mountain gorillas have relatively long clavicles, present in infancy, and a relatively long thoracic (but not lumbosacral) vertebral column.

DISCUSSION

A variety of environmental factors likely contributed to observed patterns of morphological variation among extant gorillas. We interpret the short hand and foot bones of all G. beringei as genetic adaptations to greater terrestriality in the last common ancestor of G. beringei; variation in other limb lengths to climatic adaptation, both genetic and developmental; and the larger thorax of G. b. beringei to adaptation to reduced oxygen pressure at high altitudes, again as a product of both genetic differences and environmental influences during development.

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.

Identifying and mitigating the impacts on primates of transportation and service corridors

Most primate populations are declining; 60% of species face extinction. The expansion of transportation and service corridors (T&S) (i.e., roads and railways and utility and service lines) poses a significant yet underappreciated threat. With the development of T&S corridors predicted to increase across primates' ranges, it is necessary to understand the current extent of its impacts on primates, the available options to mitigate these effectively, and recognize research and knowledge gaps. By employing a systematic search approach to identify literature that described the relationship between primates and T&S corridors, we extracted information from 327 studies published between 1980 and 2020. Our results revealed that 218 species and subspecies across 62 genera are affected, significantly more than the 92 listed by the IUCN Red List of Threatened Species. The majority of studies took place in Asia (45%), followed by mainland Africa (31%), the Neotropics (22%), and Madagascar (2%). Brazil, Indonesia, Equatorial Guinea, Vietnam, and Madagascar contained the greatest number of affected primate species. Asia featured the highest number of species affected by roads, electrical transmission lines, and pipelines and the only studies addressing the impact of rail and aerial tramways on primates. The impact of seismic lines only emerged in the literature from Africa and the Neotropics. Impacts are diverse and multifaceted, for example, animal-vehicle collisions, electrocutions, habitat loss and fragmentation, impeded movement and genetic exchange, behavioral changes, exposure to pollution, and mortality associated with hunting. Although several mitigation measures were recommended, only 41% of studies focused on their implementation, whereas only 29% evaluated their effectiveness. Finally, there was a clear bias in the species and regions benefiting from research on this topic. We recommend that government and conservation bodies recognize T&S corridors as a serious and mounting threat to primates and that further research in this area is encouraged.

The complex Y-chromosomal history of gorillas

Studies of the evolutionary relationships among gorilla populations using autosomal and mitochondrial sequences suggest that male-mediated gene flow may have been important in the past, but data on the Y-chromosomal relationships among the gorilla subspecies are limited. Here, we genotyped blood and noninvasively collected fecal samples from 12 captives and 257 wild male gorillas of known origin representing all four subspecies (Gorilla gorilla gorilla, G. g. diehli, G. beringei beringei, and G. b. graueri) at 10 Y-linked microsatellite loci resulting in 102 unique Y-haplotypes for 224 individuals. We found that western lowland gorilla (G. g. gorilla) haplotypes were consistently more diverse than any other subspecies for all measures of diversity and comprised several genetically distinct groups. However, these did not correspond to geographical proximity and some closely related haplotypes were found several hundred kilometers apart. Similarly, our broad sampling of eastern gorillas revealed that mountain (G. b. beringei) and Grauer's (G. b. graueri) gorilla Y-chromosomal haplotypes did not form distinct clusters. These observations suggest structure in the ancestral population with subsequent mixing of differentiated haplotypes by male dispersal for western lowland gorillas, and postisolation migration or incomplete lineage sorting due to short divergence times for eastern gorillas.

OUP accepted manuscript

We present Champagne, a whole-genome method for generating character matrices for phylogenomic analysis using large genomic indel events. By rigorously picking orthologous genes and locating large insertion and deletion events, Champagne delivers a character matrix that considerably reduces homoplasy compared with morphological and nucleotide-based matrices, on both established phylogenies and difficult-to-resolve nodes in the mammalian tree. Champagne provides ample evidence in the form of genomic structural variation to support incomplete lineage sorting and possible introgression in Paenungulata and human–chimp–gorilla which were previously inferred primarily through matrices composed of aligned single-nucleotide characters. Champagne also offers further evidence for Myomorpha as sister to Sciuridae and Hystricomorpha in the rodent tree. Champagne harbors distinct theoretical advantages as an automated method that produces nearly homoplasy-free character matrices on the whole-genome scale.

The effects of introgression across thousands of quantitative traits revealed by gene expression in wild tomatoes

It is now understood that introgression can serve as powerful evolutionary force, providing genetic variation that can shape the course of trait evolution. Introgression also induces a shared evolutionary history that is not captured by the species phylogeny, potentially complicating evolutionary analyses that use a species tree. Such analyses are often carried out on gene expression data across species, where the measurement of thousands of trait values allows for powerful inferences while controlling for shared phylogeny. Here, we present a Brownian motion model for quantitative trait evolution under the multispecies network coalescent framework, demonstrating that introgression can generate apparently convergent patterns of evolution when averaged across thousands of quantitative traits. We test our theoretical predictions using whole-transcriptome expression data from ovules in the wild tomato genus Solanum. Examining two sub-clades that both have evidence for post-speciation introgression, but that differ substantially in its magnitude, we find patterns of evolution that are consistent with histories of introgression in both the sign and magnitude of ovule gene expression. Additionally, in the sub-clade with a higher rate of introgression, we observe a correlation between local gene tree topology and expression similarity, implicating a role for introgressed cis-regulatory variation in generating these broad-scale patterns. Our results reveal a general role for introgression in shaping patterns of variation across many thousands of quantitative traits, and provide a framework for testing for these effects using simple model-informed predictions.

It is now known from studying large genetic datasets that species often hybridize and cross with each other over many generations – a phenomenon known as introgression. Introgression introduces new genetic variation into a population, and this variation can cause traits to be shared among the introgressing species. When researchers study the evolution of trait variation among species, this source of trait sharing is rarely accounted for. Here, we present a statistical model of the effects of introgression on trait variation. This model predicts that, when averaged across many thousands of traits, introgressing species are consistently more similar than expected from standard approaches. Researchers studying gene expression often consider the expression of many thousands of genes, making this a case where the expected effects of introgression are likely to manifest. We tested our model prediction using ovule gene expression data from the wild tomato genus Solanum, in two groups of species with evidence of historical introgression. We found that patterns of expression similarity in both groups are consistent with their histories of introgression and the predictions from our model. Our results highlight the importance of accounting for introgression as a source of trait variation among species.

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.

Exploring the potential effect of COVID-19 on an endangered great ape

The current COVID-19 pandemic has created unmeasurable damages to society at a global level, from the irreplaceable loss of life, to the massive economic losses. In addition, the disease threatens further biodiversity loss. Due to their shared physiology with humans, primates, and particularly great apes, are susceptible to the disease. However, it is still uncertain how their populations would respond in case of infection. Here, we combine stochastic population and epidemiological models to simulate the range of potential effects of COVID-19 on the probability of extinction of mountain gorillas. We find that extinction is sharply driven by increases in the basic reproductive number and that the probability of extinction is greatly exacerbated if the immunity lasts less than 6 months. These results stress the need to limit exposure of the mountain gorilla population, the park personnel and visitors, as well as the potential of vaccination campaigns to extend the immunity duration.

GC-biased gene conversion in X-chromosome palindromes conserved in human, chimpanzee, and rhesus macaque

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich as the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion since the divergence of human and rhesus macaque explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that palindromes are older than 29 million years and/or had elevated GC content at the time of their formation. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.