SARS-CoV-2 ORF3c impairs mitochondrial respiratory metabolism, oxidative stress, and autophagic flux.

Coronaviruses encode a variable number of accessory proteins that are involved in host-virus interaction, suppression of immune responses, or immune evasion. SARS-CoV-2 encodes at least twelve accessory proteins, whose roles during infection have been studied. Nevertheless, the role of the ORF3c accessory protein, an alternative open reading frame of ORF3a, has remained elusive. Herein, we show that the ORF3c protein has a mitochondrial localization and alters mitochondrial metabolism, inducing a shift from glucose to fatty acids oxidation and enhanced oxidative phosphorylation. These effects result in increased ROS production and block of the autophagic flux. In particular, ORF3c affects lysosomal acidification, blocking the normal autophagic degradation process and leading to autolysosome accumulation. We also observed different effect on autophagy for SARS-CoV-2 and batCoV RaTG13 ORF3c proteins; the 36R and 40K sites are necessary and sufficient to determine these effects.

Whole-genome sequencing of hMPXV1 in five Italian cases confirms the occurrence of the predominant epidemic lineage.

The ongoing outbreak of monkeypox virus (hMPXV1) is the largest recorded in historically nonendemic countries. Genomic surveillance has emerged as a pivotal tool to track the spread and monitor the evolution of viral pathogens. Therefore, to assess the genetic diversity of circulating hMPXV1 in northern Italy in June to July 2022, we sequenced and analyzed five complete genomes of viruses sampled from patients presenting with a typical course of hMPXV1 infection. Phylogenetic analysis confirmed that all five genomes belong to the predominant epidemic lineage (B.1). Inspection of genetic changes and comparison with the reference sequence showed the presence of 12 nucleotide substitutions. Seven are nonsynonymous mutations leading to amino acid changes in six proteins belonging to different functional classes. Moreover, 11 of these 12 nucleotide mutations involve GA>AA or TC>TT replacements, suggesting that host APOBEC3 enzymes are responsible for the generation of substitutions in circulating viruses. Finally, metagenomic analysis evidenced bacterial superinfection (Streptococcus pyogenes) in one patient. Through this study, we contributed to expand the number of complete genomes of viruses circulating in Italy and characterize them as belonging to the predominant outbreak lineage.

Evolution of the orthopoxvirus core genome.

Orthopoxviruses comprise several relevant pathogens, including the causative agent of smallpox and monkeypox virus. Analysis of orthopoxvirus genome evolution mainly focused on gene gains/losses. We instead analyzed core genes, which are conserved in all orthopoxviruses. We show that, despite their strong constraint, some genes involved in viral morphogenesis and transcription/replication were targets of pervasive positive selection, which was relatively uncommon in immunomodulatory genes. However at least three of the positively selected genes, E3L, A24R, and H3L, might have evolved in response to immune selection. Episodic positive selection was particularly common on the internal branches of the orthopox phylogeny and on the monkeypox virus lineage. The latter showed evidence of episodic positive selection at the D14L gene, which encodes a modulator of complement activation (MOPICE). Notably, two genes (B1R and A33R) targeted by episodic selection on more than one branch are involved in forms of intra-genomic conflict. Finally, we found that, in orthopoxvirus proteomes, intrinsically disordered regions (IDRs) tend to be less constrained and are common targets of positive selection. Extension of our analysis to all poxviruses showed no evidence that the IDR fraction differs with host range. Conversely, we found a strong effect of base composition, which was however not sufficient to explain IDR fraction. We thus suggest that, in poxviruses, the IDR fraction is maintained by modulating GC content to accommodate disorder-promoting codons. Overall, our data provide novel insight in orthopoxvirus evolution and provide a list of genes and sites that are expected to modulate viral phenotypes.

Simplexviruses Successfully Adapt to Their Host by Fine-Tuning Immune Responses.

Primate herpes simplex viruses are species-specific and relatively harmless to their natural hosts. However, cross-species transmission is often associated with severe disease, as exemplified by the virulence of macacine herpesvirus 1 (B virus) in humans. We performed a genome-wide scan for signals of adaptation of simplexviruses to their hominin hosts. Among core genes, we found evidence of episodic positive selection in three glycoproteins, with several selected sites located in antigenic determinants. Positively selected noncore genes were found to be involved in different immune-escape mechanisms. The herpes simplex virus (HSV)-1/HSV-2 encoded product (ICP47) of one of these genes is known to down-modulate major histocompatibility complex class I expression. This feature is not shared with B virus, which instead up-regulates Human Leukocyte Antigen (HLA)-G, an immunomodulatory molecule. By in vitro expression of different ICP47 mutants, we functionally characterized the selection signals. Results indicated that the selected sites do not represent the sole determinants of binding to the transporter associated with antigen processing (TAP). Conversely, the amino acid status at these sites was sufficient to determine HLA-G up-regulation. In fact, both HSV-1 and HSV-2 ICP47 induced HLA-G when mutated to recapitulate residues in B virus, whereas the mutated version of B virus ICP47 failed to determine HLA-G expression. These differences might contribute to the severity of B virus infection in humans. Importantly, they indicate that the evolution of ICP47 in HSV-1/HSV-2 led to the loss of an immunosuppressive effect. Thus, related simplexviruses finely tune the balance between immunosuppressive and immunostimulatory pathways to promote successful co-existence with their primate hosts.

Dating the Emergence of Human Endemic Coronaviruses.

Four endemic coronaviruses infect humans and cause mild symptoms. Because previous analyses were based on a limited number of sequences and did not control for effects that affect molecular dating, we re-assessed the timing of endemic coronavirus emergence. After controlling for recombination, selective pressure, and molecular clock model, we obtained similar tMRCA (time to the most recent common ancestor) estimates for the four coronaviruses, ranging from 72 (HCoV-229E) to 54 (HCoV-NL63) years ago. The split times of HCoV-229E and HCoV-OC43 from camel alphacoronavirus and bovine coronavirus were dated ~268 and ~99 years ago. The split times of HCoV-HKU1 and HCoV-NL63 could not be calculated, as their zoonoticic sources are unknown. To compare the timing of coronavirus emergence to that of another respiratory virus, we recorded the occurrence of influenza pandemics since 1500. Although there is no clear relationship between pandemic occurrence and human population size, the frequency of influenza pandemics seems to intensify starting around 1700, which corresponds with the initial phase of exponential increase of human population and to the emergence of HCoV-229E. The frequency of flu pandemics in the 19th century also suggests that the concurrence of HCoV-OC43 emergence and the Russian flu pandemic may be due to chance.

Homology-based classification of accessory proteins in coronavirus genomes uncovers extremely dynamic evolution of gene content.

Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super-groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike-like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion-associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross-communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.

Adaptation of the endemic coronaviruses HCoV-OC43 and HCoV-229E to the human host.

Four coronaviruses (HCoV-OC43, HCoV-HKU1, HCoV-NL63, and HCoV-229E) are endemic in human populations. All these viruses are seasonal and generate short-term immunity. Like the highly pathogenic coronaviruses, the endemic coronaviruses have zoonotic origins. Thus, understanding the evolutionary dynamics of these human viruses might provide insight into the future trajectories of SARS-CoV-2 evolution. Because the zoonotic sources of HCoV-OC43 and HCoV-229E are known, we applied a population genetics-phylogenetic approach to investigate which selective events accompanied the divergence of these viruses from the animal ones. Results indicated that positive selection drove the evolution of some accessory proteins, as well as of the membrane proteins. However, the spike proteins of both viruses and the hemagglutinin-esterase (HE) of HCoV-OC43 represented the major selection targets. Specifically, for both viruses, most positively selected sites map to the receptor-binding domains (RBDs) and are polymorphic. Molecular dating for the HCoV-229E spike protein indicated that RBD Classes I, II, III, and IV emerged 3-9 years apart. However, since the appearance of Class V (with much higher binding affinity), around 25 years ago, limited genetic diversity accumulated in the RBD. These different time intervals are not fully consistent with the hypothesis that HCoV-229E spike evolution was driven by antigenic drift. An alternative, not mutually exclusive possibility is that strains with higher affinity for the cellular receptor have out-competed strains with lower affinity. The evolution of the HCoV-OC43 spike protein was also suggested to undergo antigenic drift. However, we also found abundant signals of positive selection in HE. Whereas such signals might result from antigenic drift, as well, previous data showing co-evolution of the spike protein with HE suggest that optimization for human cell infection also drove the evolution of this virus. These data provide insight into the possible trajectories of SARS-CoV-2 evolution, especially in case the virus should become endemic.

Antigenic variation of SARS-CoV-2 in response to immune pressure.

Analysis of the bat viruses most closely related to SARS-CoV-2 indicated that the virus probably required limited adaptation to spread in humans. Nonetheless, since its introduction in human populations, SARS-CoV-2 must have been subject to the selective pressure imposed by the human immune system. We exploited the availability of a large number of high-quality SARS-CoV-2 genomes, as well as of validated epitope predictions, to show that B cell epitopes in the spike glycoprotein (S) and in the nucleocapsid protein (N) have higher diversity than nonepitope positions. Similar results were obtained for other human coronaviruses and for sarbecoviruses sampled in bats. Conversely, in the SARS-CoV-2 population, epitopes for CD4+ and CD8+ T cells were not more variable than nonepitope positions. A significant reduction in epitope variability was instead observed for some of the most immunogenic proteins (S, N, ORF8 and ORF3a). Analysis over longer evolutionary time frames indicated that this effect is not due to differential constraints. These data indicate that SARS-CoV-2 evolves to elude the host humoral immune response, whereas recognition by T cells is not actively avoided by the virus. However, we also found a trend of lower diversity of T cell epitopes for common cold coronaviruses, indicating that epitope conservation per se is not directly linked to disease severity. We suggest that conservation serves to maintain epitopes that elicit tolerizing T cell responses or induce T cells with regulatory activity.

Role of NEU3 Overexpression in the Prediction of Efficacy of EGFR-Targeted Therapies in Colon Cancer Cell Lines.

The epidermal growth factor receptor (EGFR), through the MAP kinase and PI3K-Akt-mTOR axis, plays a pivotal role in colorectal cancer (CRC) pathogenesis. The membrane-associated NEU3 sialidase interacts with and desialylates EGFR by promoting its dimerization and downstream effectors' activation. Among the targeted therapies against EGFR, the monoclonal antibody cetuximab is active only in a subgroup of patients not carrying mutations in the MAP kinase pathway. In order to better understand the EGFR-NEU3 interplay and the mechanisms of pharmacological resistance, we investigated the role of NEU3 deregulation in cetuximab-treated CRC cell lines transiently transfected with NEU3 using Western blot analysis. Our results indicate that NEU3 overexpression can enhance EGFR activation only if EGFR is overexpressed, indicating the existence of a threshold for NEU3-mediated EGFR activation. This enhancement mainly leads to the constitutive activation of the MAP kinase pathway. Consequently, we suggest that the evaluation of NEU3 expression cannot entirely substitute the evaluation of EGFR because EGFR-negative cases cannot be stimulated by NEU3. Furthermore, NEU3-mediated hyperactivation of EGFR is counterbalanced by the administration of cetuximab, hypothesizing that a combined treatment of NEU3- and EGFR-targeted therapies may represent a valid option for CRC patients, which must be investigated in the future.

Past and ongoing adaptation of human cytomegalovirus to its host.

Cytomegaloviruses (order Herpesvirales) display remarkable species-specificity as a result of long-term co-evolution with their mammalian hosts. Human cytomegalovirus (HCMV) is exquisitely adapted to our species and displays high genetic diversity. We leveraged information on inter-species divergence of primate-infecting cytomegaloviruses and intra-species diversity of clinical isolates to provide a genome-wide picture of HCMV adaptation across different time-frames. During adaptation to the human host, core viral genes were commonly targeted by positive selection. Functional characterization of adaptive mutations in the primase gene (UL70) indicated that selection favored amino acid replacements that decrease viral replication in human fibroblasts, suggesting evolution towards viral temperance. HCMV intra-species diversity was largely governed by immune system-driven selective pressure, with several adaptive variants located in antigenic domains. A significant excess of positively selected sites was also detected in the signal peptides (SPs) of viral proteins, indicating that, although they are removed from mature proteins, SPs can contribute to viral adaptation. Functional characterization of one of these SPs indicated that adaptive variants modulate the timing of cleavage by the signal peptidase and the dynamics of glycoprotein intracellular trafficking. We thus used evolutionary information to generate experimentally-testable hypotheses on the functional effect of HCMV genetic diversity and we define modulators of viral phenotypes.

Intrinsically disordered regions are abundant in simplexvirus proteomes and display signatures of positive selection.

Whereas the majority of herpesviruses co-speciated with their mammalian hosts, human herpes simplex virus 2 (HSV-2, genus Simplexvirus) most likely originated from the cross-species transmission of chimpanzee herpesvirus 1 to an ancestor of modern humans. We exploited the peculiar evolutionary history of HSV-2 to investigate the selective events that drove herpesvirus adaptation to a new host. We show that HSV-2 intrinsically disordered regions (IDRs)-that is, protein domains that do not adopt compact three-dimensional structures-are strongly enriched in positive selection signals. Analysis of viral proteomes indicated that a significantly higher portion of simplexvirus proteins is disordered compared with the proteins of other human herpesviruses. IDR abundance in simplexvirus proteomes was not a consequence of the base composition of their genomes (high G + C content). Conversely, protein function determines the IDR fraction, which is significantly higher in viral proteins that interact with human factors. We also found that the average extent of disorder in herpesvirus proteins tends to parallel that of their human interactors. These data suggest that viruses that interact with fast-evolving, disordered human proteins, in turn, evolve disordered viral interactors poised for innovation. We propose that the high IDR fraction present in simplexvirus proteomes contributes to their wider host range compared with other herpesviruses.

Evolution and Genetic Diversity of Primate Cytomegaloviruses.

Cytomegaloviruses (CMVs) infect many mammals, including humans and non-human primates (NHPs). Human cytomegalovirus (HCMV) is an important opportunistic pathogen among immunocompromised patients and represents the most common infectious cause of birth defects. HCMV possesses a large genome and very high genetic diversity. NHP-infecting CMVs share with HCMV a similar genomic organization and coding content, as well as the course of viral infection. Recent technological advances have allowed the sequencing of several HCMV strains from clinical samples and provided insight into the diversity of NHP-infecting CMVs. The emerging picture indicates that, with the exclusion of core genes (genes that have orthologs in all herpesviruses), CMV genomes are relatively plastic and diverse in terms of gene content, both at the inter- and at the intra-species level. Such variability most likely underlies the strict species-specificity of these viruses, as well as their ability to persist lifelong and with relatively little damage to their hosts. However, core genes, despite their strong conservation, also represented a target of adaptive evolution and subtle changes in their coding sequence contributed to CMV adaptation to different hosts. Indubitably, important knowledge gaps remain, the most relevant of which concerns the role of viral genetics in HCMV-associated human disease.

The influence of DCDC2 risk genetic variants on reading: Testing main and haplotypic effects.

Developmental dyslexia (DD) is a complex neurodevelopmental heritable disorder. Among DD candidate genes, DCDC2 is one of the most replicated, with rs793862, READ1 and rs793842 likely contribute to phenotypic variability in reading (dis)ability. In this study, we tested the effects of these genetic variants on DD as a categorical trait and on quantitative reading-related measures in a sample of 555 Italian nuclear families with 930 offspring, of which 687 were diagnosed with DD. We conducted both single-marker and haplotype analyses, finding that the READ1-deletion was significantly associated with reading, whereas no significant haplotype associations were found. Our findings add further evidence to support the hypothesis of a DCDC2 contribution to inter-individual variation in distinct indicators of reading (dis)ability in transparent languages (i.e., reading accuracy and speed), suggesting a potential pleiotropic effect.

The Diversity of Mammalian Hemoproteins and Microbial Heme Scavengers Is Shaped by an Arms Race for Iron Piracy.

Iron is an essential micronutrient for most living species. In mammals, hemoglobin (Hb) stores more than two thirds of the body's iron content. In the bloodstream, haptoglobin (Hp) and hemopexin (Hpx) sequester free Hb or heme. Pathogenic microorganisms usually acquire iron from their hosts and have evolved complex systems of iron piracy to circumvent nutritional immunity. Herein, we performed an evolutionary analysis of genes coding for mammalian heme-binding proteins and heme-scavengers in pathogen species. The underlying hypothesis is that these molecules are engaged in a molecular arms race. We show that positive selection drove the evolution of mammalian Hb and Hpx. Positively selected sites in Hb are located at the interaction surface with Neisseria meningitidis heme scavenger HpuA and with Staphylococcus aureus iron-regulated surface determinant B (IsdB). In turn, positively selected sites in HpuA and IsdB are located in the flexible protein regions that contact Hb. A residue in Hb (S45H) was also selected on the Caprinae branch. This site stabilizes the interaction with Trypanosoma brucei hemoglobin-haptoglobin (HbHp) receptor (TbHpHbR), a molecule that also mediates trypanosome lytic factor (TLF) entry. In TbHpHbR, positive selection drove the evolution of a variant (L210S) which allows evasion from TLF but reduces affinity for HbHp. Finally, selected sites in Hpx are located at the interaction surface with the Haemophilus influenzae hemophore HxuA, which in turn displays fast evolving sites at the Hpx-binding interface. These results shed light into host-pathogens conflicts and establish the importance of nutritional immunity as an evolutionary force.

Genetic susceptibility to infectious diseases: Current status and future perspectives from genome-wide approaches.

Genome-wide association studies (GWASs) have been widely applied to identify genetic factors that affect complex diseases or traits. Presently, the GWAS Catalog includes >2800 human studies. Of these, only a minority have investigated the susceptibility to infectious diseases or the response to therapies for the treatment or prevention of infections. Despite their limited application in the field, GWASs have provided valuable insights by pinpointing associations to both innate and adaptive immune response loci, as well as novel unexpected risk factors for infection susceptibility. Herein, we discuss some issues and caveats of GWASs for infectious diseases, we review the most recent findings ensuing from these studies, and we provide a brief summary of selected GWASs for infections in non-human mammals. We conclude that, although the general trend in the field of complex traits is to shift from GWAS to next-generation sequencing, important knowledge on infectious disease-related traits can be still gained by GWASs, especially for those conditions that have never been investigated using this approach. We suggest that future studies will benefit from the leveraging of information from the host's and pathogen's genomes, as well as from the exploration of models that incorporate heterogeneity across populations and phenotypes. Interactions within HLA genes or among HLA variants and polymorphisms located outside the major histocompatibility complex may also play an important role in shaping the susceptibility and response to invading pathogens.

Deregulation of sialidases in human normal and tumor tissues.

BACKGROUND: Aberrant sialylation is a characteristic feature associated with cancer. The four types of mammalian sialidases identified to date have been shown to behave in different manners during carcinogenesis. While NEU1, NEU2 and NEU4 have been observed to oppose malignant phenotypes, the membrane-bound sialidase NEU3 was revealed to promote cancer progression. OBJECTIVES: With the aim of improving the knowledge about sialidases deregulation in various cancer types, we investigated the amount of NEU1, NEU3 and NEU4 transcripts in paired normal and tumor tissues from 170 patients with 11 cancer types. METHODS: mRNA was extracted from patients' tissue specimens and retrotranscribed into cDNA, which was quantified by Real-Time PCR. RESULTS: We found NEU1 and NEU3 to be up regulated, while NEU4 was down regulated in most cancer types. In particular, colorectal cancer tissues showed the highest increase in NEU3 expression. Both NEU1 and NEU3 showed a strong up-regulation in ovarian cancer. CONCLUSIONS: Our data show that human sialidases are expressed at different levels in healthy tissues and are strongly deregulated in tumors. Moreover, sialidases expression in our European cohort showed significant differences from Asian populations. Some of these peculiar features open potential applications of sialidases in cancer diagnosis and therapy.

REST, a master regulator of neurogenesis, evolved under strong positive selection in humans and in non human primates.

The transcriptional repressor REST regulates many neuronal genes by binding RE1 motifs. About one third of human RE1s are recently evolved and specific to primates. As changes in the activity of a transcription factor reverberate on its downstream targets, we assessed whether REST displays fast evolutionary rates in primates. We show that REST was targeted by very strong positive selection during primate evolution. Positive selection was also evident in the human lineage, with six selected sites located in a region that surrounds a VNTR in exon 4. Analysis of expression data indicated that REST brain expression peaks during aging in humans but not in other primates. Because a REST coding variant (rs3796529) was previously associated with protection from hippocampal atrophy in elderly subjects with mild cognitive impairment (MCI), we analyzed a cohort of Alzheimer disease (AD) continuum patients. Genotyping of two coding variants (rs3796529 and rs2227902) located in the region surrounding the VNTR indicated a role for rs2227902 in modulation of hippocampal volume loss, indirectly confirming a role for REST in neuroprotection. Experimental studies will be instrumental to determine the functional effect of positively selected sites in REST and the role of REST variants in neuropreservation/neurodegeneration.

Distinct selective forces and Neanderthal introgression shaped genetic diversity at genes involved in neurodevelopmental disorders.

In addition to high intelligence, humans evolved specialized social-cognitive skills, which are specifically affected in children with autism spectrum disorder (ASD). Genes affected in ASD represent suitable candidates to study the evolution of human social cognition. We performed an evolutionary analysis on 68 genes associated to neurodevelopmental disorders; our data indicate that genetic diversity was shaped by distinct selective forces, including natural selection and introgression from archaic hominins. We discuss the possibility that segregation distortion during spermatogenesis accounts for a subset of ASD mutations. Finally, we detected modern-human-specific alleles in DYRK1A and TCF4. These variants are located within regions that display chromatin features typical of transcriptional enhancers in several brain areas, strongly suggesting a regulatory role. These SNPs thus represent candidates for association with neurodevelopmental disorders, and await experimental validation in future studies.

A common genetic variant in FOXP2 is associated with language-based learning (dis)abilities: Evidence from two Italian independent samples.

Language-based Learning Disabilities (LLDs) encompass a group of complex, comorbid, and developmentally associated deficits in communication. Language impairment and developmental dyslexia (DD) represent the most recognized forms of LLDs. Substantial genetic correlations exist between language and reading (dis)abilities. Common variants in the FOXP2 gene were consistently associated with language- and reading-related neuropsychological and neuroanatomical phenotypes. We tested the effect of a FOXP2 common variant, that is, rs6980093 (A/G), on quantitative measures of language and reading in two independent Italian samples: a population-based cohort of 699 subjects (3-11 years old) and a sample of 572 children with DD (6-18 years old). rs6980093 modulates expressive language in the general population sample, with an effect on fluency scores. In the DD sample, the variant showed an association with the accuracy in the single word reading task. rs6980093 shows distinct genetic models of association in the two cohorts, with a dominant effect of the G allele in the general population sample and heterozygote advantage in the DD cohort. We provide preliminary evidence that rs6980093 associates with language and reading (dis)abilities in two independent Italian cohorts. rs6980093 is an intronic SNP, suggesting that it (or a linked variant) modulates phenotypic association via regulation of FOXP2 expression. Because FOXP2 brain expression is finely regulated, both temporally and spatially, it is possible that the two alleles at rs6980093 differentially modulate expression levels in a developmental stage- or brain area-specific manner. This might help explaining the heterozygote advantage effect and the different genetic models in the two cohorts.

The evolutionary history of genes involved in spoken and written language: beyond FOXP2.

Humans possess a communication system based on spoken and written language. Other animals can learn vocalization by imitation, but this is not equivalent to human language. Many genes were described to be implicated in language impairment (LI) and developmental dyslexia (DD), but their evolutionary history has not been thoroughly analyzed. Herein we analyzed the evolution of ten genes involved in DD and LI. Results show that the evolutionary history of LI genes for mammals and aves was comparable in vocal-learner species and non-learners. For the human lineage, several sites showing evidence of positive selection were identified in KIAA0319 and were already present in Neanderthals and Denisovans, suggesting that any phenotypic change they entailed was shared with archaic hominins. Conversely, in FOXP2, ROBO1, ROBO2, and CNTNAP2 non-coding changes rose to high frequency after the separation from archaic hominins. These variants are promising candidates for association studies in LI and DD.