Timing and effects of template number for gene conversion of major histocompatibility complex genes in the mouse

Profiling of gene duplication patterns of sequenced teleost genomes: evidence for rapid lineage-specific genome expansion mediated by recent tandem duplications

BACKGROUND

Gene duplication has had a major impact on genome evolution. Localized (or tandem) duplication resulting from unequal crossing over and whole genome duplication are believed to be the two dominant mechanisms contributing to vertebrate genome evolution. While much scrutiny has been directed toward discerning patterns indicative of whole-genome duplication events in teleost species, less attention has been paid to the continuous nature of gene duplications and their impact on the size, gene content, functional diversity, and overall architecture of teleost genomes.

RESULTS

Here, using a Markov clustering algorithm directed approach we catalogue and analyze patterns of gene duplication in the four model teleost species with chromosomal coordinates: zebrafish, medaka, stickleback, and Tetraodon. Our analyses based on set size, duplication type, synonymous substitution rate (Ks), and gene ontology emphasize shared and lineage-specific patterns of genome evolution via gene duplication. Most strikingly, our analyses highlight the extraordinary duplication and retention rate of recent duplicates in zebrafish and their likely role in the structural and functional expansion of the zebrafish genome. We find that the zebrafish genome is remarkable in its large number of duplicated genes, small duplicate set size, biased Ks distribution toward minimal mutational divergence, and proportion of tandem and intra-chromosomal duplicates when compared with the other teleost model genomes. The observed gene duplication patterns have played significant roles in shaping the architecture of teleost genomes and appear to have contributed to the recent functional diversification and divergence of important physiological processes in zebrafish.

CONCLUSIONS

We have analyzed gene duplication patterns and duplication types among the available teleost genomes and found that a large number of genes were tandemly and intrachromosomally duplicated, suggesting their origin of independent and continuous duplication. This is particularly true for the zebrafish genome. Further analysis of the duplicated gene sets indicated that a significant portion of duplicated genes in the zebrafish genome were of recent, lineage-specific duplication events. Most strikingly, a subset of duplicated genes is enriched among the recently duplicated genes involved in immune or sensory response pathways. Such findings demonstrated the significance of continuous gene duplication as well as that of whole genome duplication in the course of genome evolution.

Evidence that replication-associated mutation alone does not explain between-chromosome differences in substitution rates

Since Haldane first noticed an excess of paternally derived mutations, it has been considered that most mutations derive from errors during germ line replication. Miyata et al. (1987) proposed that differences in the rate of neutral evolution on X, Y, and autosome can be employed to measure the extent of this male bias. This commonly applied method assumes replication to be the sole source of between-chromosome variation in substitution rates. We propose a simple test of this assumption: If true, estimates of the male bias should be independent of which two chromosomal classes are compared. Prior evidence from rodents suggested that this might not be true, but conclusions were limited by a lack of rat Y-linked sequence. We therefore sequenced two rat Y-linked bacterial artificial chromosomes and determined evolutionary rate by comparison with mouse. For estimation of rates we consider both introns and synonymous rates. Surprisingly, for both data sets the prediction of congruent estimates of alpha is strongly rejected. Indeed, some comparisons suggest a female bias with autosomes evolving faster than Y-linked sequence. We conclude that the method of Miyata et al. (1987) has the potential to provide incorrect estimates. Correcting the method requires understanding of the other causes of substitution that might differ between chromosomal classes. One possible cause is recombination-associated substitution bias for which we find some evidence. We note that if, as some suggest, this association is dominantly owing to male recombination, the high estimates of alpha seen in birds is to be expected as Z chromosomes recombine in males.

Toward an evolutionary genomics of the avian Mhc

We review recent developments in the ongoing study of the evolution of the Mhc gene family in birds, with emphasis on class II B genes and results from songbirds obtained in our laboratory. Southern blots suggest a surprising diversity in Mhc class II gene number among various songbird species (Passeriformes). We have sequenced approximately 30 kb contigs from Mhc bearing cosmid clones from two species, red-winged blackbirds (Agelaius phoeniceus) and house finches (Carpodacus mexicanus), whose demography, lifetime reproductive success, epizootics, parasitology and mate choice are among the best studied for natural populations of birds. Of three genes cloned from these species, only one appears strongly polymorphic, and one (from the house finch) is likely a pseudogene. All are similar in structure to those in chickens, albeit with introns intermediate in length between chickens and mammals. Phylogenetic analysis of available class II B peptide-binding region exons suggests that the overwhelming long-term force operating on avian genes sampled thus far has been post-speciation gene duplication and/or concerted evolution. These and other results suggest that the evolution of class II B genes in birds conforms to a mixture of several models of multigene family evolution proposed for the mammalian Mhc, incorporating ongoing homogenization, duplication and pseudogene formation. Large-scale sequencing studies in these and other species, though still in their infancy, will prove invaluable for studying the comparative structures of avian Mhcs, as well as patterns of selection, mutation and linkage disequilibrium at several scales.

Effect of gene conversion on polymorphic patterns at major histocompatibility complex loci

The pattern of polymorphisms at MHC loci was studied by computer simulations and by DNA sequence analyses, with special reference to synonymous and non-synonymous divergence among alleles. A model of overdominance plus short-term selection was simulated for a gene family with allelic and non-allelic gene conversion. It was found that the ratio of non-synonymous-to-synonymous divergence among alleles decreases as the non-allelic conversion becomes more frequent. On the other hand, the ratio of synonymous divergence at the antigen recognition site to that at non-coding region increases by non-allelic conversion. By comparing such a pattern with the observed values of mammalian genes, it is suggested that non-allelic conversion is less frequent in human genes than in genes of other mammals. A quantitative analysis on the patchwork pattern was attempted by using identity excess among nucleotide sites of a gene. Comparison of simulation results with DNA data analysis shows that the rate of allelic conversion is high.

Gene organisation determines evolution of function in the chicken MHC

Some years ago, we used our data for class I genes, proteins and peptide-binding specificities to develop the hypothesis that the chicken B-F/B-L region represents a "minimal essential MHC". In this view, the B locus contains the classical (highly expressed and polymorphic) class I alpha and class II beta multigene families, which are reduced to one or two members, with many other genes moved away or deleted from the chicken genome altogether. We found that a single dominantly expressed class I gene determines the immune response to certain infectious pathogens, due to peptide-binding specificity and cell-surface expression level. This stands in stark contrast to well-studied mammals like humans and mice, in which every haplotype is more-or-less responsive to every pathogen and vaccine, presumably due to the multigene family of MHC molecules present. In order to approach the basis for a single dominantly expressed class I molecule, we have sequenced a portion of the B complex and examined the location and polymorphism of the class I (B-F) alpha, TAP and class II (B-L) beta genes. The region is remarkably compact and simple, with many of the genes expected from the MHC of mammals absent, including LMP, class II alpha and DO genes as well as most class III region genes. However, unexpected genes were present, including tapasin and putative natural killer receptor genes. The region is also organised differently from mammals, with the TAPs in between the class I genes, the tapasin gene in between the class II (B-L) beta genes, and the C4 gene outside of the class I alpha and class II beta genes. The close proximity of TAP and class I alpha genes leads to the possibility of co-evolution, which can drive the use of a single dominantly expressed class I molecule with peptide-binding specificity like the TAP molecule. There is also a single dominantly expressed class II beta gene, but the reason for this is not yet clear. Finally, the presence of the C4 gene outside of the classical class I alpha and class II beta genes suggests the possibility that this organisation was ancestral, although a number of models of organisation and evolution are still possible, given the presence of the Rfp-Y region with non-classical class I alpha and class II beta genes as well as the presence of multigene families of B-G and rRNA genes.