Lin Sun, Chunhong Liang, Shidi Qin, Ying Zhu, Ke He. 2024: Exploring the interplay of T cell receptor-V gene copy numbers and major histocompatibility complex selection pressure in avian species: Insights into immune system evolution and reproductive investment. Avian Research, 15(1): 100204. DOI: 10.1016/j.avrs.2024.100204
Citation: Lin Sun, Chunhong Liang, Shidi Qin, Ying Zhu, Ke He. 2024: Exploring the interplay of T cell receptor-V gene copy numbers and major histocompatibility complex selection pressure in avian species: Insights into immune system evolution and reproductive investment. Avian Research, 15(1): 100204. DOI: 10.1016/j.avrs.2024.100204

Exploring the interplay of T cell receptor-V gene copy numbers and major histocompatibility complex selection pressure in avian species: Insights into immune system evolution and reproductive investment

Funds: 

the “Pioneer” and “Leading Goose” R&D Program of Zhejiang 2022C04014

Zhejiang Science and Technology Major Program on Agricultural New Variety Breeding 2021C02068-10

More Information
  • Corresponding author:

    E-mail address: yzhu@swun.edu.cn (Y. Zhu)

    E-mail address: heke@zafu.edu.cn (K. He)

  • Received Date: 27 Jan 2024
  • Rev Recd Date: 09 Aug 2024
  • Accepted Date: 12 Aug 2024
  • Publish Date: 21 Aug 2024
  • Birds, a fascinating and diverse group occupying various habitats worldwide, exhibit a wide range of life-history traits, reproductive methods, and migratory behaviors, all of which influence their immune systems. The association between major histocompatibility complex (MHC) genes and certain ecological factors in response to pathogen selection has been extensively studied; however, the role of the co-working molecule T cell receptor (TCR) remains poorly understood. This study aimed to analyze the copy numbers of TCR-V genes, the selection pressure (ω value) on MHC genes using available genomic data, and their potential ecological correlates across 93 species from 13 orders. The study was conducted using the publicly available genome data of birds. Our findings suggested that phylogeny influences the variability in TCR-V gene copy numbers and MHC selection pressure. The phylogenetic generalized least squares regression model revealed that TCR-Vαδ copy number and MHC-I selection pressure were positively associated with body mass. Clutch size was correlated with MHC selection pressure, and Migration was correlated with TCR-Vβ copy number. Further analyses revealed that the TCR-Vβ copy number was positively correlated with MHC-ⅡB selection pressure, while the TCR-Vγ copy number was negatively correlated with MHC-I peptide-binding region selection pressure. Our findings suggest that TCR-V diversity is significant in adaptive evolution and is related to species’ life-history strategies and immunological defenses and provide valuable insights into the mechanisms underlying TCR-V gene duplication and MHC selection in avian species.

  • Endosymbionts are commonly associated with insects and have interactions ranging from parasitic to mutualistic. For instance, the symbiont genus Wolbachia infects a wide range of insects and is able to manipulate the reproductive characteristics of their host by causing cytoplasmic incompatibility, parthenogenesis, male-killing and feminization (Werren et al., 2008; Saridaki and Bourtzis, 2010). Other endosymbiotic bacteria such as Buchnera aphidicola found in aphids are responsible for the synthesis of essential amino acids that are lacking in the aphids' diet of plant sap (Wilson et al., 2010). Therefore, bacterial acquisition can be critically important to an insect's fitness as it can facilitate expansion to new ecological niches (Douglas, 2009). The genome of an endosymbiont is often smaller than their free-living counterpart as it relies on its host for many functions (Wernegreen, 2015). This can result in mutations accumulating in the symbionts genome due to, e.g., a loss of a DNA repair system, genetic drift, and lack of recombination (McCutcheon and Moran, 2012), which in turn can lead to genome degeneration in the endosymbiont, ultimately resulting in symbiont extinction and replacement.

    Some of the most prominent bacterial lineages that have become established within insects belong to the class Gammaproteobacteria (Husník et al., 2011; Santos-Garcia et al., 2017), including the genus Sodalis (Gammaproteobacteria: Enterobacterales: Pectobacteriaceae). Sodalis glossinidius was first described from the Tsetse Fly (Glossina morsitans Westwood, 1851) (Dale and Maudlin, 1999), and subsequently Sodalis or Sodalis-allied species have been found in several insect groups including aphids (Aphididae), weevils (Curculionoidea), louse flies (Hippoboscidae), and sweat bees (Halictidae) (Nováková and Hypša, 2007; Wilson et al., 2010; Toju et al., 2013; Rubin et al., 2018). Recently chewing lice (Insecta: Phthiraptera) have also been added to this list (Fukatsu et al., 2007).

    Chewing lice are permeant ectoparasites of birds and mammals that have no free-living stage; therefore, they complete their whole life cycle upon their host (Marshall, 1981). Since lice are wingless, transmission between hosts generally only happens when these come into direct contact with each other. This is usually in the form of vertical transmission between parents and nestlings (Clayton and Tompkins, 1994; Lee and Clayton, 1995; Brooke, 2010), or horizontal transmission during mating (Hillgarth, 1996). However, exceptions are known, such as phoretically hitching rides between hosts on hippoboscid flies. This phoretic behavior is one way that lice may be able to come into contact with novel bird host species and could facilitate host switching (Lee et al., 2022).

    Louse endosymbionts have been known since the 1930s (Ries, 1931), but only recently have DNA sequencing and phylogenetic analysis been used to analyze these bacteria. These methods revealed that the symbiont of the Slender Pigeon Louse (Columbicola columbae Linnaeus, 1758) is allied to Sodalis glossinidius from the Tsetse Fly (Fukatsu et al., 2007). Subsequently, Sodalis-allied symbionts have been found in multiple species of Columbicola, Guimaraesiella, and Brueelia (Smith et al., 2013; D’Alessio, 2023; Sweet et al., 2023; Grossi et al., 2024), as well as some bacterial lineages that fall outside Sodalis (Smith et al., 2013). Each species of chewing lice is believed to harbor a single endosymbiotic bacterial lineage that has been maternally transmitted (Smith et al., 2013). This lineage is housed in specialized cells called bacteriocytes clustered on both sides of the abdomen in males and immature females; however, in mature females the cells were only found in ovarian tissue (Smith et al., 2013).

    Ischnoceran chewing lice primarily consume the feather barbs, which are composed of keratin and lack amino acids essential to the survival of the lice. This lack of amino acids caused Alickovic et al. (2021) to hypothesize that the endosymbionts of lice synthesize these amino acids for their hosts. Alickovic et al. (2021) examined the genome of the endosymbiont found in Columbicola wolffhuegeli (Eichler, 1952) parasitizing the Pied Imperial Pigeon (Ducula bicolor) but failed to identify many of the genes that would be expected to support metabolism of amino acids. The role of these endosymbiotic bacteria is thus still unclear.

    The phylogenies of the Sodalis-allied bacteria from Columbicola and Brueelia both have comb-like topography with short internal nodes and long terminal branches (Smith et al., 2013; D’Alessio, 2023). This may be caused by repeated recent symbiont acquisition/replacement events from a common bacterial source, and there appears to be little signal of co-speciation in these louse-symbiont associations. However, in the songbird louse genus Guimaraesiella Eichler, 1949, there is a strong signal of co-speciation between the lice and their symbionts (Grossi et al., 2024). It thus appears that some Sodalis endosymbionts have cospeciated with their louse hosts, but the factors determining whether or not co-speciation happens are unclear.

    Here we examine a fourth louse-symbiont system, that of chewing lice parasitizing shorebirds and terns (Charadriiformes). In general, each shorebird/tern species is parasitized by 1–3 species of lice in the Quadraceps-complex (in this case Lunaceps, Quadraceps, and Saemundssonia), with some species also being parasitized by the distantly related genus Carduiceps Clay & Meinertzhagen, 1939. The Quadraceps-complex includes several distinct lineages that have adapted to different microhabitats of the hosts, such as the wings and the head (Johnson et al., 2012). As multiple species of ischnoceran lice are typically easy to collect from each bird species, we ask the question: do ischnoceran lice occurring on the same host bird species share the same endosymbiotic bacteria? Additionally, does louse relatedness impact symbiont linages? That is, do closely related lice on the same or closely related bird species (e.g., two members of the Quadraceps-complex) have more similar symbiont lineages than more distally related lice parasitizing the same host species (e.g., Carduiceps and Quadraceps-complex).

    Lice were collected from birds caught on the coast of China at three locations (Dongli: 20.84° N, 110.37° E; Hebei: 20.91° N, 110.16° E; Tujiao: 20.89° N, 110.16° E) on the Leizhou Peninsula in western Guangdong and on Chiyu Island (23.32° N, 117.12° E) near Shantou in eastern Guangdong in 2020, using mist nets. Birds were caught under permit 2019-48 from the Leizhou Bureau of Natural Resources. Lice were removed by fumigation following the protocol outlined by Gustafsson et al. (2019). Birds were identified using MacKinnon and Phillipps (2000) and Arlott (2017) and taxonomy follows Clements et al. (2022). Lice were stored in absolute ethanol at −80 ℃, until DNA extraction.

    DNA extraction followed the protocol outlined by Grossi et al. (2024). For polymerase chain reactions (PCR), we used either Cytiva PureTaq Ready-To-Go beads (GE Healthcare, Vienna, Austria) or Qiagen Hot StarTaq Master Mix Kit (Qiagen, Shanghai, China), following the manufacture's protocol for 25 μL reactions. For bacteria we targeted 16S ribosomal RNA gene (16S) using primers: SodF, and R1060 (Nováková and Hypša, 2007), and thermocycle: 95 ℃ for 5 min, 35 cycles (95 ℃ for 30 s, 56 ℃ for 1 min, 72 ℃ for 1 min), 72 ℃ for 10 min. For lice we targeted cytochrome oxidase subunit 1 (COI) using primers: L6625, and H7005 (Hafner et al., 1994), and elongation factor 1α (EF-1α) using primers: EF1-For3, and EF1-CH10 (Danforth and Ji, 1998), and thermocycle: 94 ℃ for 2 min, 35 cycles (94 ℃ for 30 s, 46 ℃ for 30 s, 72 ℃ for 30 s), 72 ℃ for 7 min. PCR products were screened using gel electrophoresis and those with satisfactory bands were sent for sequencing at Tianyi Huiyuan Gene Technology, Co. Ltd. (Guangzhou, China).

    Sequences were assembled using the de novo assemble tool in Geneious Prime® v.2022.1.1, and the ends were trimmed. Sequences were aligned using MUSCLE (Edgar, 2004) with the default settings and checked manually. Three phylogenetic trees were assembled: a louse tree used to confirm identity, a bacteria tree from the lice newly sequenced in this study and a summary tree composed of bacteria sequenced in this study as well as those from other louse genera, insects and bacteria strains from GenBank (https://www.ncbi.nlm.nih.gov/genbank), for accession numbers of sequences pulled from GenBank see Appendix Table S1. Only unique louse/bacteria combination from the newly sequenced bacteria tree were carried over to the summery tree. To determine how many operational taxonomic units (OTUs) were present in the bacteria tree we used a General Mixed Yule Coalescent (Fujisawa and Barraclough, 2013) in RStudio 2023.03.1 (R Core Team, 2018) the “gmyc” function in the splits package was used to determine the number of OTUs (Joseph and Vakayil, 2022). Substitution models for each gene were evaluated in MEGA11 (Kumar et al., 2021), with the best model for lice being, CO1: HKY + G + I, EF-1α: HKY + G, and for both bacteria trees: HKY + G + I.

    The two louse genes were concatenated, and this data set was used for all downstream analyses. As an outgroup for lice, we used Campanulotes compar (Burmeister, 1838) (accession no. CO1: AF384997, and EF-1α: HQ332855) and for bacteria Vibrio cholerae (accession no. LC487865). Bayesian analysis was conducted using BEAST v.2.7.3 (Suchard et al., 2018) with the following settings: the appropriate model(s) for each gene or gene partition were selected (trees were linked for the louse data set), Gamma category count was set to four, Yule Model was used, Markov chain Monte Carlo (MCMC) was run for 1 × 108 generations, and sampled every 1000 trees. The log file produced by BEAST was examined in Tracer (Rambaut et al., 2018) to assess MCMC for convergence. TreeAnnotator v2.7.3 was then used for tree integration, we discarded the first 10,000 trees (10%) as burn-in and constructed a maximum clade credibility tree.

    After DNA extraction the gut was dissected out, and the louse exoskeleton was placed into clove oil for 30 min up to 24 h before it was slide mounted in Canada balsam and allowed to cure at room temperature for 30 days. All slides were deposited at the Institute of Zoology, Guangdong Academy of Sciences, Guangzhou, China. Specimens were examined through a Nikon Eclipse Ni microscope (Nikon Corporation, Tokyo, Japan) and identified using Hopkins and Timmermann (1954), Timmermann (1954a, 1952, 1950), Gustafsson and Olsson (2012), as well as through direct comparison with previously identified specimens in cases where no adequate illustrations or descriptions have been published.

    In total 85 individual birds, representing 11 species were examined for lice, of which 68 (80%) were infested by at least one species of louse. See Table 1 for a breakdown by bird and louse species.

    Table  1.  Bird species examined in this study, and the number of those infested with chewing lice.
    Bird species No. of birds examined Number of birds infested with lice
    Carduiceps Lunaceps Quadraceps Saemundssonia
    Calidris alpina 25 9 19
    Calidris canutus 1 1
    Calidris pygmaea 2 2
    Charadrius alexandrinus 14 2 10
    Charadrius leschenaultii 2 1
    Charadrius mongolus 3 3
    Limosa lapponica 1 1
    Onychoprion anaethetus 13 6
    Thalasseus bergii 10 6 1
    Tringa stagnatilis 10 1 1 9
    Tringa totanus 5 4
     | Show Table
    DownLoad: CSV

    All lice were found to belong to one of four genera: Carduiceps or Lunaceps (Clay & Meinertzhagen, 1939), on sandpipers (Calidrinae), and Quadraceps Clay & Meinertzhagen, 1939, or Saemundssonia Timmermann, 1936, on the other hosts. In two cases (J4451 and J4456), Lunaceps specimens were found on plovers, presumably as a result of straggling. In addition, samples from two different specimens of Spoon-billed Sandpiper (Calidris pygmaea) were found to represent different species of Lunaceps, with one (Lunaceps schismatus Gustafsson & Olsson, 2012) presumably being a straggler. The second species of Lunaceps is morphologically close to Lunaceps falcinellus Timmermann, 1954, but genetically distinct, and may represent an undescribed species. Apart from the specimens from terns, these specimens were previously reported by Grossi et al. (2023).

    The louse phylogeny (Fig. 1) was constructed from CO1, 377 bp (203 variable sites, 192 parsimony-informative sites) and EF-1α, 348 bp (99 variable sites, 74 parsimony-informative sites). In this phylogeny, the genera Lunaceps and Saemundssonia are both recovered as monophyletic with good support, but both are nested inside a paraphyletic Quadraceps; Carduiceps is placed outside the Quadraceps-complex. Within Quadraceps, species grouped together in accordance with louse host associations.

    Figure  1.  Bayesian phylogenetic tree of chewing lice infesting shorebirds and terns. Bolded taxa indicate sequences that are included in the summary phylogeny (Fig. 3), colored voucher numbers indicate lice that came from the same host individual. Only posterior probabilities greater than 0.90 were included.

    We found endosymbiotic bacteria in each of the four chewing lice genera (Table 2). The phylogeny of the endosymbiotic bacteria from the lice infesting shore birds was constructed from 16S, 890 bp (242 variable sites, 167 parsimony-informative sites) and shows two distinct clades referred to as A and B (Fig. 2). The OTU analysis found 13 distinct units, see Fig. 2.

    Table  2.  Collection and GenBank accession numbers for specimens included in this study.
    Louse species Bird host species Voucher No. Loc. Symbiont (16S) Accession No. Lice Accession No.
    CO1 Ef-1α
    Carduiceps meinertzhageni Calidris alpina J4366-02 L3 PP081401 PP082907 PP092242
    Carduiceps meinertzhageni Calidris alpina J4368-02 L3 PP081403 PP082909 PP092244
    Carduiceps meinertzhageni Calidris alpina J4458-01 L3 PP081411 PP082917 PP092252
    Carduiceps meinertzhageni Calidris alpina J4471-01 L1 PP081418 PP082924 PP092259
    Lunaceps sp. Calidris pygmaea J4442-01 L3 PP081406 PP082912 PP092247
    Lunaceps drosti Calidris canutus J4490-01 L1 PP081423 PP082929 PP092264
    Lunaceps limosella Limosa lapponica J4509-06 L1 PP081425 PP082931 PP092266
    Lunaceps schismatus Calidris pygmaea J4437-01 L3 PP081405 PP082911 PP092246
    Lunaceps schismatus Charadrius alexandrinus J4451-01 L3 PP081407 PP082913 PP092248
    Lunaceps schismatus Charadrius alexandrinus J4456-01 L3 PP081409 PP082915 PP092250
    Lunaceps schismatus Calidris alpina J4364-01 L3 PP081398 PP082904 PP092239
    Lunaceps schismatus Calidris alpina J4366-01 L3 PP081400 PP082906 PP092241
    Lunaceps schismatus Calidris alpina J4367-01 L3 PP081402 PP082908 PP092243
    Lunaceps schismatus Calidris alpina J4457-01 L3 PP081410 PP082916 PP092251
    Lunaceps schismatus Calidris alpina J4460-01 L3 PP081412 PP082918 PP092253
    Lunaceps schismatus Calidris alpina J4462-01 L3 PP081413 PP082919 PP092254
    Lunaceps schismatus Calidris alpina J4467-01 L1 PP081414 PP082920 PP092255
    Lunaceps schismatus Calidris alpina J4469-01 L3 PP081416 PP082922 PP092257
    Lunaceps schismatus Calidris alpina J4470-01 L1 PP081417 PP082923 PP092258
    Lunaceps schismatus Calidris alpina J4472-01 L3 PP081419 PP082925 PP092260
    Lunaceps schismatus Calidris alpina J4455-01 L3 PP081408 PP082914 PP092249
    Quadraceps legatus Onychoprion anaethetus J4578-01 CI PP081431 PP082937 PP092272
    Quadraceps legatus Onychoprion anaethetus J4577-04 CI PP081430 PP082936 PP092271
    Quadraceps legatus Onychoprion anaethetus J4588-02 CI PP081434 PP082940 PP092275
    Quadraceps macrocephalus Charadrius alexandrinus J4468-01 L1 PP081415 PP082921 PP092256
    Quadraceps obscurus Tringa stagnatilis J4483-01 L1 PP081421 PP082927 PP092262
    Quadraceps obtusus Tringa totanus J4489-02 L1 PP081422 PP082928 PP092263
    Quadraceps obtusus Tringa totanus J4492-02 L3 PP081424 PP082930 PP092265
    Quadraceps ptyadis Charadrius leschenaultii J4365-01 L3 PP081399 PP082905 PP092240
    Quadraceps ptyadis Charadrius mongolus J4369-01 L2 PP081404 PP082910 PP092245
    Quadraceps ptyadis Charadrius mongolus J4479-01 L1 PP081420 PP082926 PP092261
    Quadraceps sellatus Thalasseus bergii J4585-02 CI PP081433 PP082939 PP092274
    Saemundssonia laticaudata Thalasseus bergii J4583-02 CI PP081432 PP082938 PP092273
    Saemundssonia meridiana Onychoprion anaethetus J4576-01 CI PP081426 PP082932 PP092267
    Saemundssonia meridiana Onychoprion anaethetus J4576-02 CI PP081427 PP082933 PP092268
    Saemundssonia meridiana Onychoprion anaethetus J4577-01 CI PP081428 PP082934 PP092269
    Saemundssonia meridiana Onychoprion anaethetus J4577-02 CI PP081429 PP082935 PP092270
    Voucher number, the first four digits refer to the host individual and the last two digits refer to a unique louse. Thus, J4366–01 and J4366-02 are different lice from the same host individual. Location abbreviation: CI – Chiyu Island, Shantou; L1 – Dongli, Leizhou; L2 – Hebei, Leizhou; L3 – Tujiao, Leizhou, Guangdong, China.
     | Show Table
    DownLoad: CSV
    Figure  2.  Bayesian phylogenetic tree of endosymbiotic bacteria from chewing lice infesting shorebirds and terns. Bolded taxa indicate sequences that are included in the summary phylogeny (Fig. 3), colored voucher numbers indicate lice that came from the same host individual, symbols denote bird species that were infested with chewing lice with endosymbionts found in both Clades A and B. Only posterior probabilities greater than 0.90 were included. Black lines denote clades referenced to in text and grey lines indicate operational taxonomic units.

    Each louse species in this study has endosymbiotic bacteria that falls into either clade A or B, and no louse species has endosymbiotic bacteria from both clades. In two cases (Calidris alpina J4366 and Onychoprion anaethetus J4577), lice of different genera found on the same bird host individual tested positive for bacteria from the same clade. However, in the case of the lice from Thalasseus bergii (from two different host individuals J4583, J4585), the two louse species obtained tested positive for symbionts belonging to different clades. The two louse species obtained from Charadrius alexandrinus also tested positive for different symbiont clades, but in this case the Lunaceps samples are likely stragglers. Notably, each bacterial OTU was associated with exactly one louse species, but symbionts from closely related louse species did not always group together.

    Prevalence of infection for clade A: Carduiceps 80% (4/5), Lunaceps 50% (15/30), Quadraceps 11% (4/36) and Saemundssonia 63% (5/8). Prevalence of infection for clade B: Carduiceps 0%, Lunaceps 7% (2/30), Quadraceps 19% (7/36) and Saemundssonia 0%.

    One symbiont sequence from each OTU was selected for inclusion in a larger analysis of related bacteria from other louse groups as well as other insects (Fig. 3), constructed from 16S, 848 bp (316 variable sites, 220 parsimony-informative sites). Most deeper nodes of this tree are unsupported, but the division between the same two clades as in Fig. 2 is evident. Symbionts from shorebird lice are scattered throughout both clades, and few shorebird louse symbionts appear to be closely related. Notably, Sodalis symbionts from hippoboscid flies were not placed close to symbionts from Guimaraesiella spp., known to be phoretic on these flies. Clade D in this tree also include bacteria belonging to genera other than Sodalis.

    Figure  3.  Bayesian phylogenetic tree of endosymbiotic bacteria from chewing lice infesting shorebirds and terns, other genera of chewing lice, insects, and bacteria strains from GenBank (see Appendix Table S1 for accession numbers). Bolded taxa indicate newly sequenced endosymbiotic bacteria from this study. Only posterior probabilities greater than 0.90 were included. Black lines denote clades referenced to in text.

    The ischnoceran chewing lice of most shorebirds belong to the Quadraceps-complex, a morphologically diverse group of lice that includes several microhabitat specialists. Notably, whereas e.g., the head and body lice of many songbirds are distantly related (e.g., de Moya et al., 2019), those of many shorebirds are thus closely related. One notable exception to this is the small group of birds that are parasitized by lice in the genus Carduiceps, which are distantly related to the Quadraceps-complex, and possibly more closely related to lice from other host groups (Gustafsson and Olsson, 2017). This variety of relationships between lice within a relatively small group of bird hosts makes shorebird and tern lice ideal for examining the distribution patterns of their bacterial symbionts.

    The louse-bacteria systems studied to date show two almost diametrically opposite patterns: in Columbicola there appears to be little co-speciation and wide-spread replacement of bacterial strains, whereas in Guimaraesiella the symbionts appear to co-speciate extensively with their hosts. In shorebird and tern lice, we find that bacteria from each of the four louse genera examined are mixed together, with little resemblance between the louse phylogeny (Fig. 1) and that of their bacteria (Fig. 2). Put into a larger context, the bacteria of shorebird lice are scattered throughout the tree (Fig. 3), and thus potentially more reminiscent of those of Columbicola spp.

    The symbionts recovered from shorebird lice fell into two different clades. Sodalis-allied symbionts were present in all four genera examined (Carduiceps, Lunaceps, Quadraceps and Saemundssonia), whereas endosymbionts belonging to the family Enterobacteriaceae were only found in Lunaceps and Quadraceps. Potentially, this is due to the larger number of Lunaceps and Quadraceps examined, but more data is needed to evaluate whether both these clades are more widely distributed across shorebird lice. Judging from the summary tree (Fig. 3), both clades seem to be established in many different louse groups. Both of these endosymbiont groups have previously been recorded from Columbicola (Smith et al., 2013); however, only Sodalis-allied symbionts have been recorded from Brueelia and Guimaraesiella (D’Alessio, 2023; Sweet et al., 2023; Grossi et al., 2024).

    Each louse species was associated with a single bacterial OTU, even in cases where the same louse species occurs on multiple bird host species (e.g., J4365 and J4479 in Fig. 2). Thus, it does not appear that the bird hosts are the source of the endosymbiotic bacteria. However, when comparing different louse species associated with the same host species, different patterns appear. First, in the case of lice from Dunlins (Calidris alpina), each louse species was associated with a single bacterial OTU, but these were relatively closely related. Second, in the case of lice from Bridled Terns (Onychoprion anaethetus), each louse species was associated with a different bacterial symbiont lineage, and while they both were placed in Clade A in Fig. 2, they are not closely related within this clade. Finally, the two louse species from the Greater Crested Tern (Thalasseus bergii) were even more distantly related, with one symbiont lineage in each of Clades A and B (Fig. 2). Thus, host associations of the lice do not seem to predict relationships of their bacterial symbionts.

    Notably, when seen in isolation, the shorebird and tern louse symbionts appear to show a mixture of the two speciation types seen in other examined louse symbionts. The evolutionary history of the symbionts of Columbicola is characterized by short terminal branches and longer internal branches. Smith et al. (2013) ran multiple simulations and found that this typology was formed by repeated endosymbiont acquisition. Repeated replacement of Sodalis symbionts is not unique to chewing lice, but is found in many insects such as planthoppers, mealybugs, weevils and louse flies (Toju et al., 2013; Husnik and McCutcheon, 2016; Šochová et al., 2017; Michalik et al., 2021). This pattern is also seen in Clade B here, and in the two symbiont OTUs from Saemundssonia spp. In contrast, much of the rest of Clade A has shorter internal branches, more reminiscent of those found in the symbionts of Guimaraesiella spp. (Grossi et al., 2024). However, in the summary tree (Fig. 3), these differences seem to disappear, and each of the Clade A OTUs are connected to purportedly closely related species by longer internal branches.

    The bird hosts sampled for this study were all caught at similar latitudes in the same climatic zone and often in the same locality, most of which represent wintering grounds for the hosts. The exception to this is the samples from tern species, which were caught at breeding sites in a different part of Guangdong. In both cases, this leaves parts of the hosts’ ranges unsampled; however, it is unlikely that more expansive sampling would change our conclusion that bird host species does not impact the endosymbiont present in the louse. Above all, the associations between lice and their bacteria appear to be stable over very long time periods, and unlikely to change over the migration or breeding cycle of a single bird. Grossi et al. (2024) examined the endosymbiont of Guimaraesiella mcgrewi infesting four bird host species across five locations in southern China and found all G. mcgrewi harbored the same endosymbiotic strain. While there is evidence for Sodalis in the genus Columbicola to be of recent origin (∼0.4 My from free-living species [Smith et al., 2013]) these are still long-established relationships that do not change seasonally. However, it is possible that other populations of the same bird species breeding and migrating in different parts of the world may be parasitized with lice associated with other bacterial strains. More data from other host flyways are needed to address this.

    Alexandra A. Grossi: Writing – original draft, Methodology, Investigation, Formal analysis. Min Zhang: Writing – review & editing, Resources, Investigation, Funding acquisition. Fasheng Zou: Supervision, Resources, Project administration, Funding acquisition. Daniel R. Gustafsson: Writing – review & editing, Supervision, Resources, Investigation, Funding acquisition, Conceptualization.

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    We thank Chi-Yeung Choi from Duke Kunshan University, and Chunpo Tian from Shaanxi Normal University, as well as all our past colleagues for their assistance in collecting the lice used in this study.

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.avrs.2024.100206.

  • Anonymity, 2023. TCR diversity underpins immunotherapy success. Cancer Discov. 13, 1031–1032.
    Agudo, R., Carrete, M., Alcaide, M., Rico, C., Hiraldo, F., Donázar, J.A., 2012. Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird. Proc. Biol. Sci. 279, 3241–3249.
    Arstila, T.P., Casrouge, A., Baron, V., Even, J., Kanellopoulos, J., Kourilsky, P., 1999. A direct estimate of the human alphabeta T cell receptor diversity. Science 286, 958–961.
    Auld, H., MacIver, D., Klaassen, J., 2004. Heavy rainfall and waterborne disease outbreaks: the Walkerton example. J. Toxicol. Env. Heal. A 67, 1879–1887.
    Biedrzycka, A., Bielański, W., Ćmiel, A., Solarz, W., Zając, T., Migalska, M., et al., 2018. Blood parasites shape extreme major histocompatibility complex diversity in a migratory passerine. Mol. Ecol. 27, 2594–2603.
    Bush, S.E., Clayton, D.H., 2018. Anti-parasite behaviour of birds. Philos. T. Roy. Soc. B 373, 20170196.
    Chen, C., Chen, H., Zhang, Y., Thomas, H.R., Frank, M.H., He, Y., et al., 2020. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202.
    Chen, L., Qiu, Q., Jiang, Y., Wang, K., Lin, Z., Li, Z., et al., 2019. Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits. Science 364, eaav6202.
    Chen, Y., Xiang, L., Chen, P., Zhao, H., 2023. Vultures as a model for testing molecular adaptations of dietary specialization in birds. Avian Res. 14, 100128.
    Cheng, Y., Miller, M.J., Zhang, D., Xiong, Y., Hao, Y., Jia, C., et al., 2021. Parallel genomic responses to historical climate change and high elevation in East Asian songbirds. P. Natl. Acad. Sci. U.S.A. 118, e2023918118.
    Cichoń, M., Dubiec, A., 2005. Cell-mediated immunity predicts the probability of local recruitment in nestling blue tits. J. Evol. Biol. 18, 962–966.
    Davies, C.S., Worsley, S.F., Maher, K.H., Komdeur, J., Burke, T., Dugdale, H.L., et al., 2022. Immunogenetic variation shapes the gut microbiome in a natural vertebrate population. Microbiome 10, 41.
    Davis, M.M., Bjorkman, P.J., 1988. T-cell antigen receptor genes and T-cell recognition. Nature 334, 395–402.
    Dunn, P.O., Bollmer, J.L., Freeman-Gallant, C.R., Whittingham, L.A., 2013. MHC variation is related to a sexually selected ornament, survival, and parasite resistance in common yellowthroats. Evolution 67, 679–687.
    Fleischer, R., Schmid, D.W., Wasimuddin, Brändel, S.D., Rasche, A., Corman, V.M., et al., 2022. Interaction between MHC diversity and constitution, gut microbiota and Astrovirus infections in a neotropical bat. Mol. Ecol. 31, 3342–3359.
    Früh, S.P., Früh, M.A., Kaufer, B.B., Göbel, T.W., 2024. Unraveling the chicken T cell repertoire with enhanced genome annotation. Front. Immunol. 15, 1359169.
    Gurevich, A., Saveliev, V., Vyahhi, N., Tesler, G., 2013. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075.
    Hao, Y., Xiong, Y., Cheng, Y., Song, G., Jia, C., Qu, Y., et al., 2019. Comparative transcriptomics of 3 high-altitude passerine birds and their low-altitude relatives. P. Natl. Acad. Sci. U.S.A. 116, 11851–11856.
    He, K., Minias, P., Dunn, P., 2020. Long-read genome assemblies reveal extraordinary variation in the number and structure of MHC loci in birds. Genome Biol. Evol. 13, evaa270.
    Huang, X., Chen, Z., Yang, G., Xia, C., Luo, Q., Gao, X., et al., 2022. Assemblages of Plasmodium and related parasites in birds with different migration statuses. Int. J. Mol. Sci. 23, 10277.
    Irla, M., 2022. Instructive cues of thymic T cell selection. Ann. Rev. Immunol. 40, 95–119.
    Jarvis, E.D., Mirarab, S., Aberer, A.J., Li, B., Houde, P., Li, C., et al., 2014. Whole-genome analyses resolve early branches in the tree of life of modern birds. Science 346, 1320–1331.
    Knafler, G.J., Clark, J.A., Boersma, P.D., Bouzat, J.L., 2012. MHC diversity and mate choice in the magellanic penguin, Spheniscus magellanicus. J. Hered. 103, 759–768.
    Liang, C., Sun, L., Zhu, Y., Wu, J., Zhao, A., Huang, T., et al., 2024. Local chicken breeds exhibit abundant TCR-V segments but similar repertoire diversity. Dev. Comp. Immunol. 157, 105196.
    Liang, C., Sun, L., Zhu, Y., Zhao, A., Liu, H., He, K., 2023. Macroevolution of avian T cell receptor C segments using genomic data. Immunogenetics 75, 531–541.
    Liu, F., Li, J., Lin, I.Y.C., Yang, X., Ma, J., Chen, Y., et al., 2020. The genome esequencing of TCR loci in Gallus gallus revealed their distinct evolutionary features in avians. Immunohorizons 4, 33–46.
    Luzuriaga-Neira, A.R., Alvarez-Ponce, D., 2022. Rates of protein evolution across the marsupial phylogeny: heterogeneity and link to life-history traits. Genome Biol. Evol. 14, evab277.
    Merelli, I., Guffanti, A., Fabbri, M., Cocito, A., Furia, L., Grazini, U., et al., 2010. RSSsite: a reference database and prediction tool for the identification of cryptic Recombination Signal Sequences in human and murine genomes. Nucleic Acids Res. 38, W262–W267.
    Migalska, M., Sebastian, A., Radwan, J., 2019. Major histocompatibility complex class I diversity limits the repertoire of T cell receptors. P. Natl. Acad. Sci. U.S.A. 116, 5021–5026.
    Migalska, M., Węglarczyk, K., Dudek, K., Homa, J., 2024. Evolutionary trade-offs constraining the MHC gene expansion: beyond simple TCR depletion model. Front. Immunol. 14, 1240723.
    Minias, P., Gutiérrez, J.S., Dunn, P.O., 2020. Avian major histocompatibility complex copy number variation is associated with helminth richness. Biol. Lett. 16, 20200194.
    Minias, P., Pikus, E., Whittingham, L.A., Dunn, P.O., 2019. Evolution of copy number at the MHC varies across the avian tree of life. Genome Biol. Evol. 11, 17–28.
    Minias, P., Whittingham, L.A., Dunn, P.O., 2017. Coloniality and migration are related to selection on MHC genes in birds. Evolution 71, 432–441.
    O'Connor, E.A., Cornwallis, C.K., Hasselquist, D., Nilsson, J., Westerdahl, H., 2018. The evolution of immunity in relation to colonization and migration. Nat. Ecol. Evol. 2, 841–849.
    O'Connor, E.A., Hasselquist, D., Nilsson, J., Westerdahl, H., Cornwallis, C.K., 2020. Wetter climates select for higher immune gene diversity in resident, but not migratory, songbirds. Proc. Biol. Sci. 287, 20192675.
    Okanga, S., Cumming, G.S., Hockey, P.A., 2013. Avian malaria prevalence and mosquito abundance in the Western Cape, South Africa. Malaria J 12, 370.
    Orme, D., Freckleton, R.P., Thomas, G.H., Petzoldt, T., 2013. Caper: comparative analyses of phylogenetics and evolution in R. Methods Ecol. Evol. 3, 145–151.
    Pagel, M., 1997. Inferring evolutionary processes from phylogenies. Zool. Scr. 26, 331–483.
    Parra, Z.E., Miller, R.D., 2012. Comparative analysis of the chicken TCRα/δ locus. Immunogenetics 64, 641–645.
    Podos, J., Webster, M.S., 2022. Ecology and evolution of bird sounds. Curr. Biol. 32, R1100–R1104.
    Pond, S.L., Frost, S.D., Muse, S.V., 2005. HyPhy: hypothesis testing using phylogenies. Bioinformatics 21, 676–679.
    Promerová, M., Králová, T., Bryjová, A., Albrecht, T., Bryja, J., 2013. MHC class ⅡB exon 2 polymorphism in the Grey partridge (Perdix perdix) is shaped by selection, recombination and gene conversion. PLoS One 8, e69135.
    Qu, Y., Chen, C., Chen, X., Hao, Y., She, H., Wang, M., et al., 2021. The evolution of ancestral and species-specific adaptations in snowfinches at the Qinghai-Tibet Plateau. P. Natl. Acad. Sci. U.S.A. 118, e2012398118.
    Robinson, J., Halliwell, J.A., McWilliam, H., Lopez, R., Marsh, S.G., 2013. IPD–the immuno polymorphism database. Nucleic Acids Res. 41, D1234–D1240.
    Roche, B., Lebarbenchon, C., Gauthier-Clerc, M., Chang, C.M., Thomas, F., Renaud, F., et al., 2009. Water-borne transmission drives avian influenza dynamics in wild birds: the case of the 2005–2006 epidemics in the Camargue area. Infect. Genet. Evol. 9, 800–805.
    Rogers, E.J., McGuire, L., Longstaffe, F.J., Clerc, J., Kunkel, E., Fraser, E., 2022. Relating wing morphology and immune function to patterns of partial and differential bat migration using stable isotopes. J. Anim. Ecol. 91, 858–869.
    Ruiz, M.O., Chaves, L.F., Hamer, G.L., Sun, T., Brown, W.M., Walker, E.D., et al., 2010. Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in Northeast Illinois, USA. Parasit. Vectors 3, 19.
    Saino, N., Szép, T., Ambrosini, R., Romano, M., Møller, A.P., 2004. Ecological conditions during winter affect sexual selection and breeding in a migratory bird. Proc. Biol. Sci. 271, 681–686.
    Schwenke, R.A., Lazzaro, B.P., Wolfner, M.F., 2016. Reproduction-immunity trade-offs in insects. Ann. Rev. Entomol. 61, 239–256.
    Sepil, I., Lachish, S., Hinks, A.E., Sheldon, B.C., 2013. Mhc supertypes confer both qualitative and quantitative resistance to avian malaria infections in a wild bird population. Proc. Biol. Sci. 280, 20130134.
    Shah, K., Al-Haidari, A., Sun, J., Kazi, J.U., 2021. T cell receptor (TCR) signaling in health and disease. Signal Transduct. Target. Ther. 6, 412.
    Tieleman, B.I., 2018. Understanding immune function as a pace of life trait requires environmental context. Behav. Ecol. Sociobiol. 72, 55.
    Tobias, J.A., Sheard, C., Pigot, A.L., Devenish, A.J.M., Yang, J., Sayol, F., et al., 2022. AVONET: morphological, ecological and geographical data for all birds. Ecol. Lett. 25, 581–597.
    Vidović, D., Matzinger, P., 1988. Unresponsiveness to a foreign antigen can be caused by self-tolerance. Nature 336, 222–225.
    Vinkler, M., Fiddaman, S.R., Těšický, M., O'Connor, E.A., Savage, A.E., Lenz, T.L., et al., 2023. Understanding the evolution of immune genes in jawed vertebrates. J. Evol. Biol. 36, 847–873.
    Wang, Y., Wang, C., Hong, W., Tian, S., Lu, Q., Wang, B.J., et al., 2023a. Genomic analysis of Tibetan ground tits identifies molecular adaptations associated with cooperative breeding. Curr. Zool. 69, 620–630.
    Wang, Z., Wang, Y.W., Kasuga, T., Hassler, H., Lopez-Giraldez, F., Dong, C., et al., 2023b. Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa. Mol. Ecol. .
    Westerdahl, H., Waldenström, J., Hansson, B., Hasselquist, D., von Schantz, T., Bensch, S., 2005. Associations between malaria and MHC genes in a migratory songbird. Proc. Biol. Sci. 272, 1511–1518.
    Westerdahl, H., Wittzell, H., von Schantz, T., Bensch, S., 2004. MHC class I typing in a songbird with numerous loci and high polymorphism using motif-specific PCR and DGGE. Heredity 92, 534–542.
    Yang, H., Wang, Y., Jia, Z., Wang, Y., Yang, X., Wu, P., et al., 2021. Characteristics of T-cell receptor repertoire and correlation with EGFR mutations in all stages of lung cancer. Front. Oncol. 11, 537735.
    Yang, Z., Sun, Y., Ma, Y., Li, Z., Zhao, Y., Ren, L., et al., 2017. A comprehensive analysis of the germline and expressed TCR repertoire in White Peking duck. Sci. Rep. 7, 41426.
    Zhang, L., Liu, Y., Meng, G., Liang, R., Zhang, B., Xia, C., 2020. Structural and biophysical insights into the TCRαβ complex in chickens. iScience 23, 101828.
    Zou, D., Tian, S., Zhang, T., Zhuoma, N., Wu, G., Wang, M., et al., 2021. Vulture genomes reveal molecular adaptations underlying obligate scavenging and low levels of genetic diversity. Mol. Biol. Evol. 38, 3649–3663.
  • Related Articles

  • Cited by

    Periodical cited type(7)

    1. Shanrui Wu, Ru Jia, Ying Wang, et al. Prevalence, Diversity, and Virulence of Campylobacter Carried by Migratory Birds at Four Major Habitats in China. Pathogens, 2024, 13(3): 230. DOI:10.3390/pathogens13030230
    2. Zoljargal Purevdorj, Namsrai Jargal, Onolragchaa Ganbold, et al. Spatial and Temporal Variations in Waterfowl Assemblage Structures in Mongolian Lakes and the Changes Linked to the Gradient of Lake Surface Areas. Diversity, 2023, 15(3): 334. DOI:10.3390/d15030334
    3. Zhongrong Wu, Xiulin Ye, Zhongfan Kuang, et al. Positive Effects of Land Use Change on Wintering Bar-Headed Geese between 2010 and 2021. Animals, 2022, 12(22): 3142. DOI:10.3390/ani12223142
    4. MARY ANNE BISHOP, DONGPING LIU, GUOGANG ZHANG, et al. Rapid growth of the Bar-headed Goose Anser indicus wintering population in Tibet, China: 1991–2017. Bird Conservation International, 2021. DOI:10.1017/S0959270921000265
    5. Wen Wang, Fang Wang, Rongkai Hao, et al. First de novo whole genome sequencing and assembly of the bar-headed goose. PeerJ, 2020, 8: e8914. DOI:10.7717/peerj.8914
    6. Junjian Zhang, Yanbo Xie, Laixing Li, et al. Assessing site-safeguard effectiveness and habitat preferences of Bar-headed Geese (Anser indicus) at their stopover sites within the Qinghai-Tibet Plateau using GPS/GSM telemetry. Avian Research, 2020, 11(1) DOI:10.1186/s40657-020-00230-9
    7. Ruobing Zheng, Lacy Smith, Diann Prosser, et al. Investigating Home Range, Movement Pattern, and Habitat Selection of Bar-headed Geese during Breeding Season at Qinghai Lake, China. Animals, 2018, 8(10): 182. DOI:10.3390/ani8100182

    Other cited types(0)

Catalog

    Figures(3)  /  Tables(2)

    Article Metrics

    Article views (10) PDF downloads (1) Cited by(7)

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return