Adaptive Evolution of Toll-Like Receptors (TLRs) in the Family Suidae

Effect of population size and selection on Toll-like receptor diversity in populations of Galápagos mockingbirds

The interactions of evolutionary forces are difficult to analyse in free-living populations. However, when properly understood, they provide valuable insights into evolutionary biology and conservation genetics. This is particularly important for the interplay of genetic drift and natural selection in immune genes that confer resistance to disease. The Galápagos Islands are inhabited by four closely related species of mockingbirds (Mimus spp.). We used 12 different-sized populations of Galápagos mockingbirds and one population of their continental relative northern mockingbird (Mimus polyglottos) to study the effects of genetic drift on the molecular evolution of immune genes, the Toll-like receptors (TLRs: TLR1B, TLR4 and TLR15). We found that neutral genetic diversity was positively correlated with island size, indicating an important effect of genetic drift. However, for TLR1B and TLR4, there was little correlation between functional (e.g., protein) diversity and island size, and protein structural properties were largely conserved, indicating only a limited effect of genetic drift on molecular phenotype. By contrast, TLR15 was less conserved and even its putative functional polymorphism correlated with island size. The patterns observed for the three genes suggest that genetic drift does not necessarily dominate selection even in relatively small populations, but that the final outcome depends on the degree of selection constraint that is specific for each TLR locus.

African Suid Genomes Provide Insights into the Local Adaptation to Diverse African Environments

African wild suids consist of several endemic species that represent ancient members of the family Suidae and have colonized diverse habitats on the African continent. However, limited genomic resources for African wild suids hinder our understanding of their evolution and genetic diversity. In this study, we assembled high-quality genomes of a common warthog (Phacochoerus africanus), a red river hog (Potamochoerus porcus), as well as an East Asian Diannan small-ear pig (Sus scrofa). Phylogenetic analysis showed that common warthog and red river hog diverged from their common ancestor around the Miocene/Pliocene boundary, putatively predating their entry into Africa. We detected species-specific selective signals associated with sensory perception and interferon signaling pathways in common warthog and red river hog, respectively, which contributed to their local adaptation to savannah and tropical rainforest environments, respectively. The structural variation and evolving signals in genes involved in T-cell immunity, viral infection, and lymphoid development were identified in their ancestral lineage. Our results provide new insights into the evolutionary histories and divergent genetic adaptations of African suids.

CAR T Cell Locomotion in Solid Tumor Microenvironment

The promising outcomes of chimeric antigen receptor (CAR) T cell therapy in hematologic malignancies potentiates its capability in the fight against many cancers. Nevertheless, this immunotherapy modality needs significant improvements for the treatment of solid tumors. Researchers have incrementally identified limitations and constantly pursued better CAR designs. However, even if CAR T cells are armed with optimal killer functions, they must overcome and survive suppressive barriers imposed by the tumor microenvironment (TME). In this review, we will discuss in detail the important role of TME in CAR T cell trafficking and how the intrinsic barriers contribute to an immunosuppressive phenotype and cancer progression. It is of critical importance that preclinical models can closely recapitulate the in vivo TME to better predict CAR T activity. Animal models have contributed immensely to our understanding of human diseases, but the intensive care for the animals and unreliable representation of human biology suggest in vivo models cannot be the sole approach to CAR T cell therapy. On the other hand, in vitro models for CAR T cytotoxic assessment offer valuable insights to mechanistic studies at the single cell level, but they often lack in vivo complexities, inter-individual heterogeneity, or physiologically relevant spatial dimension. Understanding the advantages and limitations of preclinical models and their applications would enable more reliable prediction of better clinical outcomes.

Toll-Like Receptor Evolution: Does Temperature Matter?

Toll-like receptors (TLRs) recognize conserved pathogen-associated molecular patterns (PAMPs) and are an ancient and well-conserved group of pattern recognition receptors (PRRs). The isolation of the Antarctic continent and its unique teleost fish and microbiota prompted the present investigation into Tlr evolution. Gene homologues of tlr members in teleosts from temperate regions were present in the genome of Antarctic Nototheniidae and the non-Antarctic sister lineage Bovichtidae. Overall, in Nototheniidae apart from D. mawsoni, no major tlr gene family expansion or contraction occurred. Instead, lineage and species-specific changes in the ectodomain and LRR of Tlrs occurred, particularly in the Tlr11 superfamily that is well represented in fish. Positive selective pressure and associated sequence modifications in the TLR ectodomain and within the leucine-rich repeats (LRR), important for pathogen recognition, occurred in Tlr5, Tlr8, Tlr13, Tlr21, Tlr22, and Tlr23 presumably associated with the unique Antarctic microbiota. Exposure to lipopolysaccharide (Escherichia coli O111:B4) Gram negative bacteria did not modify tlr gene expression in N. rossii head–kidney or anterior intestine, although increased water temperature (+4°C) had a significant effect.

Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions

Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.

Hybridization with mountain hares increases the functional allelic repertoire in brown hares

Brown hares (Lepus europaeus Pallas) are able to hybridize with mountain hares (L. timidus Linnaeus) and produce fertile offspring, which results in cross-species gene flow. However, not much is known about the functional significance of this genetic introgression. Using targeted sequencing of candidate loci combined with mtDNA genotyping, we found the ancestral genetic diversity in the Finnish brown hare to be small, likely due to founder effect and range expansion, while gene flow from mountain hares constitutes an important source of functional genetic variability. Some of this variability, such as the alleles of the mountain hare thermogenin (uncoupling protein 1, UCP1), might have adaptive advantage for brown hares, whereas immunity-related MHC alleles are reciprocally exchanged and maintained via balancing selection. Our study offers a rare example where an expanding species can increase its allelic variability through hybridization with a congeneric native species, offering a route to shortcut evolutionary adaptation to the local environmental conditions.

TLR4 and TLR8 variability in Amazonian and West Indian manatee species from Brazil

Amazonian (Trichechus inunguis) and West Indian (Trichechus manatus) manatees are aquatic mammals vulnerable to extinction found in the Amazon basin and the coastal western Atlantic. Toll-like receptors (TLR) play a key role in recognizing pathogen-associated molecular patterns using leucine-rich repeats (LRRs). We described the diversity of TLR4 and TLR8 genes in these two species of manatee. Amazonian manatee showed seven SNPs in TLR4 and the eight in TLR8, while West Indian manatee shared four and six of those SNPs, respectively. In our analysis, TLR4 showed one non-conservative amino acid replacement substitution in LRR7 and LRR8, on the other hand, TLR8 was less variable and showed only conserved amino acid substitutions. Selection analysis showed that only one TLR4 site was subjected to positive selection and none in TLR8. TLR4 in manatees did not show any evidence of convergent evolution compared to species of the cetacean lineage. Differences in TLR4 and TLR8 polymorphism may be related to distinct selection by pathogens, population reduction of West Indian manatees, or an expected consequence of population expansion in Amazonian manatees. Future studies combining pathogen association and TLR polymorphism may clarify possible roles of these genes and be used for conservation purposes of manatee species.

Genetic drift and bottleneck do not influence diversity in Toll-like receptor genes at a small spatial scale in a Himalayan passerine

Genetic diversity is important for long-term viability of a population. Low genetic diversity reduces persistence and survival of populations and increases susceptibility to diseases. Comparisons of the neutral markers with functional loci such as immune genes [Toll-like receptors; TLR] can provide useful insights into evolutionary potential of a species and how the diversity of pathogens and selection pressures on their hosts are directly linked to their environment. In this study, we compare genetic diversity in neutral (eleven microsatellite loci) and adaptive (seven TLR loci) loci to determine genetic variation in a nonmigratory western Himalayan passerine, the black-throated tit (Aegithalos concinnus), distributed across an elevation gradient with varying degree of pathogen-mediated selection pressure. We further compare the diversity in TLR loci with a high-elevation sister species, the white-throated tit (Aegithalos niveogularis). Our results indicate a lack of population genetic structure in the black-throated tit and signatures of a past bottleneck. In contrast, we found high diversity in TLR loci and locus-specific (TLR7) signatures of pathogen-mediated selection, which was comparable to diversity in the white-throated tit. Levels of diversity at TLR5 locus corresponded very closely with neutral microsatellite variation. We found evidence of positive selection in TLR1LA, TLR5, and TLR7 loci highlighting the importance in pathogen recognition. Our finding demonstrates that reduction in neutral variation does not necessarily lead to reduction in functional genetic diversity and probably helps in revival of population in a widespread species.

Genetic diversity of toll-like receptor genes in the vulnerable Chinese egret (Egretta eulophotes)

Toll-like receptor (TLR) genes have recently been employed to assess genetic diversity, as they can be used to infer both demographic history and adaptation to environments with different pathogen pressure. Here, we sampled 120 individuals of the Chinese egret (Egretta eulophotes), a globally vulnerable species, from four breeding populations across China. We assessed the levels of genetic diversity, selection pressure, and population differentiation at seven TLR loci (TLR1LB, TLR2A, TLR3, TLR4, TLR5, TLR7, and TLR15). Using a variety of metrics (SNPs, heterozygosity, nucleotides, haplotypes), our analyses showed that genetic diversity was lower at 4 of the 7 TLR loci in the vulnerable Chinese egret compared to the more common little egret (Egretta garzetta). The selection test indicated TLRs, except for TLR5, were under purifying selection in TLR evolution, suggesting that low TLR genetic diversity in the Chinese egret may be caused by purifying selection. Moreover, analysis of molecular variance indicated low but significant population differentiation among four populations at all of the TLR loci in this egret. However, some comparisons based on fixation index analyses did not show significant population differentiation, and Bayesian clustering showed admixture. Our finding suggested that these four populations of the Chinese egret in China may be considered a single unit for conservation planning. These results, the new report of TLR genetic diversity in a long-distance migratory vulnerable Ardeid species, will provide fundamental TLR information for further studies on the conservation genetics of the Chinese egret and other Ardeids.

Strong selection of the TLR2 coding region among the Lagomorpha suggests an evolutionary history that differs from other mammals

Toll-like receptors (TLRs) are one of the first lines of defense against pathogens and are crucial for triggering an appropriate immune response. Among TLRs, TLR2 is functional in all vertebrates and has high ability in detecting bacterial and viral pathogen ligands. The mammals' phylogenetic tree of TLR2 showed longer branches for the Lagomorpha clade, raising the hypothesis that lagomorphs experienced an acceleration of the mutation rate. This hypothesis was confirmed by (i) Tajima's test of neutrality that revealed different evolutionary rates between lagomorphs and the remaining mammals with lagomorphs presenting higher nucleotide diversity; (ii) genetic distances were similar among lagomorphs and between lagomorphs and other mammals; and (iii) branch models reinforced the existence of an acceleration of the mutation rate in lagomorphs. These results suggest that the lagomorph TLR2 has been strongly involved in pathogen recognition, which probably caused a host-pathogen arms race that led to the observed acceleration of the mutation rate.

Widespread positive selection on cetacean TLR extracellular domain

Toll like receptors (TLRs), key members of innate immune system, can recognize a wide diversity of pathogens and initiate both innate and adaptive immune responses in vertebrate. Cetaceans must have faced new challenges of pathogens when their terrestrial relatives transitioned from the terrestrial to aquatic environment. Here, we sequenced the extracellular domain (ECD) of 10 TLRs in cetacean lineages because this region involved in the recognition of pathogens. A total of 148 sites ranging between 5-26 codons (0.01%-4.83%) were identified to be robust candidates of positive selection at the ECD of 10 TLRs. In addition, the majority (90.54%) of these positively selected codons were found to have radical amino acid changes, which strengthen the evidence of positive selection. Importantly, more radical amino acid changes in selected sites were enriched in the period of early evolutionary transition from land to semi-aquatic and from semi-aquatic to full-aquatic habitat, which might endow cetaceans with a faster adaptation to new pathogens as they transitioned into novel habitat. Interestingly, similar selective intensity was detected in both viral and non-viral TLRs in cetaceans, which was not in line with previous studies on primates and birds that reported stronger positive selection in non-viral TLRs than in viral TLRs. This result may be explained by the fact that cetaceans might have faced diversity of bacteria and viruses during its transitions from terrestrial to aquatic environment whereas both primates and birds probably being affected by only a restricted number of related viruses due to their homogeneous habitat.

The Oncopig Cancer Model: An Innovative Large Animal Translational Oncology Platform

Despite an improved understanding of cancer molecular biology, immune landscapes, and advancements in cytotoxic, biologic, and immunologic anti-cancer therapeutics, cancer remains a leading cause of death worldwide. More than 8.2 million deaths were attributed to cancer in 2012, and it is anticipated that cancer incidence will continue to rise, with 19.3 million cases expected by 2025. The development and investigation of new diagnostic modalities and innovative therapeutic tools is critical for reducing the global cancer burden. Toward this end, transitional animal models serve a crucial role in bridging the gap between fundamental diagnostic and therapeutic discoveries and human clinical trials. Such animal models offer insights into all aspects of the basic science-clinical translational cancer research continuum (screening, detection, oncogenesis, tumor biology, immunogenicity, therapeutics, and outcomes). To date, however, cancer research progress has been markedly hampered by lack of a genotypically, anatomically, and physiologically relevant large animal model. Without progressive cancer models, discoveries are hindered and cures are improbable. Herein, we describe a transgenic porcine model—the Oncopig Cancer Model (OCM)—as a next-generation large animal platform for the study of hematologic and solid tumor oncology. With mutations in key tumor suppressor and oncogenes, TP53^R167H and KRAS^G12D, the OCM recapitulates transcriptional hallmarks of human disease while also exhibiting clinically relevant histologic and genotypic tumor phenotypes. Moreover, as obesity rates increase across the global population, cancer patients commonly present clinically with multiple comorbid conditions. Due to the effects of these comorbidities on patient management, therapeutic strategies, and clinical outcomes, an ideal animal model should develop cancer on the background of representative comorbid conditions (tumor macro- and microenvironments). As observed in clinical practice, liver cirrhosis frequently precedes development of primary liver cancer or hepatocellular carcinoma. The OCM has the capacity to develop tumors in combination with such relevant comorbidities. Furthermore, studies on the tumor microenvironment demonstrate similarities between OCM and human cancer genomic landscapes. This review highlights the potential of this and other large animal platforms as transitional models to bridge the gap between basic research and clinical practice.

Polymorphisms of the Toll-Like Receptor 2 of Goats (Capra hircus) may be Associated with Somatic Cell Count in Milk

Toll-like receptor 2 (TLR2) plays an important role in recognition by the innate immune system of Gram-positive bacteria. As Gram-positive bacteria cause mastitis, we examined variations in the region of the TLR2 gene that codes for the extracellular domain. Samples of forty goats from a single dairy herd were collected, half with low SCC (≤200,000 cells/mL), and half with higher SCC. Two synonymous single nucleotide polymorphisms (SNPs) were identified: 840G > A and 1083A > G. One nonsynonymous SNP 739G > A was identified. This coded for valine or isoleucine, which have similar physiochemical properties, and was not in a region coding for a known functional domain. Surprisingly, the least square mean SCC of the heterozygous goats (146,220) was significantly lower than the SCC of homozygous GG goats (537,700; p = 0.004), although these two groups were similar in days in milk (p = 0.984), and there was no significant difference by breed (p = 0.941). Because factors other than mastitis can affect SCC and our sample sizes were limited, additional studies are needed to corroborate an association between TLR2 genotype and SCC or mastitis in goats.

Adaptive evolution of virus-sensing toll-like receptor 8 in bats

Recently, bats have gained attention as potential reservoir hosts for emerging zoonotic single-stranded (ssRNA) viruses that may prove fatal for humans and other mammals. It has been hypothesized that some features of their innate immune system may enable bats to trigger an efficient early immune response. Toll-like receptors (TLRs) represent a first line defense within the innate immune system and lie directly at the host–pathogen interface in targeting specific microbe-molecular patterns. However, the direction and strength of selection acting on TLRs are largely unknown for bats. Here, we studied the selection on viral ssRNA sensing TLR8 based on sequence data of 21 bat species. The major part (63 %) of the TLR8 gene evolved under purifying selection, likely due to functional constraints. We also found evidence for persistent positive selection acting on specific amino acid sites (7 %), especially when compared to viral TLR evolution of other mammals. All of these putatively positively selected codons were located in the ligand-binding ectodomain, some coincidenced or were in close proximity to functional sites, as suggested by the crystallographic structure of the human TLR8. This might contribute to the inter-species variation in the ability to recognize molecular patterns of viruses. TLR8 evolution within bats revealed that branches leading to ancestral and recent lineages evolved under episodic positive selection, indicating selective selection pressures in restricted bat lineages. Altogether, we found that the TLR8 displays extensive sequence variation within bats and that unique features separate them from humans and other mammals.

Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection

Background

The study of local adaptation processes is a very important research topic in the field of population genomics. There is a particular interest in the study of human populations because they underwent a process of rapid spatial expansion and faced important environmental changes that translated into changes in selective pressures. New mutations may have been selected for in the new environment and previously existing genetic variants may have become detrimental. Immune related genes may have been released from the selective pressure exerted by pathogens in the ancestral environment and new variants may have been positively selected due to pathogens present in the newly colonized habitat. Also, variants that had a selective advantage in past environments may have become deleterious in the modern world due to external stimuli including climatic, dietary and behavioral changes, which could explain the high prevalence of some polygenic diseases such as diabetes and obesity.

Results

We performed an enrichment analysis to identify gene sets enriched for signals of positive selection in humans. We used two genome scan methods, XPCLR and iHS to detect selection using a dense coverage of SNP markers combined with two gene set enrichment approaches. We identified immune related gene sets that could be involved in the protection against pathogens especially in the African population. We also identified the glycolysis & gluconeogenesis gene set, related to metabolism, which supports the thrifty genotype hypothesis invoked to explain the current high prevalence of diseases such as diabetes and obesity. Extending our analysis to the gene level, we found signals for 23 candidate genes linked to metabolic syndrome, 13 of which are new candidates for positive selection.

Conclusions

Our study provides a list of genes and gene sets associated with immunity and metabolic syndrome that are enriched for signals of positive selection in three human populations (Europeans, Africans and Asians). Our results highlight differences in the relative importance of pathogens as drivers of local adaptation in different continents and provide new insights into the evolution and high incidence of metabolic syndrome in modern human populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2783-2) contains supplementary material, which is available to authorized users.

Short-term MyD88 inhibition ameliorates cardiac graft rejection and promotes donor-specific hyporesponsiveness of skin grafts in mice

Recognition of evolutionarily conserved ligands by Toll-like receptors (TLRs) triggers signaling cascades in innate immune cells to amplify adaptive immune responses. Nearly all TLRs require MyD88 to transduce downstream signaling. MyD88 deficiency has been shown to promote the allograft acceptance in mice. However, direct evidence for therapeutic potential of MyD88 inhibitors remains lacking. Herein, we used a MyD88 inhibitor, namely ST2825, to explore its therapeutic potential and mechanisms in fully allogeneic skin and heart transplant models. Phenotypic maturation of dendritic cells stimulated by TLR ligands was alleviated by ST2825 in parallel with reduced T-cell proliferation in vitro. A short-course treatment with ST2825 significantly prolonged cardiac graft survival (mean survival time = 18.5 ± 0.92 days vs. 7.25 ± 0.46 days). ST2825-treated group had significantly reduced proinflammatory cytokines in allografts compared with control group. ST2825 combined with anti-CD154 induced long-term skin allograft acceptance in about one-third of recipients (>100 days). 'Skin-tolerant' recipients showed attenuated donor-specific IFN-γ responses, intact IL-4 responses, and compromised alloantibody responses. We conclude that MyD88 inhibitor ST2825 attenuates acute cardiac rejection and promotes donor-specific hyporesponsiveness in stringent skin transplant models. The direct evidence suggests that pharmacological inhibition of MyD88 hold promising potential for transplant rejection.