Genome-scale detection of positive selection in nine primates predicts human-virus evolutionary conflicts

NOVEL SIMPLEXVIRUS (SIMPLEXVIRUS DOLICHOTINEALPHA1) ASSOCIATED WITH FATALITY IN FOUR PATAGONIAN MARA (DOLICHOTIS PATAGONUM)

Four of seven Patagonian maras (Dolichotis patagonum) at a zoological institution developed acute neurologic signs that progressed to tetraparesis and death. All affected were young adult females (10 mon-5 yr old) that presented over 11 d. Clinical signs were rapidly progressive and unresponsive to supportive therapies. Two of the four individuals were found deceased 4 d after hospitalization. Two individuals were euthanized due to poor prognosis and decline after 6 and 8 d, respectively. Simultaneously, an additional mara developed mild and self-resolving clinical signs, including a kyphotic gait and paraparesis. On gross examination, there were widespread petechiae and ecchymoses of the skeletal muscle, myocardium, skin, pericardium, urinary bladder mucosa, and spinal cord. On histopathology, all animals had necrotizing myelitis and rhombencephalitis, with intranuclear viral inclusions in three individuals. Electron microscopy confirmed herpesviral replication and assembly complexes in neurons and oligodendrocytes. Consensus PCR performed on spinal cord, brainstem, or cerebellum revealed a novel Simplexvirus most closely related to Simplexvirus leporidalpha 4. The virus was amplified and sequenced and is referred to as Simplexvirus dolichotinealpha1. It is unknown whether this virus is endemic in Patagonian mara or whether it represents an aberrant host species. Clinicians should be aware of this virus and its potential to cause severe, rapidly progressive, life-threatening disease in this species.

Nuclear interferon-stimulated gene product maintains heterochromatin on the herpes simplex viral genome to limit lytic infection

Interferons (IFN) are expressed in and secreted from cells in response to virus infection, and they induce the expression of a variety of genes called interferon-stimulated genes (ISGs) in infected and surrounding cells to block viral infection and limit spread. The mechanisms of action of a number of cytoplasmic ISGs have been well defined, but little is known about the mechanism of action of nuclear ISGs. Constitutive levels of nuclear interferon-inducible protein 16 (IFI16) serve to induce innate signaling and epigenetic silencing of herpes simplex virus (HSV), but only when the HSV infected cell protein 0 (ICP0) E3 ligase, which promotes IFI16 degradation, is inactivated. In this study, we found that following IFN induction, the pool of IFI16 within the infected cell remains high and can restrict wild-type viral gene expression and replication due to both the induced levels of IFI16 and the IFI16-mediated repression of ICP0 levels. Restriction of viral gene expression is achieved by IFI16 promoting the maintenance of heterochromatin on the viral genome, which silences it epigenetically. These results indicate that a nuclear ISG can restrict gene expression and replication of a nuclear DNA virus by maintaining or preventing the removal of repressive heterochromatin associated with the viral genome.

FREEDA: An automated computational pipeline guides experimental testing of protein innovation

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that lead to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA that provides a simple graphical user interface requiring only a gene name; integrates widely used molecular evolution tools to detect positive selection in rodents, primates, carnivores, birds, and flies; and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 centromere proteins, we find statistical evidence of positive selection within loops and turns of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of mouse CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

Responses to predation pressure involve similar sets of genes in two divergent species of Daphnia

Species that are not closely related can express similar inducible traits, but molecular mechanisms underlying the observed responses are often unknown, nor is it known if these mechanisms are shared between such species. Here, we compared transcriptional profiles of two Daphnia species (D. mitsukuri and D. sinensis) from different subgenera, at both juvenile and adult developmental stages. Both species were exposed to the same predation threat (fish kairomones), and both showed similar induced morphological changes (reduced body length). At the early developmental stage, response to predation risk resulted in similar changes in expression levels of 23 orthologues in both species. These orthologues, involved in 107 GO categories, changed in the same direction in both species (over- or underexpressed), in comparison to non-exposed controls. Several of these orthologues were associated with DNA replication, structural constituents of cuticle or innate immune response. In both species, the differentially expressed (DE) genes on average had higher ω (dN /dS ) values than non-DE genes, suggesting that these genes had experienced greater positive selection or lower purifying selection than non-DE genes. Overall, our results suggest that similar suites of genes, responding in similar ways to predation pressure, have been retained in Daphnia for many millions of years.

Homo medicus: The transition to meat eating increased pathogen pressure and the use of pharmacological plants in Homo

The human lineage transitioned to a more carnivorous niche 2.6 mya and evolved a large body size and slower life history, which likely increased zoonotic pathogen pressure. Evidence for this increase includes increased zoonotic infections in modern hunter-gatherers and bushmeat hunters, exceptionally low stomach pH compared to other primates, and divergence in immune-related genes. These all point to change, and probably intensification, in the infectious disease environment of Homo compared to earlier hominins and other apes. At the same time, the brain, an organ in which immune responses are constrained, began to triple in size. We propose that the combination of increased zoonotic pathogen pressure and the challenges of defending a large brain and body from pathogens in a long-lived mammal, selected for intensification of the plant-based self-medication strategies already in place in apes and other primates. In support, there is evidence of medicinal plant use by hominins in the middle Paleolithic, and all cultures today have sophisticated, plant-based medical systems, add spices to food, and regularly consume psychoactive plant substances that are harmful to helminths and other pathogens. We propose that the computational challenges of discovering effective plant-based treatments, the consequent ability to consume more energy-rich animal foods, and the reduced reliance on energetically-costly immune responses helped select for increased cognitive abilities and unique exchange relationships in Homo. In the story of human evolution, which has long emphasized hunting skills, medical skills had an equal role to play.

Evolution is not Uniform Along Coding Sequences

Amino acids evolve at different speeds within protein sequences, because their functional and structural roles are different. Notably, amino acids located at the surface of proteins are known to evolve more rapidly than those in the core. In particular, amino acids at the N- and C-termini of protein sequences are likely to be more exposed than those at the core of the folded protein due to their location in the peptidic chain, and they are known to be less structured. Because of these reasons, we would expect that amino acids located at protein termini would evolve faster than residues located inside the chain. Here we test this hypothesis and found that amino acids evolve almost twice as fast at protein termini compared with those in the center, hinting at a strong topological bias along the sequence length. We further show that the distribution of solvent-accessible residues and functional domains in proteins readily explain how structural and functional constraints are weaker at their termini, leading to the observed excess of amino acid substitutions. Finally, we show that the specific evolutionary rates at protein termini may have direct consequences, notably misleading in silico methods used to infer sites under positive selection within genes. These results suggest that accounting for positional information should improve evolutionary models.

FREEDA: an automated computational pipeline guides experimental testing of protein innovation by detecting positive selection

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that leads to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA (Finder of Rapidly Evolving Exons in De novo Assemblies) that provides a simple graphical user interface requiring only a gene name, integrates widely used molecular evolution tools to detect positive selection, and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 mouse centromere proteins, we find evidence of positive selection in intrinsically disordered regions of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals^1-6, probably owing to the evolutionary pressure on males to be reproductively successful⁷. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals.

Jack of all trades: Genome assembly of Wild Jack and comparative genomics of Artocarpus

Artocarpus (Moraceae), known as breadfruits for their diverse nutritious fruits, is prized for its high-quality timber, medicinal value, and economic importance. Breadfruits are native to Southeast Asia but have been introduced to other continents. The most commonly cultivated species are Artocarpus heterophyllus (Jackfruit) and Artocarpus altilis (Breadfruit). With numerous smaller but nutritionally comparable fruits on a larger tree, Artocarpus hirsutus, also called "Wild Jack" or "Ayani", is an elusive forest species endemic to Indian Western Ghats. In this study, we sequenced and assembled the whole genome of Artocarpus hirsutus sampled from the sacred groves of Coorg, India. To decipher demographic and evolutionary history, we compared our Wild Jack genome with previously published Jackfruit and Breadfruit genomes. Demographic history reconstruction indicates a stronger effect of habitat rather than phylogeny on the population histories of these plants. Repetitive genomic regions, especially LTR Copia, strongly affected the demographic trajectory of A. heterophyllus. Upon further investigation, we found a recent lineage-specific accumulation of LTR Copia in A. heterophyllus, which had a major contribution to its larger genome size. Several genes from starch, sucrose metabolism, and plant hormone signal transduction pathways, in Artocarpus species had signatures of selection and gene family evolution. Our comparative genomic framework provides important insights by incorporating endemic species such as the Wild Jack.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

Molecular archaeology of human cognitive traits

The brains and minds of our human ancestors remain inaccessible for experimental exploration. Therefore, we reconstructed human cognitive evolution by projecting nonsynonymous/synonymous rate ratios (ω values) in mammalian phylogeny onto the anatomically modern human (AMH) brain. This atlas retraces human neurogenetic selection and allows imputation of ancestral evolution in task-related functional networks (FNs). Adaptive evolution (high ω values) is associated with excitatory neurons and synaptic function. It shifted from FNs for motor control in anthropoid ancestry (60-41 mya) to attention in ancient hominoids (26-19 mya) and hominids (19-7.4 mya). Selection in FNs for language emerged with an early hominin ancestor (7.4-1.7 mya) and was later accompanied by adaptive evolution in FNs for strategic thinking during recent (0.8 mya-present) speciation of AMHs. This pattern mirrors increasingly complex cognitive demands and suggests that co-selection for language alongside strategic thinking may have separated AMHs from their archaic Denisovan and Neanderthal relatives.

Transcriptome innovations in primates revealed by single-molecule long-read sequencing

Transcriptomic diversity greatly contributes to the fundamentals of disease, lineage-specific biology, and environmental adaptation. However, much of the actual isoform repertoire contributing to shaping primate evolution remains unknown. Here, we combined deep long- and short-read sequencing complemented with mass spectrometry proteomics in a panel of lymphoblastoid cell lines (LCLs) from human, three other great apes, and rhesus macaque, producing the largest full-length isoform catalog in primates to date. Around half of the captured isoforms are not annotated in their reference genomes, significantly expanding the gene models in primates. Furthermore, our comparative analyses unveil hundreds of transcriptomic innovations and isoform usage changes related to immune function and immunological disorders. The confluence of these evolutionary innovations with signals of positive selection and their limited impact in the proteome points to changes in alternative splicing in genes involved in immune response as an important target of recent regulatory divergence in primates.

Rapidly evolving viral motifs mostly target biophysically constrained binding pockets of host proteins

Evolutionary changes in host-virus interactions can alter the course of infection, but the biophysical and regulatory constraints that shape interface evolution remain largely unexplored. Here, we focus on viral mimicry of host-like motifs that allow binding to host domains and modulation of cellular pathways. We observe that motifs from unrelated viruses preferentially target conserved, widely expressed, and highly connected host proteins, enriched with regulatory and essential functions. The interface residues within these host domains are more conserved and bind a larger number of cellular proteins than similar motif-binding domains that are not known to interact with viruses. In contrast, rapidly evolving viral-binding human proteins form few interactions with other cellular proteins and display high tissue specificity, and their interfaces have few inter-residue contacts. Our results distinguish between conserved and rapidly evolving host-virus interfaces and show how various factors limit host capacity to evolve, allowing for efficient viral subversion of host machineries.

Centromere drive: model systems and experimental progress

Centromeres connect chromosomes and spindle microtubules to ensure faithful chromosome segregation. Paradoxically, despite this conserved function, centromeric DNA evolves rapidly and centromeric proteins show signatures of positive selection. The centromere drive hypothesis proposes that centromeric DNA can act like a selfish genetic element and drive non-Mendelian segregation during asymmetric female meiosis. Resulting fitness costs lead to genetic conflict with the rest of the genome and impose a selective pressure for centromeric proteins to adapt by suppressing the costs. Here, we describe experimental model systems for centromere drive in yellow monkeyflowers and mice, summarize key findings demonstrating centromere drive, and explain molecular mechanisms. We further discuss efforts to test if centromeric proteins are involved in suppressing drive-associated fitness costs, highlight a model for centromere drive and suppression in mice, and put forth outstanding questions for future research.

A functional map of HIV-host interactions in primary human T cells

Human Immunodeficiency Virus (HIV) relies on host molecular machinery for replication. Systematic attempts to genetically or biochemically define these host factors have yielded hundreds of candidates, but few have been functionally validated in primary cells. Here, we target 426 genes previously implicated in the HIV lifecycle through protein interaction studies for CRISPR-Cas9-mediated knock-out in primary human CD4+ T cells in order to systematically assess their functional roles in HIV replication. We achieve efficient knockout (>50% of alleles) in 364 of the targeted genes and identify 86 candidate host factors that alter HIV infection. 47 of these factors validate by multiplex gene editing in independent donors, including 23 factors with restrictive activity. Both gene editing efficiencies and HIV-1 phenotypes are highly concordant among independent donors. Importantly, over half of these factors have not been previously described to play a functional role in HIV replication, providing numerous novel avenues for understanding HIV biology. These data further suggest that host-pathogen protein-protein interaction datasets offer an enriched source of candidates for functional host factor discovery and provide an improved understanding of the mechanics of HIV replication in primary T cells.

OUP accepted manuscript

Recognition of mammalian mitochondrial promoters requires the concerted action of mitochondrial RNA polymerase (mtRNAP) and transcription initiation factors TFAM and TFB2M. In this work, we found that transcript slippage results in heterogeneity of the human mitochondrial transcripts in vivo and in vitro. This allowed us to correctly interpret the RNAseq data, identify the bona fide transcription start sites (TSS), and assign mitochondrial promoters for > 50% of mammalian species and some other vertebrates. The divergent structure of the mammalian promoters reveals previously unappreciated aspects of mtDNA evolution. The correct assignment of TSS also enabled us to establish the precise register of the DNA in the initiation complex and permitted investigation of the sequence-specific protein-DNA interactions. We determined the molecular basis of promoter recognition by mtRNAP and TFB2M, which cooperatively recognize bases near TSS in a species-specific manner. Our findings reveal a role of mitochondrial transcription machinery in mitonuclear coevolution and speciation.

Recurrent erosion of COA1/MITRAC15 exemplifies conditional gene dispensability in oxidative phosphorylation

Skeletal muscle fibers rely upon either oxidative phosphorylation or the glycolytic pathway with much less reliance on oxidative phosphorylation to achieve muscular contractions that power mechanical movements. Species with energy-intensive adaptive traits that require sudden bursts of energy have a greater dependency on glycolytic fibers. Glycolytic fibers have decreased reliance on OXPHOS and lower mitochondrial content compared to oxidative fibers. Hence, we hypothesized that gene loss might have occurred within the OXPHOS pathway in lineages that largely depend on glycolytic fibers. The protein encoded by the COA1/MITRAC15 gene with conserved orthologs found in budding yeast to humans promotes mitochondrial translation. We show that gene disrupting mutations have accumulated within the COA1 gene in the cheetah, several species of galliform birds, and rodents. The genomic region containing COA1 is a well-established evolutionary breakpoint region in mammals. Careful inspection of genome assemblies of closely related species of rodents and marsupials suggests two independent COA1 gene loss events co-occurring with chromosomal rearrangements. Besides recurrent gene loss events, we document changes in COA1 exon structure in primates and felids. The detailed evolutionary history presented in this study reveals the intricate link between skeletal muscle fiber composition and the occasional dispensability of the chaperone-like role of the COA1 gene.

Functional Diversification of Chromatin on Rapid Evolutionary Timescales

Repeat-enriched genomic regions evolve rapidly and yet support strictly conserved functions like faithful chromosome transmission and the preservation of genome integrity. The leading resolution to this paradox is that DNA repeat-packaging proteins evolve adaptively to mitigate deleterious changes in DNA repeat copy number, sequence, and organization. Exciting new research has tested this model of coevolution by engineering evolutionary mismatches between adaptively evolving chromatin proteins of one species and the DNA repeats of a close relative. Here, we review these innovative evolution-guided functional analyses. The studies demonstrate that vital, chromatin-mediated cellular processes, including transposon suppression, faithful chromosome transmission, and chromosome retention depend on species-specific versions of chromatin proteins that package species-specific DNA repeats. In many cases, the ever-evolving repeats are selfish genetic elements, raising the possibility that chromatin is a battleground of intragenomic conflict.

Comparable Number of Genes Having Experienced Positive Selection among Great Ape Species

Simple Summary

It is of great interest to quantify adaptive evolution in human lineage by studying genes under positive selection, since these genes could reveal insights into our own adaptive evolutionary history compared to our closely related species and often these genes are functionally important. We used the great apes as the subjects to detect gene-level adaptive evolution signals in all the great ape lineages and investigated the evolutionary patterns and functional relevance of these adaptive evolution signals. Even the differences in population size among these closely related great apes have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments, we found that they experienced comparable numbers of positive selection. Notably, we identified several genes that offer insights into great ape and human evolution. For example, SOD1, a gene associated with aging in humans, experienced positive selection in the common ancestor of the great ape and this positive selection may contribute to the aging evolution in great apes. Overall, an updated list of positively selected genes reported by this study not only informs us of adaptive evolution during great ape evolution, but is also helpful to the further study of non-human primate models for disease and other fields.

Abstract

Alleles that cause advantageous phenotypes with positive selection contribute to adaptive evolution. Investigations of positive selection in protein-coding genes rely on the accuracy of orthology, models, the quality of assemblies, and alignment. Here, based on the latest genome assemblies and gene annotations, we present a comparative analysis on positive selection in four great ape species and identify 211 high-confidence positively selected genes (PSGs). Even the differences in population size among these closely related great apes have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments, we found that they experienced comparable numbers of positive selection. We also uncovered that more than half of multigene families exhibited signals of positive selection, suggesting that imbalanced positive selection resulted in the functional divergence of duplicates. Moreover, at the expression level, although positive selection led to a more non-uniform pattern across tissues, the correlation between positive selection and expression patterns is diverse. Overall, this updated list of PSGs is of great significance for the further study of the phenotypic evolution in great apes.

Evolutionary mode for the functional preservation of fast-evolving Drosophila telomere capping proteins

DNA end protection is fundamental for the long-term preservation of the genome. In vertebrates the Shelterin protein complex protects telomeric DNA ends, thereby contributing to the maintenance of genome integrity. In the Drosophila genus, this function is thought to be performed by the Terminin complex, an assembly of fast-evolving subunits. Considering that DNA end protection is fundamental for successful genome replication, the accelerated evolution of Terminin subunits is counterintuitive, as conservation is supposed to maintain the assembly and concerted function of the interacting partners. This problem extends over Drosophila telomere biology and provides insight into the evolution of protein assemblies. In order to learn more about the mechanistic details of this phenomenon we have investigated the intra- and interspecies assemblies of Verrocchio and Modigliani, two Terminin subunits using in vitro assays. Based on our results and on homology-based three-dimensional models for Ver and Moi, we conclude that both proteins contain Ob-fold and contribute to the ssDNA binding of the Terminin complex. We propose that the preservation of Ver function is achieved by conservation of specific amino acids responsible for folding or localized in interacting surfaces. We also provide here the first evidence on Moi DNA binding.

Parallel pathways for recruiting effector proteins determine centromere drive and suppression

Selfish centromere DNA sequences bias their transmission to the egg in female meiosis. Evolutionary theory suggests that centromere proteins evolve to suppress costs of this "centromere drive." In hybrid mouse models with genetically different maternal and paternal centromeres, selfish centromere DNA exploits a kinetochore pathway to recruit microtubule-destabilizing proteins that act as drive effectors. We show that such functional differences are suppressed by a parallel pathway for effector recruitment by heterochromatin, which is similar between centromeres in this system. Disrupting the kinetochore pathway with a divergent allele of CENP-C reduces functional differences between centromeres, whereas disrupting heterochromatin by CENP-B deletion amplifies the differences. Molecular evolution analyses using Murinae genomes identify adaptive evolution in proteins in both pathways. We propose that centromere proteins have recurrently evolved to minimize the kinetochore pathway, which is exploited by selfish DNA, relative to the heterochromatin pathway that equalizes centromeres, while maintaining essential functions.

The antiviral state has shaped the CpG composition of the vertebrate interferome to avoid self-targeting

Antiviral defenses can sense viral RNAs and mediate their destruction. This presents a challenge for host cells since they must destroy viral RNAs while sparing the host mRNAs that encode antiviral effectors. Here, we show that highly upregulated interferon-stimulated genes (ISGs), which encode antiviral proteins, have distinctive nucleotide compositions. We propose that self-targeting by antiviral effectors has selected for ISG transcripts that occupy a less self-targeted sequence space. Following interferon (IFN) stimulation, the CpG-targeting antiviral effector zinc-finger antiviral protein (ZAP) reduces the mRNA abundance of multiple host transcripts, providing a mechanistic explanation for the repression of many (but not all) interferon-repressed genes (IRGs). Notably, IRGs tend to be relatively CpG rich. In contrast, highly upregulated ISGs tend to be strongly CpG suppressed. Thus, ZAP is an example of an effector that has not only selected compositional biases in viral genomes but also appears to have notably shaped the composition of host transcripts in the vertebrate interferome.

Our cells are poised to combat viral infection through antiviral effectors. This study proposes that as well as targeting viral RNAs, antiviral effectors sometimes target host mRNAs too; over millions of years, this has selected for compositional biases in the host’s transcriptional response to virus infection.

Convergent Loss of the Necroptosis Pathway in Disparate Mammalian Lineages Shapes Viruses Countermeasures

Programmed cell death is a vital process in the life cycle of organisms. Necroptosis, an evolutionary form of programmed necrosis, contributes to the innate immune response by killing pathogen-infected cells. This virus-host interaction pathway is organized around two components: the receptor-interacting protein kinase 3 (RIPK3), which recruits and phosphorylates the mixed lineage kinase-like protein (MLKL), inducing cellular plasma membrane rupture and cell death. Critically, the presence of necroptotic inhibitors in viral genomes validates necroptosis as an important host defense mechanism. Here, we show, counterintuitively, that in different mammalian lineages, central components of necroptosis, such as RIPK3 and MLKL, are deleted or display inactivating mutations. Frameshifts or premature stop codons are observed in all the studied species of cetaceans and leporids. In carnivores' genomes, the MLKL gene is deleted, while in a small number of species from afrotheria and rodentia premature stop codons are observed in RIPK3 and/or MLKL. Interestingly, we also found a strong correlation between the disruption of necroptosis in leporids and cetaceans and the absence of the N-terminal domain of E3-like homologs (responsible for necroptosis inhibition) in their naturally infecting poxviruses. Overall, our study provides the first comprehensive picture of the molecular evolution of necroptosis in mammals. The loss of necroptosis multiple times during mammalian evolution highlights the importance of gene/pathway loss for species adaptation and suggests that necroptosis is not required for normal mammalian development. Moreover, this study highlights a co-evolutionary relationship between poxviruses and their hosts, emphasizing the role of host adaptation in shaping virus evolution.

Evolutionary Dynamics of the Spindle Assembly Checkpoint in Eukaryotes

The tremendous diversity in eukaryotic life forms can ultimately be traced back to evolutionary modifications at the level of molecular networks. Deep understanding of these modifications will not only explain cellular diversity, but will also uncover different ways to execute similar processes and expose the evolutionary 'rules' that shape the molecular networks. Here, we review the evolutionary dynamics of the spindle assembly checkpoint (SAC), a signaling network that guards fidelity of chromosome segregation. We illustrate how the interpretation of divergent SAC systems in eukaryotic species is facilitated by combining detailed molecular knowledge of the SAC and extensive comparative genome analyses. Ultimately, expanding this to other core cellular systems and experimentally interrogating such systems in organisms from all major lineages may start outlining the routes to and eventual manifestation of the cellular diversity of eukaryotic life.

Diversification of CD1 Molecules Shapes Lipid Antigen Selectivity

Molecular studies of host-pathogen evolution have largely focused on the consequences of variation at protein-protein interaction surfaces. The potential for other microbe-associated macromolecules to promote arms race dynamics with host factors remains unclear. The cluster of differentiation 1 (CD1) family of vertebrate cell surface receptors plays a crucial role in adaptive immunity through binding and presentation of lipid antigens to T-cells. Although CD1 proteins present a variety of endogenous and microbial lipids to various T-cell types, they are less diverse within vertebrate populations than the related major histocompatibility complex (MHC) molecules. We discovered that CD1 genes exhibit a high level of divergence between simian primate species, altering predicted lipid-binding properties and T-cell receptor interactions. These findings suggest that lipid-protein conflicts have shaped CD1 genetic variation during primate evolution. Consistent with this hypothesis, multiple primate CD1 family proteins exhibit signatures of repeated positive selection at surfaces impacting antigen presentation, binding pocket morphology, and T-cell receptor accessibility. Using a molecular modeling approach, we observe that interspecies variation as well as single mutations at rapidly-evolving sites in CD1a drastically alter predicted lipid binding and structural features of the T-cell recognition surface. We further show that alterations in both endogenous and microbial lipid-binding affinities influence the ability of CD1a to undergo antigen swapping required for T-cell activation. Together these findings establish lipid-protein interactions as a critical force of host-pathogen conflict and inform potential strategies for lipid-based vaccine development.

The genome of a thorny species: comparative genomic analysis among South and North American Cactaceae

The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.

Positive Selection and Gene Expression Analyses from Salivary Glands Reveal Discrete Adaptations within the Ecologically Diverse Bat Family Phyllostomidae

The leaf-nosed bats (Phyllostomidae) are outliers among chiropterans with respect to the unusually high diversity of dietary strategies within the family. Salivary glands, owing to their functions and high ultrastructural variability among lineages, are proposed to have played an important role during the phyllostomid radiation. To identify genes underlying salivary gland functional diversification, we sequenced submandibular gland transcriptomes from phyllostomid species representative of divergent dietary strategies. From the assembled transcriptomes, we performed an array of selection tests and gene expression analyses to identify signatures of adaptation. Overall, we identified an enrichment of immunity-related gene ontology terms among 53 genes evolving under positive selection. Lineage-specific selection tests revealed several endomembrane system genes under selection in the vampire bat. Many genes that respond to insulin were under selection and differentially expressed genes pointed to modifications of amino acid synthesis pathways in plant-visitors. Results indicate salivary glands have diversified in various ways across a functional diverse clade of mammals in response to niche specializations.

Hearing Sensitivity of Primates: Recurrent and Episodic Positive Selection in Hair Cells and Stereocilia Protein-Coding Genes

The large spectrum of hearing sensitivity observed in primates results from the impact of environmental and behavioral pressures to optimize sound perception and localization. Although evidence of positive selection in auditory genes has been detected in mammals including in Hominoids, selection has never been investigated in other primates. We analyzed 123 genes highly expressed in the inner ear of 27 primate species and tested to what extent positive selection may have shaped these genes in the order Primates tree. We combined both site and branch-site tests to obtain a comprehensive picture of the positively selected genes (PSGs) involved in hearing sensitivity, and drew a detailed description of the most affected branches in the tree. We chose a conservative approach, and thus focused on confounding factors potentially affecting PSG signals (alignment, GC-biased gene conversion, duplications, heterogeneous sequencing qualities). Using site tests, we showed that around 12% of these genes are PSGs, an α selection value consistent with average human genome estimates (10-15%). Using branch-site tests, we showed that the primate tree is heterogeneously affected by positive selection, with the black snub-nosed monkey, the bushbaby, and the orangutan, being the most impacted branches. A large proportion of these genes is inclined to shape hair cells and stereocilia, which are involved in the mechanotransduction process, known to influence frequency perception. Adaptive selection, and more specifically recurrent adaptive evolution, could have acted in parallel on a set of genes (ADGRV1, USH2A, PCDH15, PTPRQ, and ATP8A2) involved in stereocilia growth and the whole complex of bundle links connecting them, in species across different habitats, including high altitude and nocturnal environments.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

Primate innate immune responses to bacterial and viral pathogens reveals an evolutionary trade-off between strength and specificity

Humans and our close evolutionary relatives respond differently to a large number of infections. Such differences are thought to result, at least in part, from interspecies differences in immune function. Here, we report on the whole-genome expression of blood leukocytes from four primate species responding to bacterial and viral stimulation. We show that apes mount a markedly stronger early transcriptional response to both viral and bacterial stimulation when compared to African and Asian monkeys. In addition, our findings suggest that apes activate a broader array of defense molecules that may be beneficial for early pathogen killing at the potential cost of increased energy expenditure and tissue damage. Our results provide insight into the evolution of immune responses in primates.

Despite their close genetic relatedness, apes and African and Asian monkeys (AAMs) differ in their susceptibility to severe bacterial and viral infections that are important causes of human disease. Such differences between humans and other primates are thought to be a result, at least in part, of interspecies differences in immune response to infection. However, because of the lack of comparative functional data across species, it remains unclear in what ways the immune systems of humans and other primates differ. Here, we report the whole-genome transcriptomic responses of ape species (human and chimpanzee) and AAMs (rhesus macaque and baboon) to bacterial and viral stimulation. We find stark differences in the responsiveness of these groups, with apes mounting a markedly stronger early transcriptional response to both viral and bacterial stimulation, altering the transcription of ∼40% more genes than AAMs. Additionally, we find that genes involved in the regulation of inflammatory and interferon responses show the most divergent early transcriptional responses across primates and that this divergence is attenuated over time. Finally, we find that relative to AAMs, apes engage a much less specific immune response to different classes of pathogens during the early hours of infection, up-regulating genes typical of anti-viral and anti-bacterial responses regardless of the nature of the stimulus. Overall, these findings suggest apes exhibit increased sensitivity to bacterial and viral immune stimulation, activating a broader array of defense molecules that may be beneficial for early pathogen killing at the potential cost of increased energy expenditure and tissue damage.

Evolutionary genetics and acclimatization in nephrology

Evolutionary processes, including mutation, migration and natural selection, have influenced the prevalence and distribution of various disorders in humans. However, despite a few well-known examples, such as the APOL1 variants - which have undergone positive genetic selection for their ability to confer resistance to Trypanosoma brucei infection but confer a higher risk of chronic kidney disease - little is known about the effects of evolutionary processes that have shaped genetic variation on kidney disease. An understanding of basic concepts in evolutionary genetics provides an opportunity to consider how findings from ancient and archaic genomes could inform our knowledge of evolution and provide insights into how population migration and genetic admixture have shaped the current distribution and landscape of human kidney-associated diseases. Differences in exposures to infectious agents, environmental toxins, dietary components and climate also have the potential to influence the evolutionary genetics of kidneys. Of note, selective pressure on loci associated with kidney disease is often from non-kidney diseases, and thus it is important to understand how the link between genome-wide selected loci and kidney disease occurs in relation to secondary nephropathies.

Retroviral Restriction Factors and Their Viral Targets: Restriction Strategies and Evolutionary Adaptations

The evolutionary conflict between retroviruses and their vertebrate hosts over millions of years has led to the emergence of cellular innate immune proteins termed restriction factors as well as their viral antagonists. Evidence accumulated in the last two decades has substantially increased our understanding of the elaborate mechanisms utilized by these restriction factors to inhibit retroviral replication, mechanisms that either directly block viral proteins or interfere with the cellular pathways hijacked by the viruses. Analyses of these complex interactions describe patterns of accelerated evolution for these restriction factors as well as the acquisition and evolution of their virus-encoded antagonists. Evidence is also mounting that many restriction factors identified for their inhibition of specific retroviruses have broader antiviral activity against additional retroviruses as well as against other viruses, and that exposure to these multiple virus challenges has shaped their adaptive evolution. In this review, we provide an overview of the restriction factors that interfere with different steps of the retroviral life cycle, describing their mechanisms of action, adaptive evolution, viral targets and the viral antagonists that evolved to counter these factors.

Richard Dickerson, Molecular Clocks, and Rates of Protein Evolution

Proteins approximately behave as molecular clocks, accumulating amino acid replacements at a more or less constant rate. Nonetheless, each protein displays a characteristic rate of evolution: whereas some proteins remain largely unaltered over large periods of time, others can rapidly accumulate amino acid replacements. An article by Richard Dickerson, published in the first issue of the Journal of Molecular Evolution (J Mol Evol 1:26-45, 1971), described the first analysis in which the rates of evolution of many proteins were compared, and the differences were interpreted in the light of their function. When comparing the sequences of fibrinopeptides, hemoglobin, and cytochrome c of different species, he observed a linear relationship between the number of amino acid replacements and divergence time. Remarkably, fibrinopeptides had evolved fast, cytochrome c had evolved slowly, and hemoglobin exhibited an intermediate rate of evolution. As the Journal of Molecular Evolution celebrates its 50th anniversary, I highlight this landmark article and reflect on its impact on the field of Molecular Evolution.

DGINN, an automated and highly-flexible pipeline for the detection of genetic innovations on protein-coding genes

Adaptive evolution has shaped major biological processes. Finding the protein-coding genes and the sites that have been subjected to adaptation during evolutionary time is a major endeavor. However, very few methods fully automate the identification of positively selected genes, and widespread sources of genetic innovations such as gene duplication and recombination are absent from most pipelines. Here, we developed DGINN, a highly-flexible and public pipeline to Detect Genetic INNovations and adaptive evolution in protein-coding genes. DGINN automates, from a gene's sequence, all steps of the evolutionary analyses necessary to detect the aforementioned innovations, including the search for homologs in databases, assignation of orthology groups, identification of duplication and recombination events, as well as detection of positive selection using five methods to increase precision and ranking of genes when a large panel is analyzed. DGINN was validated on nineteen genes with previously-characterized evolutionary histories in primates, including some engaged in host-pathogen arms-races. Our results confirm and also expand results from the literature, including novel findings on the Guanylate-binding protein family, GBPs. This establishes DGINN as an efficient tool to automatically detect genetic innovations and adaptive evolution in diverse datasets, from the user's gene of interest to a large gene list in any species range.

Relaxed constraint and functional divergence of the progesterone receptor (PGR) in the human stem-lineage

The steroid hormone progesterone, acting through the progesterone receptor (PR), a ligand-activated DNA-binding transcription factor, plays an essential role in regulating nearly every aspect of female reproductive biology. While many reproductive traits regulated by PR are conserved in mammals, Catarrhine primates evolved several derived traits including spontaneous decidualization, menstruation, and a divergent (and unknown) parturition signal, suggesting that PR may also have evolved divergent functions in Catarrhines. There is conflicting evidence, however, whether the progesterone receptor gene (PGR) was positively selected in the human lineage. Here we show that PGR evolved rapidly in the human stem-lineage (as well as other Catarrhine primates), which likely reflects an episode of relaxed selection intensity rather than positive selection. Coincident with the episode of relaxed selection intensity, ancestral sequence resurrection and functional tests indicate that the major human PR isoforms (PR-A and PR-B) evolved divergent functions in the human stem-lineage. These results suggest that the regulation of progesterone signaling by PR-A and PR-B may also have diverged in the human lineage and that non-human animal models of progesterone signaling may not faithfully recapitulate human biology.

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

Telomeres ensure chromosome length homeostasis and protection from catastrophic end-to-end chromosome fusions. All eukaryotes require this essential, strictly conserved telomere-dependent genome preservation. However, recent evolutionary analyses of mammals, plants, and flies report pervasive rapid evolution of telomere proteins. The causes of this paradoxical observation - that unconserved machinery underlies an essential, conserved function - remain enigmatic. Indeed, these fast-evolving telomere proteins bind, extend, and protect telomeric DNA, which itself evolves slowly in most systems. We hypothesize that the universally fast-evolving subtelomere - the telomere-adjacent, repetitive sequence - is a primary driver of the 'telomere paradox'. Under this model, radical sequence changes in the subtelomere perturb subtelomere-dependent, telomere functions. Compromised telomere function then spurs adaptation of telomere proteins to maintain telomere length homeostasis and protection. We propose an experimental framework that leverages both protein divergence and subtelomeric sequence divergence to test the hypothesis that subtelomere sequence evolution shapes recurrent innovation of telomere machinery.

Prisoners of war - host adaptation and its constraints on virus evolution

Recent discoveries of contemporary genotypes of hepatitis B virus and parvovirus B19 in ancient human remains demonstrate that little genetic change has occurred in these viruses over 4,500-6,000 years. Endogenous viral elements in host genomes provide separate evidence that viruses similar to many major contemporary groups circulated 100 million years ago or earlier. In this Opinion article, we argue that the extraordinary conservation of virus genome sequences is best explained by a niche-filling model in which fitness optimization is rapidly achieved in their specific hosts. Whereas short-term substitution rates reflect the accumulation of tolerated sequence changes within adapted genomes, longer-term rates increasingly resemble those of their hosts as the evolving niche moulds and effectively imprisons the virus in co-adapted virus-host relationships. Contrastingly, viruses that jump hosts undergo strong and stringent adaptive selection as they maximize their fit to their new niche. This adaptive capability may paradoxically create evolutionary stasis in long-term host relationships. While viruses can evolve and adapt rapidly, their hosts may ultimately shape their longer-term evolution.

Cytokines and chemokines involved in the defense reaction against HIV-1 and hepatitis B virus: isn't it time to use a standardized nomenclature of the involved mediators?

Discovery of new mediators of immune cell activation and interaction facilitated elucidation of the various ways of defense against infectious agents and happened some 40 years ago. Each involved group of researchers named the mediators according to their scope of investigation; often the same molecules were published at the same time with different names. To avoid confusion resulting from using different names for the same mediators and to prevent a Babylonian confusion, standardization was implemented-as in the field of metrics, music, or science including virology. For cytokines and chemokines a standard nomenclature was proposed some 10 years ago and in conclusion it should be used. In this paper the most relevant biomarkers in HIV-1 and HBV infection and their contribution during viral pathogenesis are listed.

A metaanalysis of bat phylogenetics and positive selection based on genomes and transcriptomes from 18 species

This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation.

Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

Tissue-specific genes as an underutilized resource in drug discovery

Tissue-specific genes are believed to be good drug targets due to improved safety. Here we show that this intuitive notion is not reflected in phase 1 and 2 clinical trials, despite the historic success of tissue-specific targets and their 2.3-fold overrepresentation among targets of marketed non-oncology drugs. We compare properties of tissue-specific genes and drug targets. We show that tissue-specificity of the target may also be related to efficacy of the drug. The relationship may be indirect (enrichment in Mendelian disease and PTVesc genes) or direct (elevated betweenness centrality scores for tissue-specifically produced enzymes and secreted proteins). Reduced evolutionary conservation of tissue-specific genes may represent a bottleneck for drug projects, prompting development of novel models with smaller evolutionary gap to humans. We show that the opportunities to identify tissue-specific drug targets are not exhausted and discuss potential use cases for tissue-specific genes in drug research.

Does mitochondrial DNA evolution in metazoa drive the origin of new mitochondrial proteins?

Most eukaryotic cells contain mitochondria with a genome that evolved from their α-proteobacterial ancestor. In the course of eukaryotic evolution, the mitochondrial genome underwent a dramatic reduction in size, caused by the loss and translocation of genes. This required adjustments in mitochondrial gene expression mechanisms and resulted in a complex collaborative system of mitochondrially encoded transfer RNAs and ribosomal RNAs with nuclear encoded proteins to express the mitochondrial encoded oxidative phosphorylation (OXPHOS) proteins. In this review, we examine mitochondrial gene expression from an evolutionary point of view: to what extent can we correlate changes in the mitochondrial genome in the evolutionary lineage leading to human with the origin of new nuclear encoded proteins. We dated the evolutionary origin of mitochondrial proteins that interact with mitochondrial DNA or its RNA and/or protein products in a systematic manner and compared them with documented changes in the mitochondrial DNA. We find anecdotal but accumulating evidence that metazoan RNA-interacting proteins arose in conjunction with changes of the mitochondrial DNA. We find no substantial evidence for such compensatory evolution in new OXPHOS proteins, which appear to be constrained by the ability to form supercomplexes. © 2018 IUBMB Life, 70(12):1240-1250, 2018.

Genome Sequence of Peacock Reveals the Peculiar Case of a Glittering Bird

The unique ornamental features and extreme sexual traits of Peacock have always intrigued scientists and naturalists for centuries. However, the genomic basis of these phenotypes are yet unknown. Here, we report the first genome sequence and comparative analysis of peacock with the high quality genomes of chicken, turkey, duck, flycatcher and zebra finch. Genes involved in early developmental pathways including TGF-β, BMP, and Wnt signaling, which have been shown to be involved in feather patterning, bone morphogenesis, and skeletal muscle development, revealed signs of adaptive evolution and provided useful clues on the phenotypes of peacock. Innate and adaptive immune genes involved in complement system and T-cell response also showed signs of adaptive evolution in peacock suggesting their possible role in building a robust immune system which is consistent with the predictions of the Hamilton–Zuk hypothesis. This study provides novel genomic and evolutionary insights into the molecular understanding toward the phenotypic evolution of Indian peacock.

Positive selection on human gamete-recognition genes

Coevolution of genes that encode interacting proteins expressed on the surfaces of sperm and eggs can lead to variation in reproductive compatibility between mates and reproductive isolation between members of different species. Previous studies in mice and other mammals have focused in particular on evidence for positive or diversifying selection that shapes the evolution of genes that encode sperm-binding proteins expressed in the egg coat or zona pellucida (ZP). By fitting phylogenetic models of codon evolution to data from the 1000 Genomes Project, we identified candidate sites evolving under diversifying selection in the human genes ZP3 and ZP2. We also identified one candidate site under positive selection in C4BPA, which encodes a repetitive protein similar to the mouse protein ZP3R that is expressed in the sperm head and binds to the ZP at fertilization. Results from several additional analyses that applied population genetic models to the same data were consistent with the hypothesis of selection on those candidate sites leading to coevolution of sperm- and egg-expressed genes. By contrast, we found no candidate sites under selection in a fourth gene (ZP1) that encodes an egg coat structural protein not directly involved in sperm binding. Finally, we found that two of the candidate sites (in C4BPA and ZP2) were correlated with variation in family size and birth rate among Hutterite couples, and those two candidate sites were also in linkage disequilibrium in the same Hutterite study population. All of these lines of evidence are consistent with predictions from a previously proposed hypothesis of balancing selection on epistatic interactions between C4BPA and ZP3 at fertilization that lead to the evolution of co-adapted allele pairs. Such patterns also suggest specific molecular traits that may be associated with both natural reproductive variation and clinical infertility.