Reconstructing DNA methylation maps of ancient populations.

Studying premortem DNA methylation from ancient DNA (aDNA) provides a proxy for ancient gene activity patterns, and hence valuable information on evolutionary changes in gene regulation. Due to statistical limitations, current methods to reconstruct aDNA methylation maps are constrained to high-coverage shotgun samples, which comprise a small minority of available ancient samples. Most samples are sequenced using in-situ hybridization capture sequencing which targets a predefined set of genomic positions. Here, we develop methods to reconstruct aDNA methylation maps of samples that were not sequenced using high-coverage shotgun sequencing, by way of pooling together individuals to obtain a DNA methylation map that is characteristic of a population. We show that the resulting DNA methylation maps capture meaningful biological information and allow for the detection of differential methylation across populations. We offer guidelines on how to carry out comparative studies involving ancient populations, and how to control the rate of falsely discovered differentially methylated regions. The ability to reconstruct DNA methylation maps of past populations allows for the development of a whole new frontier in paleoepigenetic research, tracing DNA methylation changes throughout human history, using data from thousands of ancient samples.

Lead in Archeological Human Bones Reflecting Historical Changes in Lead Production.

Forty years ago, in a seminal paper published in Science, Settle and Patterson used archeological and historical data to estimate the rates of worldwide lead production since the discovery of cupellation, approximately 5000 years ago. Here, we record actual lead exposure of a human population by direct measurements of the concentrations of lead in petrous bones of individuals representing approximately 12 000 years of inhabitation in Italy. This documentation of lead pollution throughout human history indicates that, remarkably, much of the estimated dynamics in lead production is replicated in human exposure. Thus, lead pollution in humans has closely followed anthropogenic lead production. This observation raises concerns that the forecasted increase in the production of lead and other metals might affect human health in the near future.

SRCP: a comprehensive pipeline for accurate annotation and quantification of circRNAs.

Here we describe a new integrative approach for accurate annotation and quantification of circRNAs named Short Read circRNA Pipeline (SRCP). Our strategy involves two steps: annotation of validated circRNAs followed by a quantification step. We show that SRCP is more sensitive than other individual pipelines and allows for more comprehensive quantification of a larger number of differentially expressed circRNAs. To facilitate the use of SRCP, we generate a comprehensive collection of validated circRNAs in five different organisms, including humans. We then utilize our approach and identify a subset of circRNAs bound to the miRNA-effector protein AGO2 in human brain samples.

LINADMIX: evaluating the effect of ancient admixture events on modern populations.

MOTIVATION: The rise in the number of genotyped ancient individuals provides an opportunity to estimate population admixture models for many populations. However, in models describing modern populations as mixtures of ancient ones, it is typically difficult to estimate the model mixing coefficients and to evaluate its fit to the data. RESULTS: We present LINADMIX, designed to tackle this problem by solving a constrained linear model when both the ancient and the modern genotypes are represented in a low-dimensional space. LINADMIX estimates the mixing coefficients and their standard errors, and computes a P-value for testing the model fit to the data. We quantified the performance of LINADMIX using an extensive set of simulated studies. We show that LINADMIX can accurately estimate admixture coefficients, and is robust to factors such as population size, genetic drift, proportion of missing data and various types of model misspecification. AVAILABILITY AND IMPLEMENTATION: LINADMIX is available as a python code at https://github.com/swidler/linadmix. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Antisense oligonucleotide-based drug development for Cystic Fibrosis patients carrying the 3849+10 kb C-to-T splicing mutation.

BACKGROUND: Antisense oligonucleotide (ASO)-based drugs for splicing modulation were recently approved for various genetic diseases with unmet need. Here we aimed to develop an ASO-based splicing modulation therapy for Cystic Fibrosis (CF) patients carrying the 3849+10 kb C-to-T splicing mutation in the CFTR gene. METHODS: We have screened, in FRT cells expressing the 3849+10 kb C-to-T splicing mutation, ~30 2'-O-Methyl-modified phosphorothioate ASOs, targeted to prevent the recognition and inclusion of a cryptic exon generated due to the mutation. The effect of highly potent ASO candidates on the splicing pattern, protein maturation and CFTR function was further analyzed in well differentiated primary human nasal and bronchial epithelial cells, derived from patients carrying at least one 3849+10 kb C-to-T allele. RESULTS: A highly potent lead ASO, efficiently delivered by free uptake, was able to significantly increase the level of correctly spliced mRNA and completely restore the CFTR function to wild type levels in cells from a homozygote patient. This ASO led to CFTR function with an average of 43% of wild type levels in cells from various heterozygote patients. Optimized efficiency of the lead ASO was further obtained with 2'-Methoxy Ethyl modification (2'MOE). CONCLUSION: The highly efficient splicing modulation and functional correction, achieved by free uptake of the selected lead ASO in various patients, demonstrate the ASO therapeutic potential benefit for CF patients carrying splicing mutations and is aimed to serve as the basis for our current clinical development.

Harnessing epigenetics to study human evolution.

Recent advances in ancient DNA extraction and high-throughput sequencing technologies enabled the high-quality sequencing of archaic genomes, including the Neanderthal and the Denisovan. While comparisons with modern humans revealed both archaic-specific and human-specific sequence changes, in the absence of gene expression information, understanding the functional implications of such genetic variations remains a major challenge. To study gene regulation in archaic humans, epigenetic research comes to our aid. DNA methylation, which is highly correlated with transcription, can be directly measured in modern samples, as well as reconstructed in ancient samples. This puts DNA methylation as a natural basis for comparative epigenetics between modern humans, archaic humans and nonhuman primates.

The Genomic History of the Bronze Age Southern Levant.

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.

Differential DNA methylation of vocal and facial anatomy genes in modern humans.

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.

Nucleotide composition affects codon usage toward the 3'-end.

The 3'-end of the coding sequence in several species is known to show specific codon usage bias. Several factors have been suggested to underlie this phenomenon, including selection against translation efficiency, selection for translation accuracy, and selection against RNA folding. All are supported by some evidence, but there is no general agreement as to which factors are the main determinants. Nor is it known how universal this phenomenon is, and whether the same factors explain it in different species. To answer these questions, we developed a measure that quantifies the codon usage bias at the gene end, and used it to compute this bias for 91 species that span the three domains of life. In addition, we characterized the codons in each species by features that allow discrimination between the different factors. Combining all these data, we were able to show that there is a universal trend to favor AT-rich codons toward the gene end. Moreover, we suggest that this trend is explained by avoidance from forming RNA secondary structures around the stop codon, which may interfere with normal translation termination.

Predicted Archaic 3D Genome Organization Reveals Genes Related to Head and Spinal Cord Separating Modern from Archaic Humans.

High coverage sequences of archaic humans enabled the reconstruction of their DNA methylation patterns. This allowed comparing gene regulation between human groups, and linking such regulatory changes to phenotypic differences. In a previous work, a detailed comparison of DNA methylation in modern humans, archaic humans, and chimpanzees revealed 873 modern human-derived differentially methylated regions (DMRs). To understand the regulatory implications of these DMRs, we defined differentially methylated genes (DMGs) as genes that harbor DMRs in their promoter or gene body. While most of the modern human-derived DMRs could be linked to DMGs, many others remained unassigned. Here, we used information on 3D genome organization to link ~70 out of the remaining 288 unassigned DMRs to genes. Combined with the previously identified DMGs, we reinforce the enrichment of these genes with vocal and facial anatomy, and additionally find significant enrichment with the spinal column, chin, hair, and scalp. These results reveal the importance of 3D genomic organization in understanding gene regulation by DNA methylation.

Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

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.

Assessing predictions of the impact of variants on splicing in CAGI5.

Precision medicine and sequence-based clinical diagnostics seek to predict disease risk or to identify causative variants from sequencing data. The Critical Assessment of Genome Interpretation (CAGI) is a community experiment consisting of genotype-phenotype prediction challenges; participants build models, undergo assessment, and share key findings. In the past, few CAGI challenges have addressed the impact of sequence variants on splicing. In CAGI5, two challenges (Vex-seq and MaPSY) involved prediction of the effect of variants, primarily single-nucleotide changes, on splicing. Although there are significant differences between these two challenges, both involved prediction of results from high-throughput exon inclusion assays. Here, we discuss the methods used to predict the impact of these variants on splicing, their performance, strengths, and weaknesses, and prospects for predicting the impact of sequence variation on splicing and disease phenotypes.

Small RNA sequences derived from pre-microRNAs in the supraspliceosome.

MicroRNAs (miRNAs) are short non-coding RNAs that negatively regulate the expression and translation of genes in healthy and diseased tissues. Herein, we characterize short RNAs from human HeLa cells found in the supraspliceosome, a nuclear dynamic machine in which pre-mRNA processing occurs. We sequenced small RNAs (<200 nt) extracted from the supraspliceosome, and identified sequences that are derived from 200 miRNAs genes. About three quarters of them are mature miRNAs, whereas the rest account for various defined regions of the pre-miRNA, and its hairpin-loop precursor. Out of these aligned sequences, 53 were undetected in cellular extract, and the abundance of additional 48 strongly differed from that in cellular extract. Notably, we describe seven abundant miRNA-derived sequences that overlap non-coding exons of their host gene. The rich collection of sequences identical to pre-miRNAs at the supraspliceosome suggests overlooked nuclear functions. Specifically, the abundant hsa-mir-99b may affect splicing of LINC01129 primary transcript through base-pairing with its exon-intron junction. Using suppression and overexpression experiments, we show that hsa-mir-7704 negatively regulates the level of the lncRNA HAGLR. We claim that in cases of extended base-pairing complementarity, such supraspliceosomal pre-miRNA sequences might have a role in transcription attenuation, maturation and processing.

Alu exaptation enriches the human transcriptome by introducing new gene ends.

In mammals, transposable elements are largely silenced, but under fortuitous circumstances may be co-opted to play a functional role. Here, we show that when Alu elements are inserted within or nearby genes in sense orientation, they may contribute to the transcriptome diversity by forming new cleavage and polyadenylation sites. We mapped these new gene ends in human onto the Alu sequence and identified three hotspots of cleavage and polyadenylation site formation. Interestingly, the native Alu sequence does not contain any canonical polyadenylation signal. We therefore studied what evolutionary processes might explain the formation of these specific hotspots of novel gene ends. We show that two of the three hotspots might have emerged from mutational processes that turned sequences that resemble polyadenylation signals into full-blown canonical signals, whereas one hotspot is tightly linked to the process of Alu insertion into the genome. Overall, Alu elements may lie behind the formation of 302 new gene end variants, affecting a total of 243 genes. Intergenic Alu elements may elongate genes by creating a downstream cleavage site, intronic Alu elements may lead to gene variants which code for truncated proteins, and 3'UTR Alu elements may result in gene variants with alternative 3'UTR.

Evidence for convergent evolution of SINE-directed Staufen-mediated mRNA decay.

Primate-specific Alu short interspersed elements (SINEs) as well as rodent-specific B and ID (B/ID) SINEs can promote Staufen-mediated decay (SMD) when present in mRNA 3'-untranslated regions (3'-UTRs). The transposable nature of SINEs, their presence in long noncoding RNAs, their interactions with Staufen, and their rapid divergence in different evolutionary lineages suggest they could have generated substantial modification of posttranscriptional gene-control networks during mammalian evolution. Some of the variation in SMD regulation produced by SINE insertion might have had a similar regulatory effect in separate mammalian lineages, leading to parallel evolution of the Staufen network by independent expansion of lineage-specific SINEs. To explore this possibility, we searched for orthologous gene pairs, each carrying a species-specific 3'-UTR SINE and each regulated by SMD, by measuring changes in mRNA abundance after individual depletion of two SMD factors, Staufen1 (STAU1) and UPF1, in both human and mouse myoblasts. We identified and confirmed orthologous gene pairs with 3'-UTR SINEs that independently function in SMD control of myoblast metabolism. Expanding to other species, we demonstrated that SINE-directed SMD likely emerged in both primate and rodent lineages >20-25 million years ago. Our work reveals a mechanism for the convergent evolution of posttranscriptional gene regulatory networks in mammals by species-specific SINE transposition and SMD.

Inferring Past Environments from Ancient Epigenomes.

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.

Identification of introns harboring functional sequence elements through positional conservation.

Many human introns carry out a function, in the sense that they are critical to maintain normal cellular activity. Their identification is fundamental to understanding cellular processes and disease. However, being noncoding elements, such functional introns are poorly predicted based on traditional approaches of sequence and structure conservation. Here, we generated a dataset of human functional introns that carry out different types of functions. We showed that functional introns share common characteristics, such as higher positional conservation along the coding sequence and reduced loss rates, regardless of their specific function. A unique property of the data is that if an intron is unknown to be functional, it still does not mean that it is indeed non-functional. We developed a probabilistic framework that explicitly accounts for this unique property, and predicts which specific human introns are functional. We show that we successfully predict function even when the algorithm is trained on introns with a different type of function. This ability has many implications in studying regulatory networks, gene regulation, the effect of mutations outside exons on human disease, and on our general understanding of intron evolution and their functional exaptation in mammals.

Gene ORGANizer: linking genes to the organs they affect.

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.