Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons.

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

The structure of songbird MHC class I reveals antigen binding that is flexible at the N-terminus and static at the C-terminus.

Long-distance migratory animals such as birds and bats have evolved to withstand selection imposed by pathogens across the globe, and pathogen richness is known to be particularly high in tropical regions. Immune genes, so-called Major Histocompatibility Complex (MHC) genes, are highly duplicated in songbirds compared to other vertebrates, and this high MHC diversity has been hypothesised to result in a unique adaptive immunity. To understand the rationale behind the evolution of the high MHC genetic diversity in songbirds, we determined the structural properties of an MHC class I protein, Acar3, from a long-distance migratory songbird, the great reed warbler Acrocephalus arundinaceus (in short: Acar). The structure of Acar3 was studied in complex with pathogen-derived antigens and shows an overall antigen presentation similar to human MHC class I. However, the peptides bound to Acar3 display an unusual conformation: Whereas the N-terminal ends of the peptides display enhanced flexibility, the conformation of their C-terminal halves is rather static. This uncommon peptide-binding mode in Acar3 is facilitated by a central Arg residue within the peptide-binding groove that fixes the backbone of the peptide at its central position, and potentially permits successful interactions between MHC class I and innate immune receptors. Our study highlights the importance of investigating the immune system of wild animals, such as birds and bats, to uncover unique immune mechanisms which may neither exist in humans nor in model organisms.

The genomic architecture of the passerine MHC region: High repeat content and contrasting evolutionary histories of single copy and tandemly duplicated MHC genes.

The major histocompatibility complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assembled and annotated the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5 Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologues among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

Seasonal dynamics of haemosporidian (Apicomplexa, Haemosporida) parasites in house sparrows Passer domesticus at four European sites: comparison between lineages and the importance of screening methods.

Infectious diseases often vary seasonally in a predictable manner, and seasonality may be responsible for geographical differences in prevalence. In temperate regions, vector-borne parasites such as malaria are expected to evolve lower virulence and a time-varying strategy to invest more in transmission when vectors are available. A previous model of seasonal variation of avian malaria described a double peak in prevalence of Plasmodium parasites in multiple hosts resulting from spring relapses and transmission to susceptible individuals in summer. However, this model was rejected by a study describing different patterns of seasonal variation of two Plasmodium spp. at the same site, with the double peak only apparent when these species were combined. Here, we assessed the seasonal variation in prevalence of haemosporidian parasites (Plasmodium, Haemoproteus and Leucocytozoon) in house sparrows (Passer domesticus) sampled across 1 year at four temperate European sites spanning a latitudinal range of 17°. We showed that parasite prevalence and diversity decreased with increasing latitude, but the parasite communities differed between sites, with only one Plasmodium lineage (P_SGS1) occurring at all sites. Moreover, the nested PCR method commonly used to detect and identify haemosporidian parasites strongly underestimated co-infections of Haemoproteus and Plasmodium, significantly biasing the pattern of seasonal variation, so additional molecular methods were used. Finally, we showed that: (i) seasonal variation in prevalence of haemosporidian parasites varied between study sites and parasite lineages/species/genera, describing further cases where the double peak model is not met; (ii) the seasonal dynamics of single lineages (P_SGS1) varied between sites; and (iii) unexpectedly, seasonality was greatest at the southernmost site, a pattern that was mostly driven by lineage H_PADOM05. Limitations of the genotyping methods and consequences of pooling (parasite lineages, sites and years) in studies of haemosporidian parasites are discussed and recommendations proposed, since these actions may obscure the patterns of prevalence and limit ecological inferences.

Parental genetic similarity and offspring performance in blue tits in relation to brood size manipulation.

In birds, as in many other taxa, higher genetic similarity of mates has long been known to reduce offspring fitness. To date, the majority of avian studies have focused on examination whether the genetic similarity of social mates predicts hatching success. Yet, increased genetic similarity of mates may also reduce offspring fitness during later life stages, including the nestling period and beyond. Here, we investigated whether parental genetic similarity influences offspring performance using data from free-living blue tits (Cyanistes caeruleus) collected across three breeding seasons. Additionally, we tested whether brood size manipulation affects the magnitude and direction of the relationship between genetic similarity of mates and offspring performance. Sixteen microsatellite markers were used to measure genetic similarity between biological parents. We found that the genetic similarity of parents negatively affects offspring immune response and this effect was independent of the experimental brood size manipulation.