Characteristics of MHC class I genes in house sparrows Passer domesticus as revealed by long cDNA transcripts and amplicon sequencing

Evaluation of Genetic Diversity and Parasite-Mediated Selection of MHC Class I Genes in Emberiza godlewskii (Passeriformes: Emberizidae)

The major histocompatibility complex (MHC) is a multi-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next-generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles, clustering into nine supertypes. Compared with other passerine birds, the number of MHC class I alleles per individual in Godlewski’s buntings is high (mean 16.1 ± 3.3, median 16). In addition, we demonstrated signatures of historical and contemporary selection on MHC genes. Reflecting historical selection, ten amino acid sites in the antigen-binding domain showed signatures of balancing selection, eight of which exhibit high amino acid polymorphism. In terms of contemporary selection, we found that specific MHC supertypes were nominally associated with the infection of two malaria parasite lineages. These findings indicate the action of historical and possibly also contemporary balancing selection and suggest negative frequency-dependent or fluctuating selection as possible selection mechanisms.

First Evidence of a Relationship Between Female Major Histocompatibility Complex Diversity and Eggshell Bacteria in House Sparrows (Passer domesticus)

Bacteria are known to exert positive and negative influences on animals’ health and fitness. Bacteria, in particular those inhabiting the skin and inner organs of vertebrates, are horizontally or vertically transmitted. Specifically, mothers of bird species can transfer bacterial strains to their offspring when the egg is passing the reproductive tract, as the eggshell rubs against the wall of the uterus. In this context, the female immune system might play an important role in influencing the vertical transmission of bacteria. Here, we investigate the relationship between the major histocompatibility complex (MHC) and cultivable eggshell bacteria originating putatively from the female urogenital tract in a captive population of house sparrows (Passer domesticus). We predict that females with a more variable MHC will transfer fewer bacteria onto the eggshells. Our results show a negative relationship between the number of functional MHC class I alleles and bacteria originating in the urinary tract and growing on a selective medium. This is the first study to find a correlation between female MHC diversity and eggshell bacteria.

Allelic diversity and patterns of selection at the major histocompatibility complex class I and II loci in a threatened shorebird, the Snowy Plover (Charadrius nivosus)

Background

Understanding the structure and variability of adaptive loci such as the major histocompatibility complex (MHC) genes is a primary research goal for evolutionary and conservation genetics. Typically, classical MHC genes show high polymorphism and are under strong balancing selection, as their products trigger the adaptive immune response in vertebrates. Here, we assess the allelic diversity and patterns of selection for MHC class I and class II loci in a threatened shorebird with highly flexible mating and parental care behaviour, the Snowy Plover (Charadrius nivosus) across its broad geographic range.

Results

We determined the allelic and nucleotide diversity for MHC class I and class II genes using samples of 250 individuals from eight breeding population of Snowy Plovers. We found 40 alleles at MHC class I and six alleles at MHC class II, with individuals carrying two to seven different alleles (mean 3.70) at MHC class I and up to two alleles (mean 1.45) at MHC class II. Diversity was higher in the peptide-binding region, which suggests balancing selection. The MHC class I locus showed stronger signatures of both positive and negative selection than the MHC class II locus. Most alleles were present in more than one population. If present, private alleles generally occurred at very low frequencies in each population, except for the private alleles of MHC class I in one island population (Puerto Rico, lineage tenuirostris).

Conclusion

Snowy Plovers exhibited an intermediate level of diversity at the MHC, similar to that reported in other Charadriiformes. The differences found in the patterns of selection between the class I and II loci are consistent with the hypothesis that different mechanisms shape the sequence evolution of MHC class I and class II genes. The rarity of private alleles across populations is consistent with high natal and breeding dispersal and the low genetic structure previously observed at neutral genetic markers in this species.

A global analysis of selection at the avian MHC

Recent advancements in sequencing technology have resulted in rapid progress in the study of the major histocompatibility complex (MHC) in non-model avian species. Here, we analyze a global dataset of avian MHC class I and class II sequences (ca. 11,000 sequences from over 250 species) to gain insight into the processes that govern macroevolution of MHC genes in birds. Analysis of substitution rates revealed striking differences in the patterns of diversifying selection between passerine and non-passerine birds. Non-passerines showed stronger selection at MHC class II, which is primarily involved in recognition of extracellular pathogens, while passerines showed stronger selection at MHC class I, which is involved in recognition of intracellular pathogens. Positions of positively selected amino-acid residues showed marked discrepancies with peptide-binding residues (PBRs) of human MHC molecules, suggesting that using a human classification of PBRs to assess selection patterns at the avian MHC may be unjustified. Finally, our analysis provided evidence that indel mutations can make a substantial contribution to adaptive variation at the avian MHC.

Characterization of Brachyspira communities from clinical cases of swine mucohaemorrhagic diarrhea through deep sequencing of the NADH oxidase (nox) gene

Swine dysentery is traditionally associated with Brachyspira hyodysenteriae, but the re-emergence of Brachyspira-associated disease in North America associated with a novel causative species, B. hampsonii, is now a concern for swine producers. The pathogenesis of Brachyspira-associated disease is not completely understood, and it is not known whether mixed infections of Brachyspira spp. are important in disease development. Deep sequencing of partial sequences of the nox gene amplified with genus-specific primers was used to detect Brachyspira spp. in 55 fecal samples from clinical cases of mucohaemorrhagic diarrhea in pigs from Western Canada that had been identified as positive for one or more Brachyspira species using established diagnostic tests. Synthetic mixtures of Brachyspira genomic DNA were included in the study to define detection limits for the technique and identify biases in detection of different species. Multiple species were detected in all clinical cases for which sufficient nox sequence data were generated (n = 47), indicating that mixed species Brachyspira infections are common, although in most cases, one species accounted for at least half of the sequences identified. In all cases, the species detected in the original diagnostic investigation of each case was also detected by nox sequencing. Results from synthetic communities indicated that the method was highly reproducible, but also indicated potential PCR bias against B. hampsonii genomovar I. Deep sequencing of the nox gene target is a suitable method for simultaneous detection of multiple Brachyspira species in clinical case material that may offer advantages over current, more targeted diagnostic approaches for investigating the significance of mixed infections in disease development.

Major Histocompatibility Complex Allele Persistence in Eurasia and America in the Genus Carduelis (Spinus) During Million Years

Introduction:

GenusCarduelis(Fringillidaefamily) includes goldfinches, siskins, redpolls, greenfinches and crossbills. Many of the species classified within this genus and other related genera have been grouped by using molecular systematics and the mitochondrial cytochrome b (mt cyt b) gene. According to this, the Eurasian siskin (C. spinus)is the only one extant direct ancestor of several North American finches; North American / South American radiations may have been originated by Eurasian siskin (or extinct relative). In the present work, we aim to perform a study of transpecies and transcontinental analyses of MHC (Major Histocompatibility Complex) Class I alleles in several genusCarduelis/Spinusspecies in order to draw evolutionary conclusions in several wild bird species belonging to the genusCarduelis / Spinus.

Materials and Methods:

Blood was taken from worldwide wild bird species. Passerine phylogeny was done after analysing mtDNA with Maximun Likelihood and Bayesian dendrograms. Major histocompatibility complex alleles were obtained by standard DNA cloning and sequencing.

Results:

We found two matches between MHC-I DNA alleles from different South American siskins at DNA level. Also, it was observed that the Eurasian siskin shares a protein with pine siskin and another with three South American siskins. Eight South American siskins species also share the same MHC protein. In addition, studied songbirds MHC class I intron 2 is longer than that ofGallus gallus.

Conclusion:

We have drawn the following conclusions: 1) We present the first direct evidence that “Minimal Essential MHC” does not exist for birds; one of its main definition characters,i.e.: small intron size does not hold for songbirds. 2) We also report that MHC genes transpecies evolution exist in birds by showing also for the first time that worldwide bird species keep the same MHC protein and DNA alleles. 3) New evidences on MHC alleles conservation from EurasianCarduelis spinus(most ancient) to South American siskins (most recent) during million years support that Eurasian siskin is the parental species for American GenusCarduelis (Spinus)species. It is uncertain whether Eurasian siskin (or extant relative) had initially an Holoartic distribution, including America.

Genes of the major histocompatibility complex highlight interactions of the innate and adaptive immune system

Background

A well-functioning immune defence is crucial for fitness, but our knowledge about the immune system and its complex interactions is still limited. Major histocompatibility complex (MHC) molecules are involved in T-cell mediated adaptive immune responses, but MHC is also highly upregulated during the initial innate immune response. The aim of our study was therefore to determine to what extent the highly polymorphic MHC is involved in interactions of the innate and adaptive immune defence and if specific functional MHC alleles (FA) or heterozygosity at the MHC are more important.

Methods

To do this we used captive house sparrows (Passer domesticus) to survey MHC diversity and immune function controlling for several environmental factors. MHC class I alleles were identified using parallel amplicon sequencing and to mirror immune function, several immunological tests that correspond to the innate and adaptive immunity were conducted.

Results

Our results reveal that MHC was linked to all immune tests, highlighting its importance for the immune defence. While all innate responses were associated with one single FA, adaptive responses (cell-mediated and humoral) were associated with several different alleles.

Discussion

We found that repeated injections of an antibody in nestlings and adults were linked to different FA and hence might affect different areas of the immune system. Also, individuals with a higher number of different FA produced a smaller secondary response, indicating a disadvantage of having numerous MHC alleles. These results demonstrate the complexity of the immune system in relation to the MHC and lay the foundation for other studies to further investigate this topic.

A quantitative and qualitative comparison of illumina MiSeq and 454 amplicon sequencing for genotyping the highly polymorphic major histocompatibility complex (MHC) in a non-model species

Background

High-throughput sequencing enables high-resolution genotyping of extremely duplicated genes. 454 amplicon sequencing (454) has become the standard technique for genotyping the major histocompatibility complex (MHC) genes in non-model organisms. However, illumina MiSeq amplicon sequencing (MiSeq), which offers a much higher read depth, is now superseding 454. The aim of this study was to quantitatively and qualitatively evaluate the performance of MiSeq in relation to 454 for genotyping MHC class I alleles using a house sparrow (Passer domesticus) dataset with pedigree information. House sparrows provide a good study system for this comparison as their MHC class I genes have been studied previously and, consequently, we had prior expectations concerning the number of alleles per individual.

Results

We found that 454 and MiSeq performed equally well in genotyping amplicons with low diversity, i.e. amplicons from individuals that had fewer than 6 alleles. Although there was a higher rate of failure in the 454 dataset in resolving amplicons with higher diversity (6–9 alleles), the same genotypes were identified by both 454 and MiSeq in 98% of cases.

Conclusions

We conclude that low diversity amplicons are equally well genotyped using either 454 or MiSeq, but the higher coverage afforded by MiSeq can lead to this approach outperforming 454 in amplicons with higher diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-017-2654-1) contains supplementary material, which is available to authorized users.

Major histocompatibility complex genes partly explain early survival in house sparrows

Environmental factors and genetic incompatibilities between parents have been suggested as important determinants for embryonic mortality and survival. The genetic set-up of the immune system, specifically the highly polymorphic major histocompatibility complex (MHC) may also influence individual resistance to infections. MHC proteins are important for an appropriate adaptive immune response and enable T-cells to separate 'self' from 'non-self'. Here we investigate the importance of MHC functional diversity for early development in birds, more specifically, if offspring survival and body mass or size depends on number of different functional MHC alleles, specific functional MHC alleles or similarity of MHC alleles in the parents. Unhatched eggs are common in clutches of many bird species. In house sparrows (Passer domesticus), embryo and nestling mortality can exceed 50%. To control for environmental factors, our study was carried out on an aviary population. We found that one specific functional MHC allele was associated with reduced nestling survival, which was additionally supported by lower body mass and a smaller tarsus when nestlings have been 6 days old. Another allele was positively associated with tarsus length at a later nestling stage (nestlings 12 days old). These results indicate that MHC alleles might influence pathogen resistance or susceptibility.

Extreme MHC class I diversity in the sedge warbler (Acrocephalus schoenobaenus); selection patterns and allelic divergence suggest that different genes have different functions

Background

Recent work suggests that gene duplications may play an important role in the evolution of immunity genes. Passerine birds, and in particular Sylvioidea warblers, have highly duplicated major histocompatibility complex (MHC) genes, which are key in immunity, compared to other vertebrates. However, reasons for this high MHC gene copy number are yet unclear. High-throughput sequencing (HTS) allows MHC genotyping even in individuals with extremely duplicated genes. This HTS data can reveal evidence of selection, which may help to unravel the putative functions of different gene copies, i.e. neofunctionalization. We performed exhaustive genotyping of MHC class I in a Sylvioidea warbler, the sedge warbler, Acrocephalus schoenobaenus, using the Illumina MiSeq technique on individuals from a wild study population.

Results

The MHC diversity in 863 genotyped individuals by far exceeds that of any other bird species described to date. A single individual could carry up to 65 different alleles, a large proportion of which are expressed (transcribed). The MHC alleles were of three different lengths differing in evidence of selection, diversity and divergence within our study population. Alleles without any deletions and alleles containing a 6 bp deletion showed characteristics of classical MHC genes, with evidence of multiple sites subject to positive selection and high sequence divergence. In contrast, alleles containing a 3 bp deletion had no sites subject to positive selection and had low divergence.

Conclusions

Our results suggest that sedge warbler MHC alleles that either have no deletion, or contain a 6 bp deletion, encode classical antigen presenting MHC molecules. In contrast, MHC alleles containing a 3 bp deletion may encode molecules with a different function. This study demonstrates that highly duplicated MHC genes can be characterised with HTS and that selection patterns can be useful for revealing neofunctionalization. Importantly, our results highlight the need to consider the putative function of different MHC genes in future studies of MHC in relation to disease resistance and fitness.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0997-9) contains supplementary material, which is available to authorized users.

Expression and phylogenetic analyses reveal paralogous lineages of putatively classical and non-classical MHC-I genes in three sparrow species (Passer)

BACKGROUND

The Major Histocompatibility Complex (MHC) plays a central role in immunity and has been given considerable attention by evolutionary ecologists due to its associations with fitness-related traits. Songbirds have unusually high numbers of MHC class I (MHC-I) genes, but it is not known whether all are expressed and equally important for immune function. Classical MHC-I genes are highly expressed, polymorphic and present peptides to T-cells whereas non-classical MHC-I genes have lower expression, are more monomorphic and do not present peptides to T-cells. To get a better understanding of the highly duplicated MHC genes in songbirds, we studied gene expression in a phylogenetic framework in three species of sparrows (house sparrow, tree sparrow and Spanish sparrow), using high-throughput sequencing. We hypothesize that sparrows could have classical and non-classical genes, as previously indicated though never tested using gene expression.

RESULTS

The phylogenetic analyses reveal two distinct types of MHC-I alleles among the three sparrow species, one with high and one with low level of polymorphism, thus resembling classical and non-classical genes, respectively. All individuals had both types of alleles, but there was copy number variation both within and among the sparrow species. However, the number of highly polymorphic alleles that were expressed did not vary between species, suggesting that the structural genomic variation is counterbalanced by conserved gene expression. Overall, 50% of the MHC-I alleles were expressed in sparrows. Expression of the highly polymorphic alleles was very variable, whereas the alleles with low polymorphism had uniformly low expression. Interestingly, within an individual only one or two alleles from the polymorphic genes were highly expressed, indicating that only a single copy of these is highly expressed.

CONCLUSIONS

Taken together, the phylogenetic reconstruction and the analyses of expression suggest that sparrows have both classical and non-classical MHC-I genes, and that the evolutionary origin of these genes predate the split of the three investigated sparrow species 7 million years ago. Because only the classical MHC-I genes are involved in antigen presentation, the function of different MHC-I genes should be considered in future ecological and evolutionary studies of MHC-I in sparrows and other songbirds.

Characterization of MHC class I in a long distance migratory wader, the Icelandic black-tailed godwit

The major histocompatibility complex (MHC) encodes proteins that are central for antigen presentation and pathogen elimination. MHC class I (MHC-I) genes have attracted a great deal of interest among researchers in ecology and evolution and have been partly characterized in a wide range of bird species. So far, the main focus has been on species within the bird orders Galliformes and Passeriformes, while Charadriiformes remain vastly underrepresented with only two species studied to date. These two Charadriiformes species exhibit striking differences in MHC-I characteristics and MHC-I diversity. We therefore set out to study a third species within Charadriiformes, the Icelandic subspecies of black-tailed godwits (Limosa limosa islandica). This subspecies is normally confined to parasite-poor environments, and we hence expected low MHC diversity. MHC-I was partially characterized first using Sanger sequencing and then using high-throughput sequencing (MiSeq) in 84 individuals. We verified 47 nucleotide alleles in open reading frame with classical MHC-I characteristics, and each individual godwit had two to seven putatively classical MHC alleles. However, in contrast to previous MHC-I data within Charadriiformes, we did not find any evidence of alleles with low sequence diversity, believed to represent non-classical MHC genes. The diversity and divergence of the godwits MHC-I genes to a large extent fell between the previous estimates within Charadriiformes. However, the MHC genes of the migratory godwits had few sites subject to positive selection, and one possible explanation could be a low exposure to pathogens.

Extensive Allelic Diversity of MHC Class I in Wild Mallard Ducks

MHC class I is critically involved in defense against viruses, and diversity from polygeny and polymorphism contributes to the breadth of the immune response and health of the population. In this article, we examine MHC class I diversity in wild mallard ducks, the natural host and reservoir of influenza A viruses. We previously showed domestic ducks predominantly use UAA, one of five MHC class I genes, but whether biased expression is also true for wild mallards is unknown. Using RT-PCR from blood, we examined expressed MHC class I alleles from 38 wild mallards (Anas platyrhynchos) and identified 61 unique alleles, typically 1 or 2 expressed alleles in each individual. To determine whether expressed alleles correspond to UAA adjacent to TAP2 as in domestic ducks, we cloned and sequenced genomic UAA-TAP2 fragments from all mallards, which matched transcripts recovered and allowed us to assign most alleles as UAA Allelic differences are primarily located in α1 and α2 domains in the residues known to interact with peptide in mammalian MHC class I, suggesting the diversity is functional. Most UAA alleles have unique residues in the cleft predicting distinct specificity; however, six alleles have an unusual conserved cleft with two cysteine residues. Residues that influence peptide-loading properties and tapasin involvement in chicken are fixed in duck alleles and suggest tapasin independence. Biased expression of one MHC class I gene may make viral escape within an individual easy, but high diversity in the population places continual pressure on the virus in the reservoir species.

Characterization of major histocompatibility complex class I loci of the lark sparrow (Chondestes grammacus) and insights into avian MHC evolution

The major histocompatibilty complex (MHC) has become increasingly important in the study of the immunocapabilities of non-model vertebrates due to its direct involvement in the immune response. The characterization of MHC class I loci in the lark sparrow (Chondestes grammacus) revealed multiple MHC class I loci with elevated genetic diversity at exon 3, evidence of differential selection between the peptide binding region (PBR) and non-PBR, and the presence of multiple pseudogenes with limited divergence. The minimum number of functional MHC class I loci was estimated at four. Sequence analysis revealed d N /d S ratios significantly less than one at non-PBR sites, indicative of negative selection, whereas PBR sites associated with antigen recognition showed ratios greater than 1 but non-significant. GenBank surveys and phylogenetic analyses of previously reported avian MHC class I sequences revealed variable signatures of evolutionary processes acting upon this gene family, including gene duplication and potential concerted evolution. An increase in the number of class I loci across species coincided with an increase in pseudogene prevalence, revealing the importance of gene duplication in the expansion of multigene families and the creation of pseudogenes.

Characterization of the house sparrow (Passer domesticus) transcriptome: a resource for molecular ecology and immunogenetics

The house sparrow (Passer domesticus) is an important model species in ecology and evolution. However, until recently, genomic resources for molecular ecological projects have been lacking in this species. Here, we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximize the yield of immune-related gene functions. After de novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3177 single-nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.