Significant Signals Of Positive Selection Research Articles

The domain exchange between 3‐domain Cry toxins contributes to their diversity and specificity

Spore-forming bacteria Bacillus thuringiensis (Bt) can infect a wide range of invertebrate hosts efficiently and specifically; therefore, it is extensively implemented in ecological biopesticide production. This phenomenon is primarily explained by the presence of highly diversified crystalliferous Cry-toxins embedded in spores; however, the exact underlying causes of these toxins' divergence remains speculative. Nonetheless, understanding these natural mechanisms could provide insights into constructing artificial insecticidal moieties, which is highly demanded due to the sporadic occurrence of resistance in pests. Three-domain Cry toxins (3-D) represent the most abundant class of Cry toxins sharing similar structural features yet are immensely diverse in sequence identity. The latter leads to the absence of thorough phylogenetic studies involving these toxins on a large-scale, and we have attempted to fill this gap. We started with gathering Cry toxin sequences from all publicly available data sources, including Bt- genome assemblies, the Genbank database, and the Bt-nomenclature. After searching for 3-D toxins followed by deduplication, we obtained 370 toxins, for which phylogenetic trees were reconstructed based on individual domains and full-length toxins. A striking dissimilarity between their topologies (which was most evident for the third domain) has allowed us to infer that the domains' evolution is somehow independent and also supported the hypothesis that the emergence of the new toxins is orchestrated by shifts of the third domain between sequences. We next ran the RDP (Recombination Detection Protocol) program, searching for recombination events. As a result, we observed more than 70 events. In other words, almost every fourth toxin is formed after the onset of the recombination event. Two-third of the domain exchanges related to the third domain; nevertheless, exchanges affected all three domains. That observation substantially extends the former hypothesis. Finally, we assessed detected events for functional significance. Not only were the events characterized by a remarkable effect on the evolutionary selection (with significant signals of positive selection), but recombination led to alterations in host specificity for several events. All things considered, we have drawn a conclusion that the domain exchange is one of the primary mechanisms underlying the emergence and specificity of Cry toxins. This observation has a potential practical implementation in using in silico experiments to construct novel artificial biopesticides.

Mitogenomics of the endemic Ethiopian rats: looking for footprints of adaptive evolution in sky islands

Organisms living in high altitude must adapt to environmental conditions with hypoxia and low temperature, e.g. by changes in the structure and function of proteins associated with oxidative phosphorylation in mitochondria. Here we analysed the signs of adaptive evolution in 27 mitogenomes of endemic Ethiopian rats (Stenocephalemys), where individual species adapted to different elevation. Significant signals of positive selection were detected in 10 of the 13 mitochondrial protein-coding genes, with a majority of functional substitutions in the NADH dehydrogenase complex. Higher frequency of positively selected sites was found in phylogenetic lineages corresponding to Afroalpine specialists.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes.

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Comparative population genomic analysis uncovers novel genomic footprints and genes associated with small body size in Chinese pony

BackgroundBody size is considered as one of the most fundamental properties of an organism. Due to intensive breeding and artificial selection throughout the domestication history, horses exhibit striking variations for heights at withers and body sizes. Debao pony (DBP), a famous Chinese horse, is known for its small body size and lives in Guangxi mountains of southern China. In this study, we employed comparative population genomics to study the genetic basis underlying the small body size of DBP breed based on the whole genome sequencing data. To detect genomic signatures of positive selection, we applied three methods based on population comparison, fixation index (FST), cross population composite likelihood ratio (XP-CLR) and nucleotide diversity (θπ), and further analyzed the results to find genomic regions under selection for body size-related traits.ResultsA number of protein-coding genes in windows with the top 1% values of FST (367 genes), XP-CLR (681 genes), and log2 (θπ ratio) (332 genes) were identified. The most significant signal of positive selection was mapped to the NELL1 gene, probably underlies the body size and development traits, and may also have been selected for short stature in the DBP population. In addition, some other loci on different chromosomes were identified to be potentially involved in the development of body size.ConclusionsResults of our study identified some positively selected genes across the horse genome, which are possibly involved in body size traits. These novel candidate genes may be useful targets for clarifying our understanding of the molecular basis of body size and as such they should be of great interest for future research into the genetic architecture of relevant traits in horse breeding program.

Adaptive Evolution of Human-Isolated H5Nx Avian Influenza A Viruses.

Avian influenza A viruses (AIVs) H5N1, first identified in 1996, are highly pathogenic in domestic poultry and continue to occasionally infect humans. In this study, we sought to identify genetic changes that occurred during their multiple invasions to humans. We evaluated all available H5Nx AIV genomes. Significant signals of positive selection were detected in 29 host-shift branches. 126 parallel evolution sites were detected on these branches, including 17 well-known sites (such as T271A, A274T, T339M, Q591K, E627K, and D701N in PB2; A134V, D154N, S223N, and R497K in HA) that play roles in allowing AIVs to cross species barriers. Our study suggests that during human infections, H5Nx viruses have experienced adaptive evolution (positive selection and convergent evolution) that allowed them to adapt to their new host environments. Analyses of adaptive evolution should be useful in identifying candidate sites that play roles in human infections, which can be tested by functional experiments.

Genetic basis of amphibian larval development along a latitudinal gradient: Gene diversity, selection and links with phenotypic variation in transcription factor C/EBP-1.

Ectotherm development rates often show adaptive divergence along climatic gradients, but the genetic basis for this variation is rarely studied. Here, we investigated the genetic basis for phenotypic variation in larval development in the moor frog Rana arvalis from five regions along a latitudinal gradient from Germany to northern Sweden. We focused on the C/EBP-1 gene, a transcription factor associated with larval development time. Allele frequencies at C/EBP-1 varied strongly among geographical regions. Overall, the distribution of alleles along the gradient was in concordance with the dual post-glacial colonization routes into Scandinavia, with a large number of alleles exclusively present along the southern colonization route. Only three of 38 alleles were shared between the routes. Analysis of contemporary selection on C/EBP-1 showed divergent selection among the regions, probably reflecting adaptation to the local environmental conditions, although this was especially strong between southern and northern regions coinciding also with lineages from different colonization routes. Overall, the C/EBP-1 gene has historically been under purifying selection, but two specific amino acid positions showed significant signals of positive selection. These positions showed divergence between southern and northern regions, and we suggest that they are functionally involved in the climatic adaptation of larval development. Using phenotypic data from a common garden experiment, we found evidence for specific C/EBP-1 alleles being correlated with larval development time, suggesting a functional role in adaptation of larval development to large-scale climatic variation.

Convergent Evolution of Human-Isolated H7N9 Avian Influenza A Viruses

Avian influenza A virus H7N9 has caused 5 epidemic waves of human infections in China since 2013. Avian influenza A viruses may face strong selection to adapt to novel conditions when establishing themselves in humans. In this study, we sought to determine whether adaptive evolution had occurred in human-isolated H7N9 viruses. We evaluated all available genomes of H7N9 avian influenza A virus. Maximum likelihood trees were separately reconstructed for all 8 genes. Signals of positive selection and convergent evolution were then detected on branches that lead to changes in host tropism (from avian to human). We found that 3 genes had significant signals of positive selection (all of them P < .05). In addition, we detected 34 sites having significant signals for parallel evolution in 8 genes (all of them P < .05), including 7 well-known sites (Q591K, E627K, and D701N in PB2 gene; R156K, V202A, and L244Q in HA; and R289K in NA) that play roles in crossing species barriers for avian influenza A viruses. Our study suggests that, during infection in humans, H7N9 viruses have undergone adaptive evolution to adapt to their new host environment and that the sites where parallel evolution occurred might play roles in crossing species barriers and respond to the new selection pressures arising from their new host environments.

Comprehensive transcriptomic analysis of Tibetan Schizothoracinae fish Gymnocypris przewalskii reveals how it adapts to a high altitude aquatic life

BackgroundUnderstanding the genetic basis of adaptation to high altitude life is of paramount importance for preserving and managing genetic diversity in highland animals. This objective has been addressed mainly in terrestrial fauna but rarely in aquatic animals. Tibetan Schizothoracinae fish is the ideal model system in evolutionary biology, carrying key insights into evolutionary genetics of speciation and adaptation at high altitude. Gymnocypris przewalskii is the newly formed Schizothoracinae fish species in the Tibetan Plateau, inhabits chronic cold, extreme saline and alkaline aquatic environment in Lake Qinghai, thus evolving the unique genomic signatures to adapt extremely severe environments.ResultsTo characterize its genomic features, we assembled de novo transcriptome of G. przewalskii from Lake Qinghai. Intriguingly, by comparative genomic analyses of G. przewalskii and 8 other fish species, we identified potential expansions in gene families related to energy metabolism, transport and developmental functions, possibly underlying the adaptation to these environmental stresses. Through comprehensive molecular evolution analyses, we found that sets of genes controlling mitochondrion, ion homoeostasis, acid-base balance and innate immunity show significant signals of positive selection. Compared to previous studies on highland fishes, we failed to identify any positively selected genes related to hypoxia response.ConclusionsOur findings provide comprehensive insights into the genetic basis of teleost fish that underlie their adaptation to extreme high altitude aquatic life on the Tibetan Plateau.

Adaptive molecular evolution of the two-pore channel 1 gene TPC1 in the karst-adapted genus Primulina (Gesneriaceae).

Limestone karst areas possess high floral diversity and endemism. The genus Primulina, which contributes to the unique calcicole flora, has high species richness and exhibit specific soil-based habitat associations that are mainly distributed on calcareous karst soils. The adaptive molecular evolutionary mechanism of the genus to karst calcium-rich environments is still not well understood. The Ca2+-permeable channel TPC1 was used in this study to test whether its gene is involved in the local adaptation of Primulina to karst high-calcium soil environments. Specific amplification and sequencing primers were designed and used to amplify the full-length coding sequences of TPC1 from cDNA of 76 Primulina species. The sequence alignment without recombination and the corresponding reconstructed phylogeny tree were used in molecular evolutionary analyses at the nucleic acid level and amino acid level, respectively. Finally, the identified sites under positive selection were labelled on the predicted secondary structure of TPC1. Seventy-six full-length coding sequences of Primulina TPC1 were obtained. The length of the sequences varied between 2220 and 2286 bp and the insertion/deletion was located at the 5' end of the sequences. No signal of substitution saturation was detected in the sequences, while significant recombination breakpoints were detected. The molecular evolutionary analyses showed that TPC1 was dominated by purifying selection and the selective pressures were not significantly different among species lineages. However, significant signals of positive selection were detected at both TPC1 codon level and amino acid level, and five sites under positive selective pressure were identified by at least three different methods. The Ca2+-permeable channel TPC1 may be involved in the local adaptation of Primulina to karst Ca2+-rich environments. Different species lineages suffered similar selective pressure associated with calcium in karst environments, and episodic diversifying selection at a few sites may play a major role in the molecular evolution of Primulina TPC1.

Uncovering Adaptation from Sequence Data: Lessons from Genome Resequencing of Four Cattle Breeds.

Detecting the molecular basis of adaptation is one of the major questions in population genetics. With the advance in sequencing technologies, nearly complete interrogation of genome-wide polymorphisms in multiple populations is becoming feasible in some species, with the expectation that it will extend quickly to new ones. Here, we investigate the advantages of sequencing for the detection of adaptive loci in multiple populations, exploiting a recently published data set in cattle (Bos taurus). We used two different approaches to detect statistically significant signals of positive selection: a within-population approach aimed at identifying hard selective sweeps and a population-differentiation approach that can capture other selection events such as soft or incomplete sweeps. We show that the two methods are complementary in that they indeed capture different kinds of selection signatures. Our study confirmed some of the well-known adaptive loci in cattle (e.g., MC1R, KIT, GHR, PLAG1, NCAPG/LCORL) and detected some new ones (e.g., ARL15, PRLR, CYP19A1, PPM1L). Compared to genome scans based on medium- or high-density SNP data, we found that sequencing offered an increased detection power and a higher resolution in the localization of selection signatures. In several cases, we could even pinpoint the underlying causal adaptive mutation or at least a very small number of possible candidates (e.g., MC1R, PLAG1). Our results on these candidates suggest that a vast majority of adaptive mutations are likely to be regulatory rather than protein-coding variants.

Molecular characteristics and evolutionary analysis of a very virulent infectious bursal disease virus.

Infectious bursal disease virus (IBDV) poses a significant threat to the poultry industry. Viral protein 2 (VP2), the major structural protein of IBDV, has been subjected to frequent mutations that have imparted tremendous genetic diversity to the virus. To determine how amino acid mutations may affect the virulence of IBDV, we built a structural model of VP2 of a very virulent strain of IBDV identified in China, vvIBDV Gx, and performed a molecular dynamics simulation of the interaction between virulence sites. The study showed that the amino acid substitutions that distinguish vvIBDV from attenuated IBDV (H253Q and T284A) favor a hydrophobic and flexible conformation of β-barrel loops in VP2, which could promote interactions between the virus and potential IBDV-specific receptors. Population sequence analysis revealed that the IBDV strains prevalent in East Asia show a significant signal of positive selection at virulence sites 253 and 284. In addition, a signal of co-evolution between sites 253 and 284 was identified. These results suggest that changes in the virulence of IBDV may result from both the interaction and the co-evolution of multiple amino acid substitutions at virulence sites.

Genome-wide analysis of cold adaptation in indigenous Siberian populations.

Following the dispersal out of Africa, where hominins evolved in warm environments for millions of years, our species has colonised different climate zones of the world, including high latitudes and cold environments. The extent to which human habitation in (sub-)Arctic regions has been enabled by cultural buffering, short-term acclimatization and genetic adaptations is not clearly understood. Present day indigenous populations of Siberia show a number of phenotypic features, such as increased basal metabolic rate, low serum lipid levels and increased blood pressure that have been attributed to adaptation to the extreme cold climate. In this study we introduce a dataset of 200 individuals from ten indigenous Siberian populations that were genotyped for 730,525 SNPs across the genome to identify genes and non-coding regions that have undergone unusually rapid allele frequency and long-range haplotype homozygosity change in the recent past. At least three distinct population clusters could be identified among the Siberians, each of which showed a number of unique signals of selection. A region on chromosome 11 (chr11:66–69 Mb) contained the largest amount of clustering of significant signals and also the strongest signals in all the different selection tests performed. We present a list of candidate cold adaption genes that showed significant signals of positive selection with our strongest signals associated with genes involved in energy regulation and metabolism (CPT1A, LRP5, THADA) and vascular smooth muscle contraction (PRKG1). By employing a new method that paints phased chromosome chunks by their ancestry we distinguish local Siberian-specific long-range haplotype signals from those introduced by admixture.

Purifying selection after episodes of recurrent adaptive diversification in fungal pathogens

Tests of selective neutrality based on comparative and population genetic data target different time scales of selection. Investigating both time scales on the same genes has rarely been done and yet can provide insights into histories of selection. In pathogens, such a comparison can elucidate whether the same genes experience recurrent positive selection across specialization events onto novel hosts and selection in the short term, e.g. by coevolution with a current host. Here, we investigated whether the genes showing signs of positive selection (based on dN/dS ratios) between pathogens specialized on different hosts also exhibited footprints of recent positive selection. We sequenced 58 DNA fragments in the model fungal plant pathogen Microbotryum lychnidis-dioicae infecting Silene latifolia. Eleven focal genes were chosen because they showed a significant signal of positive selection in a previous study and putative functions likely involved in host–parasite interactions. In addition, 47 control genes were randomly chosen among genes showing no evidence of positive selection. We used a population genetics approach to search for signatures of recent selection, controlling for population structure and demographic history. We found footprints of purifying selection in all focal genes and of recent positive selection in two of them. Signs of purifying selection were also found at the same eleven genes in a closely related species, Microbotryum silenes-dioicae, infecting Silene dioica. These results suggest that genes experiencing episodes of adaptive diversification during host shifts may subsequently be under strong functional constraint, although some can remain under positive selection in the short term.

Parallel Evolution of Auditory Genes for Echolocation in Bats and Toothed Whales

The ability of bats and toothed whales to echolocate is a remarkable case of convergent evolution. Previous genetic studies have documented parallel evolution of nucleotide sequences in Prestin and KCNQ4, both of which are associated with voltage motility during the cochlear amplification of signals. Echolocation involves complex mechanisms. The most important factors include cochlear amplification, nerve transmission, and signal re-coding. Herein, we screen three genes that play different roles in this auditory system. Cadherin 23 (Cdh23) and its ligand, protocadherin 15 (Pcdh15), are essential for bundling motility in the sensory hair. Otoferlin (Otof) responds to nerve signal transmission in the auditory inner hair cell. Signals of parallel evolution occur in all three genes in the three groups of echolocators—two groups of bats (Yangochiroptera and Rhinolophoidea) plus the dolphin. Significant signals of positive selection also occur in Cdh23 in the Rhinolophoidea and dolphin, and Pcdh15 in Yangochiroptera. In addition, adult echolocating bats have higher levels of Otof expression in the auditory cortex than do their embryos and non-echolocation bats. Cdh23 and Pcdh15 encode the upper and lower parts of tip-links, and both genes show signals of convergent evolution and positive selection in echolocators, implying that they may co-evolve to optimize cochlear amplification. Convergent evolution and expression patterns of Otof suggest the potential role of nerve and brain in echolocation. Our synthesis of gene sequence and gene expression analyses reveals that positive selection, parallel evolution, and perhaps co-evolution and gene expression affect multiple hearing genes that play different roles in audition, including voltage and bundle motility in cochlear amplification, nerve transmission, and brain function.

MHC genotypes associate with resistance to a frog-killing fungus

The emerging amphibian disease chytridiomycosis is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian populations and species differ in susceptibility to Bd, yet we know surprisingly little about the genetic basis of this natural variation. MHC loci encode peptides that initiate acquired immunity in vertebrates, making them likely candidates for determining disease susceptibility. However, MHC genes have never been characterized in the context of chytridiomycosis. Here, we performed experimental Bd infections in laboratory-reared frogs collected from five populations that show natural variation in Bd susceptibility. We found that alleles of an expressed MHC class IIB locus associate with survival following Bd infection. Across populations, MHC heterozygosity was a significant predictor of survival. Within populations, MHC heterozygotes and individuals bearing MHC allele Q had a significantly reduced risk of death, and we detected a significant signal of positive selection along the evolutionary lineage leading to allele Q. Our findings demonstrate that immunogenetic variation affects chytridiomycosis survival under controlled experimental conditions, confirming that host genetic polymorphisms contribute to chytridiomycosis resistance.

Finding candidate genes under positive selection in Non-model species: examples of genes involved in host specialization in pathogens

Numerous genes in diverse organisms have been shown to be under positive selection, especially genes involved in reproduction, adaptation to contrasting environments, hybrid inviability, and host-pathogen interactions. Looking for genes under positive selection in pathogens has been a priority in efforts to investigate coevolution dynamics and to develop vaccines or drugs. To elucidate the functions involved in host specialization, here we aimed at identifying candidate sequences that could have evolved under positive selection among closely related pathogens specialized on different hosts. For this goal, we sequenced c. 17,000-32,000 ESTs from each of four Microbotryum species, which are fungal pathogens responsible for anther smut disease on host plants in the Caryophyllaceae. Forty-two of the 372 predicted orthologous genes showed significant signal of positive selection, which represents a good number of candidate genes for further investigation. Sequencing 16 of these genes in 9 additional Microbotryum species confirmed that they have indeed been rapidly evolving in the pathogen species specialized on different hosts. The genes showing significant signals of positive selection were putatively involved in nutrient uptake from the host, secondary metabolite synthesis and secretion, respiration under stressful conditions and stress response, hyphal growth and differentiation, and regulation of expression by other genes. Many of these genes had transmembrane domains and may therefore also be involved in pathogen recognition by the host. Our approach thus revealed fruitful and should be feasible for many non-model organisms for which candidate genes for diversifying selection are needed.

Positive Selection, Relaxation, and Acceleration in the Evolution of the Human and Chimp Genome

For years evolutionary biologists have been interested in searching for the genetic bases underlying humanness. Recent efforts at a large or a complete genomic scale have been conducted to search for positively selected genes in human and in chimp. However, recently developed methods allowing for a more sensitive and controlled approach in the detection of positive selection can be employed. Here, using 13,198 genes, we have deduced the sets of genes involved in rate acceleration, positive selection, and relaxation of selective constraints in human, in chimp, and in their ancestral lineage since the divergence from murids. Significant deviations from the strict molecular clock were observed in 469 human and in 651 chimp genes. The more stringent branch-site test of positive selection detected 108 human and 577 chimp positively selected genes. An important proportion of the positively selected genes did not show a significant acceleration in rates, and similarly, many of the accelerated genes did not show significant signals of positive selection. Functional differentiation of genes under rate acceleration, positive selection, and relaxation was not statistically significant between human and chimp with the exception of terms related to G-protein coupled receptors and sensory perception. Both of these were over-represented under relaxation in human in relation to chimp. Comparing differences between derived and ancestral lineages, a more conspicuous change in trends seems to have favored positive selection in the human lineage. Since most of the positively selected genes are different under the same functional categories between these species, we suggest that the individual roles of the alternative positively selected genes may be an important factor underlying biological differences between these species.