The unique evolution of the pig LRC, a single KIR but expansion of LILR and a novel Ig receptor family

A genome assembly and transcriptome atlas of the inbred Babraham pig to illuminate porcine immunogenetic variation

The inbred Babraham pig serves as a valuable biomedical model for research due to its high level of homozygosity, including in the major histocompatibility complex (MHC) loci and likely other important immune-related gene complexes, which are generally highly diverse in outbred populations. As the ability to control for this diversity using inbred organisms is of great utility, we sought to improve this resource by generating a long-read whole genome assembly and transcriptome atlas of a Babraham pig. The genome was de novo assembled using PacBio long reads and error-corrected using Illumina short reads. Assembled contigs were then mapped to the porcine reference assembly, Sscrofa11.1, to generate chromosome-level scaffolds. The resulting TPI_Babraham_pig_v1 assembly is nearly as contiguous as Sscrofa11.1 with a contig N50 of 34.95 Mb and contig L50 of 23. The remaining sequence gaps are generally the result of poor assembly across large and highly repetitive regions such as the centromeres and tandemly duplicated gene families, including immune-related gene complexes, that often vary in gene content between haplotypes. We also further confirm homozygosity across the Babraham MHC and characterize the allele content and tissue expression of several other immune-related gene complexes, including the antibody and T cell receptor loci, the natural killer complex, and the leukocyte receptor complex. The Babraham pig genome assembly provides an alternate highly contiguous porcine genome assembly as a resource for the livestock genomics community. The assembly will also aid biomedical and veterinary research that utilizes this animal model such as when controlling for genetic variation is critical.

Comparative genomics of the Leukocyte Receptor Complex in carnivores

Background

The mammalian Leukocyte Receptor Complex (LRC) chromosomal region may contain gene families for the killer cell immunoglobulin-like receptor (KIR) and/or leukocyte immunoglobulin-like receptor (LILR) collections as well as various framing genes. This complex region is well described in humans, mice, and some domestic animals. Although single KIR genes are known in some Carnivora, their complements of LILR genes remain largely unknown due to obstacles in the assembly of regions of high homology in short-read based genomes.

Methods

As part of the analysis of felid immunogenomes, this study focuses on the search for LRC genes in reference genomes and the annotation of LILR genes in Felidae. Chromosome-level genomes based on single-molecule long-read sequencing were preferentially sought and compared to representatives of the Carnivora.

Results

Seven putatively functional LILR genes were found across the Felidae and in the Californian sea lion, four to five genes in Canidae, and four to nine genes in Mustelidae. They form two lineages, as seen in the Bovidae. The ratio of functional genes for activating LILRs to inhibitory LILRs is slightly in favor of inhibitory genes in the Felidae and the Canidae; the reverse is seen in the Californian sea lion. This ratio is even in all of the Mustelidae except the Eurasian otter, which has a predominance of activating LILRs. Various numbers of LILR pseudogenes were identified.

Conclusions

The structure of the LRC is rather conservative in felids and the other Carnivora studied. The LILR sub-region is conserved within the Felidae and has slight differences in the Canidae, but it has taken various evolutionary paths in the Mustelidae. Overall, the process of pseudogenization of LILR genes seems to be more frequent for activating receptors. Phylogenetic analysis found no direct orthologues across the Carnivora which corroborate the rapid evolution of LILRs seen in mammals.

Leukocyte immunoglobulin-like receptor subfamily B: A novel immune checkpoint molecule at the maternal-fetal interface

Due to their crucial roles in embryo implantation, maternal-fetal tolerance induction, and pregnancy progression, immune checkpoint molecules (ICMs), such as programmed cell death-1, cytotoxic T-lymphocyte antigen 4, and T cell immunoglobulin mucin 3, are considered potential targets for clinical intervention in pregnancy complications. Despite the considerable progress on these molecules, our understanding of ICMs at the maternal-fetal interface is still limited. Identification of alternative and novel ICMs and the combination of multiple ICMs is urgently needed for deeply understanding the mechanism of maternal-fetal tolerance and to discover the causes of pregnancy complications. Leukocyte immunoglobulin-like receptor subfamily B (LILRB) is a novel class of ICMs with strong negative regulatory effects on the immune response. Recent studies have revealed that LILRB is enriched in decidual immune cells and stromal cells at the maternal-fetal interface, which can modulate the biological behavior of immune cells and promote immune tolerance. In this review, we introduce the structural features, expression profiles, ligands, and orthologs of LILRB. In addition, the potential mechanisms and functions mediated by LILRB for sustaining the maternal-fetal tolerance microenvironment, remodeling the uterine spiral artery, and induction of pregnancy immune memory are summarized. We have also provided new suggestions for further understanding the roles of LILRB and potential therapeutic strategies for pregnancy-related diseases.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Genetic Association between Farrowing Rates and Swine Leukocyte Antigen Alleles or Haplotypes in Microminipigs

We have previously reported specific swine leukocyte antigen (SLA) haplotype associations with significant effects on several reproduction performance traits in a highly inbred miniature pig population of Microminipigs (MMPs). In this study, to clarify the effects on farrowing rates of SLA similarity between mating partners in the MMP population, we compared the farrowing rates as a measure of reproductive success after 1063-cumulative matings among the following three groups of mating partners: (1) completely sharing SLA class I or class II haplotypes or alleles between partners (CS), (2) only one sharing the haplotypes or alleles (OS), and (3) non-sharing the haplotypes or alleles (NS). Average farrowing rates in CS groups consisting of completely sharing SLA class II haplotypes or DRBI and DQB1 alleles were lowest in the three groups. Moreover, lower farrowing rates were indicated in mating pairs with smaller amino acid pairwise genetic distances of SLA-1, SLA-3, DRB1 and DQB1 alleles between the pairs. These results suggested that the dissimilarity of SLA class I and class II alleles between mating partners markedly improved reproductive performance; therefore, SLA alleles or haplotypes are potentially useful genetic markers for the selection of mating pairs in breeding programs and epistatic studies of reproductive traits of MMPs.

Comparative Genomics of the Major Histocompatibility Complex (MHC) of Felids

This review summarizes the current knowledge on the major histocompatibility complex (MHC) of the family Felidae. This family comprises an important domestic species, the cat, as well as a variety of free-living felids, including several endangered species. As such, the Felidae have the potential to be an informative model for studying different aspects of the biological functions of MHC genes, such as their role in disease mechanisms and adaptation to different environments, as well as the importance of genetic diversity for conservation issues in free-ranging or captive populations. Despite this potential, the current knowledge on the MHC in the family as a whole is fragmentary and based mostly on studies of the domestic cat and selected species of big cats. The overall structure of the domestic cat MHC is similar to other mammalian MHCs following the general scheme "centromere-MHC class I-MHC class III-MHC class II" with some differences in the gene contents. An unambiguously defined orthologue of the non-classical class I HLA-E gene has not been identified so far and the class II DQ and DP genes are missing or pseudogenized, respectively. A comparison with available genomes of other felids showed a generally high level of structural and sequence conservation of the MHC region. Very little and fragmentary information on in vitro and/or in vivo biological functions of felid MHC genes is available. So far, no association studies have indicated effects of MHC genetic diversity on a particular disease. No information is available on the role of MHC class I molecules in interactions with Natural Killer (NK) cell receptors or on the putative evolutionary interactions (co-evolution) of the underlying genes. A comparison of complex genomic regions encoding NK cell receptors (the Leukocyte Receptor Complex, LRC and the Natural Killer Cell Complex, NKC) in the available felid genomes showed a higher variability in the NKC compared to the LRC and the MHC regions. Studies of the genetic diversity of domestic cat populations and/or specific breeds have focused mainly on DRB genes. Not surprisingly, higher levels of MHC diversity were observed in stray cats compared to pure breeds, as evaluated by DRB sequencing as well as by MHC-linked microsatellite typing. Immunogenetic analysis in wild felids has only been performed on MHC class I and II loci in tigers, Namibian leopards and cheetahs. This information is important as part of current conservation tasks to assess the adaptive potential of endangered wild species at the human-wildlife interface, which will be essential for preserving biodiversity in a functional ecosystem.

Beware of Mycoplasma Anti-immunoglobulin Strategies

Mycoplasmas are small, genome-reduced bacteria. They are obligate parasites that can be found in a wide range of host species, including the majority of livestock animals and humans. Colonization of the host can result in a wide spectrum of outcomes. In many cases, these successful parasites are considered commensal, as they are found in the microbiota of asymptomatic carriers. Conversely, mycoplasmas can also be pathogenic, as they are associated with a range of both acute and chronic inflammatory diseases which are problematic in veterinary and human medicine. The chronicity of mycoplasma infections and the ability of these bacteria to infect even recently vaccinated individuals clearly indicate that they are able to successfully evade their host’s humoral immune response. Over the years, multiple strategies of immune evasion have been identified in mycoplasmas, with a number of them aimed at generating important antigenic diversity. More recently, mycoplasma-specific anti-immunoglobulin strategies have also been characterized. Through the expression of the immunoglobulin-binding proteins protein M or mycoplasma immunoglobulin binding (MIB), mycoplasmas have the ability to target the host’s antibodies and to prevent them from interacting with their cognate antigens. In this review, we discuss how these discoveries shed new light on the relationship between mycoplasmas and their host’s immune system. We also propose that these strategies should be taken into consideration for future studies, as they are key to our understanding of mycoplasma diseases' chronic and inflammatory nature and are probably a contributing factor to reduce vaccine efficacy.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

The mycoplasma surface proteins MIB and MIP promote the dissociation of the antibody-antigen interaction

Mycoplasma immunoglobulin binding (MIB) and mycoplasma immunoglobulin protease (MIP) are surface proteins found in the majority of mycoplasma species, acting sequentially to capture antibodies and cleave off their V_H domains. Cryo-electron microscopy structures show how MIB and MIP bind to a Fab fragment in a "hug of death" mechanism. As a result, the orientation of the V_L and V_H domains is twisted out of alignment, disrupting the antigen binding site. We also show that MIB-MIP has the ability to promote the dissociation of the antibody-antigen complex. This system is functional in cells and protects mycoplasmas from antibody-mediated agglutination. These results highlight the key role of the MIB-MIP system in immunity evasion by mycoplasmas through an unprecedented mechanism, and open exciting perspectives to use these proteins as potential tools in the antibody field.

Natural Killer Cell Receptor Genes in Camels: Another Mammalian Model

Due to production of special homodimeric heavy chain antibodies, somatic hypermutation of their T-cell receptor genes and unusually low diversity of their major histocompatibility complex genes, camels represent an important model for immunogenetic studies. Here, we analyzed genes encoding selected natural killer cell receptors with a special focus on genes encoding receptors for major histocompatibility complex (MHC) class I ligands in the two domestic camel species, Camelus dromedarius and Camelus bactrianus. Based on the dromedary genome assembly CamDro2, we characterized the genetic contents, organization, and variability of two complex genomic regions, the leukocyte receptor complex and the natural killer complex, along with the natural cytotoxicity receptor genes NCR1, NCR2, and NCR3. The genomic organization of the natural killer complex region of camels differs from cattle, the phylogenetically most closely related species. With its minimal set of KLR genes, it resembles this complex in the domestic pig. Similarly, the leukocyte receptor complex of camels is strikingly different from its cattle counterpart. With KIR pseudogenes and few LILR genes, it seems to be simpler than in the pig. The syntenies and protein sequences of the NCR1, NCR2, and NCR3 genes in the dromedary suggest that they could be human orthologues. However, only NCR1 and NCR2 have a structure of functional genes, while NCR3 appears to be a pseudogene. High sequence similarities between the two camel species as well as with the alpaca Vicugna pacos were observed. The polymorphism in all genes analyzed seems to be generally low, similar to the rest of the camel genomes. This first report on natural killer cell receptor genes in camelids adds new data to our understanding of specificities of the camel immune system and its functions, extends our genetic knowledge of the innate immune variation in dromedaries and Bactrian camels, and contributes to studies of natural killer cell receptors evolution in mammals.

The Structure, Evolution, and Gene Expression Within the Caprine Leukocyte Receptor Complex

The leukocyte receptor complex (LRC) encodes a large number of immunoglobulin (Ig)-like receptors involved in the immune response, particularly in modulating natural killer (NK) cell function. The killer cell Ig-like receptors (KIR), the leukocyte Ig-like receptors (LILR), and a recently described novel Ig-like receptor family are highly variable between species, which is consistent with rapid evolution driven by selection pressure from pathogens. Among the species studied to date, only simians (such as humans) and bovids (such as cattle and goats) have an expanded complement of KIR genes and represent an interesting model to study KIR evolution. Using recently improved genome assemblies and an assembly of bacterial artificial chromosomes, we describe the structure of the LRC, and the KIR region in particular, in goats and compare this to sheep as the assemblies allow. These species diverged from a common ancestor ~10 million years ago and from cattle ~25 million years ago. We identified conserved KIR genes common to both goats and sheep and confirm a partial sheep haplotype shared between the Rambouillet and Texel breeds. Goats and sheep have independently expanded two novel KIR subgroups, and unlike cattle or any other mammal, they do not appear to possess a functional 3DL-lineage KIR gene. Investigation of LRC gene expression using available transcriptomic data for various sheep and goat tissues largely confirmed putative gene annotation and revealed that a relatively conserved caprinae-specific KIR subgroup is expressed in macrophages. The LILR and novel Ig-like receptors were also highly expressed across a diverse range of tissues. This further step toward our understanding of the LRC receptor repertoire will help inform future studies investigating immune response variation in these species.