Early Diversification of Membrane Intrinsic Proteins (MIPs) in Eukaryotes.

Membrane intrinsic proteins (MIPs), including aquaporins (AQPs) and aquaglyceroporins (GLPs), form an ancient family of transporters for water and small solutes across biological membranes. The evolutionary history and functions of MIPs have been extensively studied in vertebrates and land plants, but their widespread presence across the eukaryotic tree of life suggests both a more complex evolutionary history and a broader set of functions than previously thought. That said, the early evolution of MIPs remains obscure. The presence of one GLP and four AQP clades across both bacteria and archaea suggests that the first eukaryotes could have possessed up to five MIPs. Here, we report on a previously unknown richness in MIP diversity across all major eukaryotic lineages, including unicellular eukaryotes, which make up the bulk of eukaryotic diversity. Three MIP clades have likely deep evolutionary origins, dating back to the last eukaryotic common ancestor (LECA), and support the presence of a complex MIP repertoire in early eukaryotes. Overall, our findings highlight the growing complexity of the reconstructed LECA genome: the dynamic evolutionary history of MIPs was set in motion when eukaryotes were in their infancy followed by radiative bursts across all main eukaryotic lineages.

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.

Feeding specialization and longer generation time are associated with relatively larger brains in bees.

Despite their miniature brains, insects exhibit substantial variation in brain size. Although the functional significance of this variation is increasingly recognized, research on whether differences in insect brain sizes are mainly the result of constraints or selective pressures has hardly been performed. Here, we address this gap by combining prospective and retrospective phylogenetic-based analyses of brain size for a major insect group, bees (superfamily Apoidea). Using a brain dataset of 93 species from North America and Europe, we found that body size was the single best predictor of brain size in bees. However, the analyses also revealed that substantial variation in brain size remained even when adjusting for body size. We consequently asked whether such variation in relative brain size might be explained by adaptive hypotheses. We found that ecologically specialized species with single generations have larger brains-relative to their body size-than generalist or multi-generation species, but we did not find an effect of sociality on relative brain size. Phylogenetic reconstruction further supported the existence of different adaptive optima for relative brain size in lineages differing in feeding specialization and reproductive strategy. Our findings shed new light on the evolution of the insect brain, highlighting the importance of ecological pressures over social factors and suggesting that these pressures are different from those previously found to influence brain evolution in other taxa.

Insights into the skin of caecilian amphibians from gene expression profiles.

BACKGROUND: Gene expression profiles can provide insights into the molecular machinery behind tissue functions and, in turn, can further our understanding of environmental responses, and developmental and evolutionary processes. During vertebrate evolution, the skin has played a crucial role, displaying a wide diversity of essential functions. To unravel the molecular basis of skin specialisations and adaptations, we compared gene expression in the skin with eight other tissues in a phylogenetically and ecologically diverse species sample of one of the most neglected vertebrate groups, the caecilian amphibians (order Gymnophiona). RESULTS: The skin of the five studied caecilian species showed a distinct gene expression profile reflecting its developmental origin and showing similarities to other epithelial tissues. We identified 59 sequences with conserved enhanced expression in the skin that might be associated with caecilian dermal specialisations. Some of the up-regulated genes shared expression patterns with human skin and potentially are involved in skin functions across vertebrates. Variation trends in gene expression were detected between mid and posterior body skin suggesting different functions between body regions. Several candidate biologically active peptides were also annotated. CONCLUSIONS: Our study provides the first atlas of differentially expressed sequences in caecilian tissues and a baseline to explore the molecular basis of the skin functions in caecilian amphibians, and more broadly in vertebrates.

What lies beneath? Molecular evolution during the radiation of caecilian amphibians.

BACKGROUND: Evolution leaves an imprint in species through genetic change. At the molecular level, evolutionary changes can be explored by studying ratios of nucleotide substitutions. The interplay among molecular evolution, derived phenotypes, and ecological ranges can provide insights into adaptive radiations. Caecilians (order Gymnophiona), probably the least known of the major lineages of vertebrates, are limbless tropical amphibians, with adults of most species burrowing in soils (fossoriality). This enigmatic order of amphibians are very distinct phenotypically from other extant amphibians and likely from the ancestor of Lissamphibia, but little to nothing is known about the molecular changes underpinning their radiation. We hypothesised that colonization of various depths of tropical soils and of freshwater habitats presented new ecological opportunities to caecilians. RESULTS: A total of 8540 candidate groups of orthologous genes from transcriptomic data of five species of caecilian amphibians and the genome of the frog Xenopus tropicalis were analysed in order to investigate the genetic machinery behind caecilian diversification. We found a total of 168 protein-coding genes with signatures of positive selection at different evolutionary times during the radiation of caecilians. The majority of these genes were related to functional elements of the cell membrane and extracellular matrix with expression in several different tissues. The first colonization of the tropical soils was connected to the largest number of protein-coding genes under positive selection in our analysis. From the results of our study, we highlighted molecular changes in genes involved in perception, reduction-oxidation processes, and aging that likely were involved in the adaptation to different soil strata. CONCLUSIONS: The genes inferred to have been under positive selection provide valuable insights into caecilian evolution, potentially underpin adaptations of caecilians to their extreme environments, and contribute to a better understanding of fossorial adaptations and molecular evolution in vertebrates.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics.

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

Multi-tissue transcriptomes of caecilian amphibians highlight incomplete knowledge of vertebrate gene families.

RNA sequencing (RNA-seq) has become one of the most powerful tools to unravel the genomic basis of biological adaptation and diversity. Although challenging, RNA-seq is particularly promising for research on non-model, secretive species that cannot be observed in nature easily and therefore remain comparatively understudied. Among such animals, the caecilians (order Gymnophiona) likely constitute the least known group of vertebrates, despite being an old and remarkably distinct lineage of amphibians. Here, we characterize multi-tissue transcriptomes for five species of caecilians that represent a broad level of diversity across the order. We identified vertebrate homologous elements of caecilian functional genes of varying tissue specificity that reveal a great number of unclassified gene families, especially for the skin. We annotated several protein domains for those unknown candidate gene families to investigate their function. We also conducted supertree analyses of a phylogenomic dataset of 1,955 candidate orthologous genes among five caecilian species and other major lineages of vertebrates, with the inferred tree being in agreement with current views of vertebrate evolution and systematics. Our study provides insights into the evolution of vertebrate protein-coding genes, and a basis for future research on the molecular elements underlying the particular biology and adaptations of caecilian amphibians.

Müller glia reactivity follows retinal injury despite the absence of the glial fibrillary acidic protein gene in Xenopus.

Intermediate filament proteins are structural components of the cellular cytoskeleton with cell-type specific expression and function. Glial fibrillary acidic protein (GFAP) is a type III intermediate filament protein and is up-regulated in glia of the nervous system in response to injury and during neurodegenerative diseases. In the retina, GFAP levels are dramatically increased in Müller glia and are thought to play a role in the extensive structural changes resulting in Müller cell hypertrophy and glial scar formation. In spite of similar changes to the morphology of Xenopus Müller cells following injury, we found that Xenopus lack a gfap gene. Other type III intermediate filament proteins were, however, significantly induced following rod photoreceptor ablation and retinal ganglion cell axotomy. The recently available X. tropicalis and X. laevis genomes indicate a small deletion most likely resulted in the loss of the gfap gene during anuran evolution. Lastly, a survey of representative species from all three extant amphibian orders including the Anura (frogs, toads), Caudata (salamanders, newts), and Gymnophiona (caecilians) suggests that deletion of the gfap locus occurred in the ancestor of all Anura after its divergence from the Caudata ancestor around 290 million years ago. Our results demonstrate that extensive changes in Müller cell morphology following retinal injury do not require GFAP in Xenopus, and other type III intermediate filament proteins may be involved in the gliotic response.