Diversity and Molecular Evolution of Antimicrobial Peptides in Caecilian Amphibians.

Antimicrobial peptides (AMPs) are key molecules in the innate immune defence of vertebrates with rapid action, broad antimicrobial spectrum, and ability to evade pathogen resistance mechanisms. To date, amphibians are the major group of vertebrates from which most AMPs have been characterised, but most studies have focused on the bioactive skin secretions of anurans (frogs and toads). In this study, we have analysed the complete genomes and/or transcriptomes of eight species of caecilian amphibians (order Gymnophiona) and characterised the diversity, molecular evolution, and antimicrobial potential of the AMP repertoire of this order of amphibians. We have identified 477 candidate AMPs within the studied caecilian genome and transcriptome datasets. These candidates are grouped into 29 AMP families, with four corresponding to peptides primarily exhibiting antimicrobial activity and 25 potentially serving as AMPs in a secondary function, either in their entirety or after cleavage. In silico prediction methods were used to identify 62 of those AMPs as peptides with promising antimicrobial activity potential. Signatures of directional selection were detected for five candidate AMPs, which may indicate adaptation to the different selective pressures imposed by evolutionary arms races with specific pathogens. These findings provide encouraging support for the expectation that caecilians, being one of the least-studied groups of vertebrates, and with ~300 million years of separate evolution, are an underexplored resource of great pharmaceutical potential that could help to contest antibiotic resistance and contribute to biomedical advance.

Multiple Instances of Adaptive Evolution in Aquaporins of Amphibious Fishes.

Aquaporins (AQPs) are a highly diverse family of transmembrane proteins involved in osmotic regulation that played an important role in the conquest of land by tetrapods. However, little is known about their possible implication in the acquisition of an amphibious lifestyle in actinopterygian fishes. Herein, we investigated the molecular evolution of AQPs in 22 amphibious actinopterygian fishes by assembling a comprehensive dataset that was used to (1) catalogue AQP paralog members and classes; (2) determine the gene family birth and death process; (3) test for positive selection in a phylogenetic framework; and (4) reconstruct structural protein models. We found evidence of adaptive evolution in 21 AQPs belonging to 5 different classes. Almost half of the tree branches and protein sites that were under positive selection were found in the AQP11 class. The detected sequence changes indicate modifications in molecular function and/or structure, which could be related to adaptation to an amphibious lifestyle. AQP11 orthologues appear to be the most promising candidates to have facilitated the processes of the water-to-land transition in amphibious fishes. Additionally, the signature of positive selection found in the AQP11b stem branch of the Gobiidae clade suggests a possible case of exaptation in this clade.

Genomic Fishing and Data Processing for Molecular Evolution Research.

Molecular evolution analyses, such as detection of adaptive/purifying selection or ancestral protein reconstruction, typically require three inputs for a target gene (or gene family) in a particular group of organisms: sequence alignment, model of evolution, and phylogenetic tree. While modern advances in high-throughput sequencing techniques have led to rapid accumulation of genomic-scale data in public repositories and databases, mining such vast amount of information often remains a challenging enterprise. Here, we describe a comprehensive, versatile workflow aimed at the preparation of genome-extracted datasets readily available for molecular evolution research. The workflow involves: (1) fishing (searching and capturing) specific gene sequences of interest from taxonomically diverse genomic data available in databases at variable levels of annotation, (2) processing and depuration of retrieved sequences, (3) production of a multiple sequence alignment, (4) selection of best-fit model of evolution, and (5) solid reconstruction of a phylogenetic tree.