Evolutionary and functional characterization of lagomorph guanylate-binding proteins: a story of gain and loss and shedding light on expression, localization and innate immunity-related functions.

Guanylate binding proteins (GBPs) are an evolutionarily ancient family of proteins that are widely distributed among eukaryotes. They belong to the dynamin superfamily of GTPases, and their expression can be partially induced by interferons (IFNs). GBPs are involved in the cell-autonomous innate immune response against bacterial, parasitic and viral infections. Evolutionary studies have shown that GBPs exhibit a pattern of gene gain and loss events, indicative for the birth-and-death model of evolution. Most species harbor large GBP gene clusters that encode multiple paralogs. Previous functional and in-depth evolutionary studies have mainly focused on murine and human GBPs. Since rabbits are another important model system for studying human diseases, we focus here on lagomorphs to broaden our understanding of the multifunctional GBP protein family by conducting evolutionary analyses and performing a molecular and functional characterization of rabbit GBPs. We observed that lagomorphs lack GBP3, 6 and 7. Furthermore, Leporidae experienced a loss of GBP2, a unique duplication of GBP5 and a massive expansion of GBP4. Gene expression analysis by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and transcriptome data revealed that leporid GBP expression varied across tissues. Overexpressed rabbit GBPs localized either uniformly and/or discretely to the cytoplasm and/or to the nucleus. Oryctolagus cuniculus (oc)GBP5L1 and rarely ocGBP5L2 were an exception, colocalizing with the trans-Golgi network (TGN). In addition, four ocGBPs were IFN-inducible and only ocGBP5L2 inhibited furin activity. In conclusion, from an evolutionary perspective, lagomorph GBPs experienced multiple gain and loss events, and the molecular and functional characteristics of ocGBP suggest a role in innate immunity.

Evolution of primate interferon-induced transmembrane proteins (IFITMs): a story of gain and loss with a differentiation into a canonical cluster and IFITM retrogenes.

Interferon-inducible transmembrane proteins (IFITMs) are a family of transmembrane proteins. The subgroup of immunity-related (IR-)IFITMs is involved in adaptive and innate immune responses, being especially active against viruses. Here, we suggest that IFITMs should be classified as (1) a canonical IFITM gene cluster, which is located on the same chromosome, and (2) IFITM retrogenes, with a random and unique location at different positions within the genome. Phylogenetic analyses of the canonical cluster revealed the existence of three novel groups of primate IFITMs (pIFITM) in the IR-IFITM clade: the prosimian pIFITMs(pro), the new world monkey pIFITMs(nwm) and the old world monkey pIFITMs(owm). Therefore, we propose a new nomenclature: IR-pIFITM1, IR-pIFITM2, IR-pIFITM3, IR-pIFITMnwm, IR-pIFITMowm, and IR-pIFITMpro. We observed divergent evolution for pIFITM5 and pIFITM10, and evidence for concerted evolution and a mechanism of birth-and-death evolution model for the IR-pIFITMs. In contrast, the IFITMs scattered throughout the genomes possessed features of retrogenes retrotransposed by class 1 transposable elements. The origin of the IFITM retrogenes correspond to more recent events. We hypothesize that the transcript of a canonical IFITM3 has been constantly retrotransposed using class 1 transposable elements resulting in the IFITM retro(pseudo)genes. The unique pattern of each species has most likely been caused by constant pseudogenization and loss of the retro(pseudo)genes. This suggests a third mechanism of evolution for the IR-IFITMs in primates, similar to the birth-and-death model of evolution, but via a transposable element mechanism, which resulted in retro(pseudo)genes.

Functional cross-species conservation of guanylate-binding proteins in innate immunity.

Guanylate binding proteins (GBPs) represent an evolutionary ancient protein family widely distributed among eukaryotes. They are interferon (IFN)-inducible guanosine triphosphatases that belong to the dynamin superfamily. GBPs are known to have a major role in the cell-autonomous innate immune response against bacterial, parasitic and viral infections and are also involved in inflammasome activation. Evolutionary studies depicted that GBPs present a pattern of gain and loss of genes in each family with several genes pseudogenized and some genes more divergent, indicative for the birth-and-death evolution process. Most species harbor large GBP gene clusters encoding multiple paralogs. Previous functional studies mainly focused on mouse and human GBPs, but more data are becoming available, broadening the understanding of this multifunctional protein family. In this review, we will provide new insights and give a broad overview about GBP evolution, conservation and their roles in all studied species, including plants, invertebrates and vertebrates, revealing how far the described features of GBPs can be transferred to other species.

The rabbit as an animal model to study innate immunity genes: Is it better than mice?

The European rabbit (Oryctolagus cuniculus) was the first animal model used to understand human diseases like rabies and syphilis. Nowadays, the rabbit is still used to study several human infectious diseases like syphilis, HIV and papillomavirus. However, due to several mainly practical reasons, it has been replaced as an animal model by mice (Mus musculus). The rabbit and mouse share a recent common ancestor and are classified in the superorder Glires which arose at approximately 82 million years ago (mya). These species diverged from the Primates' ancestor at around 92 million years ago and, as such, one expects the rabbit-human and mouse-human genetic distances to be very similar. To evaluate this hypothesis, we developed a set of tools for automatic data extraction, sequence alignment and similarity study, and a web application for visualization of the resulting data. We aligned and calculated the genetic distances for 2793 innate immune system genes from human, rabbit and mouse using sequences available in the NCBI database. The obtained results show that the rabbit-human genetic distance is lower than the mouse-human genetic distance for 88% of these genes. Furthermore, when we considered only genes with a difference in genetic distance higher than 0.05, this figure increase to 93%. These results can be explained by the increase of the mutation rates in the mouse lineage suggested by some authors and clearly show that, at least looking to the genetic distance to human genes, the European rabbit is a better model to study innate immune system genes than the mouse.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages.

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

Evolution of Guanylate Binding Protein (GBP) Genes in Muroid Rodents (Muridae and Cricetidae) Reveals an Outstanding Pattern of Gain and Loss.

Guanylate binding proteins (GBPs) are paramount in the host immunity by providing defense against invading pathogens. Multigene families related to the immune system usually show that the duplicated genes can either undergo deletion, gain new functions, or become non-functional. Here, we show that in muroids, the Gbp genes followed an unusual pattern of gain and loss of genes. Muroids present a high diversity and plasticity regarding Gbp synteny, with most species presenting two Gbp gene clusters. The phylogenetic analyses revealed seven different Gbps groups. Three of them clustered with GBP2, GBP5 and GBP6 of primates. Four new Gbp genes that appear to be exclusive to muroids were identified as Gbpa, b, c and d. A duplication event occurred in the Gbpa group in the common ancestor of Muridae and Cricetidae (~20 Mya), but both copies were deleted from the genome of Mus musculus, M. caroli and Cricetulus griseus. The Gbpb gene emerged in the ancestor of Muridae and Cricetidae and evolved independently originating Gbpb1 in Muridae, Gbpb2 and Gbpb3 in Cricetidae. Since Gbpc appears only in three species, we hypothesize that it was present in the common ancestor and deleted from most muroid genomes. The second Gbp gene cluster, Gbp6, is widespread across all muroids, indicating that this cluster emerged before the Muridae and Cricetidae radiation. An expansion of Gbp6 occurred in M. musculus and M. caroli probably to compensate the loss of Gbpa and b. Gbpd is divided in three groups and is present in most muroids suggesting that a duplication event occurred in the common ancestor of Muridae and Cricetidae. However, in Grammomys surdaster and Mus caroli, Gbpd2 is absent, and in Arvicanthis niloticus, Gbpd1 appears to have been deleted. Our results further demonstrated that primate GBP1, GBP3 and GBP7 are absent from the genome of muroids and showed that the Gbp gene annotations in muroids were incorrect. We propose a new classification based on the phylogenetic analyses and the divergence between the groups. Extrapolations to humans based on functional studies of muroid Gbps should be re-evaluated. The evolutionary analyses of muroid Gbp genes provided new insights about the evolution and function of these genes.

Evolution of the guanylate binding protein (GBP) genes: Emergence of GBP7 genes in primates and further acquisition of a unique GBP3 gene in simians.

Guanylate binding proteins (GBPs) are major players in the host immunity, providing defense against bacterial and viral invaders. Multigene families may suffer different processes of evolution. Gene families related to the immune system usually follow the birth-and-death evolution process, where duplicated genes can be deleted, gain new functions or become non-functional. We analyzed publicly available primate GBP sequences and their genomic organization and observed that GBP7 genes appear to have emerged from a duplication of GBP4 and seem to be only present in primates. Furthermore, GBP3 genes are only present in Simiiformes and probably originated from GBP1 genes. Finally, a duplication event occurred in the GBP6 in Tarsiiformes and became functional which might also explain the duplication of GBP6 in New World monkeys and Cercopithecidae. Taken together, this study provides new knowledge on the evolution of GBPs in primates and suggests that a revision of the GBPs nomenclature is necessary.

Evolution of Fc Receptor-Like Scavenger in Mammals.

Fc receptor-like (FCRL) molecules comprise a large family of receptors, homologous to the receptors for the Fc portion of immunoglobulins (FCR). Within this family, an unusual gene known to exist in mice, rats and dogs, termed FCRLS, encodes a chimeric protein with both Ig-like FCRL and type B scavenger-receptor cysteine-rich (SRCR)-like domains. In mice, FCRLS is located next to the CD5L and KIRREL1 genes. Here, we show that the curious FCRLS gene is actually present across major mammalian groups, but its annotation is generally incorrect or absent. Anchored on mouse FCRLS and FCRL2 genomic sequence alignments, phylogenetic analyses demonstrated that many mammalian sequences currently annotated as FCRL2 cluster with FCRLS, supported by a conserved genetic synteny among organisms. This analysis shows that FCRLS is present in Rodentia, some Carnivora (Canidae and Ursidae), Chiroptera, Arctiodactyla, Proboscidae, and some Primata. Thus, the FCRLS most likely originated in a eutherian mammal ancestor since it is not present in Monotremata or Marsupialia. FCRLS has a peculiar distribution pattern across mammalian lineages, being present in some species, but absent in others from the same family, as in carnivores for example. The most parsimonious hypothesis to explain this FCRLS evolution is that it was convergently lost in several independent mammalian lineages. Analyses of branch-specific nucleotide evolutionary rates, show that FCRL2 and FCRLS have similar ranges of rates across mammals, suggesting that both genes have crucial, but separate functions in the immune system. Bayesian estimates of evolutionary rates for FCRLS in mammalian lineages revealed that carnivores display the highest mutation rate after rodents. Additionally, positive diversifying selection was detected for both FCRL2 and FCRLS. Our results show that the presence of the FCRLS gene is older and more widespread across mammals than previously thought and appears to be functional, being under positive selection. Its precise physiologic role should thus be investigated.

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.

Maximum likelihood approach suggests positive selection in platelet integrin αIIbß3 in mammalian species.

Platelet integrin αIIbß3 is crucial for platelet aggregation. Although structural and functional characteristics of this protein have been extensively studied, the evolutionary pattern studies of this protein complex in mammals are scarce. Here, we addressed this question using maximum likelihood approaches to identify codons that are evolving under positive selection. Likelihood of positive selection was estimated using CODEML implemented in PAML software applied to integrin αIIbß3 derived from nucleotide sequences of 10 different mammalian species. Four codons in mature αIIb-subunit (corresponding to residues 150, 184, 193, and 370) and three codons in mature ß3-subunit (corresponding to residues 129, 440, and 444) showed signs of positive selection with posterior probabilities over 95%. The different amino acids observed for each of the positively selected residues detected showed different physicochemical properties. These results open new research avenues to understand the physiological importance of specific residues and should allow for a better understanding of the function and the different interactions of each residue within the mature protein.

The evolution of S100A7 in primates: a model of concerted and birth-and-death evolution.

The human S100A7 resides in the epidermal differentiation complex (EDC) and has been described as a key effector of innate immunity. In humans, there are five S100A7 genes located in tandem-S100A7A, S100A7P1, S100AL2, S100A7, and S100AP2. The presence of several retroelements in the S100A7A/S100A7P1 and S100A7/S100A7P2 clusters suggests that these genes were originated from a duplication around ~ 35 million years ago, during or after the divergence of Platyrrhini and Catarrhini primates. To test this hypothesis, and taking advantage of the high number of genomic sequences available in the public databases, we retrieved S100A7 gene sequences of 12 primates belonging to the Cercopithecoidea and Hominoidea (Catarrhini species). Our results support the duplication theory, with at least one gene of each cluster being identified in both Cercopithecoidea and Hominoidea species. Moreover, given the presence of an ongoing gene conversion event between S100A7 and S100A7A, a high rate of mutation in S100A7L2 and the presence of pseudogenes, we proposed a model of concerted and birth-and-death evolution to explain the evolution of S100A7 gene family. Indeed, our results suggest that S100A7L2 most likely suffered a neofunctionalization in the Catarrhini group. Being S100A7 a major protein in innate defense, we believe that our findings could open new doors in the study of this gene family in immune system.

The antiviral activity of rodent and lagomorph SERINC3 and SERINC5 is counteracted by known viral antagonists.

A first step towards the development of a human immunodeficiency virus (HIV) animal model has been the identification and surmounting of species-specific barriers encountered by HIV along its replication cycle in cells from small animals. Serine incorporator proteins 3 (SERINC3) and 5 (SERINC5) were recently identified as restriction factors that reduce HIV-1 infectivity. Here, we compared the antiviral activity of SERINC3 and SERINC5 among mice, rats and rabbits, and their susceptibility to viral counteraction to their human counterparts. In the absence of viral antagonists, rodent and lagomorph SERINC3 and SERINC5 displayed anti-HIV activity in a similar range to human controls. Vesicular stomatitis virus G protein (VSV-G) pseudotyped virions were considerably less sensitive to restriction by all SERINC3/5 orthologs. Interestingly, HIV-1 Nef, murine leukemia virus (MLV) GlycoGag and equine infectious anemia virus (EIAV) S2 counteracted the antiviral activity of all SERINC3/5 orthologs with similar efficiency. Our results demonstrate that the antiviral activity of SERINC3/5 proteins is conserved in rodents and rabbits, and can be overcome by all three previously reported viral antagonists.

An overview of the lagomorph immune system and its genetic diversity.

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Pseudogenization of CCL14 in the Ochotonidae (pika) family.

The interaction between chemokines and their receptors is crucial for inflammatory cell trafficking. CCL14 binds with high affinity to CCR5. In leporids, CCR5 underwent gene conversion with CCR2. The study of CCR5 ligands in leporid species showed that CCL8 is pseudogenized, while CCL3, CCL4 and CCL5 are functional. Here, we study the evolution of CCL14 in mammals with emphasis in the order Lagomorpha. By employing maximum likelihood methods we detected six sites under positive selection. Some of these sites are located in regions crucial for CCL14 activation and binding to receptors. Sequencing of CCL14 in Ochotona species showed that O. princeps, O. pallasi, O. alpina and O. turuchanensis have a mutation at the start codon (Met > Thr), while O. hoffmanni, O. mantchurica, O. dauurica and O. rufescens present the mammalian conserved Met. Ochotona hyperborea has the two alleles. In O. pusilla, CCL14 is a pseudogene due to a seven base pair insertion. Like CCL3, CCL4 and CCL5, CCL14 is functional in all leporids but in the Ochotonidae family it underwent a pseudogenization process. This suggests that CCL14 has an important biological role in other mammals by evolving under positive selection that has been lost in Ochotonidae (subgenera Pika and Lagotona).

An intriguing shift occurs in the novel protein phosphatase 1 binding partner, TCTEX1D4: evidence of positive selection in a pika model.

T-complex testis expressed protein 1 domain containing 4 (TCTEX1D4) contains the canonical phosphoprotein phosphatase 1 (PPP1) binding motif, composed by the amino acid sequence RVSF. We identified and validated the binding of TCTEX1D4 to PPP1 and demonstrated that indeed this protein is a novel PPP1 interacting protein. Analyses of twenty-one mammalian species available in public databases and seven Lagomorpha sequences obtained in this work showed that the PPP1 binding motif 90RVSF93 is present in all of them and is flanked by a palindromic sequence, PLGS, except in three species of pikas (Ochotona princeps, O. dauurica and O. pusilla). Furthermore, for the Ochotona species an extra glycosylation site, motif 96NLS98, and the loss of the palindromic sequence were observed. Comparison with other lagomorphs suggests that this event happened before the Ochotona radiation. The dN/dS for the sequence region comprising the PPP1 binding motif and the flanking palindrome highly supports the hypothesis that for Ochotona species this region has been evolving under positive selection. In addition, mutational screening shows that the ability of pikas TCTEX1D4 to bind to PPP1 is maintained, although the PPP1 binding motif is disrupted, and the N- and C-terminal surrounding residues are also abrogated. These observations suggest pika as an ideal model to study novel PPP1 complexes regulatory mechanisms.