Omer Gokcumen

The evolution and functional impact of human deletion variants shared with archaic hominin genomes

Allele sharing between modern and archaic hominin genomes has been variously interpreted to have originated from ancestral genetic structure or through non-African introgression from archaic hominins. However, evolution of polymorphic human deletions that are shared with archaic hominin genomes has yet to be studied. We identified 427 polymorphic human deletions that are shared with archaic hominin genomes, approximately 87% of which originated before the Human–Neandertal divergence (ancient) and only approximately 9% of which have been introgressed from Neandertals (introgressed). Recurrence, incomplete lineage sorting between human and chimp lineages, and hominid-specific insertions constitute the remaining approximately 4% of allele sharing between humans and archaic hominins. We observed that ancient deletions correspond to more than 13% of all common (>5% allele frequency) deletion variation among modern humans. Our analyses indicate that the genomic landscapes of both ancient and introgressed deletion variants were primarily shaped by purifying selection, eliminating large and exonic variants. We found 17 exonic deletions that are shared with archaic hominin genomes, including those leading to three fusion transcripts. The affected genes are involved in metabolism of external and internal compounds, growth and sperm formation, as well as susceptibility to psoriasis and Crohn’s disease. Our analyses suggest that these “exonic” deletion variants have evolved through different adaptive forces, including balancing and population-specific positive selection. Our findings reveal that genomic structural variants that are shared between humans and archaic hominin genomes are common among modern humans and can influence biomedically and evolutionarily important phenotypes.

Whole genome sequencing of Turkish genomes reveals functional private alleles and impact of genetic interactions with Europe, Asia and Africa.

BackgroundTurkey is a crossroads of major population movements throughout history and has been a hotspot of cultural interactions. Several studies have investigated the complex population history of Turkey through a limited set of genetic markers. However, to date, there have been no studies to assess the genetic variation at the whole genome level using whole genome sequencing. Here, we present whole genome sequences of 16 Turkish individuals resequenced at high coverage (32u2009×u2009-48×).ResultsWe show that the genetic variation of the contemporary Turkish population clusters with South European populations, as expected, but also shows signatures of relatively recent contribution from ancestral East Asian populations. In addition, we document a significant enrichment of non-synonymous private alleles, consistent with recent observations in European populations. A number of variants associated with skin color and total cholesterol levels show frequency differentiation between the Turkish populations and European populations. Furthermore, we have analyzed the 17q21.31 inversion polymorphism region (MAPT locus) and found increased allele frequency of 31.25% for H1/H2 inversion polymorphism when compared to European populations that show about 25% of allele frequency.ConclusionThis study provides the first map of common genetic variation from 16 western Asian individuals and thus helps fill an important geographical gap in analyzing natural human variation and human migration. Our data will help develop population-specific experimental designs for studies investigating disease associations and demographic history in Turkey.

Archaic hominin introgression in Africa contributes to functional salivary MUC7 genetic variation

Abstract One of the most abundant proteins in human saliva, mucin-7, is encoded by the MUC7 gene, which harbors copy number variable subexonic repeats (PTS-repeats) that affect the size and glycosylation potential of this protein. We recently documented the adaptive evolution of MUC7 subexonic copy number variation among primates. Yet, the evolution of MUC7 genetic variation in humans remained unexplored. Here, we found that PTS-repeat copy number variation has evolved recurrently in the human lineage, thereby generating multiple haplotypic backgrounds carrying five or six PTS-repeat copy number alleles. Contrary to previous studies, we found no associations between the copy number of PTS-repeats and protection against asthma. Instead, we revealed a significant association of MUC7 haplotypic variation with the composition of the oral microbiome. Furthermore, based on in-depth simulations, we conclude that a divergent MUC7 haplotype likely originated in an unknown African hominin population and introgressed into ancestors of modern Africans.

Atopic Dermatitis Susceptibility Variants in Filaggrin Hitchhike Hornerin Selective Sweep

Human skin has evolved rapidly, leaving evolutionary signatures in the genome. The filaggrin (FLG) gene is widely studied for its skin-barrier function in humans. The extensive genetic variation in this gene, especially common loss-of-function (LoF) mutations, has been established as primary risk factors for atopic dermatitis. To investigate the evolution of this gene, we analyzed 2,504 human genomes and genotyped the copy number variation of filaggrin repeats within FLG in 126 individuals from diverse ancestral backgrounds. We were unable to replicate a recent study claiming that LoF of FLG is adaptive in northern latitudes with lower ultraviolet light exposure. Instead, we present multiple lines of evidence suggesting that FLG genetic variation, including LoF variants, have little or no effect on fitness in modern humans. Haplotype-level scrutinization of the locus revealed signatures of a recent selective sweep in Asia, which increased the allele frequency of a haplotype group (Huxian haplogroup) in Asian populations. Functionally, we found that the Huxian haplogroup carries dozens of functional variants in FLG and hornerin (HRNR) genes, including those that are associated with atopic dermatitis susceptibility, HRNR expression levels and microbiome diversity on the skin. Our results suggest that the target of the adaptive sweep is HRNR gene function, and the functional FLG variants that involve susceptibility to atopic dermatitis, seem to hitchhike the selective sweep on HRNR. Our study presents a novel case of a locus that harbors clinically relevant common genetic variation with complex evolutionary trajectories.

The psoriasis-associated deletion of late cornified envelope genes LCE3B and LCE3C has been maintained under balancing selection since Human Denisovan divergence

BackgroundA common, 32kb deletion of LCE3B and LCE3C genes is strongly associated with psoriasis. We recently found that this deletion is ancient, predating Human-Denisovan divergence. However, it was not clear why negative selection has not removed this deletion from the population.ResultsHere, we show that the haplotype block that harbors the deletion (i) retains high allele frequency among extant and ancient human populations; (ii) harbors unusually high nucleotide variation (π, Pu2009<u20094.1u2009×u200910−3); (iii) contains an excess of intermediate frequency variants (Tajima’s D, Pu2009<u20093.9u2009×u200910−3); and (iv) has an unusually long time to coalescence to the most recent common ancestor (TSel, 0.1 quantile).ConclusionsOur results are most parsimonious with the scenario where the LCE3BC deletion has evolved under balancing selection in humans. More broadly, this is consistent with the hypothesis that a balance between autoimmunity and natural vaccination through increased exposure to pathogens maintains this deletion in humans.

Global Survey of Variation in a Human Olfactory Receptor Gene Reveals Signatures of Non-Neutral Evolution

Allelic variation at 4 loci in the human olfactory receptor gene OR7D4 is associated with perceptual variation in the sex steroid-derived odorants, androstenone, and androstadienone. Androstadienone has been linked with chemosensory identification whereas androstenone makes pork from uncastrated pigs distasteful (boar taint). In a sample of 2224 individuals from 43 populations, we identified 45 OR7D4 single nucleotide polymorphisms. Coalescent modeling of frequency-site-spectrum-based statistics identified significant deviation from neutrality in human OR7D4; individual populations with statistically significant deviations from neutrality include Gujarati, Beijing Han, Great Britain, Iberia, and Puerto Rico. Analysis of molecular variation values indicated statistically significant population differentiation driven mainly by the 4 alleles associated with androstenone perception variation; however, fixation values were low suggesting that genetic structure may not have played a strong role in creating these group divisions. We also studied OR7D4 in the genomes of extinct members of the human lineage: Altai Neandertal and Denisovan. No variants were identified in Altai but 2 were in Denisova, one of which is shared by modern humans and one of which is novel. A functional test of modern human and a synthesized mutant Denisova OR7D4 indicated no statistically significant difference in responses to androstenone between the 2 species. Our results suggest non-neutral evolution for an olfactory receptor gene.

Geographic Distribution and Adaptive Significance of Genomic Structural Variants: An Anthropological Genetics Perspective

ABSTRACT Anthropological geneticists have successfully used single-nucleotide and short tandem repeat variations across human genomes to reconstruct human history. These markers have also been used extensively to identify adaptive and phenotypic variation. The recent advent of high-throughput genomic technologies revealed an overlooked type of genomic variation: structural variants (SVs). In fact, some SVs may contribute to human adaptation in substantial and previously unexplored ways. SVs include deletions, insertions, duplications, inversions, and translocations of genomic segments that vary among individuals from the same species. SVs are much less numerous than single-nucleotide variants but account for at least seven times more variable base pairs than do single-nucleotide variants when two human genomes are compared. Moreover, recent studies have shown that SVs have higher mutation rates than singlenucleotide variants when the affected base pairs are considered, especially in certain parts of the genome. The null hypothesis for the evolution of SVs, as for single-nucleotide variants, is neutrality. Hence, drift is the primary force that shapes the current allelic distribution of most SVs. However, due to their size, a larger proportion of SVs appear to evolve under nonneutral forces (mostly purifying selection) than do single-nucleotide variants. In fact, as exemplified by several groundbreaking studies, SVs contribute to anthropologically relevant phenotypic variation and local adaptation among humans. In this review, we argue that with the advent of affordable genomic technologies, anthropological scrutiny of genomic structural variation emerges as a fertile area of inquiry to better understand human phenotypic variation. To motivate potential studies, we discuss scenarios through which structural variants (SVs) affect phenotypic variation among humans within an anthropological context. We further provide a methodological workflow in which we analyzed 1000 Genomes deletion variants and identified 16 exonic deletions that are specific to the African continent. We analyzed two of these deletion variants affecting the keratin-associated protein (KAP) cluster in a locus-specific manner. Our analysis revealed that these deletions may indeed affect phenotype and likely evolved under geography-specific positive selection. We outline all the major software and data sets for these analyses and provide the basic R and Perl codes we used for this example workflow analysis. Overall, we hope that this review will encourage and facilitate incorporation of genomic structural variation in anthropological research programs.

Variation and Functional Impact of Neanderthal Ancestry in Western Asia

Abstract Neanderthals contributed genetic material to modern humans via multiple admixture events. Initial admixture events presumably occurred in Western Asia shortly after humans migrated out of Africa. Despite being a focal point of admixture, earlier studies indicate lower Neanderthal introgression rates in some Western Asian populations as compared with other Eurasian populations. To better understand the genome-wide and phenotypic impact of Neanderthal introgression in the region, we sequenced whole genomes of nine present-day Europeans, Africans, and the Western Asian Druze at high depth, and analyzed available whole genome data from various other populations, including 16 genomes from present-day Turkey. Our results confirmed previous observations that contemporary Western Asian populations, on an average, have lower levels of Neanderthal-introgressed DNA relative to other Eurasian populations. Modern Western Asians also show comparatively high variability in Neanderthal ancestry, which may be attributed to the complex demographic history of the region. We further replicated the previously described depletion of putatively functional sequences among Neanderthal-introgressed haplotypes. Still, we find dozens of common Neanderthal-introgressed haplotypes in the Turkish sample associated with human phenotypes, including anthropometric and metabolic traits, as well as the immune response. One of these haplotypes is unusually long and harbors variants that affect the expression of members of the CCR gene family and are associated with celiac disease. Overall, our results paint a complex first picture of the genomic impact of Neanderthal introgression in the Western Asian populations.

Next-generation sequencing-based detection of germline L1-mediated transductions

BackgroundWhile active LINE-1 (L1) elements possess the ability to mobilize flanking sequences to different genomic loci through a process termed transduction influencing genomic content and structure, an approach for detecting polymorphic germline non-reference transductions in massively-parallel sequencing data has been lacking.ResultsHere we present the computational approach TIGER (Transduction Inference in GERmline genomes), enabling the discovery of non-reference L1-mediated transductions by combining L1 discovery with detection of unique insertion sequences and detailed characterization of insertion sites. We employed TIGER to characterize polymorphic transductions in fifteen genomes from non-human primate species (chimpanzee, orangutan and rhesus macaque), as well as in a human genome. We achieved high accuracy as confirmed by PCR and two single molecule DNA sequencing techniques, and uncovered differences in relative rates of transduction between primate species.ConclusionsBy enabling detection of polymorphic transductions, TIGER makes this form of relevant structural variation amenable for population and personal genome analysis.

VCFtoTree: a user-friendly tool to construct locus-specific alignments and phylogenies from thousands of anthropologically relevant genome sequences

BackgroundConstructing alignments and phylogenies for a given locus from large genome sequencing studies with relevant outgroups allow novel evolutionary and anthropological insights. However, no user-friendly tool has been developed to integrate thousands of recently available and anthropologically relevant genome sequences to construct complete sequence alignments and phylogenies.ResultsHere, we provide VCFtoTree, a user friendly tool with a graphical user interface that directly accesses online databases to download, parse and analyze genome variation data for regions of interest. Our pipeline combines popular sequence datasets and tree building algorithms with custom data parsing to generate accurate alignments and phylogenies using all the individuals from the 1000 Genomes Project, Neanderthal and Denisovan genomes, as well as reference genomes of Chimpanzee and Rhesus Macaque. It can also be applied to other phased human genomes, as well as genomes from other species. The output of our pipeline includes an alignment in FASTA format and a tree file in newick format.ConclusionVCFtoTree fulfills the increasing demand for constructing alignments and phylogenies for a given loci from thousands of available genomes. Our software provides a user friendly interface for a wider audience without prerequisite knowledge in programming. VCFtoTree can be accessed from https://github.com/duoduoo/VCFtoTree_3.0.0.