Abstract: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.

Deciphering Gorilla gorilla gorilla immunoglobulin loci in multiple genome assemblies and enrichment of IMGT resources

Deciphering of Gorilla gorilla gorilla Immunoglobulin Loci in Multiple Genome Assemblies and Enrichment of IMGT Resources

Multiple full-length homozygous IGH haplotypes from Mauritian cynomolgus macaques

Resolving haplotype variation and complex genetic architecture in the human immunoglobulin kappa chain locus in individuals of diverse ancestry

IGLoo: Profiling the Immunoglobulin Heavy chain locus in Lymphoblastoid Cell Lines with PacBio High-Fidelity Sequencing reads

Immunoglobulin Genomics in the Guinea Pig (cavia Porcellus)

A Preliminary Analysis Of The Immunoglobulin Genes In The African Elephant (Loxodonta Africana)

Genomic Screening by 454 Pyrosequencing Identifies a New Human IGHV Gene and Sixteen Other New IGHV Allelic Variants

IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data

Unveiling the Diversity of Immunoglobulin Heavy Constant Gamma (IGHG) Gene Segments in Brazilian Populations Reveals 28 Novel Alleles and Evidence of Gene Conversion and Natural Selection

Clustering-based Identification of Clonally-Related Immunoglobulin Gene Sequence Sets

A genome-wide survey of genetic variation in gorillas using reduced representation sequencing

Novel polymorphic and copy number diversity in the antibody IGH locus of South African individuals

The Inference of Phased Haplotypes for the Immunoglobulin H Chain V Region Gene Loci by Analysis of VDJ Gene Rearrangements

Complete Genomic Assembly of Mauritian Cynomolgus Macaque Killer Ig-like Receptor and Natural Killer Group 2 Haplotypes

Immunoglobulin Genomics in the Prairie Vole (microtus Ochrogaster).

Physical mapping of the giant panda immunoglobulin heavy chain constant region genes.

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

Network Inference from Glycoproteomics Data Reveals New Reactions in the IgG Glycosylation Pathway

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

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