Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

Evolutionary lineages of RT1.Ba in the Australian Rattus

In this study, the evolutionary history of the variable second exon of RT1.Ba and its adjoining intron b are compared across a number of species and subspecies of the Australian RATTUS: Three lineages are identified in the second intron across a range of Rattus species. Two of these lineages, separated by the insertion of a probable rodent short interspersed nucleotide element and by point mutations outside the indel region, are both found in each of the major clades of the endemic Australian RATTUS: This pattern of ancestral polymorphism is reflected in the adjoining exon 2 sequences, although phylogenetic constraints confirm that the clustering is not identical to that of the associated intron sequences. In addition, the coding sequences show evidence of the retention of ancestral polymorphism, with identical exon sequences found in two divergent species, and some indication of gene conversion detected for the exon sequences.     

牙鲆MHC-DAA结构及其等位基因多态性

为了探讨牙鲆MHC-DAA等位基因的多态性, 根据牙鲆 (Paralichthys olivaceus) MHC-DAA的cDNA序列设计特异引物, 扩增内含子序列。牙鲆DAA基因由4个外显子和3个内含子组成, 与大西洋鲑 (Salmo salar) DAA基因结构相同, 在内含子2中存在一个高度多态的微卫星位点(GT)n。根据获得的MHC-DAA基因组序列设计特异性引物, 在45尾牙鲆中扩增了包括完整外显子2和内含子2的长度约670 bp的DNA片段, 克隆测序后共发现30个MHC-DAA等位基因, 各等位基因的频次及其各主型中亚型的数目都不均衡。在249 bp的外显子2序列中共有55个位点发生变异, 核苷酸多样性指数为0.0887, 编码的氨基酸序列中多态变异位点31个, 其中简约信息位点30个, 单变异位点1个。非同义替代与同义替代的比率在PBR(Peptide binding region)和非PBR结合区分别为3.30和2.43, 分析证明平衡选择是牙鲆存在众多DAA等位基因的发生机制。     

Positive selection and interallelic recombination at the merozoite surface antigen-1 (MSA-1) locus of Plasmodium falciparum.

DNA sequences of alleles at the merozoite surface antigen-1 (MSA-1) gene locus of the malaria parasite Plasmodium falciparum show evidence of repeated past recombination events between alleles. These include both (1) nonreciprocal recombination events that have homogenized certain gene regions among alleles and (2) reciprocal recombination events that have combined allelic segments with divergent evolutionary histories, thereby enhancing allelic diversity. In three different gene regions, the rate of nonsynonymous nucleotide substitution significantly exceeds that of synonymous nucleotide substitution, implying that positive Darwinian selection has acted to diversify alleles at the amino acid level. The MSA-1 polymorphism seems to be quite ancient; the two major allelic types have been maintained for approximately 35 Myr.     

Polymorphism and signatures of selection in the multimammate rat DQB gene

In order to test if DQB is a good candidate marker to investigate the relationship between major histocompatibility complex genes and pathogens in natural populations of Mastomys natalensis, we assessed the polymorphism and evolutionary history of this gene. Twenty-four individuals were genotyped for exon 2 of DQB using capillary electrophoresis single-strand conformation polymorphism, cloning, and sequencing. We found 21 different alleles. Four individuals show three alleles implying a duplication event in the history of this gene. Each distinct sequence translates to give a distinct amino acid sequence and there are strong signals of positive selection on peptide binding sites. Signals of recombination were found in the sequences suggesting that recombination has played a role in generating allelic diversity. Although trans-taxon polymorphism is present at the interspecific level in DQB exon 2 sequences of Mus species, we did not find any evidence of allele sharing among Muridae genera. This indicates a temporal limit of DQB allele sharing in Muridae of less than 8 Mya. In conclusion, although DQB seems to be a good marker to investigate pathogen-driven selection, the polymorphism of gene copy number may restrict its utility in natural populations.     

In vivo processing of an intron-containing archael tRNA

In vitro studies on the processing of halobacterial tRNA introns have led to the proposal that archaeal and eukaryotic tRNA intron endonucleases have distinctly different requirements for the recognition of pre-tRNAs. Using a Haloferax volcanii in vivo expression vector we have examined the in vivo processing of modified forms of the halobacterial intron-containing tRNA_Trp gene. As observed in vitro, changes in the exon–intron boundary structure of this pre-tRNA block processing. Intron sequences, other than those at the exon–intron boundaries, are not essential for processing in vivo. We also show that conversion of the tryptophan anticodon to an opal suppressor anticodon is tolerated when the exon-intron boundary structure is maintained.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

tRNA-isoleucine-tryptophan composite gene

Transfer-RNA genes in archaea often have introns intervening between exon sequences. The structural motif at the boundary between exon and intron is the bulge-helix-bulge. Computational investigations of these boundary structures in Haloarcula marismortui lead us to propose that tRNA-isoleucine and tRNA-tryptophan genes are co-located. Precise in silico identification of the splice-sites on the bulges at the exon-intron boundaries lead us to infer that a single intron-containing composite tRNA-gene can give rise to more than one gene product.     

Sequence variation in the bovine and ovine PRNP genes

A resequencing approach was adopted to identify sequence variants in the PRNP gene that may affect susceptibility or resistance to bovine spongiform encephalopathy. The entire PRNP gene (>21 kb) was sequenced from 26 chromosomes from a group of Holstein-Friesian cows, as well as exon 3 of PRNP (>4 kb) from a further 24 chromosomes from six diverse breeds. We identified 51 variant sequences of which 42 were single nucleotide polymorphisms and nine were insertion/deletion (indel) events. The study was extended to exon 3 of the sheep PRNP gene where 23 sequence variants were observed, four of which were indels. The level of nucleotide diversity in the complete bovine PRNP gene was pi = 0.00079, which is similar to that found at the bovine T-cell receptor alpha delta joining region (pi = 0.00077), but somewhat less than that observed for the bovine leptin (pi = 0.00265). Sequence variation within exon 3 of PRNP in both cattle (pi = 0.00102) and sheep (pi = 0.00171) was greater than that for the complete PRNP gene, with sheep showing greater sequence variation in exon 3 than cattle. The level of sequence variation reported here is greater than previously thought for the bovine PRNP gene in cattle. This study highlights the contribution that recombination plays in increasing allelic diversity in this species.     

Nested Evolution of a tRNALeu(UAA) Group I Intron by both Horizontal Intron Transfer and Recombination of the Entire tRNA Locus

The origin and evolution of bacterial introns are still controversial issues. Here we present data on the distribution and evolution of a recently discovered divergent tRNA(Leu)(UAA) intron. The intron shows a higher sequence affiliation with introns in tRNA(Ile)(CAU) and tRNA(Arg)(CCU) genes in alpha- and beta-proteobacteria, respectively, than with other cyanobacterial tRNA(Leu)(UAA) group I introns. The divergent tRNA(Leu)(UAA) intron is sporadically distributed both within the Nostoc and the Microcystis radiations. The complete tRNA gene, including flanking regions and intron from Microcystis aeruginosa strain NIVA-CYA 57, was sequenced in order to elucidate the evolutionary pattern of this intron. Phylogenetic reconstruction gave statistical evidence for different phylogenies for the intron and exon sequences, supporting an evolutionary model involving horizontal intron transfer. The distribution of the tRNA gene, its flanking regions, and the introns were addressed by Southern hybridization and PCR amplification. The tRNA gene, including the flanking regions, were absent in the intronless stains but present in the intron-containing strains. This suggests that the sporadic distribution of this intron within the Microcystis genus cannot be attributed to intron mobility but rather to an instability of the entire tRNA(Leu)(UAA) intron-containing genome region. Taken together, the complete data set for the evolution of this intron can best be explained by a model involving a nested evolution of the intron, i.e., wherein the intron has been transferred horizontally (probably through a single or a few events) to a tRNA(Leu)(UAA) gene which is located within a unstable genome region.     

Molecular variation of the shuttle craft and Lim3 genes, controlling the development of the nervous system, in a natural Drosophila melanogaster population

Comparative polymorphism of the first exon and first intron of the shuttle craft (stc) and Lim3 genes and their putative regulatory 5'-flanking sequences was analyzed using 20 sequenced natural alleles. A comparison of the stc and Lim3 genes showed that the extent of polymorphism was similar in their introns and corresponded to the variation level characteristic of Drosophila melanogaster, while the putative regulatory region and first intron of the stc gene proved to be more variable than the corresponding regions of the Lim3 gene. Since the genes under study occurred on the same chromosomes isolated from one population and were close together in a region having a high recombination rate, the difference in the extent of polymorphism between the regulatory and coding regions was explained by individual characteristics of each gene. The results made it possible to assume that the extent of polymorphism of the coding gene regions is maintained by balancing selection.     

Mhc class II B gene evolution in East African cichlid fishes

 A distinctive feature of essential major histocompatibility complex (Mhc) loci is their polymorphism characterized by large genetic distances between alleles and long persistence times of allelic lineages. Since the lineages often span several successive speciations, we investigated the behavior of the Mhc alleles during or close to the speciation phase. We sequenced exon 2 of the class II B locus 4 from 232 East African cichlid fishes representing 32 related species. The divergence times of the (sub)species ranged from 6000 to 8.4 million years. Two types of evolutionary analysis were used to elucidate the pattern of exon 2 sequence divergence. First, phylogenetic methods were applied to reconstruct the most likely evolutionary pathways leading from the last common ancestor of the set to the extant sequences, and to assess the probable mechanisms involved in allelic diversification. Second, pairwise comparisons of sequences were carried out to detect differences seemingly incompatible with origin by nonparallel point mutations. The analysis revealed point mutations to be the most important mechanism behind allelic divergences, with recombination playing only an auxiliary part. Comparison of sequences from related species revealed evidence of random allelic (lineage) losses apparently associated with speciation. Sharing of identical alleles could be demonstrated between species that diverged 2 million years ago. The phylogeny of the exon was incongruent with that of the flanking introns, indicating either a high degree of convergent evolution at the peptide-binding region-encoding sites, or intron homogenization. Received: 6 December 1999 / Revised: 15 February 2000     

The complex history of a gene proposed to participate in a sexual isolation mechanism in house mice

Previous behavioral experiments showed that mouse salivary androgen-binding protein (ABP) was involved in interindividual recognition and might play a role in sexual isolation between house mouse (Mus musculus) subspecies. The pattern of evolution of Abpa, the gene for the alpha subunit of ABP, was found to be consistent with this hypothesis. Abpa apparently diverged rapidly between species and subspecies with a large excess of nonsynonymous substitutions, a lack of exon polymorphism within each of the three subspecies, and a lack of intron polymorphism in the one subspecies studied (M. musculus domesticus). Here we characterized the intron and exon sequence variations of this gene in house mouse populations from central Eurasia, a region yet unsampled and thought to be close to the cradle of the radiation of the subspecies. We also determined the intron and exon sequences in seven other species of the genus Mus. We confirmed the general pattern of rapid evolution by essentially nonsynonymous substitutions, both inter- and intraspecifically, supporting the idea that Darwinian selection has driven the evolution of this gene. We also observed a uniform intron sequence in five samples of M. musculus musculus, suggesting that a selective sweep might have occurred for that allele. In contrast to previous results, however, we found extensive intron and exon polymorphism in some house mouse populations from central Eurasia. We also found evidence for secondary admixture of the subspecies-specific alleles in regions of transition between the subspecies in central Eurasia. Furthermore, an abnormal intron phylogeny suggested that interspecific exchanges had occurred between the house mouse subspecies and three other Palearctic species. These observations appear to be at variance with the simple hypothesis that Abpa is involved in reproductive isolation. Although we do not rule out a role in recognition, the situation appears to be more complex than previously thought. Thus the selective mechanism behind the evolution of Abpa remains to be resolved, and we suggest that it may have changed during the recent colonization history of the house mouse.     

Recombination and gene conversion-like events may contribute to ABO gene diversity causing various phenotypes

We identified five different alleles, tentatively named ABO*O301, *0302, *R102, *R103, and *A110, in Japanese individuals possessing the blood group O phenotype. These alleles lack the guanine deletion at nucleotide position 261 which is shared by a majority of O alleles. Nucleotide sequence analysis revealed that *0301 and *0302 had single nonsynonymous substitutions compared with *A101 or *A102 responsible for the A1 phenotype. Analysis of intron 6 at the ABO gene by polymerase chain reaction-single-strand conformation polymorphism and direct sequencing revealed that *R102 and *R103 had chimeric sequences of A-02 and B-02, respectively, from exons 6 to 7. In the analysis of five other chimeric alleles detected in the same manner, we identified a total of four different recombination-breakpoints within or near intron 6. When 510 unrelated Japanese were examined, the frequency of the chimeric alleles generated by recombination in intron 6 or exon 7 was estimated to be 1.7%. In addition, we found that *O301, *A110, *C101, *A111, and 35% of *A102 had a unique A-B-A chimeric sequence at intron 6, presumed to originate from a gene conversion-like event. We had previously established that *A110 also had an A-O2-A chimeric sequence around nucleotide position 646 in exon 7. Thus this allele has an A-B-A-O2-A chimeric sequence from intron 6 to exon 7 probably generated by two different gene conversions. Similar patchwork sequences around nucleotide position 646 in exon 7 were observed in two other new alleles responsible for the Ax and B3 phenotypes. Thus, the site is presumably a hotspot for gene conversion. These results indicate that both recombination and gene conversion-like events play important roles in generating ABO gene diversity.     

Evolution of introns and exons of class II major histocompatibility complex genes of vertebrates

 The phylogenetic relationships and patterns of nucleotide substitution were compared for introns and exons of class II major histocompatibility complex (MHC) genes in three datasets: human DRB1, human DQA1, and cyprinid fish DAB1. In both human DRB1 and cyprinid DAB1, there was strong evidence that recombination events between alleles have occurred in such a way that intron and exon sequences of a given allele do not necessarily share the same evolutionary history. In the case of human DRB1, recombination was found to have homogenized intron 1 and intron 2 sequences relative to exon 2 sequences within lineages of alleles but not between lineages. As a result, mean divergence times of intron sequences are much more recent than those of exonic sequences. Thus, the divergence time of DRB1 introns cannot be used to date that of exons in the same alleles, and the hypothesis that most human DRB1 polymorphism is of very recent origin is not supported. Received: 5 September 1999 / Revised: 30 December 1999     

Pronounced differences of recombination activity at the sex determination locus of the honeybee, a locus under strong balancing selection

Recombination decreases the association of linked nucleotide sites and can influence levels of polymorphism in natural populations. When coupled with selection, recombination may relax potential conflict among linked genes, a concept that has played a central role in research on the evolution of recombination. The sex determination locus (SDL) of the honeybee is an informative example for exploring the combined forces of recombination, selection, and linkage on sequence evolution. Balancing selection at SDL is very strong and homozygous individuals at SDL are eliminated by worker bees. The recombination rate is increased up to four times that of the genomewide average in the region surrounding SDL. Analysis of nucleotide diversity (pi) reveals a sevenfold increase of polymorphism within the sex determination gene complementary sex determiner (csd) that rapidly declines within 45 kb to levels of genomewide estimates. Although no recombination was observed within SDL, which contains csd, analyses of heterogeneity, shared polymorphic sites, and linkage disequilibrium (LD) show that recombination has contributed to the evolution of the 5' part of some csd sequences. Gene conversion, however, has not obviously contributed to the evolution of csd sequences. The local control of recombination appears to be related to SDL function and mode of selection. The homogenizing force of recombination is reduced within SDL, which preserves allelic differences and specificity, while the increase of recombination activity around SDL relaxes conflict between SDL and linked genes.     

Genealogy-Based Methods for Inference of Historical Recombination and Gene Flow and Their Application in Saccharomyces cerevisiae

Genetic exchange between isolated populations, or introgression between species, serves as a key source of novel genetic material on which natural selection can act. While detecting historical gene flow from DNA sequence data is of much interest, many existing methods can be limited by requirements for deep population genomic sampling. In this paper, we develop a scalable genealogy-based method to detect candidate signatures of gene flow into a given population when the source of the alleles is unknown. Our method does not require sequenced samples from the source population, provided that the alleles have not reached fixation in the sampled recipient population. The method utilizes recent advances in algorithms for the efficient reconstruction of ancestral recombination graphs, which encode genealogical histories of DNA sequence data at each site, and is capable of detecting the signatures of gene flow whose footprints are of length up to single genes. Further, we employ a theoretical framework based on coalescent theory to test for statistical significance of certain recombination patterns consistent with gene flow from divergent sources. Implementing these methods for application to whole-genome sequences of environmental yeast isolates, we illustrate the power of our approach to highlight loci with unusual recombination histories. By developing innovative theory and methods to analyze signatures of gene flow from population sequence data, our work establishes a foundation for the continued study of introgression and its evolutionary relevance.     

Different splice site utilization generates diversity between the rat and human pancreatic polypeptide precursors

Pancreatic polypeptide is derived from a polyprotein precursor molecule. Although the amino acid sequences specifying the signal peptide and pancreatic polypeptide are well conserved between the rat and the human, the carboxy-terminal amino acid sequences of the precursors are highly divergent. To better understand the molecular basis of the divergence between the rat and human C-terminal peptides, we have determined the nucleotide sequence of the rat pancreatic polypeptide gene. A comparison between the primary structures of the rat and human genes reveals that the heterogeneity of the C-terminal peptides can be explained in large part by a frameshift mutation and the utilization of an alternative splice donor site in the third exon of the rat gene. As a consequence of the displaced splice site, part of the third exon of the rate gene is homologous to the sequence in the third intron of the human gene. Our results suggest that the rat and human pancreatic polypeptide genes arose from a common ancestral gene, and that differences in the C-terminal domains of the precursor reflect less strict evolutionary constraints than those imposed upon the amino-terminal domains of the precursor.