Explore the vast range of titles available.

There is an increasing interest in detecting genes, or genomic regions, that have been targeted by natural selection. The interest stems from a basic desire to learn more about evolutionary processes in humans and other organisms, and from the realization that inferences regarding selection may provide important functional information. This review provides a nonmathematical description of the issues involved in detecting selection from DNA sequences and SNP data and is intended for readers who are not familiar with population genetic theory. Particular attention is placed on issues relating to the analysis of large-scale genomic data sets.

Background A number of different statistics are used for detecting natural selection using DNA sequencing data, including statistics that are summaries of the frequency spectrum, such as Tajima’s D . These statistics are now often being applied in the analysis of Next Generation Sequencing (NGS) data. However, estimates of frequency spectra from NGS data are strongly affected by low sequencing coverage; the inherent technology dependent variation in sequencing depth causes systematic differences in the value of the statistic among genomic regions. Results We have developed an approach that accommodates the uncertainty of the data when calculating site frequency based neutrality test statistics. A salient feature of this approach is that it implicitly solves the problems of varying sequencing depth, missing data and avoids the need to infer variable sites for the analysis and thereby avoids ascertainment problems introduced by a SNP discovery process. Conclusion Using an empirical Bayes approach for fast computations, we show that this method produces results for low-coverage NGS data comparable to those achieved when the genotypes are known without uncertainty. We also validate the method in an analysis of data from the 1000 genomes project. The method is implemented in a fast framework which enables researchers to perform these neutrality tests on a genome-wide scale.

Genomic imprinting is an epigenetic phenomenon where autosomal genes display uniparental expression depending on whether they are maternally or paternally inherited. Genomic imprinting can arise from parental conflicts over resource allocation to the offspring, which could drive imprinted loci to evolve by positive selection. We investigate whether positive selection is associated with genomic imprinting in the inbreeding species Arabidopsis thaliana. Our analysis of 140 genes regulated by genomic imprinting in the A. thaliana seed endosperm demonstrates they are evolving more rapidly than expected. To investigate whether positive selection drives this evolutionary acceleration, we identified orthologs of each imprinted gene across 34 plant species and elucidated their evolutionary trajectories. Increased positive selection was sought by comparing its incidence among imprinted genes with nonimprinted controls. Strikingly, we find a statistically significant enrichment of imprinted paternally expressed genes (iPEGs) evolving under positive selection, 50.6% of the total, but no such enrichment for positive selection among imprinted maternally expressed genes (iMEGs). This suggests that maternally- and paternally expressed imprinted genes are subject to different selective pressures. Almost all positively selected amino acids were fixed across 80 sequenced A. thaliana accessions, suggestive of selective sweeps in the A. thaliana lineage. The imprinted genes under positive selection are involved in processes important for seed development including auxin biosynthesis and epigenetic regulation. Our findings support a genomic imprinting model for plants where positive selection can affect paternally expressed genes due to continued conflict with maternal sporophyte tissues, even when parental conflict is reduced in predominantly inbreeding species.

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

The emerging and ongoing outbreak of human monkeypox (hMPX) in 2022 is a serious global threat. An understanding of the evolution of the monkeypox virus (MPXV) at the single-gene level may provide clues for exploring the unique aspects of the current outbreak: rapidly expanding and sustained human-to-human transmission. For the current investigation, alleles of 156 MPXV coding genes (which account for >95% of the genomic sequence) have been gathered from roughly 1,500 isolates, including those responsible for the previous outbreaks. Using a range of molecular evolution approaches, we demonstrated that intra-species homologous recombination has a negligible effect on MPXV evolution. Despite the fact that the majority of the MPXV genes (64.10%) were subjected to negative selection at the whole gene level, 10 MPXV coding genes (MPXVgp004, 010, 012, 014, 044, 098, 138, 178, 188, and 191) were found to have a total of 15 codons or amino acid sites that are known to evolve under positive Darwinian selection. Except for MPXVgp138, almost all of these genes encode proteins that interact with the host. Of these, five ankyrin proteins (MPXVgp004, 010, 012, 178, and 188) and one Bcl-2-like protein (MPXVgp014) are involved in poxviruses’ host range determination. We discovered that the majority (80%) of positive amino acid substitutions emerged several decades ago, indicating that these sites have been under constant selection pressure and that more adaptable alleles have been circulating in the natural reservoir. This finding was also supported by the minimum spanning networks of the gene alleles. The three positive amino acid substitutions (T/A426V in MPXVgp010, A423D in MPXVgp012, and S105L in MPXVgp191) appeared in 2019 or 2022, indicating that they would be crucial for the virus’ eventual adaptation to humans. Protein modeling suggests that positive amino acid substitutions may affect protein functions in a variety of ways. Further study should focus on revealing the biological effects of positive amino acid substitutions in the genes for viral adaptation to humans, virulence, transmission, and so on. Our study advances knowledge of MPXV’s adaptive mechanism and provides insights for exploring factors that are responsible for the unique aspects of the current outbreak.

The simplest evolutionary model for catching prey by an agent (predator) is considered. The simulation is performed on the basis of a software-emulated Intel i8080 processor. Maximizing the number of catches is chosen as the objective function. This function is associated with energy dissipation and developmental time. It is shown that during Darwinian evolution, agents with an initially a random set of processor commands subsequently acquire a successful catching skill. It is found that in the process of evolution, a logarithmic relationship between astronomical and developmental times arises in agents. This result is important for the ideas available in the literature about the close connection of such concepts as time, Darwinian selection, and the maximization of entropy production.

Schizophrenia poses an evolutionary-genetic paradox because it exhibits strongly negative fitness effects and high heritability, yet it persists at a prevalence of approximately 1% across all human cultures. Recent theory has proposed a resolution: that genetic liability to schizophrenia has evolved as a secondary consequence of selection for human cognitive traits. This hypothesis predicts that genes increasing the risk of this disorder have been subject to positive selection in the evolutionary history of humans and other primates. We evaluated this prediction using tests for recent selective sweeps in human populations and maximum-likelihood tests for selection during primate evolution. Significant evidence for positive selection was evident using one or both methods for 28 of 76 genes demonstrated to mediate liability to schizophrenia, including DISC1, DTNBP1 and NRG1, which exhibit especially strong and well-replicated functional and genetic links to this disorder. Strong evidence of non-neutral, accelerated evolution was found for DISC1, particularly for exon 2, the only coding region within the schizophrenia-associated haplotype. Additionally, genes associated with schizophrenia exhibited a statistically significant enrichment in their signals of positive selection in HapMap and PAML analyses of evolution along the human lineage, when compared with a control set of genes involved in neuronal activities. The selective forces underlying adaptive evolution of these genes remain largely unknown, but these findings provide convergent evidence consistent with the hypothesis that schizophrenia represents, in part, a maladaptive by-product of adaptive changes during human evolution.

Background The relatively fast selection of symbiotic bacteria within hosts and the potential transmission of these bacteria across generations of hosts raise the question of whether interactions between host and bacteria support emergent adaptive capabilities beyond those of germ-free hosts. Results To investigate possibilities for emergent adaptations that may distinguish composite host-microbiome systems from germ-free hosts, we introduce a population genetics model of a host-microbiome system with vertical transmission of bacteria. The host and its bacteria are jointly exposed to a toxic agent, creating a toxic stress that can be alleviated by selection of resistant individuals and by secretion of a detoxification agent (“detox”). We show that toxic exposure in one generation of hosts leads to selection of resistant bacteria, which in turn, increases the toxic tolerance of the host’s offspring. Prolonged exposure to toxin over many host generations promotes anadditional form of emergent adaptation due to selection of hosts based on detox produced by their bacterial community as a whole (as opposed to properties of individual bacteria). Conclusions These findings show that interactions between pure Darwinian selections of host and its bacteria can give rise to emergent adaptive capabilities, including Lamarckian-like adaptation of the host-microbiome system. Reviewers This article was reviewed by Eugene Koonin, Yuri Wolf and Philippe Huneman.

PII-disintegrins, cysteine-rich polypeptides broadly distributed in the venoms of geographically diverse species of vipers and rattlesnakes, antagonize the adhesive functions of [beta][sub]1 and [beta][sub]3 integrin receptors. PII-disintegrins evolved in Viperidae by neofunctionalization of disintegrin-like domains of duplicated PIII-snake venom hemorrhagic metalloproteinase (SVMP) genes recruited into the venom proteome before the radiation of the advanced snakes. Minimization of the gene (loss of introns and coding regions) and the protein structures (successive loss of disulfide bonds) underpins the postduplication divergence of disintegrins. However, little is known about the underlying genetic mechanisms that have generated the structural and functional diversity among disintegrins. Phylogenetic inference and maximum likelihood-based codon substitution approaches were used to analyze the evolution of the disintegrin family. The topology of the phylogenetic tree does not parallel that of the species tree. This incongruence is consistent with that expected for a multigene family undergoing a birth-and-death process in which the appearance and disappearance of loci are being driven by selection. Cysteine and buried residues appear to be under strong purifying selection due to their role in maintaining the active conformation of disintegrins. Divergence of disintegrins is strongly influenced by positive Darwinian selection causing accelerated rate of substitution in a substantial proportion of surface-exposed disintegrin residues. Global and lineage-specific sites evolving under diversifying selection were identified. Several sites are located within the integrin-binding loop and the C-terminal tail, two regions that form a conformational functional epitope. Arginine-glycine-aspartic acid (RGD) was inferred to represent the ancestral integrin-recognition motif, which emerged from the subgroup of PIII-SVMPs bearing the RDECD sequence. The most parsimonious nucleotide substitution model required for the emergence of all known disintegrin's integrin inhibitory motifs from an ancestral RGD sequence involves a minimum of three mutations. The adaptive advantage of the emergence of motifs targeting [beta][sub]1 integrins and the role of positively selected sites located within nonfunctional disintegrin regions appear to be difficult to rationalize in the context of a predator-prey arms race. Perhaps, this represents a consequence of the neofunctionalization potential of the disintegrin domain, a feature that may underlie its recruitment into the venom proteome followed by its successful transformation into a toxin. [PUBLICATION ABSTRACT]