MHC class II pseudogene and genomic signature of a 32-kb cosmid in the house finch (Carpodacus mexicanus)

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

Using de novo genome assembly and high-throughput sequencing to characterize the MHC region in a non-model bird, the Eurasian coot

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

Trade-offs in expressed major histocompatibility complex diversity seen on a macroevolutionary scale among songbirds

To survive organisms must defend themselves against pathogens. Classical Major Histocompatibility Complex (MHC) genes play a key role in pathogen defense by encoding molecules involved in pathogen recognition. MHC gene diversity influences the variety of pathogens individuals can recognize and respond to and has consequently been a popular genetic marker for disease resistance in ecology and evolution. However, MHC diversity is predominantly estimated using genomic DNA (gDNA) with little knowledge of expressed diversity. This limits our ability to interpret the adaptive significance of variation in MHC diversity, especially in species with very many MHC genes such as songbirds. Here, we address this issue using phylogenetic comparative analyses of the number of MHC class I alleles (MHC-I diversity) in gDNA and complementary DNA (cDNA), that is, expressed alleles, across 13 songbird species. We propose three theoretical relationships that could be expected between genomic and expressed MHC-I diversity on a macroevolutionary scale and test which of these are best supported. In doing so, we show that significantly fewer MHC-I alleles than the number available are expressed, suggesting that optimal MHC-I diversity could be achieved by modulating gene expression. Understanding the relationship between genomic and expressed MHC diversity is essential for interpreting variation in MHC diversity in an evolutionary context.

Long-Read Genome Assemblies Reveal Extraordinary Variation in the Number and Structure of MHC Loci in Birds

Our knowledge of the Major Histocompatibility Complex (MHC) in birds is limited because it often consists of numerous duplicated genes within individuals that are difficult to assemble with short read sequencing technology. Long-read sequencing provides an opportunity to overcome this limitation because it allows the assembly of long regions with repetitive elements. In this study, we used genomes based on long-read sequencing to predict the number and location of MHC loci in a broad range of bird taxa. From the long-read-based genomes of 34 species, we found that there was extremely large variation in the number of MHC loci between species. Overall, there were greater numbers of both class I and II loci in passerines than nonpasserines. The highest numbers of loci (up to 193 class II loci) were found in manakins (Pipridae), which had previously not been studied at the MHC. Our results provide the first direct evidence from passerine genomes of this high level of duplication. We also found different duplication patterns between species. In some species, both MHC class I and II genes were duplicated together, whereas in most species they were duplicated independently. Our study shows that the analysis of long-read-based genomes can dramatically improve our knowledge of MHC structure, although further improvements in chromosome level assembly are needed to understand the evolutionary mechanisms producing the extraordinary interspecific variation in the architecture of the MHC region.

Avian MHC Evolution in the Era of Genomics: Phase 1.0

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.

Duplication and Variation in the Major Histocompatibility Complex Genes in Blakiston's Fish Owl, Bubo blakistoni

The major histocompatibility complex (MHC) includes many genes that are essential for the adaptive immune system, and variation in the antigen binding site (ABS) is related to resistance against pathogens. In the present study, quantitative real-time PCR indicated a larger number of MHC gene copies in the endangered population of Blakiston's fish owl (Bubo blakistoni) than in five other owl species, and massively parallel pyrosequencing detected more MHC class IIβ per individual alleles in B. blakistoni than in the other species. A chromosomal fluorescence in situ hybridization (FISH) analysis showed that the MHC class I and class IIβ loci are closely linked on a single pair of microchromosomes, indicating that the MHC genes were tandemly duplicated in a limited chromosomal region. Because B. blakistoni has twice as many MHC genes as its sister species, the tawny fish owl (Bubo flavipes), the duplication of MHC genes occurred after these species diverged by speciation. A Bayesian molecular phylogenetic analysis showed that the DAB1 and DAB2 lineages of MHC class IIβ alleles from various strigid species each formed a separate clade, indicating that the two allelic lineages preceded the radiation of Strigidae and evolved as paralogs. By contrast, the ABS sequences did not form distinct clades between DAB1 and DAB2 alleles but were intermixed, presumably due to gene conversion. Despite the low diversity of alleles per locus, B. blakistoni had many lineages of MHC class IIβ alleles. Gene duplication increases variation in the MHC genes in this species, and could have facilitated adaptation in small populations.

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.

100 million years of multigene family evolution: origin and evolution of the avian MHC class IIB

BACKGROUND

Gene duplication has led to a most remarkable adaptation involved in vertebrates' host-pathogen arms-race, the major histocompatibility complex (MHC). However, MHC duplication history is as yet poorly understood in non-mammalian vertebrates, including birds.

RESULTS

Here, we provide evidence for the evolution of two ancient avian MHC class IIB (MHCIIB) lineages by a duplication event prior to the radiation of all extant birds >100 million years ago, and document the role of concerted evolution in eroding the footprints of the avian MHCIIB duplication history.

CONCLUSIONS

Our results suggest that eroded footprints of gene duplication histories may mimic birth-death evolution and that in the avian MHC the presence of the two lineages may have been masked by elevated rates of concerted evolution in several taxa. Through the presence of a range of intermediate evolutionary stages along the homogenizing process of concerted evolution, the avian MHCIIB provides a remarkable illustration of the erosion of multigene family duplication history.

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next-generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans-species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long-term survival of the species.

Characterization of class II β chain major histocompatibility complex genes in a family of Hawaiian honeycreepers: 'amakihi (Hemignathus virens)

Hawaiian honeycreepers (Drepanidinae) have evolved in the absence of mosquitoes for over five million years. Through human activity, mosquitoes were introduced to the Hawaiian archipelago less than 200 years ago. Mosquito-vectored diseases such as avian malaria caused by Plasmodium relictum and Avipoxviruses have greatly impacted these vulnerable species. Susceptibility to these diseases is variable among and within species. Due to their function in adaptive immunity, the role of major histocompatibility complex genes (Mhc) in disease susceptibility is under investigation. In this study, we evaluate gene organization and levels of diversity of Mhc class II β chain genes (exon 2) in a captive-reared family of Hawaii 'amakihi (Hemignathus virens). A total of 233 sequences (173 bp) were obtained by PCR+1 amplification and cloning, and 5720 sequences were generated by Roche 454 pyrosequencing. We report a total of 17 alleles originating from a minimum of 14 distinct loci. We detected three linkage groups that appear to represent three distinct haplotypes. Phylogenetic analysis revealed one variable cluster resembling classical Mhc sequences (DAB) and one highly conserved, low variability cluster resembling non-classical Mhc sequences (DBB). High net evolutionary divergence values between DAB and DBB resemble that seen between chicken BLB system and YLB system genes. High amino acid identity among non-classical alleles from 12 species of passerines (DBB) and four species of Galliformes (YLB) was found, suggesting that these non-classical passerine sequences may be related to the Galliforme YLB sequences.

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine

The environment shapes host-parasite interactions, but how environmental variation affects the diversity and composition of parasite-defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC-parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous-collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra-individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency-dependent selection. We found an MHC variant associated with a ∽26% decrease in infection probability at middle elevations (1501-3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency-dependent selection and fluctuating spatial selection in the evolution of Z. capensisMHC.

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

In birds the major histocompatibility complex (MHC) organization differs both among and within orders; chickens Gallus gallus of the order Galliformes have a simple arrangement, while many songbirds of the order Passeriformes have a more complex arrangement with larger numbers of MHC class I and II genes. Chicken MHC genes are found at two independent loci, classical MHC-B and non-classical MHC-Y, whereas non-classical MHC genes are yet to be verified in passerines. Here we characterize MHC class I transcripts (α1 to α3 domain) and perform amplicon sequencing using a next-generation sequencing technique on exon 3 from house sparrow Passer domesticus (a passerine) families. Then we use phylogenetic, selection, and segregation analyses to gain a better understanding of the MHC class I organization. Trees based on the α1 and α2 domain revealed a distinct cluster with short terminal branches for transcripts with a 6-bp deletion. Interestingly, this cluster was not seen in the tree based on the α3 domain. 21 exon 3 sequences were verified in a single individual and the average numbers within an individual were nine and five for sequences with and without a 6-bp deletion, respectively. All individuals had exon 3 sequences with and without a 6-bp deletion. The sequences with a 6-bp deletion have many characteristics in common with non-classical MHC, e.g., highly conserved amino acid positions were substituted compared with the other alleles, low nucleotide diversity and just a single site was subject to positive selection. However, these alleles also have characteristics that suggest they could be classical, e.g., complete linkage and absence of a distinct cluster in a tree based on the α3 domain. Thus, we cannot determine for certain whether or not the alleles with a 6-bp deletion are non-classical based on our present data. Further analyses on segregation patterns of these alleles in combination with dating the 6-bp deletion through MHC characterization across the genus Passer may solve this matter in the future.

MHC class IIB exon 2 polymorphism in the Grey partridge (Perdix perdix) is shaped by selection, recombination and gene conversion

Among bird species, the most studied major histocompatibility complex (MHC) is the chicken MHC. Although the number of studies on MHC in free-ranging species is increasing, the knowledge on MHC variation in species closely related to chicken is required to understand the peculiarities of bird MHC evolution. Here we describe the variation of MHC class IIB (MHCIIB) exon 2 in a population of the Grey partridge (Perdix perdix), a species of high conservation concern throughout Europe and an emerging galliform model in studies of sexual selection. We found 12 alleles in 108 individuals, but in comparison to other birds surprisingly many sites show signatures of historical positive selection. Individuals displayed between two to four alleles both on genomic and complementary DNA, suggesting the presence of two functional MHCIIB loci. Recombination and gene conversion appear to be involved in generating MHCIIB diversity in the Grey partridge; two recombination breakpoints and several gene conversion events were detected. In phylogenetic analysis of galliform MHCIIB, the Grey partridge alleles do not cluster together, but are scattered through the tree instead. Thus, our results indicate that the Grey partridge MHCIIB is comparable to most other galliforms in terms of copy number and population polymorphism.

Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC

BACKGROUND

The MHC, which is regarded as the most polymorphic region in the genomes of jawed vertebrates, plays a central role in the immune system by encoding various proteins involved in the immune response. The chicken MHC-B genomic region has a highly streamlined gene content compared to mammalian MHCs. Its core region includes genes encoding Class I and Class IIB molecules but is only ~92Kb in length. Sequences of other galliform MHCs show varying degrees of similarity as that of chicken. The black grouse (Tetrao tetrix) is a wild galliform bird species which is an important model in conservation genetics and ecology. We sequenced the black grouse core MHC-B region and combined this with available data from related species (chicken, turkey, gold pheasant and quail) to perform a comparative genomics study of the galliform MHC. This kind of analysis has previously been severely hampered by the lack of genomic information on avian MHC regions, and the galliformes is still the only bird lineage where such a comparison is possible.

RESULTS

In this study, we present the complete genomic sequence of the MHC-B locus of black grouse, which is 88,390 bp long and contains 19 genes. It shows the same simplicity as, and almost perfect synteny with, the corresponding genomic region of chicken. We also use 454-transcriptome sequencing to verify expression in 17 of the black grouse MHC-B genes. Multiple sequence inversions of the TAPBP gene and TAP1-TAP2 gene block identify the recombination breakpoints near the BF and BLB genes. Some of the genes in the galliform MHC-B region also seem to have been affected by selective forces, as inferred from deviating phylogenetic signals and elevated rates of non-synonymous nucleotide substitutions.

CONCLUSIONS

We conclude that there is large synteny between the MHC-B region of the black grouse and that of other galliform birds, but that some duplications and rearrangements have occurred within this lineage. The MHC-B sequence reported here will provide a valuable resource for future studies on the evolution of the avian MHC genes and on links between immunogenetics and ecology of black grouse.

CR1 retroposons provide a new insight into the phylogeny of Phasianidae species (Aves: Galliformes)

Chicken repeat 1 (CR1) elements, a class of retroposons belonging to non-long-terminal repeats, have been recognized as powerful tools for phylogenetic studies. Here we examine the phylogenetic relationships of 11 Phasianidae species based on CR1 retroposons. Together with 19 loci reported previously, a total of 99 CR1 loci were identified from chicken genome and turkey BAC clone sequences. 75 insertion events were used to address the branching order of 11 species in Phasianidae. The topology of our tree suggests that: 1) Gallus gallus possessed a basal phylogenetic position within Phasianidae and was related to Bambusicola thoracica (BSP=100%); 2) After the split of G. gallus and B. thoracica, Arborophila rufipectus diverged from Phasianidae (BSP=100%). Nine unambiguous insertion events supported a phylogenetic position of A. rufipectus different to previous mitochondrial data suggesting a hybrid origin or an ancient introgression of A. rufipectus; and 3) 22 CR1 insertion events strongly supported the eight phasianids under investigation sharing a common ancestor. Our study has revisited the phylogenetic position of G. gallus and A. rufipectus and provided a new insight into the phylogeny of Phasianidae birds. It showed that a CR1-based methodology has a great potential to be informative within Phasianidae in resolving relationships of closely related species whose radiation and speciation have occurred very recently.

Maintenance of MHC Class IIB diversity in a recently established songbird population

We examined variation at MHC Class IIB genes in a recently established population of dark-eyed juncos (Junco hyemalis) in a coastal urban environment in southern California, USA relative to an ancestral-range population from a nearby species-typical montane environment. The founding population is estimated to have been quite small, but we predicted that variation at the major histocompatibility complex (MHC) among the founders would nevertheless be preserved owing to the high functional significance of MHC. Previous studies of MHC in songbirds have had varying degrees of success in isolating loci, as passerines show extensive MHC gene duplication. In order to compare diversity in the two populations, we employed two published approaches to sequencing MHC Class II exon 2: direct sequencing with exon-based primers, and traditional cloning and sequencing with intron-based primers. Results from both methods show that the colonist population has maintained high levels of variation. Our results also indicate varying numbers of alleles across individuals, corroborating evidence for gene duplication in songbird MHC. While future studies in songbirds may need to take a genomic approach to fully understand the structure of MHC in this lineage, our results show that it is possible to use traditional methods to reveal functional variation across populations.

Molecular cloning, polymorphism and tissue distribution of the MHC class IIB gene in the Chinese goose (Anser cygnoides)

1. The goose major histocompatibility complex (MHC) class IIB cDNA (Ancy-MHCII) was cloned by homology cloning and rapid amplification of cDNA ends by polymerase chain reaction (RACE-PCR), and the genomic structure and tissue expression were investigated. 2. Three different 5'-RACE sequences (Ancy-MHC II5'-1, Ancy-MHC II5'-2, Ancy-MHC II5'-3), one 3'-RACE sequence (Ancy-MHC II-3') and two different full length Ancy-MHC IIB cDNA sequences (Ancy-CD01, Ancy-CD02), which came from different alleles at one locus or different loci, were determined. 3. The genomic organisation is composed of 6 exons and 5 introns, with a longer intron region than that of the chicken. The alleles encode 259 and 260 amino acids in the mature protein. 4. The number of non-synonymous substitutions (dN) in the peptide-binding region of exon 2 from 8 alleles was higher than that of the synonymous substitutions (dS). 5. Tissue-specific expression of Ancy-MHC II mRNA was detected in an adult goose using RT-PCR. These results showed that Ancy-MHC II mRNA was expressed in the lung, spleen, liver, intestine, heart, kidney, pancreas, brain, skin and muscle. This is consistent with the expression of MHC class IIB in various tissues from the chicken. 6. Sequences from goose, snipe and duck clustered together when compared with known MHC class IIB sequences from the other species, significantly differing from mammals and aquatic species, indicating a pattern consistent with accepted evolutionary pathways.

Spectrum of MHC class II variability in Darwin's finches and their close relatives

The study describes >400 major histocompatibility complex (MHC) class II B exon 2 and 114 intron 2 sequences of 36 passerine bird species, 13 of which belong to the group of Darwin's finches (DFs) and the remaining 23 to close or more distant relatives of DFs in Central and South America. The data set is analyzed by a combination of judiciously selected statistical methods. The analysis reveals that reliable information concerning MHC organization, including the assignment of sequences to loci, and evolution, as well as the process of species divergence, can be obtained in the absence of genomic sequence data, if the analysis is taken several steps beyond the standard phylogenetic tree construction approach. The main findings of the present study are these: The MHC class II B region of the passerine birds is as elaborate in its organization, divergence, and genetic diversity as the MHC of the eutherian mammals, specifically the primates. Hence, the reported simplicity of the fowl MHC is an oddity. With the help of appropriate markers, the divergence of the MHC genes can be traced deep in the phylogeny of the bird taxa. Transspecies polymorphism is rampant at many of the bird MHC loci. In this respect, the DFs behave as if they were a single, genetically undifferentiated population. There is thus far no indication of alleles that could be considered species, genus, or even DF group specific. The implication of these findings is that DFs are in the midst of adaptive radiations, in which morphological differentiation into species is running ahead of genetic differentiation in genetic systems such as the MHC or the mitochondrial DNA. The radiations are so young that there has not been enough time to sort out polymorphisms at most of the loci among the morphologically differentiating species. These findings parallel those on Lake Victoria haplochromine fishes. Several of the DF MHC allelic lineages can be traced back to the MHC genes of the species Tiaris obscura, which we identified previously as the closest extant relative of DFs in continental America.

Reduced MHC and neutral variation in the Galápagos hawk, an island endemic

BACKGROUND

Genes at the major histocompatibility complex (MHC) are known for high levels of polymorphism maintained by balancing selection. In small or bottlenecked populations, however, genetic drift may be strong enough to overwhelm the effect of balancing selection, resulting in reduced MHC variability. In this study we investigated MHC evolution in two recently diverged bird species: the endemic Galápagos hawk (Buteo galapagoensis), which occurs in small, isolated island populations, and its widespread mainland relative, the Swainson's hawk (B. swainsoni).

RESULTS

We amplified at least two MHC class II B gene copies in each species. We recovered only three different sequences from 32 Galápagos hawks, while we amplified 20 unique sequences in 20 Swainson's hawks. Most of the sequences clustered into two groups in a phylogenetic network, with one group likely representing pseudogenes or nonclassical loci. Neutral genetic diversity at 17 microsatellite loci was also reduced in the Galápagos hawk compared to the Swainson's hawk.

CONCLUSIONS

The corresponding loss in neutral diversity suggests that the reduced variability present at Galápagos hawk MHC class II B genes compared to the Swainson's hawk is primarily due to a founder event followed by ongoing genetic drift in small populations. However, purifying selection could also explain the low number of MHC alleles present. This lack of variation at genes involved in the adaptive immune response could be cause for concern should novel diseases reach the archipelago.

Characterization and evolution of MHC class II B genes in Ardeid birds

Major histocompatibility complex (MHC) is a multi-gene family that is very suitable to investigate a wide range of open questions in evolutionary ecology. In this study, we characterized two expressed MHC class II B genes (DAB1 and DAB2) in the Grey Heron (Aves: Ardea cinerea). We further developed the primer pairs to amplify and sequence two MHC class II B loci in ten ardeid birds. Phylogenetic analysis revealed that different parts of the genes showed different evolutionary patterns. The exon 2 sequences tended to cluster two gene-specific lineages. In each lineage, exon 2 sequences from several species showed closer relationships than sequences within species, and two shared identical alleles were found between species (Egretta sacra and Nycticorax nycticorax; Egretta garzetta and Bubulcus ibis), supporting the hypothesis of trans-species polymorphism. In contrast, the species-specific intron 2 plus partial exon 3 tree suggested that DAB1 and DAB2 were subject to concerted evolution. GENECONV analyses showed the gene exchange played an important role in the ardeid MHC evolution.

Genetic mapping of the major histocompatibility complex in the zebra finch (Taeniopygia guttata)

Genes of the major histocompatibility complex (MHC) have received much attention in immunology, genetics, and ecology because they are highly polymorphic and play important roles in parasite resistance and mate choice. Until recently, the MHC of passerine birds was not well-described. However, the genome sequencing of the zebra finch (Taeniopygia guttata) has partially redressed this gap in our knowledge of avian MHC genes. Here, we contribute further to the understanding of the zebra finch MHC organization by mapping SNPs within or close to known MHC genes in the zebra finch genome. MHC class I and IIB genes were both mapped to zebra finch chromosome 16, and there was no evidence that MHC class I genes are located on chromosome 22 (as suggested by the genome assembly). We confirm the location in the MHC region on chromosome 16 for several other genes (BRD2, FLOT1, TRIM7.2, GNB2L1, and CSNK2B). Two of these (CSNK2B and FLOT1) have not previously been mapped in any other bird species. In line with previous results, we also find that orthologs to the immune-related genes B-NK and CLEC2D, which are part of the MHC region in chicken, are situated on zebra finch chromosome Z and not among other MHC genes in the zebra finch.

454 sequencing reveals extreme complexity of the class II Major Histocompatibility Complex in the collared flycatcher

Background

Because of their functional significance, the Major Histocompatibility Complex (MHC) class I and II genes have been the subject of continuous interest in the fields of ecology, evolution and conservation. In some vertebrate groups MHC consists of multiple loci with similar alleles; therefore, the multiple loci must be genotyped simultaneously. In such complex systems, understanding of the evolutionary patterns and their causes has been limited due to challenges posed by genotyping.

Results

Here we used 454 amplicon sequencing to characterize MHC class IIB exon 2 variation in the collared flycatcher, an important organism in evolutionary and immuno-ecological studies. On the basis of over 152,000 sequencing reads we identified 194 putative alleles in 237 individuals. We found an extreme complexity of the MHC class IIB in the collared flycatchers, with our estimates pointing to the presence of at least nine expressed loci and a large, though difficult to estimate precisely, number of pseudogene loci. Many similar alleles occurred in the pseudogenes indicating either a series of recent duplications or extensive concerted evolution. The expressed alleles showed unambiguous signals of historical selection and the occurrence of apparent interlocus exchange of alleles. Placing the collared flycatcher's MHC sequences in the context of passerine diversity revealed transspecific MHC class II evolution within the Muscicapidae family.

Conclusions

454 amplicon sequencing is an effective tool for advancing our understanding of the MHC class II structure and evolutionary patterns in Passeriformes. We found a highly dynamic pattern of evolution of MHC class IIB genes with strong signals of selection and pronounced sequence divergence in expressed genes, in contrast to the apparent sequence homogenization in pseudogenes. We show that next generation sequencing offers a universal, affordable method for the characterization and, in perspective, genotyping of MHC systems of virtually any complexity.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide-binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.

Gene duplication and fragmentation in the zebra finch major histocompatibility complex

BACKGROUND

Due to its high polymorphism and importance for disease resistance, the major histocompatibility complex (MHC) has been an important focus of many vertebrate genome projects. Avian MHC organization is of particular interest because the chicken Gallus gallus, the avian species with the best characterized MHC, possesses a highly streamlined minimal essential MHC, which is linked to resistance against specific pathogens. It remains unclear the extent to which this organization describes the situation in other birds and whether it represents a derived or ancestral condition. The sequencing of the zebra finch Taeniopygia guttata genome, in combination with targeted bacterial artificial chromosome (BAC) sequencing, has allowed us to characterize an MHC from a highly divergent and diverse avian lineage, the passerines.

RESULTS

The zebra finch MHC exhibits a complex structure and history involving gene duplication and fragmentation. The zebra finch MHC includes multiple Class I and Class II genes, some of which appear to be pseudogenes, and spans a much more extensive genomic region than the chicken MHC, as evidenced by the presence of MHC genes on each of seven BACs spanning 739 kb. Cytogenetic (FISH) evidence and the genome assembly itself place core MHC genes on as many as four chromosomes with TAP and Class I genes mapping to different chromosomes. MHC Class II regions are further characterized by high endogenous retroviral content. Lastly, we find strong evidence of selection acting on sites within passerine MHC Class I and Class II genes.

CONCLUSION

The zebra finch MHC differs markedly from that of the chicken, the only other bird species with a complete genome sequence. The apparent lack of synteny between TAP and the expressed MHC Class I locus is in fact reminiscent of a pattern seen in some mammalian lineages and may represent convergent evolution. Our analyses of the zebra finch MHC suggest a complex history involving chromosomal fission, gene duplication and translocation in the history of the MHC in birds, and highlight striking differences in MHC structure and organization among avian lineages.

Trans-species polymorphism of the Mhc class II DRB-like gene in banded penguins (genus Spheniscus)

The Major Histocompatibility Complex (Mhc) class II DRB locus of vertebrates is highly polymorphic and some alleles may be shared between closely related species as a result of balancing selection in association with resistance to parasites. In this study, we developed a new set of PCR primers to amplify, clone, and sequence overlapping portions of the Mhc class II DRB-like gene from the 5'UTR end to intron 3, including exons 1, 2, and 3 and introns 1 and 2 in four species (20 Humboldt, six African, five Magellanic, and three Galapagos penguins) of penguin from the genus Spheniscus (Sphe). Analysis of gene sequence variation by the neighbor-joining method of 21 Sphe sequences and 20 previously published sequences from four other penguin species revealed overlapping clades within the Sphe species, but species-specific clades for the other penguin species. The overlap of the DRB-like gene sequence variants between the four Sphe species suggests that, despite their allopatric distribution, the Sphe species are closely related and that some shared DRB1 alleles may have undergone a trans-species inheritance because of balancing selection and/or recent rapid speciation. The new primers and PCR assays that we have developed for the identification of the DRB1 DNA and protein sequence variations appear to be useful for the characterization of the molecular evolution of the gene in closely related Penguin species and might be helpful for the assessment of the genetic health and the management of the conservation and captivity of these endangered species.

Structure and evolution of a new avian MHC class II B gene in a sub-Antarctic seabird, the thin-billed prion (Procellariiformes: Pachyptila belcheri)

The major histocompatibility complex encodes molecules that present foreign peptides to T cells of the immune system. The peptide binding region (PBR) of these molecules is among the most polymorphic regions found in vertebrate taxa. Genomic cloning approaches are improving our understanding of the evolution of this multigene family in nonmodel avian groups. By building a cosmid library, a new MHC class II B gene, Pabe-DAB1, was isolated and characterized at the genomic level in a sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri). Pabe-DAB1 exhibits the hallmark structural features of functional MHC class II loci. Direct sequencing of the PBR encoding exon in a panel of prions revealed significantly higher levels of genetic diversity compared to two noncoding neutral loci, with most alleles differing by at least one replacement substitution in the peptide binding codons. We estimated evolutionary dynamics for Pabe-DAB1 using a variety of Bayesian and other approaches. Evidence for balancing selection comes from a spatially variable ratio of nonsynonymous-to-synonymous substitutions (mean d (N)/d (S) = 2.87) in the PBR, with sites predicted to be functionally relevant exhibiting the highest omega values. We estimate the population recombination rate to be approximately 0.3 per site per generation, indicating an important role for recombination in generating polymorphism at this locus. Pabe-DAB1 is among the few avian class II loci characterized at the genomic level and with a known intron-exon structure, a feature that greatly facilitated the amplification and sequencing of a single MHC locus in this species.

An MHC component to kin recognition and mate choice in birds: predictions, progress, and prospects

The major histocompatibility complex (Mhc) has been identified as a locus influencing disease resistance, mate choice, and kin recognition in mammals and fish. However, it is unclear whether the mechanisms by which Mhc genes influence behavior in mammals are applicable to other nonmammalian vertebrates such as birds. We review the biology of Mhc genes with particular reference to their relevance to avian mating and social systems. New genomics approaches recently have been applied to the Mhcs of chickens, quail, and several icons of avian behavioral ecology, including red-winged blackbirds (Agelaius phoeniceus) and house finches (Carpodacus mexicanus). The predominance of concerted evolution at avian Mhc loci makes such methods attractive for providing access to this complicated multigene family. Although some biological processes influenced by Mhc in mammals are physiologically implausible for birds, Mhc could influence cues that form well-known bases for mate choice in birds by influencing the health and vigor of individuals. The tight associations of Mhc variation and disease resistance in chickens raise hope that finding associations of Mhc genes, disease, and mate choice in natural populations of birds will be as fruitful as in mammalian systems.

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus)

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.

Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey

During the last decade, the major histocompatibility complex (MHC) has received much attention in the fields of evolutionary and conservation biology because of its potential implications in many biological processes. New insights into the gene structure and evolution of MHC genes can be gained through study of additional lineages of birds not yet investigated at the genomic level. In this study, we characterized MHC class II B genes in five families of birds of prey (Accipitridae, Pandionidae, Strigidae, Tytonidae, and Falconidae). Using PCR approaches, we isolated genomic MHC sequences up to 1300 bp spanning exons 1 to 3 in 26 representatives of each raptor lineage, finding no stop codons or frameshift mutations in any coding region. A survey of diversity across the entirety of exon 2 in the lesser kestrel Falco naumanni reported 26 alleles in 21 individuals. Bayesian analysis revealed 21 positively selected amino acid sites, which suggests that the MHC genes described here are functional and probably expressed. Finally, through interlocus comparisons and phylogenetic analysis, we also discuss genetic evidence for concerted and transspecies evolution in the raptor MHC.

Molecular evolutionary genomics of birds

Insight into the molecular evolution of birds has been offered by the steady accumulation of avian DNA sequence data, recently culminating in the first draft sequence of an avian genome, that of chicken. By studying avian molecular evolution we can learn about adaptations and phenotypic evolution in birds, and also gain an understanding of the similarities and differences between mammalian and avian genomes. In both these lineages, there is pronounced isochore structure with highly variable GC content. However, while mammalian isochores are decaying, they are maintained in the chicken lineage, which is consistent with a biased gene conversion model where the high and variable recombination rate of birds reinforces heterogeneity in GC. In Galliformes, GC is positively correlated with the rate of nucleotide substitution; the mean neutral mutation rate is 0.12-0.15% at each site per million years but this estimate comes with significant local variation in the rate of mutation. Comparative genomics reveals lower d(N)/d(S) ratios on micro- compared to macrochromosomes, possibly due to population genetic effects or a non-random distribution of genes with respect to chromosome size. A non-random genomic distribution is shown by genes with sex-biased expression, with male-biased genes over-represented and female-biased genes under-represented on the Z chromosome. A strong effect of selection is evident on the non-recombining W chromosome with high d(N)/d(S) ratios and limited polymorphism. Nucleotide diversity in chicken is estimated at 4-5 x 10(-3) which might be seen as surprisingly high given presumed bottlenecks during domestication, but is lower than that recently observed in several natural populations of other species. Several important aspects of the molecular evolutionary process of birds remain to be understood and it can be anticipated that the upcoming genome sequence of a second bird species, the zebra finch, as well as the integration of data on gene expression, shall further advance our knowledge of avian evolution.

Low MHC variation in the endangered Galápagos penguin (Spheniscus mendiculus)

The major histocompatibility complex (MHC) is one of the most polymorphic regions of the genome, likely due to balancing selection acting to maintain alleles over time. Lack of MHC variability has been attributed to factors such as genetic drift in small populations and relaxed selection pressure. The Galápagos penguin (Spheniscus mendiculus), endemic to the Galápagos Islands, is the only penguin that occurs on the equator. It relies upon cold, nutrient-rich upwellings and experiences severe population declines when ocean temperatures rise during El Niño events. These bottlenecks, occurring in an already small population, have likely resulted in reduced genetic diversity in this species. In this study, we used MHC class II exon 2 sequence data from a DRB1-like gene to characterize the amount of genetic variation at the MHC in 30 Galápagos penguins, as well as one Magellanic penguin (S. magellanicus) and two king penguins (Aptenodytes patagonicus), and compared it to that in five other penguin species for which published data exist. We found that the Galápagos penguin had the lowest MHC diversity (as measured by number of polymorphic sites and average divergence among alleles) of the eight penguin species studied. A phylogenetic analysis showed that Galápagos penguin MHC sequences are most closely related to Humboldt penguin (Spheniscus humboldti) sequences, its putative sister species based on other loci. An excess of non-synonymous mutations and a pattern of trans-specific evolution in the neighbor-joining tree suggest that selection is acting on the penguin MHC.

Evolutionary genetics of Carpodacus mexicanus, a recently colonized host of a bacterial pathogen, Mycoplasma gallisepticum

We present molecular data documenting how introduction to the eastern United States and an epizootic involving a bacterial pathogen has affected the genetic diversity of house finches, a cardueline songbird. Population bottlenecks during introduction can cause loss of genetic variation and may negatively affect a population's ability to adapt to novel stressors such as disease. Although a genome-wide survey using Amplified Fragment Length Polymorphism (AFLP) markers suggests little loss of genetic diversity in introduced populations, an epizootic of bacterial Mycoplasma has nonetheless caused dramatic declines in the eastern US population. Sequence analysis of a candidate gene for pathogen resistance in the Major Histocompatibity Complex (MHC) in pre- and post-epizootic population samples reveals allele frequency shifts since introduction of the pathogen, but similar shifts are also observed in control populations not exposed to the bacteria, and in a neutral non-coding locus. Expression studies using a novel subtractive hybridization approach indicate decreased expression of the class II MHC locus upon exposure to Mycoplasma, a pattern also seen in MHC class I loci in mice infected with cytomegalovirus and consistent with manipulation of the finch immune system by Mycoplasma. These results will be further expanded using experimental studies as well as examination of evolution of the pathogen genome itself.

Computing distribution of scale independent motifs in biological sequences

The use of Chaos Game Representation (CGR) or its generalization, Universal Sequence Maps (USM), to describe the distribution of biological sequences has been found objectionable because of the fractal structure of that coordinate system. Consequently, the investigation of distribution of symbolic motifs at multiple scales is hampered by an inexact association between distance and sequence dissimilarity. A solution to this problem could unleash the use of iterative maps as phase-state representation of sequences where its statistical properties can be conveniently investigated. In this study a family of kernel density functions is described that accommodates the fractal nature of iterative function representations of symbolic sequences and, consequently, enables the exact investigation of sequence motifs of arbitrary lengths in that scale-independent representation. Furthermore, the proposed kernel density includes both Markovian succession and currently used alignment-free sequence dissimilarity metrics as special solutions. Therefore, the fractal kernel described is in fact a generalization that provides a common framework for a diverse suite of sequence analysis techniques.

A comparison of variability and population structure for major histocompatibility complex and microsatellite loci in California coastal steelhead (Oncorhynchus mykiss Walbaum)

The major histocompatibility complex (MHC) contains genes integral to immune response in vertebrates. MHC genes have been shown to be under selection in a number of vertebrate taxa, making them intriguing for population genetic studies. We have conducted a survey of genetic variation in an MHC class II gene for steelhead trout from 24 sites in coastal California and compared this variation to that observed at 16 presumably neutral microsatellite loci. A high amount of allelic variation was observed at the MHC when compared to previously published studies on other Pacific salmonids. Elevated nonsynonymous substitutions, relative to synonymous substitutions, were detected at the MHC gene, indicating the signature of historical balancing selection. The MHC data were tested for correlations to and deviations from the patterns found with the microsatellite data. Estimates of allelic richness for the MHC gene and for the microsatellites were positively correlated, as were estimates of population differentiation (F(ST)). An analysis for F(ST) outliers indicates that the MHC locus has an elevated F(ST) relative to the neutral expectation, although a significant result was found for only one particular geographical subgroup. Relatively uniform allele frequency distributions were detected in four populations, although this finding may be partially due to recent population bottlenecks. These results indicate that, at the scale studied here, drift and migration play a major role in the observed geographical variability of MHC genes in steelhead, and that contemporary selection is relatively weak and difficult to detect.

Patterns of variation in MHC class II beta loci of the little greenbul (Andropadus virens) with comments on MHC evolution in birds

We have isolated major histocompatibility complex (MHC) class II beta loci from the little greenbul (Andropadus virens), an African songbird. We utilized preexisting information about conserved regions of the avian MHC to design primers to amplify a pool of sequences representing multiple loci. From this pool, a unique locus spanning 1109 bp that we designate as Anvi-DAB1 was cloned and sequenced. We designed locus-specific primers based on this sequence information and amplified six alleles from seven individuals. Compared to other A. virens MHC sequences obtained from genomic DNA or cDNA, the variability of sequences from Anvi-DAB1 was low and the ratio of nonsynonymous to synonymous substitution was much less than one, suggesting that Anvi-DAB1 may either be a pseudogene or a nonclassical MHC locus. Phylogenetic analysis revealed that the Anvi-DAB1 locus was highly divergent when compared with other passerine or A. virens genomic or transcribed MHC sequences. The use of conserved MHC primers followed by analysis of cloned sequences allows rapid isolation of MHC loci from exotic species and avoids laborious large-scale cloning and sequencing.

A cDNA macroarray approach to parasite-induced gene expression changes in a songbird host: genetic response of house finches to experimental infection by Mycoplasma gallisepticum

In 1994, the bacterial parasite Mycoplasma gallisepticum expanded its host range and swept through populations of a novel host--eastern US populations of the house finch (Carpodacus mexicanus). This epizootic caused a dramatic decline in finch population numbers, has been shown to have caused strong selection on house finch morphology, and presumably caused evolutionary change at the molecular level as finches evolved enhanced resistance. As a first step toward identifying finch genes that respond to infection by Mycoplasma and which may have experienced natural selection by this parasite, we used suppression subtractive hybridization (SSH) and cDNA macroarray approaches to identify differentially expressed genes regulated by the Mycoplasma parasite. Two subtractive cDNA libraries consisting of 16,512 clones were developed from spleen using an experimentally uninfected bird as the 'tester' and an infected bird as 'driver', and vice versa. Two hundred and twenty cDNA clones corresponding 34 genes with known vertebrate homologues and a large number of novel transcripts were found to be qualitatively up- or down-regulated genes by high-density filter hybridization. These gene expression changes were further confirmed by a high throughout reverse Northern blot approach and in specific cases by targeted Northern analysis. blast searches show that heat shock protein (HSP) 90, MHC II-associated invariant chain (CD74), T-cell immunoglobulin mucin 1 (TIM1), as well as numerous novel expressed genes not found in the databases were up- or down-regulated by the host in response to this parasite. Our results and macroarray resources provide a foundation for molecular co-evolutionary studies of the Mycoplasma parasite and its recently colonized avian host.

A comparison of variation between a MHC pseudogene and microsatellite loci of the little greenbul (Andropadus virens)

Background

We investigated genetic variation of a major histcompatibility complex (MHC) pseudogene (Anvi-DAB1) in the little greenbul (Andropadus virens) from four localities in Cameroon and one in Ivory Coast, West Africa. Previous microsatellite and mitochondrial DNA analyses had revealed little or no genetic differentiation among Cameroon localities but significant differentiation between localities in Cameroon and Ivory Coast.

Results

Levels of genetic variation, heterozygosity, and allelic diversity were high for the MHC pseudogene in Cameroon. Nucleotide diversity of the MHC pseudogene in Cameroon and Ivory Coast was comparable to levels observed in other avian species that have been studied for variation in nuclear genes. An excess of rare variants for the MHC pseudogene was found in the Cameroon population, but this excess was not statistically significant. Pairwise measures of population differentiation revealed high divergence between Cameroon and Ivory Coast for microsatellites and the MHC locus, although for the latter distance measures were much higher than the comparable microsatellite distances.

Conclusion

We provide the first ever comparison of variation in a putative MHC pseudogene to variation in neutral loci in a passerine bird. Our results are consistence with the action of neutral processes on the pseudogene and suggest they can provide an independent perspective on demographic history and population substructure.

Phylogenetics of modern birds in the era of genomics

In the 14 years since the first higher-level bird phylogenies based on DNA sequence data, avian phylogenetics has witnessed the advent and maturation of the genomics era, the completion of the chicken genome and a suite of technologies that promise to add considerably to the agenda of avian phylogenetics. In this review, we summarize current approaches and data characteristics of recent higher-level bird studies and suggest a number of as yet untested molecular and analytical approaches for the unfolding tree of life for birds. A variety of comparative genomics strategies, including adoption of objective quality scores for sequence data, analysis of contiguous DNA sequences provided by large-insert genomic libraries, and the systematic use of retroposon insertions and other rare genomic changes all promise an integrated phylogenetics that is solidly grounded in genome evolution. The avian genome is an excellent testing ground for such approaches because of the more balanced representation of single-copy and repetitive DNA regions than in mammals. Although comparative genomics has a number of obvious uses in avian phylogenetics, its application to large numbers of taxa poses a number of methodological and infrastructural challenges, and can be greatly facilitated by a 'community genomics' approach in which the modest sequencing throughputs of single PI laboratories are pooled to produce larger, complementary datasets. Although the polymerase chain reaction era of avian phylogenetics is far from complete, the comparative genomics era-with its ability to vastly increase the number and type of molecular characters and to provide a genomic context for these characters-will usher in a host of new perspectives and opportunities for integrating genome evolution and avian phylogenetics.

MHC class I typing in a songbird with numerous loci and high polymorphism using motif-specific PCR and DGGE

The major histocompatibility complex (MHC) has a central role in the specific immune defence of vertebrates. Exon 3 of MHC class I genes encodes the domain that binds and presents peptides from pathogens that trigger immune reactions. Here we develop a fast population screening method for detecting genetic variation in the MHC class I genes of birds. We found evidence of at least 15 exon 3 sequences in the investigated great reed warbler individual. The organisation of the great reed warbler MHC class I genes suggested that a locus-specific screening protocol is impractical due to the high similarity between alleles across loci, including the introns flanking exon 3. Therefore, we used motif-specific PCR to amplify two subsets of alleles (exon 3 sequences) that were separated with by DGGE. The motif-specific primers amplify a substantial proportion of the transcribed class I alleles (2-12 alleles per individual) from as many as six class I loci. Although not exhaustive, this gives a reliable estimate of the class I variation. The method is highly repeatable and more sensitive in detecting genetic variation than the RFLP method. The motif-specific primers also allow us to avoid screening pseudogenes. In our study population of great reed warblers, we found a high level of genetic variation in MHC class I, and no less than 234 DGGE genotypes were detected among 248 screened individuals.