Putative Positive Selection Research Articles

Adaptive Evolution of Two Distinct Adaptive Haplotypes of Neanderthal Origin at the Immunoglobulin Heavy-chain Locus in East Asian and European Populations.

Immunoglobulins (Igs) have a crucial role in humoral immunity. Two recent studies have reported a high-frequency Neanderthal-introgressed haplotype throughout Eurasia and a high-frequency Neanderthal-introgressed haplotype specific to southern East Asia at the immunoglobulin heavy-chain (IGH) gene locus on chromosome 14q32.33. Surprisingly, we found the previously reported high-frequency Neanderthal-introgressed haplotype does not exist throughout Eurasia. Instead, our study identified two distinct high-frequency haplotypes of putative Neanderthal origin in East Asia and Europe, although they shared introgressed alleles. Notably, the alleles of putative Neanderthal origin reduced the expression of IGHG1 and increased the expression of IGHG2 and IGHG3 in various tissues. These putatively introgressed alleles also affected the production of IgG1 upon antigen stimulation and increased the risk of systemic lupus erythematosus. Additionally, the greatest genetic differentiation across the whole genome between southern and northern East Asians was observed for the East Asian haplotype of putative Neanderthal origin. The frequency decreased from southern to northern East Asia and correlated positively with the genome-wide proportion of southern East Asian ancestry, indicating that this putative positive selection likely occurred in the common ancestor of southern East Asian populations before the admixture with northern East Asian populations.

Novel Gene Rearrangement Pattern in Pachycrepoideus vindemmiae Mitochondrial Genome: New Gene Order in Pteromalidae (Hymenoptera: Chalcidoidea).

The mitochondrial genomes of Muscidifurax similadanacus, M. sinesensilla, Nasonia vitripennis, and Pachycrepoideus vindemmiae were sequenced to better understand the structural evolution of Pteromalidae mitogenomes. These newly sequenced mitogenomes all contained 37 genes. Nucleotide composition was AT-biased and the majority of the protein-coding genes exhibited a negative AT skew. All 13 protein-coding genes (PCGs) initiated with the standard start codon of ATN, excepted for nad1 of N. vitripennis, which started with TTG, and terminated with a typical stop codon TAA/TAG or an incomplete stop codon T. All transfer RNA (tRNA) genes were predicted to fold into the typical clover-leaf secondary structures, except for trnS1, which lacks the DHU arm in all species. In P. vindemmiae, trnR and trnQ lack the DHU arm and TΨC arm, respectively. Although most genes evolved under a strong purifying selection, the Ka/Ks value of the atp8 gene of P. vindemmiae was greater than 1, indicating putative positive selection. A novel transposition of trnR in P. vindemmiae was revealed, which was the first of this kind to be reported in Pteromalidae. Two kinds of datasets (PCG12 and AA) and two inference methods (maximum likelihood and Bayesian inference) were used to reconstruct a phylogenetic hypothesis for the newly sequenced mitogenomes of Pteromalidae and those deposited in GenBank. The topologies obtained recovered the monophyly of the three subfamilies included. Pachyneurinae and Pteromalinae were recovered as sister families, and both appeared sister to Sycophaginae. The pairwise breakpoint distances of mitogenome rearrangements were estimated to infer phylogeny among pteromalid species. The topology obtained was not totally congruent with those reconstructed using the ML and BI methods.

Archaic introgression contributed to the pre-agriculture adaptation of vitamin B1 metabolism in East Asia.

Thiamine (vitamin B1) is an essential micronutrient. Genes involved in thiamine metabolisms, such as SLC19A2, SLC35F3, and SLC35F4, were assumed to be underlying positive selection in East Asians, but the detailed mechanism remains unknown. Here, we analyzed genome data of 3,823 individuals representing 223 global populations and identified the adaptive haplotypes at thiamine genes. Interestingly, the putative adaptive haplotype at SLC35F4 was of Neanderthal ancestry, while that at SLC35F3 was also likely of archaic origins. Leveraging new methods and available ancient DNA data, we further demonstrated that the beneficial haplotypes reached a high frequency at least 10,000 years ago and are maintained persistently in present-day East Asians. We argue that pathogens, rather than agriculture developed ∼10,000 years ago in East Asia, were likely the initial driving force of the putative positive selection. Notably, the first American people did not carry the putative adaptive haplotype at SLC35F4.

Worldwide Population Structure of the Coffee Rust Fungus Hemileia vastatrix Is Strongly Shaped by Local Adaptation and Breeding History.

The devastating disease coffee leaf rust, caused by Hemileia vastatrix, has been a major constraint to worldwide coffee production. Recently, H.vastatrix populations were shown to be structured into three divergent genetic lineages with marked host specialization (C1, C2, and C3). However, there is yet no overall understanding of the population dynamics and adaptation of the most widespread and epidemiological relevant H.vastatrix group (C3). We used restriction site-associated DNA sequencing to generate 13,804 single nucleotide polymorphisms (SNPs) across a worldwide collection of 99 H. vastatrix isolates. Phylogenetic analyses uncovered a well-supported structuring within C3, with three main subgroups (SGs; SGI, SGII, and SGIII), which seem to reflect the historical distribution, breeding, and exchange of coffee cultivars. SGI shows a ladder-like diversification pattern and occurs across all four continents sampled, SGII is mainly restricted to Africa, and SGIII is observed only in Timor, revealing a higher genetic differentiation. Outlier and association tests globally identified 112 SNPs under putative positive selection, which impacted population structure. In particular, 29 overlapping SNPs per se seemed to have an extremely strong effect on H. vastatrix population divergence. We also found exclusive and fixed alleles associated with the SGs supporting local adaptation. Functional annotation revealed that transposable elements may play a role in host adaptation. Our study provides a higher-resolution perspective on the evolutionary history of H. vastatrix on cultivated coffee, showing its strong ability to adapt and the strength of the selective force imposed by coffee hosts, which should be taken into account when designing strategies for pathogen dissemination control and selective breeding.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

Evolutionary forces at work in partitiviruses.

The family Partitiviridae consists of dsRNA viruses with genome separated into two segments and encoding replicase and capsid protein only. We examined the nucleotide diversity expressed as the ratio dN/dS of nonsynonymous and synonymous substitutions, which has been calculated for 12 representative viruses of all five genera of partitiviruses. We can state that strong purifying selection works on both the RdRp and CP genes and propose that putative positive selection occurs also on the RdRp genes in two viruses. Among the 95 evaluated viruses, wherein both segments had been sequenced, 8 viruses in betapartitiviruses and 9 in alphapartitiviruses were identified as reassortment candidates because they differ extremely in their CP identity even as they are related in terms of RdRp. Furthermore, there are indications that reassortants are present among isolates of different viruses.

Persistence and continuous evolution of the human respiratory syncytial virus in northern Taiwan for two decades

The study aimed to characterize the molecular epidemiology, phylogenetic relationship, and population dynamics of the G protein gene in clinical respiratory syncytial virus (RSV) strains isolated from northern Taiwan. We analyzed a total of 160 and 116 G protein gene sequences of RSV-A and RSV-B representative strains, respectively, from 804 clinical viral stocks collected between July 2000 and June 2016. Population dynamic patterns of the RSV G protein gene were analyzed using Bayesian inference through the Markov chain Monte Carlo framework. A phylogenetic analysis revealed that RSV-A from Taiwan could be categorized into GA2, GA5, and GA7 lineages. GA2 of RSV-A could be further divided into NA1, NA2, NA4, and ON1 clades. These RSV-A lineages has been replaced over time, whereas RSV-B strains from Taiwan continually evolved from a single lineage with significant time-dependent waves. Four putative positive selection sites were observed in both RSV-A and RSV-B. The Bayesian skyline plot revealed that the local population dynamics of RSV were associated with lineage displacement events. Both circulating subtypes and population dynamics represented a unique local pattern. Our results affirm the necessity of continuing molecular surveillance of RSV to attain a more comprehensive understanding of epidemics.

Comparative transcriptomics in three Passerida species provides insights into the evolution of avian mitochondrial complex I

Recent studies have shown that mitochondria play a crucial role in cellular energy production through the oxidative phosphorylation (OXPHOS) system. Complex I (NADH:ubiquinone oxidoreductase), the first and largest enzyme complex of the OXPHOS system, includes both nuclear- and mitochondrial-encoded proteins. However, the patterns of natural selection and phylogenetic implications of complex I in birds still remain unclear. In this study, we combined transcriptomic and phylogenetic analyses to comprehensively determine the evolution of avian complex I. The transcriptomes of three Passerida species (Leiothrix lutea, Spodiopsar sericeus, and Passer montanus) were obtained using the Illumina HiSeq™ 2500 system. More than 192,000,000 clean reads were assembled in a total of 828,267 transcripts. Evolutionary selection analysis suggested that six genes of the core subunits in avian complex I may have undergone putative positive selection. Notably, we found that the mean dN/dS (ω) ratio for mitochondrial genes of core subunits was significantly lower than that for nuclear genes of non-core subunits within complex I. The constructed maximum-parsimony, maximum-likelihood, and Bayesian inference phylogenetic trees were based on 44 complex I genes. We verified that the family Paridae (represented by Parus major and Pseudopodoces humilis) was clustered with Musicicapoidea. Our results provide new insights into the evolution of avian mitochondrial complex I.

Temporal genetic patterns of diversity and structure evidence chaotic genetic patchiness in a spiny lobster.

Population structure of many marine organisms is spatially patchy and varies within and between years, a phenomenon defined as chaotic genetic patchiness. This results from the combination of planktonic larval dispersal and environmental stochasticity. Additionally, in species with bi-partite life, postsettlement selection can magnify these genetic differences. The high fecundity (up to 500,000 eggs annually) and protracted larval duration (12-24months) and dispersal of the southern rock lobster, Jasus edwardsii, make it a good test species for chaotic genetic patchiness and selection during early benthic life. Here, we used double digest restriction site-associated DNA sequencing (ddRADseq) to investigate chaotic genetic patchiness and postsettlement selection in this species. We assessed differences in genetic structure and diversity of recently settled pueruli across four settlement years and between two sites in southeast Australia separated by approximately 1,000km. Postsettlement selection was investigated by identifying loci under putative positive selection between recently settled pueruli and postpueruli and quantifying differences in the magnitude and strength of the selection at each year and site. Genetic differences within and among sites through time in neutral SNP markers indicated chaotic genetic patchiness. Recently settled puerulus at the southernmost site exhibited lower genetic diversity during years of low puerulus catches, further supporting this hypothesis. Finally, analyses of outlier SNPs detected fluctuations in the magnitude and strength of the markers putatively under positive selection over space and time. One locus under putative positive selection was consistent at both locations during the same years, suggesting the existence of weak postsettlement selection.

Persistence of Lineage IV Peste-des-petits ruminants virus within Israel since 1993: An evolutionary perspective.

Peste-des-petits ruminants (PPR) is one of the most important infectious diseases of domesticated small ruminants. From the initial identification in 1942 in West Africa, PPR virus (PPRV) has spread throughout much of the developing world. PPRV is now considered endemic throughout Africa, with the notable exception of South Africa, the Middle-East and Israel, as well as South-, East-, and Central Asia. Despite this widespread dispersal, the evolution and transmission of PPRV in endemic populations is not well understood. This understanding will be critical in the planning of rational measures to eradicate PPRV by the planned time as defined by the FAO and OIE. To further advance the understanding of the evolution of PPRV the full genome sequence of 18 viruses isolated from Israel from consecutive years between 1997–2014 were generated. This data set is unique and crucial for the understanding of the evolution of PPRV, as it represents the first set of full-length sequence data available from consecutive years from a single geographic location. Analysis of these full genome sequences shows 96.2–99.9% nucleotide conservation across the Israel isolates and further demonstrates the strong purifying selection pressures on PPRV within Israel and globally. Four amino acid substitutions indicative of putative positive selection were additionally identified within the Israel isolates. The mean substitution rate per site per year was estimated to be 9.22 x 10−4 (95% HPD 6.206 x 10−4–1.26 x 10−3). Using Bayesian and phylogenetic analyses we further demonstrate that the PPRV isolates from Israel belongs to linage IV and form a single strong regional cluster within all other lineage IV viruses circulating worldwide implying a single incursion into Israel.

Flux Control in Glycolysis Varies Across the Tree of Life.

Biochemical thought posits that rate-limiting steps (defined here as points of flux control) are strongly selected as points of pathway regulation and control and are thus expected to be evolutionarily conserved. Conversely, population genetic thought based upon the concepts of mutation-selection-drift balance at the pathway level might suggest variation in flux controlling steps over evolutionary time. Glycolysis, as one of the most conserved and best characterized pathways, was studied to evaluate its evolutionary conservation. The flux controlling step in glycolysis was found to vary over the tree of life. Further, phylogenetic analysis suggested at least 60 events of gene duplication and additional events of putative positive selection that might alter pathway kinetic properties. Together, these results suggest that even with presumed largely negative selection on pathway output on glycolysis, the co-evolutionary process under the hood is dynamic.

Transcriptome Sequencing of Lima Bean (Phaseolus lunatus) to Identify Putative Positive Selection in Phaseolus and Legumes.

The identification of genes under positive selection is a central goal of evolutionary biology. Many legume species, including Phaseolus vulgaris (common bean) and Phaseolus lunatus (lima bean), have important ecological and economic value. In this study, we sequenced and assembled the transcriptome of one Phaseolus species, lima bean. A comparison with the genomes of six other legume species, including the common bean, Medicago, lotus, soybean, chickpea, and pigeonpea, revealed 15 and 4 orthologous groups with signatures of positive selection among the two Phaseolus species and among the seven legume species, respectively. Characterization of these positively selected genes using Non redundant (nr) annotation, gene ontology (GO) classification, GO term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that these genes are mostly involved in thylakoids, photosynthesis and metabolism. This study identified genes that may be related to the divergence of the Phaseolus and legume species. These detected genes are particularly good candidates for subsequent functional studies.

Teratozoospermia: spotlight on the main genetic actors in the human

Male infertility affects >20 million men worldwide and represents a major health concern. Although multifactorial, male infertility has a strong genetic basis which has so far not been extensively studied. Recent studies of consanguineous families and of small cohorts of phenotypically homogeneous patients have however allowed the identification of a number of autosomal recessive causes of teratozoospermia. Homozygous mutations of aurora kinase C (AURKC) were first described to be responsible for most cases of macrozoospermia. Other genes defects have later been identified in spermatogenesis associated 16 (SPATA16) and dpy-19-like 2 (DPY19L2) in patients with globozoospermia and more recently in dynein, axonemal, heavy chain 1 (DNAH1) in a heterogeneous group of patients presenting with flagellar abnormalities previously described as dysplasia of the fibrous sheath or short/stump tail syndromes, which we propose to call multiple morphological abnormalities of the flagella (MMAF). A comprehensive review of the scientific literature available in PubMed/Medline was conducted for studies on human genetics, experimental models and physiopathology related to teratozoospermia in particular globozoospermia, large headed spermatozoa and flagellar abnormalities. The search included all articles with an English abstract available online before September 2014. Molecular studies of numerous unrelated patients with globozoospermia and large-headed spermatozoa confirmed that mutations in DPY19L2 and AURKC are mainly responsible for their respective pathological phenotype. In globozoospermia, the deletion of the totality of the DPY19L2 gene represents ∼ 81% of the pathological alleles but point mutations affecting the protein function have also been described. In macrozoospermia only two recurrent mutations were identified in AURKC, accounting for almost all the pathological alleles, raising the possibility of a putative positive selection of heterozygous individuals. The recent identification of DNAH1 mutations in a proportion of patients with MMAF is promising but emphasizes that this phenotype is genetically heterogeneous. Moreover, the identification of mutations in a dynein strengthens the emerging point of view that MMAF may be a phenotypic variation of the classical forms of primary ciliary dyskinesia. Based on data from human and animal models, the MMAF phenotype seems to be favored by defects directly or indirectly affecting the central pair of axonemal microtubules of the sperm flagella. The studies described here provide valuable information regarding the genetic and molecular defects causing infertility, to improve our understanding of the physiopathology of teratozoospermia while giving a detailed characterization of specific features of spermatogenesis. Furthermore, these findings have a significant influence on the diagnostic strategy for teratozoospermic patients allowing the clinician to provide the patient with informed genetic counseling, to adopt the best course of treatment and to develop personalized medicine directly targeting the defective gene products.

Adaptive Evolution as a Predictor of Species-Specific Innate Immune Response.

It has been proposed that positive selection may be associated with protein functional change. For example, human and macaque have different outcomes to HIV infection and it has been shown that residues under positive selection in the macaque TRIM5α receptor locate to the region known to influence species-specific response to HIV. In general, however, the relationship between sequence and function has proven difficult to fully elucidate, and it is the role of large-scale studies to help bridge this gap in our understanding by revealing major patterns in the data that correlate genotype with function or phenotype. In this study, we investigate the level of species-specific positive selection in innate immune genes from human and mouse. In total, we analyzed 456 innate immune genes using codon-based models of evolution, comparing human, mouse, and 19 other vertebrate species to identify putative species-specific positive selection. Then we used population genomic data from the recently completed Neanderthal genome project, the 1000 human genomes project, and the 17 laboratory mouse genomes project to determine whether the residues that were putatively positively selected are fixed or variable in these populations. We find evidence of species-specific positive selection on both the human and the mouse branches and we show that the classes of genes under positive selection cluster by function and by interaction. Data from this study provide us with targets to test the relationship between positive selection and protein function and ultimately to test the relationship between positive selection and discordant phenotypes.

Functionality of Intergenic Transcription: An Evolutionary Comparison

Although a large proportion of human transcription occurs outside the boundaries of known genes, the functional significance of this transcription remains unknown. We have compared the expression patterns of known genes as well as intergenic transcripts within the ENCODE regions between humans and chimpanzees in brain, heart, testis, and lymphoblastoid cell lines. We find that intergenic transcripts show patterns of tissue-specific conservation of their expression, which are comparable to exonic transcripts of known genes. This suggests that intergenic transcripts are subject to functional constraints that restrict their rate of evolutionary change as well as putative positive selection to an extent comparable to that of classical protein-coding genes. In brain and testis, we find that part of this intergenic transcription is caused by widespread use of alternative promoters. Further, we find that about half of the expression differences between humans and chimpanzees are due to intergenic transcripts.

Reevaluation of amino acid variability of the human immunodeficiency virus type 1 gp120 envelope glycoprotein and prediction of new discontinuous epitopes.

To elucidate the evolutionary mechanisms of the human immunodeficiency virus type 1 gp120 envelope glycoprotein at the single-site level, the degree of amino acid variation and the numbers of synonymous and nonsynonymous substitutions were examined in 186 nucleotide sequences for gp120 (subtype B). Analyses of amino acid variabilities showed that the level of variability was very different from site to site in both conserved (C1 to C5) and variable (V1 to V5) regions previously assigned. To examine the relative importance of positive and negative selection for each amino acid position, the numbers of synonymous and nonsynonymous substitutions that occurred at each codon position were estimated by taking phylogenetic relationships into account. Among the 414 codon positions examined, we identified 33 positions where nonsynonymous substitutions were significantly predominant. These positions where positive selection may be operating, which we call putative positive selection (PS) sites, were found not only in the variable loops but also in the conserved regions (C1 to C4). In particular, we found seven PS sites at the surface positions of the alpha-helix (positions 335 to 347 in the C3 region) in the opposite face for CD4 binding. Furthermore, two PS sites in the C2 region and four PS sites in the C4 region were detected in the same face of the protein. The PS sites found in the C2, C3, and C4 regions were separated in the amino acid sequence but close together in the three-dimensional structure. This observation suggests the existence of discontinuous epitopes in the protein's surface including this alpha-helix, although the antigenicity of this area has not been reported yet.

A B-cell receptor-specific selection step governs immature to mature B cell differentiation.

Seventy percent of peripheral immature conventional (B2) B cells fail to develop into mature B cells. The nature of this cell loss has not been characterized; the process that governs which immature B cells develop into long-lived peripheral B cells could be either stochastic or selective. Here, we demonstrate that this step is in fact selective, in that the fate of an immature B cell is highly dependent on its Ig receptor specificity. A significant skewing of the B cell receptor repertoire occurs by the time cells enter the mature B cell fraction, which indicates that there is selection of only a minority of immature B cells to become mature B cells. Because only a few heavy-light chain pairs are enhanced of the diverse available repertoire, we favor the idea that selection is positive for these few heavy-light chain pairs rather than negative against nearly all others. Because most immature B cells are lost at this transition, this putative positive selection event is likely to be a major force shaping the mature B cell receptor repertoire available for adaptive immune responses.