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.     

Evolutionary origin and diversification of the mammalian CD1 antigen genes.

CD1 antigens are cell-surface glycoproteins which have a molecular structure which is similar (consisting of extracellular domains alpha 1, alpha 2, and alpha 3, a transmembrane portion, and a cytoplasmic tail) to that of class I MHC molecules. Phylogenetic analysis of mammalian CD1 DNA sequences revealed that these genes are more closely related to the class I major histocompatibility complex (MHC) than to the class II MHC and that mammalian genes are more closely related to avian class I MHC genes than they are to mammalian class I MHC genes. The CD1 genes form a multigene family with different numbers of genes in different species (five in human, eight in rabbit, and two in mouse). Known CD1 genes are grouped into the following three families, on the basis of evolutionary relationship: (1) the human HCD1B gene and a partial sequence from the domestic rabbit, (2) the human HCD1A and HCD1C genes, and (3) the human HCD1D and HCD1E genes plus the two mouse genes and a sequence from the cottontail rabbit. The alpha 1 and alpha 2 domains of CD1 are much less conserved at the amino acid level than are the corresponding domains of class I MHC molecules, but the alpha 3 domain of CD1 seems to be still more conserved than the well-conserved alpha 3 domain of class I MHC molecules. Furthermore, in the human CD1 gene family, interlocus exon exchange has homogenized alpha 3 domains of all CD1 genes except HCD1C.     

Isolation and characterization of major histocompatibility complex (MHC) class II B genes in the Barn owl (Aves: Tyto alba)

We isolated major histocompatibility complex class II B (MHCIIB) genes in the Barn owl (Tyto alba). A PCR-based approach combined with primer walking on genomic and complementary DNA as well as Southern blot analyses revealed the presence of two MHCIIB genes, both being expressed in spleen, liver, and blood. Characteristic structural features of MHCIIB genes as well as their expression and high non-synonymous substitution rates in the region involved in antigen binding suggest that both genes are functional. MHC organization in the Barn owl is simple compared to passerine species that show multiple duplications, and resembles the minimal essential MHC of chicken.     

Critical review of NGS analyses for de novo genotyping multigene families

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next‐generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra‐deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the ‘degree of change’ (DOC) method avoids using arbitrary cut‐off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third‐generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second‐generation sequencing approaches.     

Sequence of a Complete Chicken BG Haplotype Shows Dynamic Expansion and Contraction of Two Gene Lineages with Particular Expression Patterns

Many genes important in immunity are found as multigene families. The butyrophilin genes are members of the B7 family, playing diverse roles in co-regulation and perhaps in antigen presentation. In humans, a fixed number of butyrophilin genes are found in and around the major histocompatibility complex (MHC), and show striking association with particular autoimmune diseases. In chickens, BG genes encode homologues with somewhat different domain organisation. Only a few BG genes have been characterised, one involved in actin-myosin interaction in the intestinal brush border, and another implicated in resistance to viral diseases. We characterise all BG genes in B12 chickens, finding a multigene family organised as tandem repeats in the BG region outside the MHC, a single gene in the MHC (the BF-BL region), and another single gene on a different chromosome. There is a precise cell and tissue expression for each gene, but overall there are two kinds, those expressed by haemopoietic cells and those expressed in tissues (presumably non-haemopoietic cells), correlating with two different kinds of promoters and 5′ untranslated regions (5′UTR). However, the multigene family in the BG region contains many hybrid genes, suggesting recombination and/or deletion as major evolutionary forces. We identify BG genes in the chicken whole genome shotgun sequence, as well as by comparison to other haplotypes by fibre fluorescence in situ hybridisation, confirming dynamic expansion and contraction within the BG region. Thus, the BG genes in chickens are undergoing much more rapid evolution compared to their homologues in mammals, for reasons yet to be understood.     

Genomic analysis of human multigene families using chromosome-specific vectorette PCR.

We report a technique for the rapid determination of genomic structure of individual members of human interspersed multigene families which circumvents the requirement for genomic clone isolation. In this approach, vectorette libraries were constructed from human/rodent somatic cell hybrid DNA harbouring single members of the gene family. Using these libraries as PCR templates with nested gene-specific primers in combination with a common vectorette primer resulted in the amplification of gene-specific products suitable for the subsequent determination of intron/exon structure. We have applied this technique to characterise members of two gene families.     

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host–pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1′665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen‐mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host–pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.     

Birth-and-death evolution of the internalin multigene family in Listeria

The birth-and-death model of multigene family evolution describes patterns of gene origination, diversification and loss within multigene families. Since it was first developed in the 1990s, the model has been found to characterize a large number of eukaryotic multigene families. In this paper, we report for the first time a bacterial multigene family that undergoes birth-and-death evolution. By analyzing the evolutionary relationships among internalins, a relatively large and diverse family of genes associated with key virulence functions in Listeria, we demonstrate the importance of birth-and-death evolution in the diversification of this important bacterial pathogen. We also detected two instances of lateral gene transfer within the internalins, but the estimated frequency would have been much higher had it not been analyzed within the context of birth-and-death evolutionary dynamics and a phenomenon that we term "paralog-sorting", which involves the unequal transmittal of gene duplicates during or subsequent to the speciation process. As such, in addition to providing the first demonstration of birth-and-death evolution within a bacterial multigene family, our results indicate that the extent of lateral transfer in bacterial multigene families should be re-examined in the light of birth-and-death evolution.     

Gene conversions within the skeletal myosin multigene family.

Comparisons of the nucleotide sequences of the light meromyosin (LMM) region of developmentally regulated fast chicken myosin heavy chain (MHC) isoforms indicates that chicken MHC isoforms are more similar to each other than to MHC isoforms in other species. The sequence data provide evidence that gene conversion events have occurred recently among the isoforms. An embryonic (Cemb1) isoform and neonatal isoform have the most extensive regions of sequence identity. Similar gene conversion events are present in the rat alpha- and beta-cardiac MHCs, but were not obvious in the LMM of developmentally regulated fast human MHC isoforms. The data suggest that gene conversion events can play a significant role in the evolution of the MHC multigene families and that concerted evolution of the chicken multigene family occurred after the divergence of mammals and avians.     

Genetic Variation of the Major Histocompatibility Complex (MHC Class II B Gene) in the Threatened Hume’s Pheasant,Syrmaticus humiae

Major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates and encode molecules that play a crucial role in pathogen resistance. As a result of their diversity, they have received much attention in the fields of evolutionary and conservation biology. Here, we described the genetic variation of MHC class II B (MHCIIB) exon 2 in a wild population of Hume’s pheasant (Syrmaticus humiae), which has suffered a dramatic decline in population over the last three decades across its ranges in the face of heavy exploitation and habitat loss. Twenty-four distinct alleles were found in 73 S. humiae specimens. We found seven shared alleles among four geographical groups as well as six rare MHCIIB alleles. Most individuals displayed between one to five alleles, suggesting that there are at least three MHCIIB loci of the Hume’s pheasant. The d _N ⁄ d _S ratio at putative antigen-binding sites (ABS) was significantly greater than one, indicating balancing selection is acting on MHCIIB exon 2. Additionally, recombination and gene conversion contributed to generating MHCIIB diversity in the Hume’s pheasant. One to three recombination events and seventy-five significant gene conversion events were observed within the Hume’s pheasant MHCIIB loci. The phylogenetic tree and network analysis revealed that the Hume’s pheasant alleles do not cluster together, but are scattered through the tree or network indicating a trans-species evolutionary mode. These findings revealed the evolution of the Hume’s pheasant MHC after suffering extreme habitat fragmentation.     

Theoretical population genetics of repeated genes forming a multigene family

The evolution of repeated genes forming a multigene family in a finite population is studied with special reference to the probability of gene identity, i.e., the identity probability of two gene units chosen from the gene family. This quantity is called clonality and is defined as the sum of squares of the frequencies of gene lineages in the family. The multigene family is undergoing continuous unequal somatic crossing over, ordinary interchromosomal crossing over, mutation and random frequency drift. Two measures of clonality are used: clonality within one chromosome and that between two different chromosomes. The equilibrium properties of the means, the variances and the covariance of the two measures of clonality are investigated by using the diffusion equation method under the assumption of constant number of gene units in the multigene family. Some models of natural selection based on clonality are considered. The possible significance of the variance and covariance of clonality among the chromosomes on the adaptive differentiation of gene families such as those producing antibodies is discussed.     

Blood parasites shape extreme major histocompatibility complex diversity in a migratory passerine

Pathogens are one of the main forces driving the evolution and maintenance of the highly polymorphic genes of the vertebrate major histocompatibility complex (MHC). Although MHC proteins are crucial in pathogen recognition, it is still poorly understood how pathogen‐mediated selection promotes and maintains MHC diversity, and especially so in host species with highly duplicated MHC genes. Sedge warblers (Acrocephalus schoenobaenus) have highly duplicated MHC genes, and using data from high‐throughput MHC genotyping, we were able to investigate to what extent avian malaria parasites explain temporal MHC class I supertype fluctuations in a long‐term study population. We investigated infection status and infection intensities of two different strains of Haemoproteus, that is avian malaria parasites that are known to have significant fitness consequences in sedge warblers. We found that prevalence of avian malaria in carriers of specific MHC class I supertypes was a significant predictor of their frequency changes between years. This finding suggests that avian malaria infections partly drive the temporal fluctuations of the MHC class I supertypes. Furthermore, we found that individuals with a large number of different supertypes had higher resistance to avian malaria, but there was no evidence for an optimal MHC class I diversity. Thus, the two studied malaria parasite strains appear to select for a high MHC class I supertype diversity. Such selection may explain the maintenance of the extremely high number of MHC class I gene copies in sedge warblers and possibly also in other passerines where avian malaria is a common disease.     

Perils of paralogy: using HSP70 genes for inferring organismal phylogenies

Conserved genes have found their way into the mainstream of molecular systematics. Many of these genes are members of multigene families. A difficulty with using single genes of multigene families for phylogenetic inference is that genes from one species may be paralogous to those from another taxon. We focus attention on this problem using heat shock 70 (HSP70) genes. Using polymerase chain reaction techniques with genomic DNA, we isolated and sequenced 123 distinct sequences from 12 species of sharks. Phylogenetic analysis indicated that the sequences cluster with constituitively expressed cytoplasmic heat shock-like genes. Three highly divergent gene clades were sampled. A number of similar sequences were sampled from each species within each distinct gene clade. Comparison of published species trees with an HSP70 gene tree inferred using Bayesian phylogenetic analysis revealed several cases of gene duplication and differential sorting of gene lineages within this group of sharks. Gene tree parsimony based on the objective criteria of duplication and losses showed that previously published hypotheses of species relationships and two novel hypothesis based on Bayesian phylogenetics were concordant with the history of HSP70 gene duplication and loss. By contrast, two published hypotheses based on morphological data were not significantly different from the null hypothesis of a random association between species relatedness and the HSP70 gene tree. These results suggest that gene tree parsimony using data from multigene families can be used for inferring species relationships or testing published alternative hypotheses. More importantly, the results suggest that systematic studies relying on phylogenetic inferences from HSP70 genes may by plagued by unrecognized paralogy of sampled genes. Our results underscore the distinction between gene and species trees and highlight an underappreciated source of discordance between gene trees and organismal phylogeny, i.e., unrecognized paralogy of sampled genes.     

The Cohesive Population Genetics of Molecular Drive

The long-term population genetics of multigene families is influenced by several biased and unbiased mechanisms of nonreciprocal exchanges (gene conversion, unequal exchanges, transposition) between member genes, often distributed on several chromosomes. These mechanisms cause fluctuations in the copy number of variant genes in an individual and lead to a gradual replacement of an original family of n genes (A) in N number of individuals by a variant gene (a). The process for spreading a variant gene through a family and through a population is called molecular drive. Consideration of the known slow rates of nonreciprocal exchanges predicts that the population variance in the copy number of gene a per individual is small at any given generation during molecular drive. Genotypes at a given generation are expected only to range over a small section of all possible genotypes from one extreme (n number of A) to the other (n number of a). A theory is developed for estimating the size of the population variance by using the concept of identity coefficients. In particular, the variance in the course of spreading of a single mutant gene of a multigene family was investigated in detail, and the theory of identity coefficients at the state of steady decay of genetic variability proved to be useful. Monte Carlo simulations and numerical analysis based on realistic rates of exchange in families of known size reveal the correctness of the theoretical prediction and also assess the effect of bias in turnover. The population dynamics of molecular drive in gradually increasing the mean copy number of a variant gene without the generation of a large variance (population cohesion) is of significance regarding potential interactions between natural selection and molecular drive.     

Characterization of cDNA and genomic sequences corresponding to an embryonic myosin heavy chain

We report here the isolation and characterization of cDNA and genomic sequences corresponding to a rat embryonic myosin heavy chain (MHC) protein. This gene, which is present as a single copy in the rat genome, comprises about 25 kilobase pairs of DNA and contains approximately 80% intronic sequences. The embryonic MHC gene belongs to a highly conserved multigene family, and exhibits a high degree of nucleotide and amino acid sequence conservation with other sarcomeric MHC genes from nematode to man. S1 nuclease mapping experiments using cDNA and genomic probes show that this MHC gene is transiently expressed during skeletal muscle development. Its mRNA is detected in fetal skeletal muscle during early development and persists up to 2 weeks after birth with the overlapping expression of neonatal and adult skeletal MHC mRNAs. However, this MHC is not expressed in the adult skeletal muscle with the exception of extraocular muscle fibers. The transient expression during muscle development of the isoform produced by this gene and its sequential replacement by other MHCs raises interesting questions about the mechanism controlling MHC isozyme transitions and the physiological significance of the individual MHCs in muscle fibers.     

Using Illumina next generation sequencing technologies to sequence multigene families in de novo species

The advent of Next Generation Sequencing Technology (NGST) has revolutionized molecular biology research, allowing for rapid gene/genome sequencing from a multitude of diverse species. As high throughput sequencing becomes more accessible, more efficient workflows must be developed to deal with the amounts of data produced and better assemble the genomes of de novo lineages. We combine traditional laboratory methods with Illumina NGST to amplify and sequence the largest mammalian multigene family, the Olfactory Receptor gene family, for species with and without a reference genome. We develop novel assembly methods to annotate and filter these data, which can be utilized for any gene family or any species. We find no significant difference between the ratio of genes within their respective gene families of our data compared with available genomic data. Using simulated data we explore the limitations of short‐read sequence data and our assembly in recovering this gene family. We highlight the benefits and shortcomings of these methods. Compared with data generated from traditional polymerase chain reaction, cloning and Sanger sequencing methodologies, sequence data generated using our pipeline increases yield and sequencing efficiency without reducing the number of unique genes amplified. A cloning step is not required, therefore shortening data generation time. The novel downstream methodologies and workflows described provide a tool to be utilized by many fields of biology, to access and analyze the vast quantities of data generated. By combining laboratory and in silico methods, we provide a means of extracting genomic information for multigene families without complete genome sequencing.     

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.     

Poxvirus genome evolution by gene gain and loss

The poxviruses (Poxviridae) are a family of viruses with double-stranded DNA genomes and substantial numbers (often >200) of genes per genome. We studied the patterns of gene gain and loss over the evolutionary history of 17 poxvirus complete genomes. A phylogeny based on gene family presence/absence showed good agreement with families based on concatenated amino acid sequences of conserved single-copy genes. Gene duplications in poxviruses were often lineage specific, and the most extensively duplicated viral gene families were found in only a few of the genomes analyzed. A total of 34 gene families were found to include a member in at least one of the poxvirus genomes analyzed and at least one animal genome; in 16 (47%) of these families, there was evidence of recent horizontal gene transfer (HGT) from host to virus. Gene families with evidence of HGT included several involved in host immune defense mechanisms (the MHC class I, interleukin-10, interleukin-24, interleukin-18, the interferon gamma receptor, and tumor necrosis factor receptor II) and others (glutaredoxin and glutathione peroxidase) involved in resistance of cells to oxidative stress. Thus "capture" of host genes by HGT has been a recurrent feature of poxvirus evolution and has played an important role in adapting the virus to survive host antiviral defense mechanisms.     

Complex evolution of vitellogenin genes in salmonid fishes

Vitellogenins (Vtg) are usually encoded by small multigene families containing up to six genes. With 20 tandemly arranged genes, the rainbow trout ( Oncorhynchus mykiss) is an exception to this rule. PCR amplification, cloning and sequence analysis of Vtg genes in other salmonid species revealed the existence of two paralogous gene clusters, designated Vtg-A and Vtg-B. Southern hybridization showed that the number of genes varies from 2 to 30 copies from one species to another, as well as between the two gene clusters. All Coregonus, Thymallus, Salmo and Salvelinus species studied have both gene clusters, while Oncorhynchus species possess only the Vtg-A locus. Phylogenetic trees constructed from Vtg sequences revealed conflicting nodes with the consensus tree based on morphological and anatomical data. Vtg sequences support the grouping ( Salmo, ( Salvelinus, Oncorhynchus)) instead of the accepted consensus ( Salvelinus, ( Salmo, Oncorhynchus)). Structural data on gene organization also support the contention that Salvelinus and Oncorhynchus are sister taxa. Evolutionary implications for the Vtg gene clusters in salmonids are discussed.