Alessandra Mozzi

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.

OASes and STING: Adaptive Evolution in Concert

OAS (2′–5′-oligoadenylate synthases) proteins and cyclic GMP–AMP synthase (cGAS, gene symbol: MB21D1) patrol the cytoplasm for the presence of foreign nucleic acids. Upon binding to double-stranded RNA or double-stranded DNA, OAS proteins and cGAS produce nucleotide second messengers to activate RNase L and STING (stimulator of interferon genes, gene symbol: TMEM173), respectively; this leads to the initiation of antiviral responses. We analyzed the evolutionary history of the MB21D1–TMEM173 and OAS–RNASEL axes in primates and bats and found evidence of widespread positive selection in both orders. In TMEM173, residue 230, a major determinant of response to natural ligands and to mimetic drugs (e.g., DMXAA), was positively selected in Primates and Chiroptera. In both orders, selection also targeted an α-helix/loop element in RNase L that modulates the enzyme preference for single-stranded RNA versus stem loops. Analysis of positively selected sites in OAS1, OAS2, and MB21D1 revealed parallel evolution, with the corresponding residues being selected in different genes. As this cannot result from gene conversion, these data suggest that selective pressure acting on OAS and MB21D1 genes is related to nucleic acid recognition and to the specific mechanism of enzyme activation, which requires a conformational change. Finally, a population genetics-phylogenetics analysis in humans, chimpanzees, and gorillas detected several positively selected sites in most genes. Data herein shed light into species-specific differences in infection susceptibility and in response to synthetic compounds, with relevance for the design of synthetic compounds as vaccine adjuvants.

Extensive Positive Selection Drives the Evolution of Nonstructural Proteins in Lineage C Betacoronaviruses

ABSTRACT Middle East respiratory syndrome-related coronavirus (MERS-CoV) spreads to humans via zoonotic transmission from camels. MERS-CoV belongs to lineage C of betacoronaviruses (betaCoVs), which also includes viruses isolated from bats and hedgehogs. A large portion of the betaCoV genome consists of two open reading frames (ORF1a and ORF1b) that are translated into polyproteins. These are cleaved by viral proteases to generate 16 nonstructural proteins (nsp1 to nsp16) which compose the viral replication-transcription complex. We investigated the evolution of ORF1a and ORF1b in lineage C betaCoVs. Results indicated widespread positive selection, acting mostly on ORF1a. The proportion of positively selected sites in ORF1a was much higher than that previously reported for the surface-exposed spike protein. Selected sites were unevenly distributed, with nsp3 representing the preferential target. Several pairs of coevolving sites were also detected, possibly indicating epistatic interactions; most of these were located in nsp3. Adaptive evolution at nsp3 is ongoing in MERS-CoV strains, and two selected sites (G720 and R911) were detected in the protease domain. While position 720 is variable in camel-derived viruses, suggesting that the selective event does not represent a specific adaptation to humans, the R911C substitution was observed only in human-derived MERS-CoV isolates, including the viral strain responsible for the recent South Korean outbreak. It will be extremely important to assess whether these changes affect host range or other viral phenotypes. More generally, data herein indicate that CoV nsp3 represents a major selection target and that nsp3 sequencing should be envisaged in monitoring programs and field surveys. IMPORTANCE Both severe acute respiratory syndrome coronavirus (SARS-CoV) and MERS-CoV originated in bats and spread to humans via an intermediate host. This clearly highlights the potential for coronavirus host shifting and the relevance of understanding the molecular events underlying the adaptation to new host species. We investigated the evolution of ORF1a and ORF1b in lineage C betaCoVs and in 87 sequenced MERS-CoV isolates. Results indicated widespread positive selection, stronger in ORF1a than in ORF1b. Several selected sites were found to be located in functionally relevant protein regions, and some of them corresponded to functional mutations in other coronaviruses. The proportion of selected sites we identified in ORF1a is much higher than that for the surface-exposed spike protein. This observation suggests that adaptive evolution in ORF1a might contribute to host shifts or immune evasion. Data herein also indicate that genetic diversity at nonstructural proteins should be taken into account when antiviral compounds are developed.

Natural Selection at the Brush-Border: Adaptations to Carbohydrate Diets in Humans and Other Mammals

Dietary shifts can drive molecular evolution in mammals and a major transition in human history, the agricultural revolution, favored carbohydrate consumption. We investigated the evolutionary history of nine genes encoding brush-border proteins involved in carbohydrate digestion/absorption. Results indicated widespread adaptive evolution in mammals, with several branches experiencing episodic selection, particularly strong in bats. Many positively selected sites map to functional protein regions (e.g., within glucosidase catalytic crevices), with parallel evolution at SI (sucrase-isomaltase) and MGAM (maltase-glucoamylase). In human populations, five genes were targeted by positive selection acting on noncoding variants within regulatory elements. Analysis of ancient DNA samples indicated that most derived alleles were already present in the Paleolithic. Positively selected variants at SLC2A5 (fructose transporter) were an exception and possibly spread following the domestication of specific fruit crops. We conclude that agriculture determined no major selective event at carbohydrate metabolism genes in humans, with implications for susceptibility to metabolic disorders.

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 mammalian complement system as an epitome of host-pathogen genetic conflicts.

The complement system is an innate immunity effector mechanism; its action is antagonized by a wide array of pathogens and complement evasion determines the virulence of several infections. We investigated the evolutionary history of the complement system and of bacterial‐encoded complement‐interacting proteins. Complement components targeted by several pathogens evolved under strong selective pressure in primates, with selection acting on residues at the contact interface with microbial/viral proteins. Positively selected sites in CFH and C4BPA account for the human specificity of gonococcal infection. Bacterial interactors, evolved adaptively as well, with selected sites located at interaction surfaces with primate complement proteins. These results epitomize the expectation under a genetic conflict scenario whereby the hosts and the pathogens genes evolve within binding avoidance‐binding seeking dynamics. In silico mutagenesis and protein–protein docking analyses supported this by showing that positively selected sites, both in the hosts and in the pathogens interacting partner, modulate binding.

Diverse selective regimes shape genetic diversity at ADAR genes and at their coding targets

A-to-I RNA editing operated by ADAR enzymes is extremely common in mammals. Several editing events in coding regions have pivotal physiological roles and affect protein sequence (recoding events) or function. We analyzed the evolutionary history of the 3 ADAR family genes and of their coding targets. Evolutionary analysis indicated that ADAR evolved adaptively in primates, with the strongest selection in the unique N-terminal domain of the interferon-inducible isoform. Positively selected residues in the human lineage were also detected in the ADAR deaminase domain and in the RNA binding domains of ADARB1 and ADARB2. During the recent history of human populations distinct variants in the 3 genes increased in frequency as a result of local selective pressures. Most selected variants are located within regulatory regions and some are in linkage disequilibrium with eQTLs in monocytes. Finally, analysis of conservation scores of coding editing sites indicated that editing events are counter-selected within regions that are poorly tolerant to change. Nevertheless, a minority of recoding events occurs at highly conserved positions and possibly represents the functional fraction. These events are enriched in pathways related to HIV-1 infection and to epidermis/hair development. Thus, both ADAR genes and their targets evolved under variable selective regimes, including purifying and positive selection. Pressures related to immune response likely represented major drivers of evolution for ADAR genes. As for their coding targets, we suggest that most editing events are slightly deleterious, although a minority may be beneficial and contribute to antiviral response and skin homeostasis.

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.

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.

Albuminoid Genes: Evolving at the Interface of Dispensability and Selection

The albuminoid gene family comprises vitamin D-binding protein (GC), alpha-fetoprotein (AFP), afamin (AFM), and albumin (ALB). Albumin is the most abundant human serum protein, and, as the other family members, acts as a transporter of endogenous and exogenous substances including thyroxine, fatty acids, and drugs. Instead, the major cargo of GC is 25-hydroxyvitamin D. We performed an evolutionary study of albuminoid genes and we show that ALB evolved adaptively in mammals. Most positively selected sites are located within albumin-binding sites for fatty acids and thyroxine, as well as at the contact surface with neonatal Fc receptor. Positive selection was also detected for residues forming the prostaglandin-binding pocket. Adaptation to hibernation/torpor might explain the signatures of episodic positive selection we detected for few mammalian lineages. Application of a population genetics–phylogenetics approach showed that purifying selection represented a major force acting on albuminoid genes in both humans and chimpanzees, with the strongest constraint observed for human GC. Population genetic analysis revealed that GC was also the target of locally exerted selective pressure, which drove the frequency increase of different haplotypes in distinct human populations. A search for known variants that modulate GC and 25-hydroxyvitamin D concentrations revealed linkage disequilibrium with positively selected variants, although European and Asian major GC haplotypes carry alleles with reported opposite effect on GC concentration. Data herein indicate that albumin, an extremely abundant housekeeping protein, was the target of pervasive and episodic selection in mammals, whereas GC represented a selection target during the recent evolution of human populations.