Evolutionary genomics of immunoglobulin-encoding Loci in vertebrates

Ruminant livestock TR V(D)J genes and CDR3 repertoire

Ruminant livestock exhibit certain immune characteristics that make them valuable models for studying T cell receptor diversity and immune responses. This resistance is attributed to their well-developed immune system, comprising both innate and adaptive components. In this review, we delve into the intricate workings of the immune system of ruminant livestock, focusing on innate immunity and adaptive immunity. Specifically, we discuss the TR V(D)J genes (including TRB, TRG, and TRA/D chain) and the characteristics of the complementary determining region 3 (CDR3) repertoire in bovine and ovine species, shedding light on the diversity and functionality of the T-cell receptor(TCR) repertoire in these species. Understanding the distinct features of these germline genes and CDR3 repertoires is essential for unraveling the complexities of immune responses in ruminant livestock. Lastly, we outline future prospects in this field, emphasizing the importance of further research to enhance our understanding of ruminant livestock immunity and its potential applications in disease management, vaccine development, and breeding strategies.

Molecular structure and functional responses of IgM, IgT and IgD to Flavobacterium covae and Streptococcus iniae infection in Asian seabass (Lates calcarifer Bloch, 1790)

The Asian seabass (Lates calcarifer) faces significant disease threats, which are exacerbated by intensive farming practices and environmental changes. Therefore, understanding its immune system is crucial. The current study presents a comprehensive analysis of immune-related genes in Asian seabass peripheral blood leukocytes (PBLs) using Iso-seq technology, identifying 16 key pathways associated with 7857 immune-related genes, comprising 634 unique immune-related genes. The research marks the first comprehensive report on the entire immunoglobulin repertoire in Asian seabass, revealing specific characteristics of immunoglobulin heavy chain constant region transcripts, including IgM (Cμ, ighm), IgT (Cτ, ight), and IgD (Cδ, ighd). The study confirms the presence of membrane-bound form, ighmmb, ightmb, ighdmb of IgM, IgT and IgD and secreted form, ighmsc and ightsc of IgM and IgT, respectively, with similar structural patterns and conserved features in amino acids across immunoglobulin molecules, including cysteine residues crucial for structural integrity observed in other teleost species. In response to bacterial infections by Flavobacterium covae (formerly F. columnare genomovar II) and Streptococcus iniae, both secreted and membrane-bound forms of IgM (ighmmb and ighmsc) and IgT (ightmb and ightsc) show significant expression, indicating their roles in systemic and mucosal immunity. The expression of membrane-bound form IgD gene, ighdmb, predominantly exhibits targeted upregulation in PBLs, suggesting a regulatory role in B cell-mediated immunity. The findings underscore the dynamic and tissue-specific expression of immunoglobulin repertoires, ighmmb, ighmsc, ightmb, ightsc and ighdmb in Asian seabass, indicating a sophisticated immune response to bacterial pathogens. These findings have practical implications for fish aquaculture, and disease control strategies, serving as a valuable resource for advancing research in Asian seabass immunology.

Bat humoral immunity and its role in viral pathogenesis, transmission, and zoonosis

Bats harbor viruses that can cause severe disease and death in humans including filoviruses (e.g., Ebola virus), henipaviruses (e.g., Hendra virus), and coronaviruses (e.g., SARS-CoV). Bats often tolerate these viruses without noticeable adverse immunological effects or succumbing to disease. Previous studies have largely focused on the role of the bat's innate immune response to control viral pathogenesis, but little is known about bat adaptive immunity. A key component of adaptive immunity is the humoral response, comprised of antibodies that can specifically recognize viral antigens with high affinity. The antibody genes within the 1,400 known bat species are highly diverse, and these genetic differences help shape fundamental aspects of the antibody repertoire, including starting diversity and viral antigen recognition. Whether antibodies in bats protect, mediate viral clearance, and prevent transmission within bat populations is poorly defined. Furthermore, it is unclear how neutralizing activity and Fc-mediated effector functions contribute to bat immunity. Although bats have canonical Fc genes (e.g., mu, gamma, alpha, and epsilon), the copy number and sequences of their Fc genes differ from those of humans and mice. The function of bat antibodies targeting viral antigens has been speculated based on sequencing data and polyclonal sera, but functional and biochemical data of monoclonal antibodies are lacking. In this review, we summarize current knowledge of bat humoral immunity, including variation between species, their potential protective role(s) against viral transmission and replication, and address how these antibodies may contribute to population dynamics within bats communities. A deeper understanding of bat adaptive immunity will provide insight into immune control of transmission and replication for emerging viruses with the potential for zoonotic spillover.

An Atlas of Grass Carp IgM+ B Cells in Homeostasis and Bacterial Infection Helps to Reveal the Unique Heterogeneity of B Cells in Early Vertebrates

Teleost B cells are primitive lymphocytes with both innate and adaptive immune functions. However, the heterogeneity and differentiation trajectory of teleost B cells remain largely unknown. In this study, the landscape of grass carp IgM+ (gcIgM+) B cells was revealed by single-cell RNA sequencing. The results showed that gcIgM+ B cells mainly comprise six populations: (im)mature B cells, innate B cells, proliferating B cells, plasma cells, CD22+ cells, and CD34+ cells, among which innate B cells and proliferating B cells were uncommon B cell subsets with, to our knowledge, new characteristics. Remarkably, three functional IgMs were discovered in grass carp, and a significant percentage of gcIgM+ B cells, especially plasma cells, expressed multiple Igμ genes (Igμ1, Igμ2, and/or Igμ3). More importantly, through single-cell sorting combined with Sanger sequencing, we found that distinct VHDJH recombination patterns of Igμ genes were present in single IgM+ B cells, indicating that individual teleost B cells might produce multiple Abs by coexpressing rearranged IgM subclass genes. Moreover, the percentage of IgM1highIgM2highIgM3high plasma cells increased significantly after bacterial infection, suggesting that individual plasma cells might tend to produce multiple IgMs to resist the infection in teleost fish. In summary, to our knowledge, this study not only helps to uncover the unique heterogeneity of B cells in early vertebrates but also provided significant new evidence supporting the recently proposed "one cell-multiple Abs" paradigm, challenging the classical rule of "one cell-one Ab."

Evolution of the Immunoglobulin Isotypes-Variations of Biophysical Properties among Animal Classes

The adaptive immune system arose around 500 million years ago in jawed fish, and, since then, it has mediated the immune defense against pathogens in all vertebrates. Antibodies play a central role in the immune reaction, recognizing and attacking external invaders. During the evolutionary process, several immunoglobulin isotypes emerged, each having a characteristic structural organization and dedicated function. In this work, we investigate the evolution of the immunoglobulin isotypes, in order to highlight the relevant features that were preserved over time and the parts that, instead, mutated. The residues that are coupled in the evolution process are often involved in intra- or interdomain interactions, meaning that they are fundamental to maintaining the immunoglobulin fold and to ensuring interactions with other domains. The explosive growth of available sequences allows us to point out the evolutionary conserved residues and compare the biophysical properties among different animal classes and isotypes. Our study offers a general overview of the evolution of immunoglobulin isotypes and advances the knowledge of their characteristic biophysical properties, as a first step in guiding protein design from evolution.

Evolutive emergence and divergence of an Ig regulatory node: An environmental sensor getting cues from the aryl hydrocarbon receptor?

One gene, the immunoglobulin heavy chain (IgH) gene, is responsible for the expression of all the different antibody isotypes. Transcriptional regulation of the IgH gene is complex and involves several regulatory elements including a large element at the 3' end of the IgH gene locus (3'RR). Animal models have demonstrated an essential role of the 3'RR in the ability of B cells to express high affinity antibodies and to express different antibody classes. Additionally, environmental chemicals such as aryl hydrocarbon receptor (AhR) ligands modulate mouse 3'RR activity that mirrors the effects of these chemicals on antibody production and immunocompetence in mouse models. Although first discovered as a mediator of the toxicity induced by the high affinity ligand 2,3,7,8-tetracholordibenzo-p-dioxin (dioxin), understanding of the AhR has expanded to a physiological role in preserving homeostasis and maintaining immunocompetence. We posit that the AhR also plays a role in human antibody production and that the 3'RR is not only an IgH regulatory node but also an environmental sensor receiving signals through intrinsic and extrinsic pathways, including the AhR. This review will 1) highlight the emerging role of the AhR as a key transducer between environmental signals and altered immune function; 2) examine the current state of knowledge regarding IgH gene regulation and the role of the AhR in modulation of Ig production; 3) describe the evolution of the IgH gene that resulted in species and population differences; and 4) explore the evidence supporting the environmental sensing capacity of the 3'RR and the AhR as a transducer of these cues. This review will also underscore the need for studies focused on human models due to the premise that understanding genetic differences in the human population and the signaling pathways that converge at the 3'RR will provide valuable insight into individual sensitivities to environmental factors and antibody-mediated disease conditions, including emerging infections such as SARS-CoV-2.

Gene prediction in the immunoglobulin loci

The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) is critical for studies of the adaptive immune system. However, since the standard gene finding algorithms are not suitable for predicting IG genes, they have been semi-manually annotated in very few species. We developed the IGDetective algorithm for predicting IG genes and applied it to species with the assembled IG loci. IGDetective generated the first large collection of IG genes across many species and enabled their evolutionary analysis, including the analysis of the "bat IG diversity" hypothesis. This analysis revealed extremely conserved V genes in evolutionary distant species indicating that these genes may be subjected to the same selective pressure, e.g., pressure driven by common pathogens. IGDetective also revealed extremely diverged V genes and a new family of evolutionary conserved V genes in bats with unusual noncanonical cysteines. Moreover, in difference from all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Since cysteines form disulfide bonds, we hypothesize that these cysteine-rich V genes might generate antibodies with noncanonical conformations and could potentially form a unique part of the immune repertoire in bats. We also analyzed the diversity landscape of the recombination signal sequences and revealed their features that trigger the high/low usage of the IG genes.

IgM Deficiency in Autoimmune Blistering Mucocutaneous Diseases Following Various Treatments: Long Term Follow-Up and Relevant Observations

IgM deficiency has been reported in patients with many autoimmune diseases treated with Rituximab (RTX). It has not been studied, in detail, in autoimmune mucocutaneous blistering diseases (AIMBD). Our objectives were: (i) Examine the dynamics of IgM levels in patients with and without RTX. (ii) Influence of reduced serum IgM levels on clinical and laboratory parameters. (iii) Explore the possible molecular and cellular basis for reduced serum IgM levels. This retrospective study that was conducted in a single-center from 2000 to 2020. Serial IgM levels were studied in 348 patients with five AIMBD (pemphigus vulgaris, pemphigus foliaceus, bullous pemphigoid, mucous membrane pemphigoid, and ocular cicatricial pemphigoid) and found decreased in 55 patients treated with RTX, IVIG, and conventional immunosuppressive therapy (CIST). Hence the incidence of decreased serum IgM is low. The incidence of decreased IgM in patients treated with RTX was 19.6%, in patients treated with IVIG and CIST, it was 52.8% amongst the 55 patients. IgM levels in the post-RTX group were statistically significantly different from the IVIG group (p<0.018) and CIST group (p<0.001). There were no statistically significant differences between the groups in other clinical and laboratory measures. Decreased serum IgM did not affect depletion or repopulation of CD19+ B cells. Patients in the three groups achieved clinical and serological remission, in spite of decreased IgM levels. Decrease in IgM was isolated, since IgG and IgA were normal throughout the study period. Decreased IgM persisted at the same level, while the patients were in clinical remission, for several years. In spite of persistent decreased IgM levels, the patients did not develop infections, tumors, other autoimmune diseases, or warrant hospitalization. Studies on IgM deficiency in knockout mice provided valuable insights. There is no universally accepted mechanism that defines decreased IgM levels in AIMBD. The data is complex, multifactorial, sometimes contradictory, and not well understood. Nonetheless, data in this study provides novel information that enhances our understanding of the biology of IgM in health and disease.

The evolutionary puzzle solution for the origins of the partial loss of the Cτ2 exon in notothenioid fishes

Cryonotothenioidea is the main group of fishes that thrive in the extremely cold Antarctic environment, thanks to the acquisition of peculiar morphological, physiological and molecular adaptations. We have previously disclosed that IgM, the main immunoglobulin isotype in teleosts, display typical cold-adapted features. Recently, we have analyzed the gene encoding the heavy chain constant region (CH) of the IgT isotype from the Antarctic teleost Trematomus bernacchii (family Nototheniidae), characterized by the near-complete deletion of the CH2 domain. Here, we aimed to track the loss of the CH2 domain along notothenioid phylogeny and to identify its ancestral origins. To this end, we obtained the IgT gene sequences from several species belonging to the Antarctic families Nototheniidae, Bathydraconidae and Artedidraconidae. All species display a CH2 remnant of variable size, encoded by a short Cτ2 exon, which retains functional splicing sites and therefore is included in the mature transcript. We also considered representative species from the three non-Antarctic families: Eleginopsioidea (Eleginops maclovinus), Pseudaphritioidea (Pseudaphritis urvillii) and Bovichtidae (Bovichtus diacanthus and Cottoperca gobio). Even though only E. maclovinus, the sister taxa of Cryonotothenioidea, shared the partial loss of Cτ2, the other non-Antarctic notothenioid species displayed early molecular signatures of this event. These results shed light on the evolutionary path that underlies the origins of this remarkable gene structural modification.

On Origin and Evolution of the Antibody Molecule

Simple Summary

Like many other molecules playing vital functions in animals, the antibody molecule possesses a complex structure with distinctive features. The structure of the basic unit, i.e., the immunoglobulin domain of very ancient origin is substantially simple. However, high complexity resides in the types and numbers of the domains composing the whole molecule. The emergence of the antibody molecule during evolution overturned the effectiveness of the organisms’ defense system. The particular organization of the coding genes, the mechanisms generating antibody diversity, and the plasticity of the overall protein structure, attest to an extraordinary successful evolutionary history. Here, we attempt to trace, across the evolutionary scale, the very early origins of the most significant features characterizing the structure of the antibody molecule and of the molecular mechanisms underlying its major role in recognizing an almost unlimited number of pathogens.

Abstract

The vertebrate immune system provides a powerful defense because of the ability to potentially recognize an unlimited number of pathogens. The antibody molecule, also termed immunoglobulin (Ig) is one of the major mediators of the immune response. It is built up from two types of Ig domains: the variable domain, which provides the capability to recognize and bind a potentially infinite range of foreign substances, and the constant domains, which exert the effector functions. In the last 20 years, advances in our understanding of the molecular mechanisms and structural features of antibody in mammals and in a variety of other organisms have uncovered the underlying principles and complexity of this fundamental molecule. One notable evolutionary topic is the origin and evolution of antibody. Many aspects have been clearly stated, but some others remain limited or obscure. By considering a wide range of prokaryotic and eukaryotic organisms through a literature survey about the topic, we have provided an integrated view of the emergence of antibodies in evolution and underlined the very ancient origins.

Programmed DNA Elimination in Vertebrates

Over the last few decades, an increasing number of vertebrate taxa have been identified that undergo programmed genome rearrangement, or programmed DNA loss, during development. In these organisms, the genome of germ cells is often reproducibly different from the genome of all other cells within the body. Although we clearly have not identified all vertebrate taxa that undergo programmed genome loss, the list of species known to undergo loss now represents ∼10% of vertebrate species, including several basally diverging lineages. Recent studies have shed new light on the targets and mechanisms of DNA loss and their association with canonical modes of DNA silencing. Ultimately, expansion of these studies into a larger collection of taxa will aid in reconstructing patterns of shared/independent ancestry of programmed DNA loss in the vertebrate lineage, as well as more recent evolutionary events that have shaped the structure and content of eliminated DNA.

Analysis of structural variants in four African cichlids highlights an association with developmental and immune related genes

BACKGROUND

East African lake cichlids are one of the most impressive examples of an adaptive radiation. Independently in Lake Victoria, Tanganyika, and Malawi, several hundreds of species arose within the last 10 million to 100,000 years. Whereas most analyses in cichlids focused on nucleotide substitutions across species to investigate the genetic bases of this explosive radiation, to date, no study has investigated the contribution of structural variants (SVs) in the evolution of adaptive traits across the three Great Lakes of East Africa.

RESULTS

Here, we annotate and characterize the repertoires and evolutionary potential of different SV classes (deletion, duplication, inversion, insertions and translocations) in four cichlid species: Haplochromis burtoni, Metriaclima zebra, Neolamprologus brichardi and Pundamilia nyererei. We investigate the patterns of gain and loss evolution for each SV type, enabling the identification of lineage specific events. Both deletions and inversions show a significant overlap with SINE elements, while inversions additionally show a limited, but significant association with DNA transposons. Inverted regions are enriched for genes regulating behaviour, or involved in skeletal and visual system development. We also find that duplicated regions show enrichment for genes associated with "antigen processing and presentation" and other immune related categories. Our pipeline and results were further tested by PCR validation of selected deletions and inversions, which confirmed respectively 7 out of 10 and 6 out of 9 events.

CONCLUSIONS

Altogether, we provide the first comprehensive overview of rearrangement evolution in East African cichlids, and some important insights into their likely contribution to adaptation.

A comprehensive analysis of the genomic organization, expression and phylogeny of immunoglobulin light chain genes in pigeon (Columba livia)

Previous studies on immunoglobulin light chain (IgL) genes in avian species are limited to Galloanseres, and few studies have investigated IgL genes in Neoaves, which includes most living birds. Based on published genome data, we demonstrate that the pigeon (Columba livia) IgL locus spans approximately 24 kb of DNA and contains twenty Vλ segments located upstream of a single pair of Jλ-Cλ. Among the identified Vλ gene segments, four segments are structurally intact and all four segments are able to recombine with Jλ. Moreover, the four functional Vλ segments are preferentially utilized in VλJλ recombination. Phylogenetic analysis suggests that the presence of the four functional Vλ segments in pigeon was likely generated by gene duplication that occurred after the divergence of pigeon and other birds. Our study provides insight into IgL gene evolution and evolutionary diversity of Ig genes in birds.

Submegabase copy number variations arise during cerebral cortical neurogenesis as revealed by single-cell whole-genome sequencing

Reports of copy number variations (CNVs) within single human brain cells have been limited to megabase-scale alterations in relatively few cells, leaving unclear when CNVs first arise and whether their generation is regulated. Answering these questions has been limited by an absence of experimentation with model organisms that allow developmental assessments infeasible with human samples. Here, we identify the existence and developmental dynamics of cerebral cortical CNVs in mouse, showing that their prevalence increases through midneurogenesis. Our improved sequencing approach also allowed characterization of previously undocumented neural CNVs below 1 Mb in size, comprising half of all alterations. These data demonstrate the existence of myriad CNVs, which genomically diversify neural cells before incorporation into the mature organization of the brain.

Somatic copy number variations (CNVs) exist in the brain, but their genesis, prevalence, forms, and biological impact remain unclear, even within experimentally tractable animal models. We combined a transposase-based amplification (TbA) methodology for single-cell whole-genome sequencing with a bioinformatic approach for filtering unreliable CNVs (FUnC), developed from machine learning trained on lymphocyte V(D)J recombination. TbA–FUnC offered superior genomic coverage and removed >90% of false-positive CNV calls, allowing extensive examination of submegabase CNVs from over 500 cells throughout the neurogenic period of cerebral cortical development in Mus musculus. Thousands of previously undocumented CNVs were identified. Half were less than 1 Mb in size, with deletions 4× more common than amplification events, and were randomly distributed throughout the genome. However, CNV prevalence during embryonic cortical development was nonrandom, peaking at midneurogenesis with levels triple those found at younger ages before falling to intermediate quantities. These data identify pervasive small and large CNVs as early contributors to neural genomic mosaicism, producing genomically diverse cellular building blocks that form the highly organized, mature brain.

Molecular characterization and expression analysis of three subclasses of IgT in rainbow trout (Oncorhynchus mykiss)

As the teleost specific immunoglobulin, IgT plays important roles in systemic and mucosal immunity. In the current study, in rainbow trout, we have cloned the heavy chain (Igτ) genes of a secretory form of IgT2 as well as the membrane and secretory forms of a third IgT subclass, termed IgT3. Conserved cysteine and tryptophan residues that are crucial for the folding of the immunoglobulin domain as well as hydrophobic and hydrophilic residues within CART motif were identified in all IgT subclasses. Through analysis of the rainbow trout genome assembly, Igτ3 gene was found localized upstream of Igτ1 gene, while Igτ2 gene situated on another scaffold. At the transcriptional level, Igτ1 was mainly expressed in both systemic and mucosal lymphoid tissues, while Igτ2 was largely expressed in systemic lymphoid organs. After LPS and poly (I:C) treatment, Igτ1 and Igτ2 genes exhibited different expression profiles. Interestingly the transcriptional level of Igτ3 was negligible, although its protein product could be identified in trout serum. Importantly, a previously reported monoclonal antibody directed against trout IgT1 was able to recognize IgT2 and IgT3. These data demonstrate that there exist three subclasses of IgT in rainbow trout, and that their heavy chain genes display different expression patterns during stimulation. Overall, our data reflect the diversity and complexity of immunoglobulin in trout, thus provide a better understanding of the IgT system in the immune response of teleost fish.

Involvement of cancer-derived IgG in the proliferation, migration and invasion of bladder cancer cells

It is widely accepted that immunoglobulin (Ig), the classical immune molecule, is extensively expressed in many cell types other than B-cells (non-B-IgG), including some malignant cells. The expression of Ig in malignant cells has been associated with a poor prognosis. In the present study, immunohistochemical analysis detected strong positive staining of IgG in three bladder cancer cell lines, the cancer cells in 77 bladder cancer patient samples and the cells in 3 cystitis glandularis tissue samples, while negative staining was observed in 4 specimens of normal transitional epithelial tissues. Importantly, functional transcripts of IgG with unique VHDJH rearrangement patterns were also found in bladder cancer cells. The knockdown of IgG in bladder cancer cell lines using small interfering RNA significantly inhibited the proliferation, migration and invasion of the cells. Notably, high IgG expression, as determined by immunostaining, was significantly correlated with a high histological grade and recurrence. The results of the present study suggested that IgG expression is involved in the malignant biological behavior and poor prognosis of bladder cancer. Therefore, IgG may serve as a novel target for bladder cancer therapy.

Significant Differences in Physicochemical Properties of Human Immunoglobulin Kappa and Lambda CDR3 Regions

Antibody variable regions are composed of a heavy and a light chain, and in humans, there are two light chain isotypes: kappa and lambda. Despite their importance in receptor editing, the light chain is often overlooked in the antibody literature, with the focus being on the heavy chain complementarity-determining region (CDR)-H3 region. In this paper, we set out to investigate the physicochemical and structural differences between human kappa and lambda light chain CDR regions. We constructed a dataset containing over 29,000 light chain variable region sequences from IgM-transcribing, newly formed B cells isolated from human bone marrow and peripheral blood. We also used a published human naïve dataset to investigate the CDR-H3 properties of heavy chains paired with kappa and lambda light chains and probed the Protein Data Bank to investigate the structural differences between kappa and lambda antibody CDR regions. We found that kappa and lambda light chains have very different CDR physicochemical and structural properties, whereas the heavy chains with which they are paired do not differ significantly. We also observed that the mean CDR3 N nucleotide addition in the kappa, lambda, and heavy chain gene rearrangements are correlated within donors but can differ between donors. This indicates that terminal deoxynucleotidyl transferase may work with differing efficiencies between different people but the same efficiency in the different classes of immunoglobulin chain within one person. We have observed large differences in the physicochemical and structural properties of kappa and lambda light chain CDR regions. This may reflect different roles in the humoral immune response.

Unique Features of Fish Immune Repertoires: Particularities of Adaptive Immunity Within the Largest Group of Vertebrates

Fishes (i.e., teleost fishes) are the largest group of vertebrates. Although their immune system is based on the fundamental receptors, pathways, and cell types found in all groups of vertebrates, fishes show a diversity of particular features that challenge some classical concepts of immunology. In this chapter, we discuss the particularities of fish immune repertoires from a comparative perspective. We examine how allelic exclusion can be achieved when multiple Ig loci are present, how isotypic diversity and functional specificity impact clonal complexity, how loss of the MHC class II molecules affects the cooperation between T and B cells, and how deep sequencing technologies bring new insights about somatic hypermutation in the absence of germinal centers. The unique coexistence of two distinct B-cell lineages respectively specialized in systemic and mucosal responses is also discussed. Finally, we try to show that the diverse adaptations of immune repertoires in teleosts can help in understanding how somatic adaptive mechanisms of immunity evolved in parallel in different lineages across vertebrates.

Antibody-specific model of amino acid substitution for immunological inferences from alignments of antibody sequences

Antibodies are glycoproteins produced by the immune system as a dynamically adaptive line of defense against invading pathogens. Very elegant and specific mutational mechanisms allow B lymphocytes to produce a large and diversified repertoire of antibodies, which is modified and enhanced throughout all adulthood. One of these mechanisms is somatic hypermutation, which stochastically mutates nucleotides in the antibody genes, forming new sequences with different properties and, eventually, higher affinity and selectivity to the pathogenic target. As somatic hypermutation involves fast mutation of antibody sequences, this process can be described using a Markov substitution model of molecular evolution. Here, using large sets of antibody sequences from mice and humans, we infer an empirical amino acid substitution model AB, which is specific to antibody sequences. Compared with existing general amino acid models, we show that the AB model provides significantly better description for the somatic evolution of mice and human antibody sequences, as demonstrated on large next generation sequencing (NGS) antibody data. General amino acid models are reflective of conservation at the protein level due to functional constraints, with most frequent amino acids exchanges taking place between residues with the same or similar physicochemical properties. In contrast, within the variable part of antibody sequences we observed an elevated frequency of exchanges between amino acids with distinct physicochemical properties. This is indicative of a sui generis mutational mechanism, specific to antibody somatic hypermutation. We illustrate this property of antibody sequences by a comparative analysis of the network modularity implied by the AB model and general amino acid substitution models. We recommend using the new model for computational studies of antibody sequence maturation, including inference of alignments and phylogenetic trees describing antibody somatic hypermutation in large NGS data sets. The AB model is implemented in the open-source software CodonPhyML (http://sourceforge.net/projects/codonphyml) and can be downloaded and supplied by the user to ProGraphMSA (http://sourceforge.net/projects/prographmsa) or other alignment and phylogeny reconstruction programs that allow for user-defined substitution models.

The immunoglobulins of cold-blooded vertebrates

Although lymphocyte-like cells secreting somatically-recombining receptors have been identified in the jawless fishes (hagfish and lamprey), the cartilaginous fishes (sharks, skates, rays and chimaera) are the most phylogenetically distant group relative to mammals in which bona fide immunoglobulins (Igs) have been found. Studies of the antibodies and humoral immune responses of cartilaginous fishes and other cold-blooded vertebrates (bony fishes, amphibians and reptiles) are not only revealing information about the emergence and roles of the different Ig heavy and light chain isotypes, but also the evolution of specialised adaptive features such as isotype switching, somatic hypermutation and affinity maturation. It is becoming increasingly apparent that while the adaptive immune response in these vertebrate lineages arose a long time ago, it is most definitely not primitive and has evolved to become complex and sophisticated. This review will summarise what is currently known about the immunoglobulins of cold-blooded vertebrates and highlight the differences, and commonalities, between these and more “conventional” mammalian species.

Genome complexity in the coelacanth is reflected in its adaptive immune system

We have analyzed the available genome and transcriptome resources from the coelacanth in order to characterize genes involved in adaptive immunity. Two highly distinctive IgW-encoding loci have been identified that exhibit a unique genomic organization, including a multiplicity of tandemly repeated constant region exons. The overall organization of the IgW loci precludes typical heavy chain class switching. A locus encoding IgM could not be identified either computationally or by using several different experimental strategies. Four distinct sets of genes encoding Ig light chains were identified. This includes a variant sigma-type Ig light chain previously identified only in cartilaginous fishes and which is now provisionally denoted sigma-2. Genes encoding α/β and γ/δ T-cell receptors, and CD3, CD4, and CD8 co-receptors also were characterized. Ig heavy chain variable region genes and TCR components are interspersed within the TCR α/δ locus; this organization previously was reported only in tetrapods and raises questions regarding evolution and functional cooption of genes encoding variable regions. The composition, organization and syntenic conservation of the major histocompatibility complex locus have been characterized. We also identified large numbers of genes encoding cytokines and their receptors, and other genes associated with adaptive immunity. In terms of sequence identity and organization, the adaptive immune genes of the coelacanth more closely resemble orthologous genes in tetrapods than those in teleost fishes, consistent with current phylogenomic interpretations. Overall, the work reported described herein highlights the complexity inherent in the coelacanth genome and provides a rich catalog of immune genes for future investigations.