Ana Lemos de Matos

Evolution and divergence of the mammalian SAMD9 / SAMD9L gene family

BackgroundThe physiological functions of the human Sterile Alpha Motif Domain-containing 9 (SAMD9) gene and its chromosomally adjacent paralogue, SAMD9-like (SAMD9L), currently remain unknown. However, the direct links between the deleterious mutations or deletions in these two genes and several human disorders, such as inherited inflammatory calcified tumors and acute myeloid leukemia, suggest their biological importance. SAMD9 and SAMD9L have also recently been shown to play key roles in the innate immune responses to stimuli such as viral infection. We were particularly interested in understanding the mammalian evolutionary history of these two genes. The phylogeny of SAMD9 and SAMD9L genes was reconstructed using the Maximum Likelihood method. Furthermore, six different methods were applied to detect SAMD9 and SAMD9L codons under selective pressure: the site-specific model M8 implemented in the codeml program in PAML software and five methods available on the Datamonkey web server, including the Single Likelihood Ancestor Counting method, the Fixed Effect Likelihood method, the Random Effect Likelihood method, the Mixed Effects Model of Evolution method and the Fast Unbiased Bayesian AppRoximation method. Additionally, the house mouse (Mus musculus) genome has lost the SAMD9 gene, while keeping SAMD9L intact, prompting us to investigate whether this loss is a unique event during evolution.ResultsOur evolutionary analyses suggest that SAMD9 and SAMD9L arose through an ancestral gene duplication event after the divergence of Marsupialia from Placentalia. Additionally, selection analyses demonstrated that both genes have been subjected to positive evolutionary selection. The absence of either SAMD9 or SAMD9L genes from some mammalian species supports a partial functional redundancy between the two genes.ConclusionsTo the best of our knowledge, this work is the first study on the evolutionary history of mammalian SAMD9 and SAMD9L genes. We conclude that evolutionary selective pressure has acted on both of these two genes since their divergence, suggesting their importance in multiple cellular processes, such as the immune responses to viral pathogens.

Positive Evolutionary Selection On the RIG-I-Like Receptor Genes in Mammals

The mammalian RIG-I-like receptors, RIG-I, MDA5 and LGP2, are a family of DExD/H box RNA helicases responsible for the cytoplasmic detection of viral RNA. These receptors detect a variety of RNA viruses, or DNA viruses that express unusual RNA species, many of which are responsible for a great number of severe and lethal diseases. Host innate sentinel proteins involved in pathogen recognition must rapidly evolve in a dynamic arms race with pathogens, and thus are subjected to long-term positive selection pressures to avoid potential infections. Using six codon-based Maximum Likelihood methods, we were able to identify specific codons under positive selection in each of these three genes. The highest number of positively selected codons was detected in MDA5, but a great percentage of these codons were located outside of the currently defined protein domains for MDA5, which likely reflects the imposition of both functional and structural constraints. Additionally, our results support LGP2 as being the least prone to evolutionary change, since the lowest number of codons under selection was observed for this gene. On the other hand, the preponderance of positively selected codons for RIG-I were detected in known protein functional domains, suggesting that pressure has been imposed by the vast number of viruses that are recognized by this RNA helicase. Furthermore, the RIG-I repressor domain, the region responsible for recognizing and binding to its RNA substrates, exhibited the strongest evidence of selective pressures. Branch-site analyses were performed and several species branches on the three receptor gene trees showed evidence of episodic positive selection. In conclusion, by looking for evidence of positive evolutionary selection on mammalian RIG-I-like receptor genes, we propose that a multitude of viruses have crafted the receptors biological function in host defense, specifically for the RIG-I gene, contributing to the innate species-specific resistance/susceptibility to diverse viral pathogens.

Evolution of viral sensing RIG-I-like receptor genes in Leporidae genera Oryctolagus, Sylvilagus, and Lepus

One of the most severe European rabbit (Oryctolagus cuniculus) pathogens is myxoma virus (MYXV), a rabbit-specific leporipoxvirus that causes the highly lethal disease myxomatosis. Other leporid genera, Sylvilagus and Lepus, encompass species with variable susceptibilities to MYXV, but these do not develop the lethal form of the disease. The protective role of the retinoic acid-inducible gene-I (RIG-I/DDX58) in sensing MYXV in nonpermissive human myeloid cells prompted the study of the RIG-I-like receptor (RLR) family evolution in the three leporid genera. This viral-sensor family also includes the melanoma differentiation-associated factor 5 (MDA5/IFIH1), and the laboratory of genetics and physiology 2 (LGP2/DHX58). Considering specifically the MYXV susceptible host (European rabbit) and one of the virus natural long-term hosts (Sylvilagus bachmani, brush rabbit), the amino acid differences of positively selected sites in RIG-I between the two species were located in the protein region responsible for viral RNA recognition and binding, the repressor domain. Such differences might play a determinant role in how MYXV is sensed. When looking for episodic selection on MDA5 and LGP2 of the eastern cottontail (Sylvilagus floridanus), we also uncovered evidence of selective pressures that might be exerted by a species-specific leporipoxvirus, the Shope fibroma virus. Finally, a putative alternative splicing case was identified in Oryctolagus and Lepus MDA5 isoforms, corresponding to the deletion of one specific exon. This study provided the first insights into the evolution of the leporid RLR gene family that helps illuminate the origins of the species-specific innate responses to pathogens and more specifically to MYXV.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

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 the MCP-2/CCL8 chemokine gene in European rabbit (genus Oryctolagus), but not in species of Cottontail rabbit (Sylvilagus) and Hare (Lepus).

BackgroundRecent studies in human have highlighted the importance of the monocyte chemotactic proteins (MCP) in leukocyte trafficking and their effects in inflammatory processes, tumor progression, and HIV-1 infection. In European rabbit (Oryctolagus cuniculus) one of the prime MCP targets, the chemokine receptor CCR5 underwent a unique structural alteration. Until now, no homologue of MCP-2/CCL8a, MCP-3/CCL7 or MCP-4/CCL13 genes have been reported for this species. This is interesting, because at least the first two genes are expressed in most, if not all, mammals studied, and appear to be implicated in a variety of important chemokine ligand-receptor interactions. By assessing the Rabbit Whole Genome Sequence (WGS) data we have searched for orthologs of the mammalian genes of the MCP-Eotaxin cluster.ResultsWe have localized the orthologs of these chemokine genes in the genome of European rabbit and compared them to those of leporid genera which do (i.e. Oryctolagus and Bunolagus) or do not share the CCR5 alteration with European rabbit (i.e. Lepus and Sylvilagus). Of the Rabbit orthologs of the CCL8, CCL7, and CCL13 genes only the last two were potentially functional, although showing some structural anomalies at the protein level. The ortholog of MCP-2/CCL8 appeared to be pseudogenized by deleterious nucleotide substitutions affecting exon1 and exon2. By analyzing both genomic and cDNA products, these studies were extended to wild specimens of four genera of the Leporidae family: Oryctolagus, Bunolagus, Lepus, and Sylvilagus. It appeared that the anomalies of the MCP-3/CCL7 and MCP-4/CCL13 proteins are shared among the different species of leporids. In contrast, whereas MCP-2/CCL8 was pseudogenized in every studied specimen of the Oryctolagus - Bunolagus lineage, this gene was intact in species of the Lepus - Sylvilagus lineage, and was, at least in Lepus, correctly transcribed.ConclusionThe biological function of a gene was often revealed in situations of dysfunction or gene loss. Infections with Myxoma virus (MYXV) tend to be fatal in European rabbit (genus Oryctolagus), while being harmless in Hares (genus Lepus) and benign in Cottontail rabbit (genus Sylvilagus), the natural hosts of the virus. This communication should stimulate research on a possible role of MCP-2/CCL8 in poxvirus related pathogenicity.

Study of Sylvilagus rabbit TRIM5α species-specific domain: how ancient endoviruses could have shaped the antiviral repertoire in Lagomorpha

BackgroundSince the first report of the antiretroviral restriction factor TRIM5α in primates, several orthologs in other mammals have been described. Recent studies suggest that leporid retroviruses like RELIK, the first reported endogenous lentivirus ever, may have imposed positive selection in TRIM5α orthologs of the European rabbit and European brown hare. Considering that RELIK must already have been present in a common ancestor of the leporid genera Lepus, Sylvilagus and Oryctolagus, we extended the study of evolutionary patterns of TRIM5α to other members of the Leporidae family, particularly to the genus Sylvilagus. Therefore, we obtained the TRIM5α nucleotide sequences of additional subspecies and species of the three leporid genera. We also compared lagomorph TRIM5α deduced protein sequences and established TRIM5α gene and TRIM5α protein phylogenies.ResultsThe deduced protein sequence of Iberian hare TRIM5α was 89% identical to European rabbit TRIM5α, although high divergence was observed at the PRYSPRY v1 region between rabbit and the identified alleles from this hare species (allele 1: 50% divergence; allele 2: 53% divergence). A high identity was expected between the Sylvilagus and Oryctolagus TRIM5α proteins and, in fact, the Sylvilagus TRIM5α was 91% identical to the Oryctolagus protein. Nevertheless, the PRYSPRY v1 region was only 50% similar between these genera. Selection analysis of Lagomorpha TRIM5α proteins identified 25 positively-selected codons, 11 of which are located in the PRYSPRY v1 region, responsible for species specific differences in viral capsid recognition.ConclusionsBy extending Lagomorpha TRIM5α studies to an additional genus known to bear RELIK, we verified that the divergent species-specific pattern observed between the Oryctolagus and Lepus PRYSPRY-domains is also present in Sylvilagus TRIM5α. This work is one of the first known studies that compare the evolution of the antiretroviral restriction factor TRIM5α in different mammalian groups, Lagomorpha and Primates.

Genetic characterization of interleukins (IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18) with relevant biological roles in lagomorphs

ILs, as essential innate immune modulators, are involved in an array of biological processes. In the European rabbit (Oryctolagus cuniculus) IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18 have been implicated in inflammatory processes and in the immune response against rabbit hemorrhagic disease virus and myxoma virus infections. In this study we characterized these ILs in six Lagomorpha species (European rabbit, pygmy rabbit, two cottontail rabbit species, European brown hare and American pika). Overall, these ILs are conserved between lagomorphs, including in their exon/intron structure. Most differences were observed between leporids and American pika. Indeed, when comparing both, some relevant differences were observed in American pika, such as the location of the stop codon in IL-1α and IL-2, the existence of a different transcript in IL8 and the number of cysteine residues in IL-1β. Changes at N-glycosylation motifs were also detected in IL-1, IL-10, IL-12B and IL-15. IL-1α is the protein that presents the highest evolutionary distances, which is in contrast to IL-12A where the distances between lagomorphs are the lowest. For all these ILs, sequences of human and European rabbit are more closely related than between human and mouse or European rabbit and mouse.

Ex Vivo Oncolytic Virotherapy with Myxoma Virus Arms Multiple Allogeneic Bone Marrow Transplant Leukocytes to Enhance Graft versus Tumor

Allogeneic stem cell transplant-derived T cells have the potential to seek and eliminate sites of residual cancer that escaped primary therapy. Oncolytic myxoma virus (MYXV) exhibits potent anti-cancer efficacy against human cancers like multiple myeloma (MM) and can arm transplant-derived T cells to become more effective cancer killers in vitro and in an immunodeficient xenotransplant murine model. Here, we tested ex vivo MYXV virotherapy against residual murine MM in immunocompetent mice using an allogeneic mouse-mouse model. In contrast to all human MM cell lines previously tested, the murine MM cell line tested here was highly resistant to direct MYXV infection and oncolysis in vitro. Despite this in vitro resistance, we found that ex vivo MYXV-armed allogeneic bone marrow (BM) transplantation dramatically ablated pre-seeded residual MM in vivo. Unexpectedly, we show that both neutrophils and activated T cells from the donor function as virus-armed carrier cells, and MYXV-preloaded cells enhanced MM killing. Our results demonstrate a novel therapeutic paradigm for residual cancer, in which multiple classes of allotransplant leukocytes can be armed by MYXV ex vivo to enhance the graft-versus-tumor effects.

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.