Detection Of Positive Selection Research Articles

Tree sequences as a general-purpose tool for population genetic inference.

As population genetics data increases in size new methods have been developed to store genetic information in efficient ways, such as tree sequences. These data structures are computationally and storage efficient, but are not interchangeable with existing data structures used for many population genetic inference methodologies such as the use of convolutional neural networks (CNNs) applied to population genetic alignments. To better utilize these new data structures we propose and implement a graph convolutional network (GCN) to directly learn from tree sequence topology and node data, allowing for the use of neural network applications without an intermediate step of converting tree sequences to population genetic alignment format. We then compare our approach to standard CNN approaches on a set of previously defined benchmarking tasks including recombination rate estimation, positive selection detection, introgression detection, and demographic model parameter inference. We show that tree sequences can be directly learned from using a GCN approach and can be used to perform well on these common population genetics inference tasks with accuracies roughly matching or even exceeding that of a CNN-based method. As tree sequences become more widely used in population genetics research we foresee developments and optimizations of this work to provide a foundation for population genetics inference moving forward.

Tree sequences as a general-purpose tool for population genetic inference.

As population genetics data increases in size new methods have been developed to store genetic information in efficient ways, such as tree sequences. These data structures are computationally and storage efficient, but are not interchangeable with existing data structures used for many population genetic inference methodologies such as the use of convolutional neural networks (CNNs) applied to population genetic alignments. To better utilize these new data structures we propose and implement a graph convolutional network (GCN) to directly learn from tree sequence topology and node data, allowing for the use of neural network applications without an intermediate step of converting tree sequences to population genetic alignment format. We then compare our approach to standard CNN approaches on a set of previously defined benchmarking tasks including recombination rate estimation, positive selection detection, introgression detection, and demographic model parameter inference. We show that tree sequences can be directly learned from using a GCN approach and can be used to perform well on these common population genetics inference tasks with accuracies roughly matching or even exceeding that of a CNN-based method. As tree sequences become more widely used in population genetics research we foresee developments and optimizations of this work to provide a foundation for population genetics inference moving forward.

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution

BackgroundSchistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected.MethodsWe obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima’s D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes.ResultsOur study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host.ConclusionsOur study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.Graphical

Genetic characterization, structural analysis, and detection of positive selection in small heat shock proteins of Cypriniformes and Clupeiformes.

In the current investigation, a total of 42 full-length, non-redundant small heat shock proteins (sHsp) were detected in Cyprinus carpio, Labeo rohita, Danio rerio, Salmo salar, Oncorhynchus mykiss, and Clupea harengus. The sHsp genes were classified into three groups based on phylogenetic analysis. All the sHsps were shown to have higher aliphatic index values, which is an indication that these proteins are more thermally stable. The hydrophilic nature of sHsps was deduced from the fact that all fish species had negative GRAVY scores. In all of the representative fish species, sHsp genes were assigned to distinct chromosomes in an inconsistent and unequal manner. Segmental duplications are the main events that have contributed to the expansion of the sHsp genes in all species. We were also able to determine the selective pressure that was placed on particular codons and discovered several significant coding sites within the coding region of sHsps. Eventually, diversifying positive selection was found to be connected with evolutionary changes in sHsp proteins, which showed that gene evolution controlled the fish adaption event in response to environmental conditions. Clarification of the links between sHsps and environmental stress in fish will be achieved through rigorous genomic comparison, which will also yield substantial new insights.

Two distinct Fusarium graminearum populations colonized European wheat in the past two decades.

Fusarium graminearum is the main causal agent of Fusarium head blight (FHB) disease in wheat in Europe. To reveal population structure and to pinpoint genetic targets of selection we studied genomes of 96 strains of F. graminearum using population genomics. Bayesian and phylogenomic analyses indicated that the F. graminearum emergence in Europe could be linked to two independently evolving populations termed here as East European (EE) and West European (WE) population. The EE strains are primarily prevalent in Eastern Europe, but to a lesser extent also in western and southern areas. In contrast, the WE population appears to be endemic to Western Europe. Both populations evolved in response to population-specific selection forces, resulting in distinct localized adaptations that allowed them to migrate into their environmental niche. The detection of positive selection in genes with protein/zinc ion binding domains, transcription factors and in genes encoding proteins involved in transmembrane transport highlights their important role in driving evolutionary novelty that allow F. graminearum to increase adaptation to the host and/or environment. F. graminearum also maintained distinct sets of accessory genes showing population-specific conservation. Among them, genes involved in host invasion and virulence such as those encoding proteins with high homology to tannase/feruloyl esterase and genes encoding proteins with functions related to oxidation-reduction were mostly found in the WE population. Our findings shed light on genetic features related to microevolutionary divergence of F. graminearum and reveal relevant genes for further functional research aiming at better control of this pathogen.

Tracking the Selective Pressure Profile and Gene Flow of SARS-CoV-2 Delta Variant in Italy from April to October 2021 and Frequencies of Key Mutations from Three Representative Italian Regions.

The SARS-CoV-2 Delta variant of concern (VOC) was often associated with serious clinical course of the COVID-19 disease. Herein, we investigated the selective pressure, gene flow and evaluation on the frequencies of mutations causing amino acid substitutions in the Delta variant in three Italian regions. A total of 1500 SARS-CoV-2 Delta genomes, collected in Italy from April to October 2021 were investigated, including a subset of 596 from three Italian regions. The selective pressure and the frequency of amino acid substitutions and the prediction of their possible impact on the stability of the proteins were investigated. Delta variant dataset, in this study, identified 68 sites under positive selection: 16 in the spike (23.5%), 11 in nsp2 (16.2%) and 10 in nsp12 (14.7%) genes. Three of the positive sites in the spike were located in the receptor-binding domain (RBD). In Delta genomes from the three regions, 6 changes were identified as very common (>83.7%), 4 as common (>64.0%), 21 at low frequency (2.1%-25.0%) and 29 rare (≤2.0%). The detection of positive selection on key mutations may represent a model to identify recurrent signature mutations of the virus.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments.

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

HaploBlocks: Efficient Detection of Positive Selection in Large Population Genomic Datasets.

Genomic regions under positive selection harbor variation linked for example to adaptation. Most tools for detecting positively selected variants have computational resource requirements rendering them impractical on population genomic datasets with hundreds of thousands of individuals or more. We have developed and implemented an efficient haplotype-based approach able to scan large datasets and accurately detect positive selection. We achieve this by combining a pattern matching approach based on the positional Burrows-Wheeler transform with model-based inference which only requires the evaluation of closed-form expressions. We evaluate our approach with simulations, and find it to be both sensitive and specific. The computational resource requirements quantified using UK Biobank data indicate that our implementation is scalable to population genomic datasets with millions of individuals. Our approach may serve as an algorithmic blueprint for the era of "big data" genomics: a combinatorial core coupled with statistical inference in closed form.

Detection of positive selection acting on protein surfaces at the whole-genome scale in the human malaria parasite Plasmodium falciparum.

The host-parasite evolutionary arms race is a fundamental process with medical implications. During this process, the host develops parasite resistance, and the parasite develops host immune evasion strategies. Thus, this process accelerates relevant protein evolution. This study test hypothesizes that proteins subject to sequence evolution structural constraints play a crucial role and that these constraints hinder the modification of such proteins in this process. These hypotheses were tested using Plasmodium falciparum model and evaluated protein structures predicted for the entire proteome by the AlphaFold method. Based on dN/dS test results and P. falciparum and P. reichenowi comparisons, the presented approach identified proteins subject to purifying selection acting on the whole sequence and buried residues (dN < dS) and positive selection on nonburied residues. Of the 26 proteins, some known antigens (ring-exported protein 3, RAP protein, erythrocyte binding antigen-140, and protein P47) targeted by the host immune system are promising vaccine candidates. The set also contained 11 enzymes, including FIKK kinase, which modifies host proteins. This set was compared with genes for which the dN/dS test suggested that positive selection acts on the whole gene (i.e., dN > dS). The present study found that such genes encode enzymes and antigenic vaccine candidates less frequently than genes for which evolution is not subject to selection constraints and positive selection acts on only exposed residues. The analysis was repeated comparing P. falciparum with P. alderi, which is more distantly related. The study discusses the potential implications of the presented methodology for rational vaccine design and the parasitology and evolutionary biology fields.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts.

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

Loss and fixation of strongly favoured new variants: Understanding and extending Haldane’s result via the Wright–Fisher model

The implications of Haldane’s analysis for the fixation of beneficial alleles lies at the heart of much of ‘population genetic thinking’ and underlies many approaches that have been tailored to the detection of positive selection. Within the framework of a branching process, Haldane gave an approximation for the probability that fixation ultimately occurs when the selective advantage of a beneficial allele is small (≪1). Here, we make no use of branching processes. Rather, we work solely within a finite-population Wright–Fisher framework. We use this framework to analyse where Haldane’s result applies, and extend Haldane’s analysis. In particular, we present results for: (i) the domain of applicability of Haldane’s analysis; (ii) the probability that loss occurs up to a given time; (iii) the probabilities that loss and fixation ultimately occur; (iv) an analytic approximation associated with the probability of loss and fixation ultimately occurring; (v) quantification of the crossover from weak to strong selection; (vi) determination of the number of invasive alleles that have a significant probability (>0.95) of invading a novel population. We note that the results obtained for (ii), (iii) and (iv) hold for an arbitrary initial number of mutations, and for selection that can be arbitrarily strong. Our results have fundamental implications for population and conservation genetics, and open up new avenues to identify traces of historically beneficial alleles through comparative genomics.

ImpMKT: the imputed McDonald and Kreitman test, a straightforward correction that significantly increases the evidence of positive selection of the McDonald and Kreitman test at the gene level.

The McDonald and Kreitman test is one of the most powerful and widely used methods to detect and quantify recurrent natural selection in DNA sequence data. One of its main limitations is the underestimation of positive selection due to the presence of slightly deleterious variants segregating at low frequencies. Although several approaches have been developed to overcome this limitation, most of them work on gene pooled analyses. Here, we present the imputed McDonald and Kreitman test (impMKT), a new straightforward approach for the detection of positive selection and other selection components of the distribution of fitness effects at the gene level. We compare imputed McDonald and Kreitman test with other widely used McDonald and Kreitman test approaches considering both simulated and empirical data. By applying imputed McDonald and Kreitman test to humans and Drosophila data at the gene level, we substantially increase the statistical evidence of positive selection with respect to previous approaches (e.g. by 50% and 157% compared with the McDonald and Kreitman test in Drosophila and humans, respectively). Finally, we review the minimum number of genes required to obtain a reliable estimation of the proportion of adaptive substitution (α) in gene pooled analyses by using the imputed McDonald and Kreitman test compared with other McDonald and Kreitman test implementations. Because of its simplicity and increased power to detect recurrent positive selection on genes, we propose the imputed McDonald and Kreitman test as the first straightforward approach for testing specific evolutionary hypotheses at the gene level. The software implementation and population genomics data are available at the web-server imkt.uab.cat.

A chromosome-level Populus qiongdaoensis genome assembly provides insights into tropical adaptation and a cryptic turnover of sex determination.

Populus species have long been used as model organisms to study the adaptability of trees and the evolution of sex chromosomes. As a species belonging to the section Populus and limited to tropical areas, the P. qiongdaoensis genome contains important information for tropical poplar studies and protection. Here, we report a chromosome-level genome assembly and annotation of a female P. qiongdaoensis. Gene family clustering, positive selection detection and historical reconstruction of population dynamics revealed the tropical adaptation of P. qiongdaoensis, and showed convergent evolution with another tropical poplar, P. ilicifolia, at the molecular level, especially on some functional genes (e.g., PIF3 and PIL1). In addition, we also identified a ZW sex determination system on chromosome 19 of P. qiongdaoensis, and inferred that it seems to have a similar sex determination mechanism to other poplars, controlled by a type-A cytokinin response regulator (RR) gene. However, comparison and phylogenetic analysis of the sex determination regions confirmed a cryptic sex turnover event in the section Populus, which may be caused by the translocation and duplication of the RR gene driven by Helitron-like transposable elements. Our study provides new insights into the environmental adaptation and sex chromosome evolution of poplars, and emphasizes the importance of using long read sequencing in ecological and evolutionary inferences of plants.

DGINN, an automated and highly-flexible pipeline for the detection of genetic innovations on protein-coding genes.

Adaptive evolution has shaped major biological processes. Finding the protein-coding genes and the sites that have been subjected to adaptation during evolutionary time is a major endeavor. However, very few methods fully automate the identification of positively selected genes, and widespread sources of genetic innovations such as gene duplication and recombination are absent from most pipelines. Here, we developed DGINN, a highly-flexible and public pipeline to Detect Genetic INNovations and adaptive evolution in protein-coding genes. DGINN automates, from a gene's sequence, all steps of the evolutionary analyses necessary to detect the aforementioned innovations, including the search for homologs in databases, assignation of orthology groups, identification of duplication and recombination events, as well as detection of positive selection using five methods to increase precision and ranking of genes when a large panel is analyzed. DGINN was validated on nineteen genes with previously-characterized evolutionary histories in primates, including some engaged in host-pathogen arms-races. Our results confirm and also expand results from the literature, including novel findings on the Guanylate-binding protein family, GBPs. This establishes DGINN as an efficient tool to automatically detect genetic innovations and adaptive evolution in diverse datasets, from the user's gene of interest to a large gene list in any species range.

Bioinformatic Detection of Positive Selection Pressure in Plant Pathogens: The Neutral Theory of Molecular Sequence Evolution in Action.

The genomes of plant pathogenic fungi and oomycetes are often exposed to strong positive selection pressure. During speciation, shifts in host range and preference can lead to major adaptive changes. Furthermore, evolution of total host resistance to most isolates can force rapid evolutionary changes in host-specific pathogens. Crop pathogens are subjected to particularly intense selective pressures from monocultures and fungicides. Detection of the footprints of positive selection in plant pathogen genomes is a worthwhile endeavor as it aids understanding of the fundamental biology of these important organisms. There are two main classes of test for detection of positively selected alleles. Tests based on the ratio of non-synonymous to synonymous substitutions per site detect the footprints of multiple fixation events between divergent lineages. Thus, they are well-suited to the study of ancient adaptation events spanning speciations. On the other hand, tests that scan genomes for local fluctuations in allelic diversity within populations are suitable for detection of recent positive selection in populations. In this review, I briefly describe some of the more widely used tests of positive selection and the theory underlying them. I then discuss various examples of their application to plant pathogen genomes, emphasizing the types of genes that are associated with signatures of positive selection. I conclude with a discussion of the practicality of such tests for identification of pathogen genes of interest and the important features of pathogen ecology that must be taken into account for accurate interpretation.

Coevolved Positions Represent Key Functional Properties in the Trypsin-Like Serine Proteases Protein Family

Trypsin-like serine proteases are a group of homologous enzymes which exert multiple roles in both vertebrate and invertebrate organisms. Key properties of these enzymes include their activation from an inactive zymogen form to their active form by cleavage of residues in their N-terminus, the presence of a conserved catalytic triad of residues, and the existence of different patterns of substrate selectivity for residue cleavage between the various members of this protein family. In this article, we apply the decomposition of residue coevolution networks computational method to find sets of residues related to some of these key properties, especially to zymogen activation. Positive selection detection, normal modes analysis, and the calculation of thermal couplings between the bovine trypsinogen and bovine trypsin structures residues yielded further information for understanding the zymogen activation process and highlighted the importance of some of the coevolved set residues during these transitions.

RAiSD-X

Detecting traces of positive selection in genomes carries theoretical significance and has practical applications from shedding light on the forces that drive adaptive evolution to the design of more effective drug treatments. The size of genomic datasets currently grows at an unprecedented pace, fueled by continuous advances in DNA sequencing technologies, leading to ever-increasing compute and memory requirements for meaningful genomic analyses. The majority of existing methods for positive selection detection either are not designed to handle whole genomes or scale poorly with the sample size; they inevitably resort to a runtime versus accuracy tradeoff, raising an alarming concern for the feasibility of future large-scale scans. To this end, we present RAiSD-X, a high-performance system that relies on a decoupled access-execute processing paradigm for efficient FPGA acceleration and couples a novel, to our knowledge, sliding-window algorithm for the recently introduced μ statistic with a mutation-driven hashing technique to rapidly detect patterns in the data. RAiSD-X achieves up to three orders of magnitude faster processing than widely used software implementations, and more importantly, it can exhaustively scan thousands of human chromosomes in minutes, yielding a scalable full-system solution for future studies of positive selection in species of flora and fauna.

Diet and Adaptive Evolution of Alanine-Glyoxylate Aminotransferase Mitochondrial Targeting in Birds.

Adaptations to different diets represent a hallmark of animal diversity. The diets of birds are highly variable, making them an excellent model system for studying adaptive evolution driven by dietary changes. To test whether molecular adaptations to diet have occurred during the evolution of birds, we examined a dietary enzyme alanine-glyoxylate aminotransferase (AGT), which tends to target mitochondria in carnivorous mammals, peroxisomes in herbivorous mammals, and both mitochondria and peroxisomes in omnivorous mammals. A total of 31 bird species were examined in this study, which included representatives of most major avian lineages. Of these, 29 have an intact mitochondrial targeting sequence (MTS) of AGT. This finding is in stark contrast to mammals, which showed a number of independent losses of the MTS. Our cell-based functional assays revealed that the efficiency of AGT mitochondrial targeting was greatly reduced in unrelated lineages of granivorous birds, yet it tended to be high in insectivorous and carnivorous lineages. Furthermore, we found that proportions of animal tissue in avian diets were positively correlated with mitochondrial targeting efficiencies that were experimentally determined, but not with those that were computationally predicted. Adaptive evolution of AGT mitochondrial targeting in birds was further supported by the detection of positive selection on MTS regions. Our study contributes to the understanding of how diet drives molecular adaptations in animals, and suggests that caution must be taken when computationally predicting protein subcellular targeting.

Cluster expansion of apolipoprotein D (ApoD) genes in teleost fishes

BackgroundGene and genome duplication play important roles in the evolution of gene function. Compared to individual duplicated genes, gene clusters attract particular attention considering their frequent associations with innovation and adaptation. Here, we report for the first time the expansion of the apolipoprotein D (ApoD) ligand-transporter genes in a cluster manner specific to teleost fishes.ResultsBased on comparative genomic and transcriptomic analyses, protein 3D structure comparison, positive selection detection and breakpoints detection, the single ApoD gene in the ancestor expanded into two clusters following a dynamic evolutionary pattern in teleost fishes. Orthologous genes show conserved expression patterns, whereas lineage-specific duplicated genes show tissue-specific expression patterns and even evolve new gene expression profiles. Positive selection occurred in branches before and after gene duplication, especially for lineage-specific duplicated genes. Cluster analyses based on protein 3D structure comparisons, especially comparisons of the four loops at the opening side, show gene duplication-segregating patterns. Duplicated ApoD genes are predicted to be associated with forkhead transcription factors and MAPK genes. ApoD clusters are located next to the breakpoints of genome rearrangements.ConclusionsHere, we report the expansion of ApoD genes specific to teleost fishes in a cluster manner for the first time. Neofunctionalization and subfunctionalization were observed at both the protein and expression levels after duplication. Evidence from different aspects—i.e., abnormal expression-induced disease in humans, fish-specific expansion, predicted associations with forkhead transcription factors and MAPK genes, specific expression patterns in tissues related to sexual selection and adaptation, duplicated genes under positive selection and their location next to the breakpoints of genome rearrangements—suggests the potentially advantageous roles of ApoD genes in teleost fishes. The cluster expansion of ApoD genes specific to teleost fishes provides thus an ideal evo-devo model for studying gene duplication, cluster maintenance and new gene function emergence.

Evidence of positive selection suggests possible role of aquaporins in the water-to-land transition of mudskippers

Aquaporins are integral membrane proteins that exchange water and small solutes. They played an important role in the colonisation of terrestrial environments by tetrapod ancestors via the appearance of three exclusive paralogs. Like early tetrapods, mudskippers represent an independent case of amphibious lifestyle evolution that is unparalleled by other extant fish groups. Given this lifestyle parallelism and that aquaporins were relevant for tetrapod terrestrialisation, this study examines the aquaporins in mudskippers to investigate whether similar changes in aquaporins could have possibly occurred during their water-to-land transition. We have catalogued aquaporin genes in four mudskipper genomes and studied their diversity and molecular evolution (including detection of positive selection) in a broad phylogenetic context of vertebrates. Our genomic screening returned 55 aquaporin genes for mudskippers (none of them constituting novel paralogs) that can be assigned to 10 different known classes. We detected signatures of positive selection in AQP10a and AQP11b in mudskippers (both the entire clade and the clade containing the most terrestrial species, implying different evolutionary times). This suggests possible alteration of the molecular function of such paralogs caused by changes at specific protein sequence positions, some of them located in relatively close proximity to parts of the molecule involved in pore formation and substrate selectivity. Given the importance of aquaporins for osmotic regulation in fishes, it might be possible that these selective changes (perhaps allowing permeability to new solutes) could have played a role during the adaptation of mudskippers to an amphibious lifestyle.