ABSTRACT
Predicting phenotypes from genomic data is a key goal in genetics, but for most complex phenotypes, predictions are hampered by incomplete genotype-to-phenotype mapping. Here, we describe a more attainable approach than quantitative predictions, which is aimed at qualitatively predicting phenotypic differences. Despite incomplete genotype-to-phenotype mapping, we show that it is relatively easy to determine which of two individuals has a greater phenotypic value. This question is central in many scenarios, e.g., comparing disease risk between individuals, the yield of crop strains, or the anatomy of extinct vs extant species. To evaluate prediction accuracy, i.e., the probability that the individual with the greater predicted phenotype indeed has a greater phenotypic value, we developed an estimator of the ratio between known and unknown effects on the phenotype. We evaluated prediction accuracy using human data from tens of thousands of individuals from either the same family or the same population, as well as data from different species. We found that, in many cases, even when only a small fraction of the loci affecting a phenotype is known, the individual with the greater phenotypic value can be identified with over 90% accuracy. Our approach also circumvents some of the limitations in transferring genetic association results across populations. Overall, we introduce an approach that enables accurate predictions of key information on phenotypes - the direction of phenotypic difference - and suggest that more phenotypic information can be extracted from genomic data than previously appreciated.
ABSTRACT
Measuring allele-specific expression in interspecies hybrids is a powerful way to detect cis-regulatory changes underlying adaptation. However, it remains difficult to identify genes most likely to explain species-specific traits. Here, we outline a simple strategy that leverages population-scale allele-specific RNA-seq data to identify genes that show constrained cis-regulation within species yet show divergence between species. Applying this strategy to data from human-chimpanzee hybrid cortical organoids, we identify signatures of lineage-specific selection on genes related to saccharide metabolism, neurodegeneration, and primary cilia. We also highlight cis-regulatory divergence in CUX1 and EDNRB that may shape the trajectory of human brain development.
Subject(s)
Evolution, Molecular , Organoids , Humans , Alleles , Pan troglodytes , AnimalsABSTRACT
Mucus protects the epithelial cells of the digestive and respiratory tracts from pathogens and other hazards. Progress in determining the molecular mechanisms of mucus barrier function has been limited by the lack of high-resolution structural information on mucins, the giant, secreted, gel-forming glycoproteins that are the major constituents of mucus. Here, we report how mucin structures we determined enabled the discovery of an unanticipated protective role of mucus: managing the toxic transition metal copper. Using two juxtaposed copper binding sites, one for Cu2+ and the other for Cu1+, the intestinal mucin, MUC2, prevents copper toxicity by blocking futile redox cycling and the squandering of dietary antioxidants, while nevertheless permitting uptake of this important trace metal into cells. These findings emphasize the value of molecular structure in advancing mucosal biology, while introducing mucins, produced in massive quantities to guard extensive mucosal surfaces, as extracellular copper chaperones.
Subject(s)
Copper , Mucins , Mucins/metabolism , Mucin-2 , Copper/analysis , Copper/metabolism , Intestines , Mucus/metabolism , Intestinal Mucosa/metabolismABSTRACT
The Neanderthal and Denisovan genomes enabled the discovery of sequences that differ between modern and archaic humans, the majority of which are noncoding. However, our understanding of the regulatory consequences of these differences remains limited, in part due to the decay of regulatory marks in ancient samples. Here, we used a massively parallel reporter assay in embryonic stem cells, neural progenitor cells, and bone osteoblasts to investigate the regulatory effects of the 14,042 single-nucleotide modern human-specific variants. Overall, 1791 (13%) of sequences containing these variants showed active regulatory activity, and 407 (23%) of these drove differential expression between human groups. Differentially active sequences were associated with divergent transcription factor binding motifs, and with genes enriched for vocal tract and brain anatomy and function. This work provides insight into the regulatory function of variants that emerged along the modern human lineage and the recent evolution of human gene expression.
Subject(s)
Embryonic Stem Cells/metabolism , Gene Expression Regulation , Osteoblasts/metabolism , Polymorphism, Single Nucleotide , Stem Cells/metabolism , Transcription Factors/metabolism , Genome, Human , HumansABSTRACT
Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for nongenetic confounding factors. We differentiated these cells into cranial neural crest cells, the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific sixfold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.
Subject(s)
Chimera/genetics , Ellis-Van Creveld Syndrome/genetics , Intercellular Signaling Peptides and Proteins/genetics , Neural Crest/metabolism , Pan troglodytes/genetics , Skull/metabolism , Animals , Biological Evolution , Cell Differentiation , Chimera/metabolism , Ellis-Van Creveld Syndrome/metabolism , Ellis-Van Creveld Syndrome/pathology , Female , Gene Expression , Genotype , Humans , Induced Pluripotent Stem Cells/cytology , Induced Pluripotent Stem Cells/metabolism , Intercellular Signaling Peptides and Proteins/deficiency , Male , Mice , Mice, Knockout , Neural Crest/pathology , Pan troglodytes/anatomy & histology , Pan troglodytes/metabolism , Phenotype , Signal Transduction , Skull/anatomy & histology , Species Specificity , TetraploidyABSTRACT
We report genome-wide DNA data for 73 individuals from five archaeological sites across the Bronze and Iron Ages Southern Levant. These individuals, who share the "Canaanite" material culture, can be modeled as descending from two sources: (1) earlier local Neolithic populations and (2) populations related to the Chalcolithic Zagros or the Bronze Age Caucasus. The non-local contribution increased over time, as evinced by three outliers who can be modeled as descendants of recent migrants. We show evidence that different "Canaanite" groups genetically resemble each other more than other populations. We find that Levant-related modern populations typically have substantial ancestry coming from populations related to the Chalcolithic Zagros and the Bronze Age Southern Levant. These groups also harbor ancestry from sources we cannot fully model with the available data, highlighting the critical role of post-Bronze-Age migrations into the region over the past 3,000 years.
Subject(s)
DNA, Ancient/analysis , Ethnicity/genetics , Gene Flow/genetics , Archaeology/methods , DNA, Mitochondrial/genetics , Ethnicity/history , Gene Flow/physiology , Genetic Variation/genetics , Genetics, Population/methods , Genome, Human/genetics , Genomics/methods , Haplotypes , History, Ancient , Human Migration/history , Humans , Mediterranean Region , Middle East , Sequence Analysis, DNAABSTRACT
Changes in potential regulatory elements are thought to be key drivers of phenotypic divergence. However, identifying changes to regulatory elements that underlie human-specific traits has proven very challenging. Here, we use 63 reconstructed and experimentally measured DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that likely emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes associated with face and vocal tract anatomy went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-associated genes (SOX9, ACAN, COL2A1, NFIX and XYLT1). We propose that these repression patterns appeared after the split from Neanderthals and Denisovans, and that they might have played a key role in shaping the modern human face and vocal tract.
Subject(s)
DNA Methylation , DNA, Ancient , Face/anatomy & histology , Phenotype , Phonation/genetics , Adult , Aged , Animals , Cells, Cultured , Child , Chondrocytes , Evolution, Molecular , Female , Gene Regulatory Networks , Genetic Speciation , Humans , Larynx/anatomy & histology , Male , Middle Aged , Neanderthals/genetics , Pan troglodytes/genetics , Primary Cell Culture , Tongue/anatomy & histology , Vocal Cords/anatomy & histology , Vocalization, AnimalABSTRACT
Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.
Subject(s)
DNA Methylation/genetics , Neanderthals/anatomy & histology , Neanderthals/genetics , Pan troglodytes/anatomy & histology , Pan troglodytes/genetics , Phenotype , Animals , Base Sequence , Databases, Genetic , Extinction, Biological , Fossils , Genome, Human/genetics , Humans , Polymorphism, Single Nucleotide/genetics , Skeleton , SkullABSTRACT
Analyzing the conditions in which past individuals lived is key to understanding the environments and cultural transitions to which humans had to adapt. Here, we suggest a methodology to probe into past environments, using reconstructed premortem DNA methylation maps of ancient individuals. We review a large body of research showing that differential DNA methylation is associated with changes in various external and internal factors, and propose that loci whose DNA methylation level is environmentally responsive could serve as markers to infer about ancient daily life, diseases, nutrition, exposure to toxins, and more. We demonstrate this approach by showing that hunger-related DNA methylation changes are found in ancient hunter-gatherers. The strategy we present here opens a window to reconstruct previously inaccessible aspects of the lives of past individuals.
Subject(s)
DNA, Ancient/analysis , Epigenomics/methods , Animals , DNA/genetics , DNA Methylation/genetics , DNA Methylation/physiology , Environment , Epigenesis, Genetic/genetics , Hominidae/genetics , HumansABSTRACT
One of the biggest challenges in studying how genes work is understanding their effect on the physiology and anatomy of the body. Existing tools try to address this using indirect features, such as expression levels and biochemical pathways. Here, we present Gene ORGANizer (geneorganizer.huji.ac.il), a phenotype-based tool that directly links human genes to the body parts they affect. It is built upon an exhaustive curated database that links >7000 genes to â¼150 anatomical parts using >150 000 gene-organ associations. The tool offers user-friendly platforms to analyze the anatomical effects of individual genes, and identify trends within groups of genes. We demonstrate how Gene ORGANizer can be used to make new discoveries, showing that chromosome X is enriched with genes affecting facial features, that positive selection targets genes with more constrained phenotypic effects, and more. We expect Gene ORGANizer to be useful in a variety of evolutionary, medical and molecular studies aimed at understanding the phenotypic effects of genes.
Subject(s)
Databases, Genetic , Genetic Association Studies , Genome, Human , Software , Chromosomes, Human, X/chemistry , Female , Genomic Imprinting , Genotype , Humans , Male , Metabolic Networks and Pathways/genetics , Organ Specificity/genetics , PhenotypeABSTRACT
A recent study conducted the first genome-wide scan for selection in Inuit from Greenland using single nucleotide polymorphism chip data. Here, we report that selection in the region with the second most extreme signal of positive selection in Greenlandic Inuit favored a deeply divergent haplotype that is closely related to the sequence in the Denisovan genome, and was likely introgressed from an archaic population. The region contains two genes, WARS2 and TBX15, and has previously been associated with adipose tissue differentiation and body-fat distribution in humans. We show that the adaptively introgressed allele has been under selection in a much larger geographic region than just Greenland. Furthermore, it is associated with changes in expression of WARS2 and TBX15 in multiple tissues including the adrenal gland and subcutaneous adipose tissue, and with regional DNA methylation changes in TBX15.
Subject(s)
Adaptation, Biological/genetics , Inuit/genetics , T-Box Domain Proteins/genetics , Adipose Tissue/physiology , Alleles , Animals , DNA Methylation , DNA, Ancient , Greenland , Haplotypes , Humans , Models, Genetic , Neanderthals , Polymorphism, Single Nucleotide , Selection, Genetic , Sequence Analysis, DNA/methodsABSTRACT
Pluripotent self-renewing embryonic stem cells (ESCs) have been the focus of a growing number of high-throughput experiments, revealing the genome-wide locations of hundreds of transcription factors and histone modifications. While most of these datasets were used in a specific context, all datasets combined offer a comprehensive view of chromatin characteristics and regulatory elements that govern cell states. Here, using hundreds of datasets in ESCs, we generated colocalization maps of chromatin proteins and modifications, and built a discovery pipeline for regulatory proteins of gene families. By comparing genome-wide binding data with over-expression and knockdown analysis of hundreds of genes, we discovered that the pluripotency-related factor NR5A2 separates mitochondrial from cytosolic ribosomal genes, regulating their expression. We further show that genes with a common chromatin profile are enriched for distinct Gene Ontology (GO) categories. Our approach can be generalized to reveal common regulators of any gene group; discover novel gene families, and identify common genomic elements based on shared chromatin features.
Subject(s)
Chromatin/genetics , Human Embryonic Stem Cells/cytology , Mouse Embryonic Stem Cells/cytology , Receptors, Cytoplasmic and Nuclear/metabolism , Animals , Binding Sites , Cell Differentiation/genetics , Histones/genetics , Humans , Mice , Mitochondrial Ribosomes/metabolism , Ribosomal Proteins/genetics , Transcription Factors/metabolismABSTRACT
Recent years have witnessed the rise of ancient DNA (aDNA) technology, allowing comparative genomics to be carried out at unprecedented time resolution. While it is relatively straightforward to use aDNA to identify recent genomic changes, it is much less clear how to utilize it to study changes in epigenetic regulation. Here we review recent works demonstrating that highly degraded aDNA still contains sufficient information to allow reconstruction of epigenetic signals, including DNA methylation and nucleosome positioning maps. We discuss challenges arising from the tissue specificity of epigenetics, and show how some of them might in fact turn into advantages. Finally, we introduce a method to infer methylation states in tissues that do not tend to be preserved over time.
Subject(s)
Epigenesis, Genetic , Hominidae/genetics , Animals , DNA/genetics , DNA MethylationABSTRACT
Ancient DNA sequencing has recently provided high-coverage archaic human genomes. However, the evolution of epigenetic regulation along the human lineage remains largely unexplored. We reconstructed the full DNA methylation maps of the Neandertal and the Denisovan by harnessing the natural degradation processes of methylated and unmethylated cytosines. Comparing these ancient methylation maps to those of present-day humans, we identified ~2000 differentially methylated regions (DMRs). Particularly, we found substantial methylation changes in the HOXD cluster that may explain anatomical differences between archaic and present-day humans. Additionally, we found that DMRs are significantly more likely to be associated with diseases. This study provides insight into the epigenetic landscape of our closest evolutionary relatives and opens a window to explore the epigenomes of extinct species.
Subject(s)
DNA Methylation , Epigenesis, Genetic , Evolution, Molecular , Genome, Human , Neanderthals/genetics , Animals , HumansABSTRACT
Histones, the building blocks of eukaryotic chromatin, are essential for genome packaging, function and regulation. However, little is known about their transcriptional regulation. Here we conducted a comprehensive computational analysis, based on chromatin immunoprecipitation-sequencing and -microarray analysis (ChIP-seq and ChIP-chip) data of over 50 transcription factors and histone modifications in mouse embryonic stem cells. Enrichment scores supported by gene expression data from gene knockout studies identified E2f1 and E2f4 as master regulators of histone genes, CTCF and Zfx as repressors of core and linker histones, respectively, and Smad1, Smad2, YY1 and Ep300 as restricted or cell type-specific regulators. We propose that histone gene regulation is substantially more complex than previously thought, and that a combination of factors orchestrate histone gene regulation, from strict synchronization with S phase to targeted regulation of specific histone subtypes.