An ancestral miR-1304 allele present in Neanderthals regulates genes involved in enamel formation and could explain dental differences with modern humans

Next-generation sequencing profiling of miRNAs in individuals with 22q11.2 deletion syndrome revealed altered expression of miR-185-5p

BACKGROUND

The 22q11.2 deletion syndrome (22q11.2DS) is a microdeletion syndrome with highly variable phenotypic manifestations, even though most patients present the typical 3 Mb microdeletion, usually affecting the same ~ 106 genes. One of the genes affected by this deletion is DGCR8, which plays a crucial role in miRNA biogenesis. Therefore, the haploinsufficiency of DGCR8 due to this microdeletion can alter the modulation of the expression of several miRNAs involved in a range of biological processes.

RESULTS

In this study, we used next-generation sequencing to evaluate the miRNAs profiles in the peripheral blood of 12 individuals with typical 22q11DS compared to 12 healthy matched controls. We used the DESeq2 package for differential gene expression analysis and the DIANA-miTED dataset to verify the expression of differentially expressed miRNAs in other tissues. We used miRWalk to predict the target genes of differentially expressed miRNAs. Here, we described two differentially expressed miRNAs in patients compared to controls: hsa-miR-1304-3p, located outside the 22q11.2 region, upregulated in patients, and hsa-miR-185-5p, located in the 22q11.2 region, which showed downregulation. Expression of miR-185-5p is observed in tissues frequently affected in patients with 22q11DS, and previous studies have reported its downregulation in individuals with 22q11DS. hsa-miR-1304-3p has low expression in blood and, thus, needs more validation, though using a sensitive technology allowed us to identify differences in expression between patients and controls.

CONCLUSIONS

Thus, lower expression of miR-185-5p can be related to the 22q11.2 deletion and DGCR8 haploinsufficiency, leading to phenotypic consequences in 22q11.2DS patients, while higher expression of hsa-miR-1304-3p might be related to individual genomic variances due to the heterogeneous background of the Brazilian population.

Denisovan introgression has shaped the immune system of present-day Papuans

Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.

Denisovan and Neanderthal archaic introgression differentially impacted the genetics of complex traits in modern populations

BACKGROUND

Introgression from extinct Neanderthal and Denisovan human species has been shown to contribute to the genetic pool of modern human populations and their phenotypic spectrum. Evidence of how Neanderthal introgression shaped the genetics of human traits and diseases has been extensively studied in populations of European descent, with signatures of admixture reported for instance in genes associated with pigmentation, immunity, and metabolic traits. However, limited information is currently available about the impact of archaic introgression on other ancestry groups. Additionally, to date, no study has been conducted with respect to the impact of Denisovan introgression on the health and disease of modern populations. Here, we compare the way evolutionary pressures shaped the genetics of complex traits in East Asian and European populations, and provide evidence of the impact of Denisovan introgression on the health of East Asian and Central/South Asian populations.

RESULTS

Leveraging genome-wide association statistics from the Biobank Japan and UK Biobank, we assessed whether Denisovan and Neanderthal introgression together with other evolutionary genomic signatures were enriched for the heritability of physiological and pathological conditions in populations of East Asian and European descent. In EAS, Denisovan-introgressed loci were enriched for coronary artery disease heritability (1.69-fold enrichment, p=0.003). No enrichment for archaic introgression was observed in EUR. We also performed a phenome-wide association study of Denisovan and Neanderthal alleles in six ancestry groups available in the UK Biobank. In EAS, the Denisovan-introgressed SNP rs62391664 in the major histocompatibility complex region was associated with albumin/globulin ratio (beta=-0.17, p=3.57×10^-7). Neanderthal-introgressed alleles were associated with psychiatric and cognitive traits in EAS (e.g., "No Bipolar or Depression"-rs79043717 beta=-1.5, p=1.1×10^-7), and with blood biomarkers (e.g., alkaline phosphatase-rs11244089 beta=0.1, p=3.69×10^-116) and red hair color (rs60733936 beta=-0.86, p=4.49×10^-165) in EUR. In the other ancestry groups, Neanderthal alleles were associated with several traits, also including the use of certain medications (e.g., Central/South East Asia: indapamide - rs732632 beta=-2.38, p=5.22×10^-7).

CONCLUSIONS

Our study provides novel evidence regarding the impact of archaic introgression on the genetics of complex traits in worldwide populations, highlighting the specific contribution of Denisovan introgression in EAS populations.

Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.

Similarities and differences in the dental tissue proportions of the deciduous and permanent canines of Early and Middle Pleistocene human populations

The two- and three-dimensional assessment of dental tissues has become routine in human taxonomic studies throughout the years. Nonetheless, most of our knowledge of the variability of the enamel and dentine dimensions of the human evolutionary lineage comes from the study of permanent dentition, and particularly from molars. This leads to a biased view of the variability of these features. Due to their early formation and rapid development, the deciduous teeth allow more simplified inferences regarding the processes involved in the dental tissue development of each group. Therefore, their study could be very valuable in dental palaeohistology. In this research, we have explored the dental tissue proportions of the deciduous canines belonging to some human samples of the Early and Middle Pleistocene. The purpose of this was to discuss the meaning of the similarities and differences observed in their histological pattern, as well as to evaluate the degree of covariance with that observed in the permanent dentition of these populations. Our results show that, although there are some similarities in the dental tissue proportions between the deciduous and permanent canines of the study samples, the two dental classes do not provide a similar or comparable pictures of the dental tissue pattern present in the dentition of fossil hominins. Future works on the dental tissue patterns of the anterior and posterior dentition, including deciduous teeth, of fossil samples, may help to shed light on this hypothesis.

rs12416605:C>T in MIR938 associates with gastric cancer through affecting the regulation of the CXCL12 chemokine gene

Background

MicroRNAs are small regulatory RNAs with important roles in carcinogenesis. Genetic variants in these regulatory molecules may contribute to disease. We aim to identify allelic variants in microRNAs as susceptibility factors to gastric cancer using association studies and functional approaches.

Methods

Twenty‐one single nucleotide variants potentially functional, because of their location in either the seed, mature or precursor region of 22 microRNAs, were selected for association studies. Genetic association with gastric cancer in 365 cases and 1,284 matched controls (European Prospective Investigation into Cancer and Nutrition Cohort) was analysed using logistic regression. MicroRNA overexpression, transcriptome analysis, and target gene validation experiments were performed for functional studies.

Results

rs3746444:T>C, in the seed of MIR499A and mature MIR499B, associated with the cardia adenocarcinoma location; rs12416605:C>T, in the seed of MIR938, associated with the diffuse subtype; and rs2114358:T>C, in the precursor MIR1206, associated with the noncardia phenotype. In all cases, the association was inverse, indicating a protective affect against gastric cancer of the three minor allelic variants. MIR499 rs3746444:T>C and MIR1206 rs2114358:T>C are reported to affect the expression of these miRNAs, but the effect of MIR938 rs12416605:C>T is unknown yet. Functional approaches showed that the expression of MIR938 is affected by rs12416605:C>T and revealed that MIR938 could regulate a subset of cancer‐related genes in an allele‐specific fashion. Furthermore, we demonstrated that CXCL12, a chemokine participating in gastric cancer metastasis, is specifically regulated by only one of the rs12416605:C>T alleles.

Conclusion

rs12416605 appears to be involved in gastric cancer by affecting the regulatory function of MIR938 on genes related to this cancer type, particularly on CXCL12 posttranscriptional regulation.

Impact and Evolutionary Determinants of Neanderthal Introgression on Transcriptional and Post-Transcriptional Regulation

Archaic admixture is increasingly recognized as an important source of diversity in modern humans, and Neanderthal haplotypes cover 1%-3% of the genome of present-day Eurasians. Recent work has shown that archaic introgression has contributed to human phenotypic diversity, mostly through the regulation of gene expression. Yet the mechanisms through which archaic variants alter gene expression and the forces driving the introgression landscape at regulatory regions remain elusive. Here, we explored the impact of archaic introgression on transcriptional and post-transcriptional regulation. We focused on promoters and enhancers across 127 different tissues as well as on microRNA (miRNA)-mediated regulation. Although miRNAs themselves harbor few archaic variants, we found that some of these variants may have a strong impact on miRNA-mediated gene regulation. Enhancers were by far the regulatory elements most affected by archaic introgression: up to one-third of the tissues we tested presented significant enrichments. Specifically, we found strong enrichments of archaic variants in adipose-related tissues and primary T cells, even after accounting for various genomic and evolutionary confounders such as recombination rate and background selection. Interestingly, we identified signatures of adaptive introgression at enhancers of some key regulators of adipogenesis, raising the interesting hypothesis of a possible adaptation of early Eurasians to colder climates. Collectively, this study sheds new light on the mechanisms through which archaic admixture has impacted gene regulation in Eurasians and, more generally, increases our understanding of the contribution of Neanderthals to the regulation of acquired immunity and adipose homeostasis in modern humans.

Motifome comparison between modern human, Neanderthal and Denisovan

BACKGROUND

The availability of the genomes of two archaic humans, Neanderthal and Denisovan, and that of modern humans provides researchers an opportunity to investigate genetic differences between these three subspecies on a genome-wide scale. Here we describe an algorithm that predicts statistically significant motifs based on the difference between a given motif's actual and expected distributions. The algorithm was previously applied to plants but was modified for this work.

RESULTS

The result of applying the algorithm to the human, Neanderthal, and Denisovan genomes is a catalog of potential regulatory motifs in these three human subspecies. We examined the distributions of these motifs in genetic elements including human retroviruses, human accelerated regions, and human accelerated conserved noncoding sequences regions. Differences in these distributions could be the origin of differences in phenotype between the three subspecies. Twenty significant motifs common to all three genomes were found; thirty-three were found in endogenous retroviruses in Neanderthal and Denisovan. Ten of these motifs mapped to the 22 bp core of MiR-1304. The core of this genetic element regulates the ENAM and AMTN genes, which take part in odontogenesis and whose 3' UTRs contained significant motifs. The introns of 20 genes were found to contain a large number of significant motifs, which were also overrepresented in 49 human accelerated regions. These genes include NAV2, SorCS2, TRAPPC9, GRID1, PRDM16, CAMTA1, and ASIC which are all involved in neuroregulation. Further analysis of these genes using the GO database indicates that many are associated with neurodevelopment. Also, varying numbers of significant motifs were found to occur in regions of the Neanderthal and Denisovan genomes that are missing from the human genome, suggesting further functional differences between modern and archaic humans.

CONCLUSION

Although Neanderthal and Denisovan are now extinct, detailed examination of elements from their genomes can shed light on possible phenotypic and cognitive differences between these two archaic human subspecies and modern humans. Genetic similarities and differences between these three subspecies and other fossil hominids would also be of interest.

Computational identification of mutually homologous Zika virus miRNAs that target microcephaly genes

Background Zika virus (ZIKV) has been associated with a variety of neuropathologies, including microcephaly. We hypothesize that ZIKV genes activate host microRNAs (miRNAs) causing dysfunctional development of human fetal brains.

Objectives/methods A bioinformatics search for miRNA genome-wide binding sites in the prototypic ZIKV (strain MR766) was undertaken to hunt for miRNAs with significant similarities with MCPH genetic sequences responsible for inducing MCHP in human fetal brains.

Results Six ZIKV miRNAs were found to share mutual homology with 12 MCPH genetic sequences responsible for inducing MCPH. Noteworthy was miR-1304, which expressed 100% identity to six different MCPH genes.

Conclusions We suggest that following infection of fetal neurons ZIKV may modulate the action of various miRNAs, and miR-1304 in particular, resulting in microcephaly.

Zika virus: An emerging challenge for obstetrics and gynecology

Microcephaly is a rare birth defect, however, the re-emerging mosquito and sexual transmitted flavivirus, Zika virus (ZIKV), had changed the situation and caused an urgent challenge for the obstetrics and gynecology. This review will brief summarize the epidemiology and virology of ZIKV. And compared the animal models that had developed for the ZIKV infections. These animal models will be benefit for the development of vaccines and anti-ZIKV drugs. Furthermore, the genes that are involved in the causation of microcephaly were also summarized. Finally, the Wnt signal is critical for the brain development as well as innate immune response. Based on previous literatures, we proposed that ZIKV-induced microcephaly might result from the influence of Wnt/β-catenin signaling pathway through the regulation of miRNA-34.

Neanderthal and Denisova tooth protein variants in present-day humans

Environment parameters, diet and genetic factors interact to shape tooth morphostructure. In the human lineage, archaic and modern hominins show differences in dental traits, including enamel thickness, but variability also exists among living populations. Several polymorphisms, in particular in the non-collagenous extracellular matrix proteins of the tooth hard tissues, like enamelin, are involved in dental structure variation and defects and may be associated with dental disorders or susceptibility to caries. To gain insights into the relationships between tooth protein polymorphisms and dental structural morphology and defects, we searched for non-synonymous polymorphisms in tooth proteins from Neanderthal and Denisova hominins. The objective was to identify archaic-specific missense variants that may explain the dental morphostructural variability between extinct and modern humans, and to explore their putative impact on present-day dental phenotypes. Thirteen non-collagenous extracellular matrix proteins specific to hard dental tissues have been selected, searched in the publicly available sequence databases of Neanderthal and Denisova individuals and compared with modern human genome data. A total of 16 non-synonymous polymorphisms were identified in 6 proteins (ameloblastin, amelotin, cementum protein 1, dentin matrix acidic phosphoprotein 1, enamelin and matrix Gla protein). Most of them are encoded by dentin and enamel genes located on chromosome 4, previously reported to show signs of archaic introgression within Africa. Among the variants shared with modern humans, two are ancestral (common with apes) and one is the derived enamelin major variant, T648I (rs7671281), associated with a thinner enamel and specific to the Homo lineage. All the others are specific to Neanderthals and Denisova, and are found at a very low frequency in modern Africans or East and South Asians, suggesting that they may be related to particular dental traits or disease susceptibility in these populations. This modern regional distribution of archaic dental polymorphisms may reflect persistence of archaic variants in some populations and may contribute in part to the geographic dental variations described in modern humans.

MicroRNA-1304 suppresses human non-small cell lung cancer cell growth in vitro by targeting heme oxygenase-1

Previous studies have shown that microRNA-1304 (miR-1304) is dysregulated in certain types of cancers, including non-small cell lung cancer (NSCLC), and might be involved in tumor survival and/or growth. In this study we investigated the direct target of miR-1304 and its function in NSCLC in vitro. Human lung adenocarcinoma cell lines (A549 and NCI-H1975) were studied. The cell proliferation and survival were investigated via cell counting, MTT and colony-formation assays. Cell apoptosis and cell cycle were examined using annexin V-PE/7-AAD and PI staining assays, respectively. The dual-luciferase reporter assay was used to verify post-transcriptional regulation of heme oxygenase-1 (HO-1) by miR-1304. CRISPR/Cas9 was used to deplete endogenous miR-1304. Overexpression of MiR-1304 significantly decreased the number and viability of NSCLC cells and colony formation, and induced cell apoptosis and G₀/G₁ phase cell cycle arrest. HO-1 was demonstrated to be a direct target of miR-1304 in NSCLC cells. Restoration of HO-1 expression by hemin (20 μmol/L) abolished the inhibition of miR-1304 on cell growth and rescued miR-1304-induced apoptosis in A549 cells. Suppression of endogenous miR-1304 with anti-1304 significantly increased HO-1 expression and promoted cell growth and survival in A549 cells. In 17 human NSCLC tissue samples, miR-1304 expression was significantly decreased, while HO-1 expression was significantly increased as compared to normal lung tissues. MicroRNA-1304 is a tumor suppressor and HO-1 is its direct target in NSCLC. The results suggest the potential for miR-1304 as a therapeutic target for NSCLC.

MicroRNA Genetic Variation: From Population Analysis to Functional Implications of Three Allele Variants Associated with Cancer

Nucleotide variants in microRNA regions have been associated with disease; nevertheless, few studies still have addressed the allele-dependent effect of these changes. We studied microRNA genetic variation in human populations and found that while low-frequency variants accumulate indistinctly in microRNA regions, the mature and seed regions tend to be depleted of high-frequency variants, probably as a result of purifying selection. Comparison of pairwise population fixation indexes among regions showed that the seed had higher population fixation indexes than the other regions, suggesting the existence of local adaptation in the seed region. We further performed functional studies of three microRNA variants associated with cancer (rs2910164:C > G in MIR146A, rs11614913:C > T in MIR196A2, and rs3746444:A > G in both MIR499A and MIR499B). We found differences in the expression between alleles and in the regulation of several genes involved in cancer, such as TP53, KIT, CDH1, CLH, and TERT, which may result in changes in regulatory networks related to tumorigenesis. Furthermore, luciferase-based assays showed that MIR499A could be regulating the cadherin CDH1 and the cell adhesion molecule CLH1 in an allele-dependent fashion. A better understanding of the effect of microRNA variants associated with disease could be key in our way to a more personalized medicine.

Functional Implications of Human-Specific Changes in Great Ape microRNAs

microRNAs are crucial post-transcriptional regulators of gene expression involved in a wide range of biological processes. Although microRNAs are highly conserved among species, the functional implications of existing lineage-specific changes and their role in determining differences between humans and other great apes have not been specifically addressed. We analyzed the recent evolutionary history of 1,595 human microRNAs by looking at their intra- and inter-species variation in great apes using high-coverage sequenced genomes of 82 individuals including gorillas, orangutans, bonobos, chimpanzees and humans. We explored the strength of purifying selection among microRNA regions and found that the seed and mature regions are under similar and stronger constraint than the precursor region. We further constructed a comprehensive catalogue of microRNA species-specific nucleotide substitutions among great apes and, for the first time, investigated the biological relevance that human-specific changes in microRNAs may have had in great ape evolution. Expression and functional analyses of four microRNAs (miR-299-3p, miR-503-3p, miR-508-3p and miR-541-3p) revealed that lineage-specific nucleotide substitutions and changes in the length of these microRNAs alter their expression as well as the repertoires of target genes and regulatory networks. We suggest that the studied molecular changes could have modified crucial microRNA functions shaping phenotypes that, ultimately, became human-specific. Our work provides a frame to study the impact that regulatory changes may have in the recent evolution of our species.

Ancient genomics

The past decade has witnessed a revolution in ancient DNA (aDNA) research. Although the field's focus was previously limited to mitochondrial DNA and a few nuclear markers, whole genome sequences from the deep past can now be retrieved. This breakthrough is tightly connected to the massive sequence throughput of next generation sequencing platforms and the ability to target short and degraded DNA molecules. Many ancient specimens previously unsuitable for DNA analyses because of extensive degradation can now successfully be used as source materials. Additionally, the analytical power obtained by increasing the number of sequence reads to billions effectively means that contamination issues that have haunted aDNA research for decades, particularly in human studies, can now be efficiently and confidently quantified. At present, whole genomes have been sequenced from ancient anatomically modern humans, archaic hominins, ancient pathogens and megafaunal species. Those have revealed important functional and phenotypic information, as well as unexpected adaptation, migration and admixture patterns. As such, the field of aDNA has entered the new era of genomics and has provided valuable information when testing specific hypotheses related to the past.

Major transitions in human evolution revisited: a tribute to ancient DNA

The origin and diversification of modern humans have been characterized by major evolutionary transitions and demographic changes. Patterns of genetic variation within modern populations can help with reconstructing this ∼200 thousand year-long population history. However, by combining this information with genomic data from ancient remains, one can now directly access our evolutionary past and reveal our population history in much greater detail. This review outlines the main recent achievements in ancient DNA research and illustrates how the field recently moved from the polymerase chain reaction (PCR) amplification of short mitochondrial fragments to whole-genome sequencing and thereby revisited our own history. Ancient DNA research has revealed the routes that our ancestors took when colonizing the planet, whom they admixed with, how they domesticated plant and animal species, how they genetically responded to changes in lifestyle, and also, which pathogens decimated their populations. These approaches promise to soon solve many pending controversies about our own origins that are indecipherable from modern patterns of genetic variation alone, and therefore provide an extremely powerful toolkit for a new generation of molecular anthropologists.

IsomiR expression profiles in human lymphoblastoid cell lines exhibit population and gender dependencies

For many years it was believed that each mature microRNA (miRNA) existed as a single entity with fixed endpoints and a 'static' and unchangeable primary sequence. However, recent evidence suggests that mature miRNAs are more 'dynamic' and that each miRNA precursor arm gives rise to multiple isoforms, the isomiRs. Here we report on our identification of numerous and abundant isomiRs in the lymphoblastoid cell lines (LCLs) of 452 men and women from five different population groups. Unexpectedly, we find that these isomiRs exhibit an expression profile that is population-dependent and gender-dependent. This is important as it indicates that the LCLs of each gender/population combination have their own unique collection of mature miRNA transcripts. Moreover, each identified isomiR has its own characteristic abundance that remains consistent across biological replicates indicating that these are not degradation products. The primary sequences of identified isomiRs differ from the known miRBase miRNA either at their 5´-endpoint (leads to a different 'seed' sequence and suggests a different targetome), their 3´-endpoint, or both simultaneously. Our analysis of Argonaute PAR-CLIP data from LCLs supports the association of many of these newly identified isomiRs with the Argonaute silencing complex and thus their functional roles through participation in the RNA interference pathway.

The role of gene regulatory factors in the evolutionary history of humans

Deciphering the molecular basis of how modern human phenotypes have evolved is one of the most fascinating challenges in biology. Here, we will focus on the roles of gene regulatory factors (GRFs), in particular transcription factors (TFs) and long non-coding RNAs (lncRNAs) during human evolution. We will present examples of TFs and lncRNAs that have changed or show signs of positive selection in humans compared to chimpanzees, in modern humans compared to archaic humans, or within modern human populations. On the basis of current knowledge about the functions of these GRF genes, we speculate that they have been involved in speciation as well as in shaping phenotypes such as brain functions, skeletal morphology, and metabolic processes.

Study of Modern Human Evolution via Comparative Analysis with the Neanderthal Genome

Many other human species appeared in evolution in the last 6 million years that have not been able to survive to modern times and are broadly known as archaic humans, as opposed to the extant modern humans. It has always been considered fascinating to compare the modern human genome with that of archaic humans to identify modern human-specific sequence variants and figure out those that made modern humans different from their predecessors or cousin species. Neanderthals are the latest humans to become extinct, and many factors made them the best representatives of archaic humans. Even though a number of comparisons have been made sporadically between Neanderthals and modern humans, mostly following a candidate gene approach, the major breakthrough took place with the sequencing of the Neanderthal genome. The initial genome-wide comparison, based on the first draft of the Neanderthal genome, has generated some interesting inferences regarding variations in functional elements that are not shared by the two species and the debated admixture question. However, there are certain other genetic elements that were not included or included at a smaller scale in those studies, and they should be compared comprehensively to better understand the molecular make-up of modern humans and their phenotypic characteristics. Besides briefly discussing the important outcomes of the comparative analyses made so far between modern humans and Neanderthals, we propose that future comparative studies may include retrotransposons, pseudogenes, and conserved non-coding regions, all of which might have played significant roles during the evolution of modern humans.

Longstanding dental pathology in Neandertals from El Sidrón (Asturias, Spain) with a probable familial basis

Two Neandertal specimens from El Sidrón, northern Spain, show evidence of retained left mandibular deciduous canines. These individuals share the same mitochondrial (mtDNA) haplotype, indicating they are maternally related and suggesting a potential heritable basis for these dental anomalies. Radiographs and medical CT scans provide evidence of further, more extensive dental pathology in one of these specimens. An anomalous deciduous canine crown morphology that developed before birth subsequently suffered a fracture of the crown exposing the pulp sometime after eruption into functional occlusion. This led to death of the tooth, periapical granuloma formation and arrested deciduous canine root growth at an estimated age of 2.5 years. At some point the underlying permanent canine tooth became horizontally displaced and came to lie low in the trabecular bone of the mandibular corpus. A dentigerous cyst then developed around the crown. Anterior growth displacement of the mandible continued around the stationary permanent canine, leaving it posteriorly positioned in the mandibular corpus by the end of the growth period beneath the third permanent molar roots, which, in turn, suggests a largely horizontal growth vector. Subsequent longstanding repeated infections of the expanding cyst cavity are evidenced by bouts of bone deposition and resorption of the boundary walls of the cyst cavity. This resulted in the establishment of two permanent bony drainage sinuses, one through the buccal plate of the alveolar bone anteriorly, immediately beneath the infected deciduous canine root, and the other through the buccal plate anterior to the mesial root of the first permanent molar. It is probable that this complicated temporal sequence of dental pathologies had an initial heritable trigger that progressed in an unusually complex way in one of these individuals. During life, this individual may have been largely unaware of this ongoing pathology.

Archaic human genomics

For much of the 20th century, the predominant view of human evolutionary history was derived from the fossil record. Homo erectus was seen arising in Africa from an earlier member of the genus and then spreading throughout the Old World and into the Oceania. A regional continuity model of anagenetic change from H. erectus via various intermediate archaic species into the modern humans in each of the regions inhabited by H. erectus was labeled the multiregional model of human evolution (MRE). A contrasting model positing a single origin, in Africa, of anatomically modern H. sapiens with some populations later migrating out of Africa and replacing the local archaic populations throughout the world with complete replacement became known as the recent African origin (RAO) model. Proponents of both models used different interpretations of the fossil record to bolster their views for decades. In the 1980s, molecular genetic techniques began providing evidence from modern human variation that allowed not only the different models of modern human origins to be tested but also the exploration demographic history and the types of selection that different regions of the genome and even specific traits had undergone. The majority of researchers interpreted these data as strongly supporting the RAO model, especially analyses of mitochondrial DNA (mtDNA). Extrapolating backward from modern patterns of variation and using various calibration points and substitution rates, a consensus arose that saw modern humans evolving from an African population around 200,000 years ago. Much later, around 50,000 years ago, a subset of this population migrated out of Africa replacing Neanderthals in Europe and western Asia as well as archaics in eastern Asia and Oceania. mtDNA sequences from more than two-dozen Neanderthals and early modern humans re-enforced this consensus. In 2010, however, the complete draft genomes of Neanderthals and of heretofore unknown hominins from Siberia, called Denisovans, demonstrated gene flow between these archaic human species and modern Eurasians but not sub-Saharan Africans. Although the levels of gene flow may be very limited, this unexpected finding does not fit well with either the RAO model or MRE model. More thorough sampling of modern human diversity, additional fossil discoveries, and the sequencing of additional hominin fossils are necessary to throw light onto our origins and our history.