Colon cancer associated genes exhibit signatures of positive selection at functionally significant positions

Evolutionary Origin of Germline Pathogenic MUTYH Variations in Modern Humans

MUTYH plays an essential role in preventing oxidation-caused DNA damage. Pathogenic germline variations in MUTYH damage its function, causing intestinal polyposis and colorectal cancer. Determination of the evolutionary origin of the variation is essential to understanding the etiological relationship between MUTYH variation and cancer development. In this study, we analyzed the origins of pathogenic germline variants in human MUTYH. Using a phylogenic approach, we searched pathogenic MUTYH variants in modern humans in the MUTYH of 99 vertebrates across eight clades. We did not find pathogenic variants shared between modern humans and the non-human vertebrates following the evolutionary tree, ruling out the possibility of cross-species conservation as the origin of human pathogenic variants in MUTYH. We then searched the variants in the MUTYH of 5031 ancient humans and extinct Neanderthals and Denisovans. We identified 24 pathogenic variants in 42 ancient humans dated between 30,570 and 480 years before present (BP), and three pathogenic variants in Neanderthals dated between 65,000 and 38,310 years BP. Data from our study revealed that human pathogenic MUTYH variants mostly arose in recent human history and were partially inherited from Neanderthals.

Divergent evolution and purifying selection of the Type 2 diabetes gene sequences in Drosophila: a phylogenomic study

The recently developed phylogenomic approach provides a unique way to identify disease risk or protective allele in any organism. While risk alleles evolve mostly under purifying selection, protective alleles are evolving either under balancing or positive selection. Owing to insufficient information, authors employed the phylogenomic approach to detect the nature of selection acting on type 2 diabetes (T2D) genes in Drosophila genus using various models of CODEML utility of PAML. The obtained result revealed that T2D gene sequences are evolving under purifying selection. However, only a few sites in membrane proteins encoded via CG8051, ZnT35C, and kar, are significantly evolving under positive selection under specific scenarios, which might be because of positive or adaptive evolution in response to changing niche, diet or other factors. In the near future, this information will be highly useful in the field of evolutionary medicine and the drug discovery process.

Resolving Peto's paradox: Modeling the potential effects of size-related metabolic changes, and of the evolution of immune policing and cancer suppression

The intrinsic risk of cancer increases with body size and longevity; however, big long-lived species do not exhibit this increase, a contradiction named Peto's paradox. Five hypotheses potentially resolving this paradox were modeled using the multistage model of carcinogenesis. The five hypotheses were based on (1) intrinsic changes in metabolic rate with body size; adaptive increase in immune policing of (2) cancer cells or (3) cells with driver mutations; or adaptive increase in cancer suppression via (4) decreased somatic mutation rate, or (5) increased genetic control. Parameter changes needed to stabilize cancer risk in three types of cancer were estimated for tissues scaled from mouse size and longevity to human and blue whale levels. The metabolic rate hypothesis alone was rejected due to a conflict between the required interspecific effect with the observed intraspecific effect of size on cancer risk, but some metabolic change was optionally incorporated in the other models. Necessary parameter changes in immune policing and somatic mutation rate far exceeded values observed; however, natural selection increasing the genetic suppression of cancer was generally consistent with data. Such adaptive increases in genetic control of cancers in large and/or long-lived animals raise the possibility that nonmodel animals will reveal novel anticancer mechanisms.

Positive Selection of Squalene Synthase in Cucurbitaceae Plants

Triterpenoid saponins are secondary metabolites synthesized through isoprenoid pathways in plants. Cucurbitaceae represent an important plant family in which many species contain cucurbitacins as secondary metabolites synthesized through isoprenoid and triterpenoid pathways. Squalene synthase (SQS) is required for the biosynthesis of isoprenoids, but the forces driving the evolution of SQS remain undetermined. In this study, 10 SQS cDNA sequences cloned from 10 species of Cucurbitaceae and 49 sequences of SQS downloaded from GenBank and UniProt databases were analyzed in a phylogenetic framework to identify the evolutionary forces for functional divergence. Through phylogenetic construction and positive selection analysis, we found that SQS sequences are under positive selection. The sites of positive selection map to functional and transmembrane domains. 180L, 189S, 194S, 196S, 265I, 289P, 389P, 390T, 407S, 408A, 410R, and 414N were identified as sites of positive selection that are important during terpenoid synthesis and map to transmembrane domains. 196S and 407S are phosphorylated and influence SQS catalysis and triterpenoid accumulation. These results reveal that positive selection is an important evolutionary force for SQS in plants. This provides new information into the molecular evolution of SQS within the Cucurbitaceae family.

New insights on unspecific peroxygenases: superfamily reclassification and evolution

BACKGROUND

Unspecific peroxygenases (UPO) (EC 1.11.2.1) represent an intriguing oxidoreductase sub-subclass of heme proteins with peroxygenase and peroxidase activity. With over 300 identified substrates, UPOs catalyze numerous oxidations including 1- or 2- electron oxygenation, selective oxyfunctionalizations, which make them most significant in organic syntheses and potentially attractive as industrial biocatalysts. There are very few UPOs available with distinct properties, notably, MroUPO which shows behavior ranging between UPO and another heme-thiolate peroxidase, called Chloroperoxidase (CPO). It prompted us to search for more UPOs in fungal kingdom which led us to studying their relationship with CPO.

RESULTS

In this study, we searched for novel UPOs in more than 800 fungal genomes and found 113 putative UPO-encoding sequences distributed in 35 different fungal species (or strains), amongst which single sequence per species were subjected to phylogeny study along with CPOs. Our phylogenetic study show that the UPOs are distributed in Basidiomycota and Ascomycota phyla of fungi. The sequence analysis helped to classify the UPOs into five distinct subfamilies: classic AaeUPO and four new subfamilies with potential new traits. We have also shown that each of these five subfamilies (supported by) have their own signature motifs. Surprisingly, some of the CPOs appeared to be a type of UPOs indicating that they were previously identified incorrectly. Selection pressure was observed on important motifs in UPOs which could have driven their functional divergence. Furthermore, the sites having different evolutionary rates caused by the functional divergence were also identified on some motifs along with the other relevant amino acid residues. Finally, we predicted critical amino acids responsible for the functional divergence in the UPOs and identified some sequence differences among UPOs, CPOs, and MroUPO to predict it's ranging behavior.

CONCLUSION

This study discovers new UPOs, provides a glimpse of their evolution from CPOs, and presents new insight on their functional divergence. We present a new classification of UPOs and shed new light on its phylogenetics. These different UPOs may exhibit a wide range of characteristics and specificities which may help in various fields of synthetic chemistry and industrial biocatalysts, and may as well lead to an advancement towards the understanding of physiological role of UPOs in fungi.

Comprehensive Analysis of the Chitinase Family Genes in Tomato (Solanum lycopersicum)

Chitinase catalyzes the hydrolysis of chitin β-1,4 linkages. However, plants cannot produce chitin, suggesting that plant chitinases do not have the same function as animals. This study investigated the chitinase gene family in tomato and divided into eight groups via phylogenetic analyses with Arabidopsis and rice members. Conserved gene structures and motif arrangements indicated their functional relevance with each group. These genes were nonrandomly distributed across the tomato chromosomes, and tandem duplication contributed to the expansion of this gene family. Synteny analysis also established orthology relationships and functional linkages between Arabidopsis and tomato chitinase genes. Several positive selection sites were identified, which may contribute to the functional divergence of the protein family in evolution. In addition, differential expression profiles of the tomato chitinase genes were also investigated at some developmental stages, or under different biotic and abiotic stresses. Finally, functional network analysis found 124 physical or functional interactions, implying the diversity of physiological functions of the family proteins. These results provide a foundation for the exploration of the chitinase genes in plants and will offer some insights for further functional studies.

Molecular Evolution of the Vacuolar Iron Transporter (VIT) Family Genes in 14 Plant Species

The vacuolar iron transporter (VIT) proteins are involved in the storage and transport of iron. However, the evolution of this gene family in plants is unknown. In this study, I first identified 114 VIT genes in 14 plant species and classified these genes into seven groups by phylogenetic analysis. Conserved gene organization and motif distribution implied conserved function in each group. I also found that tandem duplication, segmental duplication and transposition contributed to the expansion of this gene family. Additionally, several positive selection sites were identified. Divergent expression patterns of soybean VIT genes were further investigated in different development stages and under iron stress. Functional network analysis exhibited 211 physical or functional interactions. The results will provide the basis for further functional studies of the VIT genes in plants.

Positive selection and functional divergence of farnesyl pyrophosphate synthase genes in plants

Background

Farnesyl pyrophosphate synthase (FPS) belongs to the short-chain prenyltransferase family, and it performs a conserved and essential role in the terpenoid biosynthesis pathway. However, its classification, evolutionary history, and the forces driving the evolution of FPS genes in plants remain poorly understood.

Results

Phylogeny and positive selection analysis was used to identify the evolutionary forces that led to the functional divergence of FPS in plants, and recombinant detection was undertaken using the Genetic Algorithm for Recombination Detection (GARD) method. The dataset included 68 FPS variation pattern sequences (2 gymnosperms, 10 monocotyledons, 54 dicotyledons, and 2 outgroups). This study revealed that the FPS gene was under positive selection in plants. No recombinant within the FPS gene was found. Therefore, it was inferred that the positive selection of FPS had not been influenced by a recombinant episode. The positively selected sites were mainly located in the catalytic center and functional areas, which indicated that the 98S and 234D were important positively selected sites for plant FPS in the terpenoid biosynthesis pathway. They were located in the FPS conserved domain of the catalytic site. We inferred that the diversification of FPS genes was associated with functional divergence and could be driven by positive selection.

Conclusions

It was clear that protein sequence evolution via positive selection was able to drive adaptive diversification in plant FPS proteins. This study provides information on the classification and positive selection of plant FPS genes, and the results could be useful for further research on the regulation of triterpenoid biosynthesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-017-0081-4) contains supplementary material, which is available to authorized users.

Positive selection in cathelicidin host defense peptides: adaptation to exogenous pathogens or endogenous receptors?

The cause of adaptive protein evolution includes internal (for example, co-evolution of ligand-receptor pairs) and external (for example, adaptation to different ecological niches) mechanisms. Host defense peptides (HDPs) are a class of vertebrate-specific cationic antimicrobial peptides evolving under positive selection. Besides their antibiotic activity, HDPs also exert an effect on multiple host immune cells, thus providing an ideal model to study selective agents driving their evolution. On the basis of a combination of computational and experimental approaches, we studied the evolution of LL-37-type HDPs in mammals, the mature peptide of cathelicidin CAP18 (herein termed CAP18-MP) and investigated the driving force behind the evolution. Using codon-substitution maximum likelihood models, we analyzed CAP18-MPs in 40 species belonging to nine mammalian Orders and identified four positively selected sites (PSSs) that all are located on two terminal unordered regions of CAP18-MPs. Grafting the two positively selected regions of human or whale CAP18-MP to the α-helical scaffold of a rabbit homolog (substituting its corresponding parts) led to no alterations in antibacterial activity, spectrum and action mode. Likewise, further deletion of the two terminal regions did not alter its functional features. Evolutionary conservation analysis of mammalian FPR2, a receptor known to interact with the C-terminal positively selected region of LL-37, revealed high evolutionary variability in its ligand-binding extracellular loop domains, matching sequence diversity of the unordered regions in CAP18-MPs. This is the first report describing that the signature of positive selection of cathelicidins is not associated with their direct bactericidal activity, but rather with the evolutionary variability of their endogenous receptors.

Phylogenetic analyses reveal molecular signatures associated with functional divergence among Subtilisin like Serine Proteases are linked to lifestyle transitions in Hypocreales

Background

Subtilisin-like serine proteases or Subtilases in fungi are important for penetration and colonization of host. In Hypocreales, these proteins share several properties with other fungal, bacterial, plant and mammalian homologs. However, adoption of specific roles in entomopathogenesis may be governed by attainment of unique biochemical and structural features during the evolutionary course. Due to such functional shifts Subtilases coded by different family members of Hypocreales acquire distinct features according to respective hosts and lifestyle. We conducted phylogenetic and DIVERGE analyses and identified important protein residues that putatively assign functional specificity to Subtilases in fungal families/species under the order Hypocreales.

Results

A total of 161 Subtilases coded by 10 species from five different families under the fungal order Hypocreales was included in the analysis. Based on the presence of conserved domains, the Subtilase genes were divided into three subfamilies, Subtilisin (S08.005), Proteinase K (S08.054) and Serine-carboxyl peptidases (S53.001). These subfamilies were investigated for phylogenetic associations, protein residues under positive selection and functional divergence among paralogous clades. The observations were co-related with the life-styles of the fungal families/species. Phylogenetic and Divergence analyses of Subtilisin (S08.005) and Proteinase K (S08.054) families of proteins revealed that the paralogous clades were clear-cut representation of familial origin of the protein sequences. We observed divergence between the paralogous clades of plant-pathogenic fungi (Nectriaceae), insect-pathogenic fungi (Cordycipitaceae/Clavicipitaceae) and nematophagous fungi (Ophiocordycipitaceae). In addition, Subtilase genes from the nematode-parasitic fungus Purpureocillium lilacinum made a unique cluster which putatively indicated that the fungus might have developed distinctive mechanisms for nematode-pathogenesis. Our evolutionary genetics analysis revealed evidence of positive selection on the Subtilisin (S08.005) and Proteinase K (S08.054) protein sequences of the entomopathogenic and nematophagous species belonging to Cordycipitaceae, Clavicipitaceae and Ophiocordycipitaceae families of Hypocreales.

Conclusions

Our study provided new insights into the evolution of Subtilisin like serine proteases in Hypocreales, a fungal order largely consisting of biological control species. Subtilisin (S08.005) and Proteinase K (S08.054) proteins seemed to play important roles during life style modifications among different families and species of Hypocreales. Protein residues found significant in functional divergence analysis in the present study may provide support for protein engineering in future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0793-y) contains supplementary material, which is available to authorized users.

Positive Selection Linked with Generation of Novel Mammalian Dentition Patterns

A diverse group of genes are involved in the tooth development of mammals. Several studies, focused mainly on mice and rats, have provided a detailed depiction of the processes coordinating tooth formation and shape. Here we surveyed 236 tooth-associated genes in 39 mammalian genomes and tested for signatures of selection to assess patterns of molecular adaptation in genes regulating mammalian dentition. Of the 236 genes, 31 (∼13.1%) showed strong signatures of positive selection that may be responsible for the phenotypic diversity observed in mammalian dentition. Mammalian-specific tooth-associated genes had accelerated mutation rates compared with older genes found across all vertebrates. More recently evolved genes had fewer interactions (either genetic or physical), were associated with fewer Gene Ontology terms and had faster evolutionary rates compared with older genes. The introns of these positively selected genes also exhibited accelerated evolutionary rates, which may reflect additional adaptive pressure in the intronic regions that are associated with regulatory processes that influence tooth-gene networks. The positively selected genes were mainly involved in processes like mineralization and structural organization of tooth specific tissues such as enamel and dentin. Of the 236 analyzed genes, 12 mammalian-specific genes (younger genes) provided insights on diversification of mammalian teeth as they have higher evolutionary rates and exhibit different expression profiles compared with older genes. Our results suggest that the evolution and development of mammalian dentition occurred in part through positive selection acting on genes that previously had other functions.

Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection

Background

The study of local adaptation processes is a very important research topic in the field of population genomics. There is a particular interest in the study of human populations because they underwent a process of rapid spatial expansion and faced important environmental changes that translated into changes in selective pressures. New mutations may have been selected for in the new environment and previously existing genetic variants may have become detrimental. Immune related genes may have been released from the selective pressure exerted by pathogens in the ancestral environment and new variants may have been positively selected due to pathogens present in the newly colonized habitat. Also, variants that had a selective advantage in past environments may have become deleterious in the modern world due to external stimuli including climatic, dietary and behavioral changes, which could explain the high prevalence of some polygenic diseases such as diabetes and obesity.

Results

We performed an enrichment analysis to identify gene sets enriched for signals of positive selection in humans. We used two genome scan methods, XPCLR and iHS to detect selection using a dense coverage of SNP markers combined with two gene set enrichment approaches. We identified immune related gene sets that could be involved in the protection against pathogens especially in the African population. We also identified the glycolysis & gluconeogenesis gene set, related to metabolism, which supports the thrifty genotype hypothesis invoked to explain the current high prevalence of diseases such as diabetes and obesity. Extending our analysis to the gene level, we found signals for 23 candidate genes linked to metabolic syndrome, 13 of which are new candidates for positive selection.

Conclusions

Our study provides a list of genes and gene sets associated with immunity and metabolic syndrome that are enriched for signals of positive selection in three human populations (Europeans, Africans and Asians). Our results highlight differences in the relative importance of pathogens as drivers of local adaptation in different continents and provide new insights into the evolution and high incidence of metabolic syndrome in modern human populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2783-2) contains supplementary material, which is available to authorized users.

Target-Driven Positive Selection at Hot Spots of Scorpion Toxins Uncovers Their Potential in Design of Insecticides

Positive selection sites (PSSs), a class of amino acid sites with an excess of nonsynonymous to synonymous substitutions, are indicators of adaptive molecular evolution and have been detected in many protein families involved in a diversity of biological processes by statistical approaches. However, few studies are conducted to evaluate their functional significance and the driving force behind the evolution (i.e., agent of selection). Scorpion α-toxins are a class of multigene family of peptide neurotoxins affecting voltage-gated Na(+ )(Nav) channels, whose members exhibit differential potency and preference for insect and mammalian Nav channels. In this study, we undertook a systematical molecular dissection of nearly all the PSSs newly characterized in the Mesobuthus α-toxin family and a two-residue insertion ((19)AlaPhe(20)) located within a positively selected loop via mutational analysis of α-like MeuNaTxα-5, one member affecting both insect and mammalian Nav channels. This allows to identify hot-spot residues on its functional face involved in interaction with the receptor site of Nav channels, which comprises two PSSs (Ile(40) and Leu(41)) and the small insertion, both located on two spatially separated functional loops. Mutations at these hot-spots resulted in a remarkably decreased anti-mammalian activity in MeuNaTxα-5 with partially impaired or enhanced insecticide activity, suggesting the potential of PSSs in designing promising candidate insecticides from scorpion α-like toxins. Based on an experiment-guided toxin-channel complex model and high evolutionary variability in the receptor site of predators and prey of scorpions, we provide new evidence for target-driven adaptive evolution of scorpion toxins to deal with their targets' diversity.

GATA transcription factors in vertebrates: evolutionary, structural and functional interplay

GATA transcription factors perform conserved and essential roles during animal development, including germ-layer specification, hematopoiesis, and cardiogenesis. The evolutionary history and the changes in selection pressures following duplication of the six GATA family members in vertebrates have not been completely understood. Recently, we explored multiple databases to find GATAs in different vertebrate species. Using these sequences, we have performed molecular phylogenetic analyses using Maximum Likelihood and Bayesian methods, and statistical tests of tree topologies, to ascertain the phylogenetic relationship and selection pressures among GATA proteins. Seventy-one full-length cDNA sequences from 24 vertebrate species were extracted from multiple databases. By phylogenetic analyses, we investigated the origin, conservation, and evolution of the GATAs. Six GATA genes in vertebrates might be formed by gene duplication. The inferred evolutionary transitions that separate members which belong to different gene clusters correlated with changes in functional properties. Selection analysis and protein structure analysis were combined to explain Darwinian selection in GATA sequences and these changes brought putative biological significance. 26 positive selection sites were detected in this process. This study reveals the evolutionary history of vertebrate GATA paralogous and positively selected sites likely relevant for the distinct functional properties of the paralogs. It provides a new perspective for understanding the origin and evolution and biological functions of GATAs, which will help to uncover the GATAs' biological roles, evolution and their relationship with associated diseases; in addition, other complex multidomain families and also larger superfamilies can be investigated in a similar way.

Common genetic variants regulating ADD3 gene expression alter biliary atresia risk

BACKGROUND & AIMS

Biliary atresia (BA) is a rare and most severe cholestatic disease in neonates, but the pathogenic mechanisms are unknown. Through a previous genome wide association study (GWAS) on Han Chinese, we discovered association of the 10q24.2 region encompassing ADD3 and XPNPEP1 genes, which was replicated in Chinese and Thai populations. This study aims to fully characterize the genetic architecture at 10q24.2 and to reveal the link between the genetic variants and BA.

METHODS

We genotyped 107 single nucleotide polymorphisms (SNPs) in 10q24.2 in 339 Han Chinese patients and 401 matched controls using Sequenom. Exhaustive follow-up studies of the association signals were performed.

RESULTS

The combined BA-association p-value of the GWAS SNP (rs17095355) achieved 6.06×10(-10). Further, we revealed the common risk haplotype encompassing 5 tagging-SNPs, capturing the risk-predisposing alleles in 10q24.2 [p=5.32×10(-11); odds ratio, OR: 2.38; confidence interval, CI: (2.14-2.62)]. Through Sanger sequencing, no deleterious rare variants (RVs) residing in the risk haplotype were found, dismissing the theory of "synthetic" association. Moreover, in bioinformatics and in vivo genotype-expression investigations, the BA-associated potentially regulatory SNPs correlated with ADD3 gene expression (n=36; p=0.0030). Remarkably, the risk haplotype frequency coincides with BA incidences in the population, and, positive selection (favoring the derived alleles that arose from mutations) was evident at the ADD3 locus, suggesting a possible role for the BA-associated common variants in shaping the general population diversity.

CONCLUSIONS

Common genetic variants in 10q24.2 can alter BA risk by regulating ADD3 expression levels in the liver, and may exert an effect on disease epidemiology and on the general population.

Molecular adaptation of telomere associated genes in mammals

BACKGROUND

Placental mammals display a huge range of life history traits, including size, longevity, metabolic rate and germ line generation time. Although a number of general trends have been proposed between these traits, there are exceptions that warrant further investigation. Species such as naked mole rat, human and certain bat species all exhibit extreme longevity with respect to body size. It has long been established that telomeres and telomere maintenance have a clear role in ageing but it has not yet been established whether there is evidence for adaptation in telomere maintenance proteins that could account for increased longevity in these species.

RESULTS

Here we carry out a molecular investigation of selective pressure variation, specifically focusing on telomere associated genes across placental mammals. In general we observe a large number of instances of positive selection acting on telomere genes. Although these signatures of selection overall are not significantly correlated with either longevity or body size we do identify positive selection in the microbat species Myotis lucifugus in functionally important regions of the telomere maintenance genes DKC1 and TERT, and in naked mole rat in the DNA repair gene BRCA1.

CONCLUSION

These results demonstrate the multifarious selective pressures acting across the mammal phylogeny driving lineage-specific adaptations of telomere associated genes. Our results show that regardless of the longevity of a species, these proteins have evolved under positive selection thereby removing increased longevity as the single selective force driving this rapid rate of evolution. However, evidence of molecular adaptations specific to naked mole rat and Myotis lucifugus highlight functionally significant regions in genes that may alter the way in which telomeres are regulated and maintained in these longer-lived species.

Natural selection and functional diversification of the epidermal growth factor receptor EGFR family in vertebrates

BACKGROUND

Genes that have been subject to adaptive evolution can produce varying degrees of pathology or differing symptomatology. ErbB family receptor activation will initiate a number of downstream signaling pathways, such as mitogen-activated protein kinase (MAPK), activator of transcription (STAT), the modulation of calcium channels, and so on, all of which lead to aggressive tumor behavior. However, the evolutionary mechanisms operating in the retention of ErbB family genes and the changes in selection pressures are not clear.

RESULTS

Sixty-two full-length cDNA sequences from 27 vertebrate species were extracted from the UniProt protein database, NCBI's GenBank and the Ensembl database. The result of phylogenetic analysis showed that the four ErbB family members in vertebrates might be formed by gene duplication. In order to determine the mode of evolution in vertebrates, selection analysis and functional divergence analysis were combined to explain the relationship of the site-specific evolution and functional divergence in the vertebrate ErbB family. Our results indicate that the acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of ErbBs, and positive selections were detected in the ErbB family.

CONCLUSION

An evolutional phylogeny of 27 vertebrates was presented in our study; the tree showed that the genes have evolved through duplications followed by purifying selection, except for seven sites, which evolved by positive selection. There was one common site with positive selection and functional divergence. In the process of functional differentiation evolving through gene duplication, relaxed selection may play an important part.