Evolution of CCL11: genetic characterization in lagomorphs and evidence of positive and purifying selection in mammals

Adaptive selection in the evolution of programmed cell death-1 and its ligands in vertebrates

Programmed cell death-1 (PD-1) and its ligands, particularly PD-L1 and PD-L2, are the most important proteins responsible for signaling T-cell inhibition and arbitrating immune homeostasis and tolerance mechanisms. However, the adaptive evolution of these genes is poorly understood. In this study, we aligned protein-coding genes from vertebrate species to evaluate positive selection constraints and evolution in the PD1, PD-L1 and PD-L2 genes conserved across up to 166 vertebrate species, with an average of 55 species per gene. We determined that although the positive selection was obvious, an average of 5.3% of codons underwent positive selection in the three genes across vertebrate lineages, and increased positive selection pressure was detected in both the Ig-like domains and transmembrane domains of the proteins. Moreover, the PD1, PD-L1 and PD-L2 genes were highly expressed in almost all tissues of the selected species indicating a distinct expression pattern in different tissues among most species. Our study reveals that adaptive selection plays a key role in the evolution of PD1 and its ligands in the majority of vertebrate species, which is in agreement with the contribution of these residues to the mechanisms of pathogen identification and coevolution in the complexity and novelties of vertebrate immune systems.

Not so unique to Primates: The independent adaptive evolution of TRIM5 in Lagomorpha lineage

The plethora of restriction factors with the ability to inhibit the replication of retroviruses have been widely studied and genetic hallmarks of evolutionary selective pressures in Primates have been well documented. One example is the tripartite motif-containing protein 5 alpha (TRIM5α), a cytoplasmic factor that restricts retroviral infection in a species-specific fashion. In Lagomorphs, similarly to what has been observed in Primates, the specificity of TRIM5 restriction has been assigned to the PRYSPRY domain. In this study, we present the first insight of an intra-genus variability within the Lagomorpha TRIM5 PRYSPRY domain. Remarkably, and considering just the 32 residue-long v1 region of this domain, the deduced amino acid sequences of Daurian pika (Ochotona dauurica) and steppe pika (O. pusilla) evidenced a high divergence when compared to the remaining Ochotona species, presenting values of 44% and 66% of amino acid differences, respectively. The same evolutionary pattern was also observed when comparing the v1 region of two Sylvilagus species members (47% divergence). However, and unexpectedly, the PRYSPRY domain of Lepus species exhibited a great conservation. Our results show a high level of variation in the PRYSPRY domain of Lagomorpha species that belong to the same genus. This suggests that, throughout evolution, the Lagomorpha TRIM5 should have been influenced by constant selective pressures, likely as a result of multiple different retroviral infections.

Evolution of CCL16 in Glires (Rodentia and Lagomorpha) shows an unusual random pseudogenization pattern

Background

The C-C motif chemokine ligand 16 (CCL16) is a potent pro-inflammatory chemokine and a chemoattractant for monocytes and lymphocytes. In normal plasma, it is present at high concentrations and elicits its effects on cells by interacting with cell surface chemokine receptors. In the European rabbit and in rodents such as mouse, rat and guinea pig, CCL16 was identified as a pseudogene, while in the thirteen-lined ground squirrel it appears to be potentially functional. To gain insight into the evolution of this gene in the superorder Glires (rodents and lagomorphs), we amplified the CCL16 gene from eleven Leporidae and seven Ochotonidae species.

Results

We compared our sequences with CCL16 sequences of twelve rodent species retrieved from public databases. The data show that for all leporid species studied CCL16 is a pseudogene. This is primarily due to mutations at the canonical Cys Cys motif, creating either premature stop codons, or disrupting amino acid replacements. In the Mexican cottontail, CCL16 is pseudogenized due to a frameshift deletion. Additionally, in the exon 1 (signal peptide), there are frameshift deletions present in all leporids studied. In contrast, in Ochotona species, CCL16 is potentially functional, except for an allele in Hoffmann’s pika. In rodents, CCL16 is functional in a number of species, but patterns of pseudogenization similar to those observed in lagomorphs also exist.

Conclusions

Our results suggest that while functional in the Glires ancestor, CCL16 underwent pseudogenization in some species. This process occurred stochastically or in specific lineages at different moments in the evolution of Glires. These observations suggest that the CCL16 had different evolutionary constrains in the Glires group that could be associated with the CCL16 biological function.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1390-7) contains supplementary material, which is available to authorized users.

Evolutionary studies on the betaretrovirus RERV-H in the Leporidae family reveal an endogenization in the ancestor of Oryctolagus, Bunolagus and Pentalagus at 9 million years ago

RERV-H was first identified in human tissues and mistaken for a human exogenous retrovirus. However, the integration sites carried by this virus showed that it was instead a European rabbit (Oryctolagus cuniculus) endogenous retrovirus. The first clones retrieved from European rabbit samples represented defective proviruses, although estimation of proviral copy numbers found in the European rabbit genome ranged from hundreds to thousands. Screening for the presence of RERV-H showed the absence of the virus in two other lagomorphs, pika (Ochotona) and hares (Lepus), which diverged from rabbits about 35 and 12 million years ago, respectively. Using a PCR-based approach, samples of seven different Lagomorph genera were tested for the presence of RERV-H. It was possible to amplify a proviral fragment corresponding to RNaseH from Oryctolagus, Bunolagus and Pentalagus genomic samples. The amplification of proviral DNA in species other than Oryctolagus revealed that this virus was endogenized in their common ancestor, roughly 9 million years ago. Using the European rabbit genome sequence OryCun2.0, it was possible to find multiple copies spread throughout the genome and several complete proviral genomes were retrieved. Some copies contained full open reading frames for all viral components. The lack of a complete genome in the other Lagomorph species did not allow further analyses of the provirus, although more deleterious mutations were found in Bunolagus and Pentalagus than in Oryctolagus RNaseH-amplified sequences. To what extent RERV-H and other endogenous viruses might have had an impact on the rabbit genome and its immune system remains elusive.

Mechanistic insights into molecular evolution of species-specific differential glycosaminoglycan binding surfaces in growth-related oncogene chemokines

Chemokines are chemotactic cytokines involved in leucocyte trafficking to infected tissue. Growth-related oncogene (GRO) chemokines namely CXCL1, CXCL2 and CXCL3 are neutrophil activating chemokines sharing a conserved three-dimensional structure, but encompassing functional diversity due to gene duplication and evolutionary events. However, the evolutionary mechanisms including selection pressures involved in diversification of GRO genes have not yet been characterized. Here, we performed comprehensive evolutionary analysis of GRO genes among different mammalian species. Phylogenetic analysis illustrated a species-specific evolution pattern. Selection analysis evidenced that these genes have undergone concerted evolution. Seventeen positively selected sites were obtained, although the majority of the protein is under purifying selection. Interestingly, these positively selected sites are more concentrated on the C-terminal/glycosaminoglycan (GAG) binding and dimerization segment compared to receptor binding domain. Substitution rate analysis confirmed the C-terminal domain of GRO genes as the highest substituted segment. Further, structural analysis established that the nucleotide alterations in the GAG binding domain are the source of surface charge modulation, thus generating the differential GAG binding surfaces and multiple binding sites as per evolutionary pressure, although the helical surface is primordial for GAG binding. Indeed, such variable electrostatic surfaces are crucial to regulate chemokine gradient formation during a host's defence against pathogens and also explain the significance of chemokine promiscuity.

The remnant of the European rabbit (Oryctolagus cuniculus) IgD gene

Although IgD first appeared, along with IgM, in the cartilaginous fishes and has been retained throughout subsequent vertebrate evolution, it has been lost in a diverse group of vertebrate species. We previously showed that, unlike vertebrates that express IgD, the rabbit lacks an IgD (Cδ) gene within 13.5 kb downstream of the IgM gene. We report here that, by conducting BLAST searches of rabbit Ig heavy chain genomic DNA with known mammalian IgD exons, we identified the remnant of the rabbit Cδ gene approximately 21 kb downstream of the IgM gene. The remnant Cδ locus lacks the δCH1 and hinge exons, but contains truncated δCH2 and δCH3 exons, as well as largely intact, but non-functional, secretory and transmembrane exons. In addition, we report that the Cδ gene probably became non-functional in leporids at least prior to the divergence of rabbits and hares ~12 million years ago.