Positive Diversifying Selection Research Articles

Data mining reveals tissue-specific expression and host lineage-associated forms of Apis mellifera filamentous virus.

Apis mellifera filamentous virus (AmFV) is a large double-stranded DNA virus of uncertain phylogenetic position that infects honey bees (Apis mellifera). Little is known about AmFV evolution or molecular aspects of infection. Accurate annotation of open-reading frames (ORFs) is challenged by weak homology to other known viruses. This study was undertaken to evaluate ORFs (including coding-frame conservation, codon bias, and purifying selection), quantify genetic variation within AmFV, identify host characteristics that covary with infection rate, and examine viral expression patterns in different tissues. Short-read data were accessed from the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI). Sequence reads were downloaded from accessions meeting search criteria and scanned for kmers representative of AmFV genomic sequence. Samples with kmer counts above specified thresholds were downloaded in full for mapping to reference sequences and de novo assembly. At least three distinct evolutionary lineages of AmFV exist. Clade 1 predominates in Europe but in the Americas and Africa it is replaced by the other clades as infection level increases in hosts. Only clade 3 was found at high relative abundance in hosts with African ancestry, whereas all clades achieved high relative abundance in bees of non-African ancestry. In Europe and Africa, clade 2 was generally detected only in low-level infections but was locally dominant in some North American samples. The geographic distribution of clade 3 was consistent with an introduction to the Americas with 'Africanized' honey bees in the 1950s. Localized genomic regions of very high nucleotide divergence in individual isolates suggest recombination with additional, as-yet unidentified AmFV lineages. A set of 155 high-confidence ORFs was annotated based on evolutionary conservation in six AmFV genome sequences representative of the three clades. Pairwise protein-level identity averaged 94.6% across ORFs (range 77.1-100%), which generally exhibited low evolutionary rates and moderate to strong codon bias. However, no robust example of positive diversifying selection on coding sequence was found in these alignments. Most of the genome was detected in RNA short-read alignments. Transcriptome assembly often yielded contigs in excess of 50 kb and containing ORFs in both orientations, and the termini of long transcripts were associated with tandem repeats. Lower levels of AmFV RNA were detected in brain tissue compared to abdominal tissue, and a distinct set of ORFs had minimal to no detectable expression in brain tissue. A scan of DNA accessions from the parasitic mite Varroa destructor was inconclusive with respect to replication in that species. Collectively, these results expand our understanding of this enigmatic virus, revealing transcriptional complexity and co-evolutionary associations with host lineage.

Genetic diversity and natural selection of apical membrane antigen-1 (ama-1) in Cameroonian Plasmodium falciparum isolates

Antigenic variation associated with genetic diversity in global Plasmodium falciparum apical membrane antigen-1 (PfAMA-1) is a major impediment to designing an effective malaria vaccine. Here, we report the first study on genetic diversity and natural selection of the Pfama-1 gene in P. falciparum isolates from Cameroon. A total of 328 P. falciparum positive samples collected during 2016 and 2019 from five localities of Cameroon were analysed. The ectodomain coding fragment of Pfama-1 gene was amplified for polymorphism profiling and natural selection analysis. A total of 108 distinct haplotypes were found in 203 P. falciparum isolates with considerable nucleotide diversity (π = 0.016) and haplotype diversity (Hd = 0.976). Most amino acid substitutions detected were scattered in ectodomain-I and few specific mutations viz P145L, K148Q, K462I, L463F, N471K, S482L, E537G, K546R and I547F were seen only in Cameroonian isolates. A tendency of natural selection towards positive diversifying selection was observed (Taj-D = 2.058). Five positively selected codon sites (P145L, S283L, Q308E/K, P330S and I547F) were identified, which overlapped with predicted B-cell epitopes and red blood cell (RBC) binding sites, suggesting their potential implication in host immune pressure and parasite-RBC binding complex modulation. The Cameroonian P. falciparum populations indicated a moderate level of genetic differentiation when compared with global sequences, with few exceptions from Vietnam and Venezuela. Our findings provide baseline data on existing Pfama-1 gene polymorphisms in Cameroonian field isolates, which will be useful information for malaria vaccine design.

Serological relationships among strains of grapevine leafroll-associated virus 4 reflect the evolutive behavior of its coat protein gene

This research studied serological relationships and genetic diversity of Argentinean isolates of grapevine leafroll-associated virus-4 (GLRaV-4). Phylogenetic analysis of coat protein (CP) sequences from 19 local isolates revealed clustering with the previously described GLRaV-4 strain 5, strain 6, and strain 9 groups. Evolutionary sequence analysis of the obtained and database-available sequences showed evidence of recombination events. Additionally, both CP N- and C-terminal regions appeared to be under purifying selection, but the N-terminal region presented seven sites under positive selection, with a dN/dS ratio 5-fold greater than that of the C-terminal region. Serological reactivity against monoclonal antibodies supports a higher occurrence probability for linear epitopes in the N-terminal region, as inferred by the sequence analysis. The obtained results reflect an unusual evolutionary behavior of the CP that, together with protein serological reactivity, suggests biological significance of the observed variability. Highlights: The occurrence of GLRaV-4 strains -5, -6 and -9 was recorded in Argentina. The coat protein sequence showed a highly variable N-terminal region among the strains. An evolutionary analysis inferred the occurrence of sites under positive diversifying selection pressure. The most variable region of the coat protein presented the highest probability of occurrence of linear epitopes.

The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae)

Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts’ nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.

Antiviral function and viral antagonism of the rapidly evolving dynein activating adaptor NINL.

Viruses interact with the intracellular transport machinery to promote viral replication. Such host-virus interactions can drive host gene adaptation, leaving signatures of pathogen-driven evolution in host genomes. Here, we leverage these genetic signatures to identify the dynein activating adaptor, ninein-like (NINL), as a critical component in the antiviral innate immune response and as a target of viral antagonism. Unique among genes encoding components of active dynein complexes, NINL has evolved under recurrent positive (diversifying) selection, particularly in its carboxy-terminal cargo-binding region. Consistent with a role for NINL in host immunity, we demonstrate that NINL knockout cells exhibit an impaired response to interferon, resulting in increased permissiveness to viral replication. Moreover, we show that proteases encoded by diverse picornaviruses and coronaviruses cleave and disrupt NINL function in a host- and virus-specific manner. Our work reveals the importance of NINL in the antiviral response and the utility of using signatures of host-virus genetic conflicts to uncover new components of antiviral immunity and targets of viral antagonism.

Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas.

BackgroundThe fathead minnow (Pimephales promelas) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation.MethodsHere we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950.ResultsWe recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow (P. notatus) as sister to fathead minnow, with the slim minnow (P. tenellus) and bullhead minnow (P. vigilax) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (FST) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.

Evolution of the recombination regulator PRDM9 in minke whales

BackgroundPRDM9 is a key regulator of meiotic recombination in most metazoans, responsible for reshuffling parental genomes. During meiosis, the PRDM9 protein recognizes and binds specific target motifs via its array of C2H2 zinc-fingers encoded by a rapidly evolving minisatellite. The gene coding for PRDM9 is the only speciation gene identified in vertebrates to date and shows high variation, particularly in the DNA-recognizing positions of the zinc-finger array, within and between species. Across all vertebrate genomes studied for PRDM9 evolution, only one genome lacks variability between repeat types – that of the North Pacific minke whale. This study aims to understand the evolution and diversity of Prdm9 in minke whales, which display the most unusual genome reference allele of Prdm9 so far discovered in mammals.ResultsMinke whales possess all the features characteristic of PRDM9-directed recombination, including complete KRAB, SSXRD and SET domains and a rapidly evolving array of C2H2-type-Zincfingers (ZnF) with evidence of rapid evolution, particularly at DNA-recognizing positions that evolve under positive diversifying selection. Seventeen novel PRDM9 variants were identified within the Antarctic minke whale species, plus a single distinct PRDM9 variant in Common minke whales – shared across North Atlantic and North Pacific minke whale subspecies boundaries.ConclusionThe PRDM9 ZnF array evolves rapidly, in minke whales, with at least one DNA-recognizing position under positive selection. Extensive PRDM9 diversity is observed, particularly in the Antarctic in minke whales. Common minke whales shared a specific Prdm9 allele across subspecies boundaries, suggesting incomplete speciation by the mechanisms associated with PRDM9 hybrid sterility.

Genetic Diversity and Selection of Plasmodium vivax Apical Membrane Antigen-1 in China-Myanmar Border of Yunnan Province, China, 2009-2016.

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is an important vaccine candidate for vivax malaria. However, antigenic variation within PvAMA-1 is a major obstacle to the design of a global protective malaria vaccine. In this study, we analyzed the genetic polymorphism and selection of the PvAMA-1 gene from 152 P. vivax isolates from imported cases to China, collected in the China–Myanmar border (CMB) area in Yunnan Province (YP) during 2009–2011 (n = 71) and 2014–2016 (n = 81), in comparison with PvAMA-1 gene information from Myanmar (n = 73), collected from public data. The overall nucleotide diversity of the PvAMA-1 gene from the 152 YP isolates was 0.007 with 76 haplotypes identified (Hd = 0.958). Results from the population structure suggested three groups among the YP and Myanmar isolates with optimized clusters value of K = 7. In addition, YP (2014–2016) isolates generally lacked some K components that were commonly found in YP (2009–2011) and Myanmar. Meanwhile, PvAMA-1 domain I is found to be the dominant target of positive diversifying selection and most mutation loci were found in this domain. The mutation frequencies of D107N/A, R112K/T, K120R, E145A, E277K, and R438H in PvAMA-1 were more than 70% in the YP isolates. In conclusion, high genetic diversity and positive selection were found in the PvAMA-1 gene from YP isolates, which are significant findings for the design and development of PvAMA-1-based malaria vaccine.

Functional genomics of abiotic environmental adaptation in lacertid lizards and other vertebrates.

Understanding the genomic basis of adaptation to different abiotic environments is important in the context of climate change and resulting short-term environmental fluctuations. Using functional and comparative genomics approaches, we here investigated whether signatures of genomic adaptation to a set of environmental parameters are concentrated in specific subsets of genes and functions in lacertid lizards and other vertebrates. We first identify 200 genes with signatures of positive diversifying selection from transcriptomes of 24 species of lacertid lizards and demonstrate their involvement in physiological and morphological adaptations to climate. To understand how functionally similar these genes are to previously predicted candidate functions for climate adaptation and to compare them with other vertebrate species, we then performed a meta-analysis of 1,100 genes under selection obtained from -omics studies in vertebrate species adapted to different abiotic factors. We found that the vertebrate gene set formed a tightly connected interactome, which was to 23% enriched in previously predicted functions of adaptation to climate, and to a large part (18%) involved in organismal stress response. We found a much higher degree of identical genes being repeatedly selected among different animal groups (43.6%), and of functional similarity and post-translational modifications than expected by chance, and no clear functional division between genes used for ectotherm and endotherm physiological strategies. In total, 171 out of 200 genes of Lacertidae were part of this network. These results highlight an important role of a comparatively small set of genes and their functions in environmental adaptation and narrow the set of candidate pathways and markers to be used in future research on adaptation and stress response related to climate change.

Genetic Diversity Analysis of Surface-Related Antigen (SRA) in Plasmodium falciparum Imported From Africa to China

Plasmodium falciparum surface-related antigen (SRA) is located on the surfaces of gametocyte and merozoite and has the structural and functional characteristics of potential targets for multistage vaccine development. However, little information is available regarding the genetic polymorphism of pfsra. To determine the extent of genetic variation about P. falciparum by characterizing the sra sequence, 74 P. falciparum samples were collected from migrant workers who returned to China from 12 countries of Africa between 2015 and 2019. The full length of the sra gene was amplified and sequenced. The average pairwise nucleotide diversities (π) of P. falciparum sra gene was 0.00132, and the haplotype diversity (Hd) was 0.770. The average number of nucleotide differences (k) for pfsra was 3.049. The ratio of non-synonymous (dN) to synonymous (dS) substitutions across sites (dN/dS) was 1.365. Amino acid substitutions of P. falciparum SRA could be categorized into 35 unique amino acid variants. Neutrality tests showed that the polymorphism of PfSRA was maintained by positive diversifying selection, which indicated its role as a potential target of protective immune responses and a vaccine candidate. Overall, the ability of the N-terminal of PfSRA antibodies to evoke inhibition of merozoite invasion of erythrocytes and conserved amino acid at low genetic diversity suggest that the N-terminal of PfSRA could be evaluated as a vaccine candidate against P. falciparum infection.

Weak selection on synonymous codons substantially inflates dN/dS estimates in bacteria

Synonymous codon substitutions are not always selectively neutral as revealed by several types of analyses, including studies of codon usage patterns among genes. We analyzed codon usage in 13 bacterial genomes sampled from across a large order of bacteria, Enterobacterales, and identified presumptively neutral and selected classes of synonymous substitutions. To estimate substitution rates, given a neutral/selected classification of synonymous substitutions, we developed a flexible [Formula: see text] substitution model that allows multiple classes of synonymous substitutions. Under this multiclass synonymous substitution (MSS) model, the denominator of [Formula: see text] includes only the strictly neutral class of synonymous substitutions. On average, the value of [Formula: see text] under the MSS model was 80% of that under the standard codon model in which all synonymous substitutions are assumed to be neutral. The indication is that conventional [Formula: see text] analyses overestimate these values and thus overestimate the frequency of positive diversifying selection and underestimate the strength of purifying selection. To quantify the strength of selection necessary to explain this reduction, we developed a model of selected compensatory codon substitutions. The reduction in synonymous substitution rate, and thus the contribution that selection makes to codon bias variation among genes, can be adequately explained by very weak selection, with a mean product of population size and selection coefficient, [Formula: see text].

Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen.

Virus host shifts are generally associated with novel adaptations to exploit the cells of the new host species optimally. Surprisingly, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has apparently required little to no significant adaptation to humans since the start of the Coronavirus Disease 2019 (COVID-19) pandemic and to October 2020. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus the early SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects evidence for significant positive episodic diversifying selection acting at the base of the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in these ancestral bat hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor about 1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. While an undiscovered "facilitating" intermediate species cannot be discounted, collectively, our results support the progenitor of SARS-CoV-2 being capable of efficient human-human transmission as a consequence of its adaptive evolutionary history in bats, not humans, which created a relatively generalist virus.

Evolution of Fc Receptor-Like Scavenger in Mammals.

Fc receptor-like (FCRL) molecules comprise a large family of receptors, homologous to the receptors for the Fc portion of immunoglobulins (FCR). Within this family, an unusual gene known to exist in mice, rats and dogs, termed FCRLS, encodes a chimeric protein with both Ig-like FCRL and type B scavenger-receptor cysteine-rich (SRCR)-like domains. In mice, FCRLS is located next to the CD5L and KIRREL1 genes. Here, we show that the curious FCRLS gene is actually present across major mammalian groups, but its annotation is generally incorrect or absent. Anchored on mouse FCRLS and FCRL2 genomic sequence alignments, phylogenetic analyses demonstrated that many mammalian sequences currently annotated as FCRL2 cluster with FCRLS, supported by a conserved genetic synteny among organisms. This analysis shows that FCRLS is present in Rodentia, some Carnivora (Canidae and Ursidae), Chiroptera, Arctiodactyla, Proboscidae, and some Primata. Thus, the FCRLS most likely originated in a eutherian mammal ancestor since it is not present in Monotremata or Marsupialia. FCRLS has a peculiar distribution pattern across mammalian lineages, being present in some species, but absent in others from the same family, as in carnivores for example. The most parsimonious hypothesis to explain this FCRLS evolution is that it was convergently lost in several independent mammalian lineages. Analyses of branch-specific nucleotide evolutionary rates, show that FCRL2 and FCRLS have similar ranges of rates across mammals, suggesting that both genes have crucial, but separate functions in the immune system. Bayesian estimates of evolutionary rates for FCRLS in mammalian lineages revealed that carnivores display the highest mutation rate after rodents. Additionally, positive diversifying selection was detected for both FCRL2 and FCRLS. Our results show that the presence of the FCRLS gene is older and more widespread across mammals than previously thought and appears to be functional, being under positive selection. Its precise physiologic role should thus be investigated.

Expansion and Accelerated Evolution of 9-Exon Odorant Receptors in Polistes Paper Wasps.

Independent origins of sociality in bees and ants are associated with independent expansions of particular odorant receptor (OR) gene subfamilies. In ants, one clade within the OR gene family, the 9-exon subfamily, has dramatically expanded. These receptors detect cuticular hydrocarbons (CHCs), key social signaling molecules in insects. It is unclear to what extent 9-exon OR subfamily expansion is associated with the independent evolution of sociality across Hymenoptera, warranting studies of taxa with independently derived social behavior. Here, we describe OR gene family evolution in the northern paper wasp, Polistes fuscatus, and compare it to four additional paper wasp species spanning ∼40 million years of evolutionary divergence. We find 200 putatively functional OR genes in P. fuscatus, matching predictions from neuroanatomy, and more than half of these are in the 9-exon subfamily. Most OR gene expansions are tandemly arrayed at orthologous loci in Polistes genomes, and microsynteny analysis shows species-specific gain and loss of 9-exon ORs within tandem arrays. There is evidence of episodic positive diversifying selection shaping ORs in expanded subfamilies. Values of omega (dN/dS) are higher among 9-exon ORs compared to other OR subfamilies. Within the Polistes OR gene tree, branches in the 9-exon OR clade experience relaxed negative (relaxed purifying) selection relative to other branches in the tree. Patterns of OR evolution within Polistes are consistent with 9-exon OR function in CHC perception by combinatorial coding, with both natural selection and neutral drift contributing to interspecies differences in gene copy number and sequence.

Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

Identifying Candidate Genetic Markers of CDV Cross-Species Pathogenicity in African Lions.

Canine distemper virus (CDV) is a multi-host pathogen with variable clinical outcomes of infection across and within species. We used whole-genome sequencing (WGS) to search for viral markers correlated with clinical distemper in African lions. To identify candidate markers, we first documented single-nucleotide polymorphisms (SNPs) differentiating CDV strains associated with different clinical outcomes in lions in East Africa. We then conducted evolutionary analyses on WGS from all global CDV lineages to identify loci subject to selection. SNPs that both differentiated East African strains and were under selection were mapped to a phylogenetic tree representing global CDV diversity to assess if candidate markers correlated with documented outbreaks of clinical distemper in lions (n = 3). Of 54 SNPs differentiating East African strains, ten were under positive or episodic diversifying selection and 20 occurred in the clinical strain despite strong purifying selection at those loci. Candidate markers were in functional domains of the RNP complex (n = 19), the matrix protein (n = 4), on CDV glycoproteins (n = 5), and on the V protein (n = 1). We found mutations at two loci in common between sequences from three CDV outbreaks of clinical distemper in African lions; one in the signaling lymphocytic activation molecule receptor (SLAM)-binding region of the hemagglutinin protein and another in the catalytic center of phosphodiester bond formation on the large polymerase protein. These results suggest convergent evolution at these sites may have a functional role in clinical distemper outbreaks in African lions and uncover potential novel barriers to pathogenicity in this species.

Adaptive Molecular Evolution of AKT3 Gene for Positive Diversifying Selection in Mammals.

The V-Akt Murine Thymoma Viral Oncogene Homolog 3 (AKT3) gene is of the serine/threonine-protein kinase family and influences the production of milk fats and cholesterol by acting on the sterol administrative area restricting protein (SREBP). The AKT3 gene is highly preserved in animals, and during lactation in cattle, its expression increases. The AKT3 gene is expressed in the digestive system, mammary gland, and immune cells. A phylogenetic investigation was performed to clarify the evolutionary role of AKT3, by maximum probability. The AKT3 gene sequence data of various mammalian species was evident even with animals undergoing breeding selection. From 39 mammalian species studied, there was a signal of positive diversifying selection with Hominidae at 13Q, 16G, 23R, 24P, 121P, 294K, 327V, 376L, 397K, 445T, and 471F among other codon sites of the AKT3 gene. These sites were codes for amino acids such as arginine, proline, lysine, and leucine indicating major roles for the function of immunological proteins, and in particular, the study highlighted the importance of changes in gene expression of AKT3 on immunity.

Population-specific, recent positive directional selection suggests adaptation of human male reproductive genes to different environmental conditions

BackgroundRecent human transcriptomic analyses revealed a very large number of testis-enriched genes, many of which are involved in spermatogenesis. This comprehensive transcriptomic data lead us to the question whether positive selection was a decisive force influencing the evolution and variability of testis-enriched genes in humans. We used two methodological approaches to detect different levels of positive selection, namely episodic positive diversifying selection (i.e., past selection) in the human lineage within primate phylogeny, potentially driven by sperm competition, and recent positive directional selection in contemporary human populations, which would indicate adaptation to different environments.ResultsIn the human lineage (after correction for multiple testing) we found that only the gene TULP2, for which no functional data are yet available, is subject to episodic positive diversifying selection. Using less stringent statistical criteria (uncorrected p-values), also the gene SPATA16, which has a pivotal role in male fertility and for which episodes of adaptive evolution have been suggested, also displays a putative signal of diversifying selection in the human branch. At the same time, we found evidence for recent positive directional selection acting on several human testis-enriched genes (MORC1, SLC9B1, ROPN1L, DMRT1, PLCZ1, RNF17, FAM71D and WBP2NL) that play important roles in human spermatogenesis and fertilization. Most of these genes are population-specifically under positive selection.ConclusionEpisodic diversifying selection, possibly driven by sperm competition, was not an important force driving the evolution of testis-enriched genes in the human lineage. Population-specific, recent positive directional selection suggests an adaptation of male reproductive genes to different environmental conditions. Positive selection acts on eQTLS and sQTLs, indicating selective effects on important gene regulatory functions. In particular, the transcriptional diversity regulated by sQTLs in testis-enriched genes may be important for spermatocytes to respond to environmental and physiological stress.

Origin and Evolution of Two Independently Duplicated Genes Encoding UDP- Glucose: Glycoprotein Glucosyltransferases in Caenorhabditis and Vertebrates.

UDP- glucose: glycoprotein glucosyltransferase (UGGT) is a protein that operates as the gatekeeper for the endoplasmic reticulum (ER) quality control mechanism of glycoprotein folding. It is known that vertebrates and Caenorhabditis genomes harbor two uggt gene copies that exhibit differences in their properties.Bayesian phylogenetic inference based on 195 UGGT and UGGT-like protein sequences of an ample spectrum of eukaryotic species showed that uggt genes went through independent duplications in Caenorhabditis and vertebrates. In both lineages, the catalytic domain of the duplicated genes was subjected to a strong purifying selective pressure, while the recognition domain was subjected to episodic positive diversifying selection. Selective relaxation in the recognition domain was more pronounced in Caenorhabditis uggt-b than in vertebrates uggt-2. Structural bioinformatics analysis revealed that Caenorhabditis UGGT-b protein lacks essential sequences proposed to be involved in the recognition of unfolded proteins. When we assayed glucosyltrasferase activity of a chimeric protein composed by Caenorhabditis uggt-b recognition domain fused to S. pombe catalytic domain expressed in yeast, no activity was detected.The present results support the conservation of the UGGT activity in the catalytic domain and a putative divergent function of the recognition domain for the UGGT2 protein in vertebrates, which would have gone through a specialization process. In Caenorhabditis, uggt-b evolved under different constraints compared to uggt-a which, by means of a putative neofunctionalization process, resulted in a non-redundant paralog. The non-canonical function of uggt-b in the worm lineage highlights the need to take precautions before generalizing gene functions in model organisms.

Episodic positive diversifying selection on key immune system genes in major avian lineages.

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.