Polymorphism across an exon-intron boundary in an avian Mhc class II B gene.

Twenty-three sequence haplotypes spanning the boundary of the second exon and intron of a red-winged blackbird Mhc class II B gene, Agph-DAB1, are presented. The polymorphism of the exon segment is distributed in two divergent allelic lineages which appear to be maintained by balancing selection. The silent nucleotide diversity of the exon (pi = 0.101) is more than five times that of the intron (pi = 0.018) and decays rapidly across the exon-intron boundary. Additionally, genealogical reconstruction indicates that divergence from a common ancestor in the exon sample is over four times that of the intron. The intron sequences reveal a pattern of polymorphism which is characteristic of directional selection, rather than a pattern expected from linkage to a balanced polymorphism. These results suggest that the evolutionary histories of these two adjacent regions have been disassociated by recombination or gene conversion. The estimated population recombination parameter between the exon and the intron is sufficiently high (4NeC = 8.545) to explain the homogenization of intron sequences. Compatibility analyses estimate that these events primarily occur from the exon-intron boundary to about 20-30 bases into the intron. Additionally, the observation that divergent exon alleles share identical intron sequence supports the conclusion of disassociation of exon and intron evolutionary histories by recombination.     

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

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.     

Evolution of the ABO blood group gene in Japanese macaque

We determined 5 sequences of Japanese macaque ABO blood group gene exon 7 (ca. 0.5 kb) and 2 sequences for exon 5 and intron 6 (ca. 1.7 kb). We compared those data with published sequences of other Old World monkey species, and the results suggest that alleles A and B were polymorphic in the ancestral species of macaques, and that B type allele evolved independently in macaque and baboon lineages.     

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     

Full-length sequence analysis of the HLA-DRB1 locus suggests a recent origin of alleles

The HLA region harbors some of the most polymorphic loci in the human genome. Among them is the class II locus HLA-DRB1, with more than 400 known alleles. The age of the polymorphism and the rate at which new alleles are generated at HLA loci has caused much controversy over the years. Previous studies have mostly been restricted to the 270 base pairs that constitute the second exon and represent the most variable part of the gene. Here, we investigate the evolutionary history of the HLA-DRB1 locus on the basis of an analysis of 15 genomic full-length alleles (10-15 kb). In addition, the variation in 49 complete coding sequences and 322 exon 2 sequences were analyzed. When excluding exon 2 from the analysis, the diversity at the synonymous sites was found to be similar to the intron diversity. The overall diversity in noncoding region was also similar to the genome average. The DRB1*03 lineage has been found in human, chimpanzee, bonobo, gorilla, and orangutan. An ancestral "proto HLA-DRB1*03 lineage" appeared to have diverged in the last 5 million years into the human-specific lineages *08, *11, *13, and *14. With exception to exon 2, both the coding- and the noncoding diversity suggests a recent origin (<1 million years ago) for most of the alleles at the HLA-DRB1 locus. Sites encoding for amino acids involved in antigen binding [antigen recognizing sites (ARS)] appear to have a more ancient origin. Taken together, the recent origin of most alleles, the high diversity between allelic lineages, and the ancient origin of sequence motifs in exon 2, is consistent with a relatively rapid generation of novel alleles by gene conversion like events.     

Evolution and trans-species polymorphism of MHC class IIbeta genes in cyprinid fish

The polymorphism of DAB genes encoding MHC IIbeta was investigated in 11 cyprinid species from central Europe. The species belonged to four subfamilies: Cyprininae, Tincinae, Gobioninae and Leuciscinae. Two paralogous groups of sequences, DAB1 and DAB3, were recognised according to the similarity of their nucleotide and amino-acid sequences and from phylogenetic analyses using either partial exon 2 or partial exon 3 sequences. A high allelic variability among species was found for exon 2, indicating extensive MHC polymorphism. Time divergence estimation supports the separation of DAB1 and DAB3 groups predating the separation into fish subfamilies, and a cyprinid origin of the DAB genes. Phylogenetic trees using exon 2 support the hypothesis of trans-species polymorphism, which appears to be limited to the subfamily level, i.e. the presence of sequences from different species in the same allelic group was more often recognised within subfamilies Cyprininae and Leuciscinae than between them. Phylogenetic trees using exon 3 reflect the phylogenetic patterns previously found for Cyprinidae systematics. Specific nucleotides and amino-acids in exon 3 that separate both subfamilies, as well as the species within the Cyprininae subfamily were observed. A lack of segregation in leuciscin species was recognised and the alleles of different leuciscin species tend to share similar motifs in exon 3. This could be explained by the ancient and complicated dispersion history of Cyprininae and the radiation of Leuciscinae. The effects of selective pressures were investigated: (1) within species, (2) among lineages, and (3) among sites. From intraspecific analyses, exon 2 sequences were identified as the targets of diversifying selection, whilst the evolution of exon 3 seems to be under the influence of purifying selection. The analyses among lineages indicate positive selection in many branches when using exon 2, therefore confirming trans-species polymorphism, whilst the DAB lineages of exon 3 are potentially submitted to purifying selection to some extent. Moreover, our results suggest the secondary acquisition of function of DAB1 group after duplication. The analyses among sites reveal that exon 2 exhibits sites under positive selection mostly corresponding to the putative PBR sites involved in the alpha-helix structure of the protein.     

Population genetics of speciation in nonmodel organisms: I. Ancestral polymorphism in mangroves

The level of DNA polymorphism in the ancestral species at the time of speciation can be estimated using DNA sequences from many loci sampled from 2 or more extant species. The comparison between ancestral and extant polymorphism can be informative about the population genetics of speciation. In this study, we collected and analyzed DNA sequences of approximately 60 genes from 4 species of Sonneratia, a common genus of mangroves on the Indo-Pacific coasts. We found that the 3 ancestral species were comparable to each other in terms of level of polymorphism. However, the ancestral species at the time of speciation were substantially more polymorphic than the extant geographical populations. This ancestral polymorphism is in fact larger than, or at least equal to, the level of polymorphism of the entire species across extant geographical populations. The observations are not fully compatible with speciation by strict allopatry. We suggest that, at the time of speciation, the ancestral species consisted of interconnected but strongly divided geographical populations. This population structure would give rise to high level of polymorphism across species range. This approach of studying the speciation history by genomic means should be applicable to nonmodel organisms.     

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     

Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon

Polymorphism of the nucleotide sequences encoding 149 amino acids of linked major histocompatibility complex (Mhc) class II 131 and 132 peptides, and of the intervening intron (548–773 base pairs), was examined within and among seven Pacific salmon (Oncorhynchus) species. Levels of nucleotide diversity were higher for theB1 sequence than forB2 or the intron in comparisons both within and between species. For the codons of the peptide binding region of the BI sequence, the level of nonsynonymous nucleotide substitution (d_N) exceeded the level of synonymous substitution (d_S) by a factor of ten for within-species comparisons, and by a factor of four for between-species comparisons. The excess of d_N indicates that balancing selection maintains diversity at this salmonidMhc class II locus, as is common forMhc loci in other vertebrates. Levels of nucleotide diversity for both the exon and intron sequences were greater among than within species, and there were numerous species-specific nucleotides present in both the coding and noncoding regions. Thus, neighbor-joining analysis of both the intron and exon regions provided phylogenies in which the sequences clustered strongly by species. There was little evidence of shared ancestral (trans-species) polymorphism in the exon phylogeny, and the intron phylogeny depicted standard relationships among the Pacific salmon species. The lack of shared allelicB1 lineages in these closely related species may result from severe bottlenecks that occurred during speciation or during the ice ages that glaciated the rim of the north Pacific Ocean approximately every 100 000 years in the Pleistocene.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers U34692-U34720     

Coevolution of genomic intron number and splice sites

Spliceosomal intron numbers and boundary sequences vary dramatically in eukaryotes. We found a striking correspondence between low intron number and strong sequence conservation of 5' splice sites (5'ss) across eukaryotic genomes. The phylogenetic pattern suggests that ancestral 5'ss were relatively weakly conserved, but that some lineages independently underwent both major intron loss and 5'ss strengthening. It seems that eukaryotic ancestors had relatively large intron numbers and 'weak' 5'ss, a pattern associated with frequent alternative splicing in modern organisms.     

Comparison of Bombyx mori and Helicoverpa armigera cytoplasmic actin genes provides clues to the evolution of actin genes in insects

The cytoplasmic actin genes BmA3 and BmA4 of Bombyx mori were found clustered in a single genomic clone in the same orientation. As a similar clustering of the two cytoplasmic actin genes Ha3a and Ha3b also occurs in another lepidopteran, Helicoverpa armigera, we analyzed the sequence of the pair of genes from each species. Due to the high conservation of cytoplasmic actins, the coding sequence of the four genes was easily aligned, allowing the detection of similarities in noncoding exon and intron sequences as well as in flanking sequences. All four genes exhibited a conserved intron inserted in codon 117, an original position not encountered in other species. It can thus be postulated that all of these genes derived from a common ancestral gene carrying this intron after a single event of insertion. The comparison of the four genes revealed that the genes of B. mori and H. armigera are related in two different ways: the coding sequence and the intron that interrupts it are more similar between paralogous genes within each species than between orthologous genes of the two species. In contrast, the other (noncoding) regions exhibited the greatest similarity between a gene of one species and a gene of the other species, defining two pairs of orthologous genes, BmA3 and HaA3a on one hand and BmA4 and HaA3b on the other. However, in each species, the very high similarities of the coding sequence and of the single intron that interrupts it strongly suggest that gene conversion events have homogenized this part of the sequence. As the divergence of the B. mori genes was higher than that of the H. armigera genes, we postulated that the gene conversion occurred earlier in the B. mori lineage. This leads us to hypothesize that gene conversion could also be responsible for the original transfer of the common intron to the second gene copy before the divergence of the B. mori and H. armigera lineages.     

Mitochondrial and chloroplast DNA-based phylogeny of Pelargonium (Geraniaceae)

Overall phylogenetic relationships within the genus Pelargonium (Geraniaceae) were inferred based on DNA sequences from mitochondrial(mt)-encoded nad1 b/c exons and from chloroplast(cp)-encoded trnL (UAA) 5' exon-trnF (GAA) exon regions using two species of Geranium and Sarcocaulon vanderetiae as outgroups. The group II intron between nad1 exons b and c was found to be absent from the Pelargonium, Geranium, and Sarcocaulon sequences presented here as well as from Erodium, which is the first recorded loss of this intron in angiosperms. Separate phylogenetic analyses of the mtDNA and cpDNA data sets produced largely congruent topologies, indicating linkage between mitochondrial and chloroplast genome inheritance. Simultaneous analysis of the combined data sets yielded a well-resolved topology with high clade support exhibiting a basic split into small and large chromosome species, the first group containing two lineages and the latter three. One large chromosome lineage (x = 11) comprises species from sections Myrrhidium and Chorisma and is sister to a lineage comprising P. mutans (x = 11) and species from section Jenkinsonia (x = 9). Sister to these two lineages is a lineage comprising species from sections Ciconium (x = 9) and Subsucculentia (x = 10). Cladistic evaluation of this pattern suggests that x = 11 is the ancestral basic chromosome number for the genus.     

Evolutionary origins of retroposon lineages of Mhc class II Ab alleles

Major histocompatibility complex (Mhc) class II Ab genes have evolved into three distinct lineages. While lineage 2 alleles differ from lineage 1 alleles by the insertion of a retroposon in intron 2, the basis for the extremely large intron 2 in lineage 3 alleles has heretofore been undetermined. In this report, we demonstrate by nucleotide sequencing that the genomic sequences of prototypic alleles from all three lineages diverge significantly and that lineage 3 is derived from lineage 2 by two insertional events in intron 2. One insert, composed of a member of B1 short interspersed repetitive elements (SINEs), occurs 508 base pairs (bp) 3 of exon 2, and the other, 1141 bp 3 of exon 2 within the retroposon that distinguishes lineage 2 from lineage 1. To assess the evolutionary stability of these lineages and the extent of ancestral polymorphisms of Ab within Mus species, we extended our restriction site polymorphism analysis to include 86 alleles from 120 independently derived H2 haplotypes from 12 separate species and subspecies of Mus. A phylogenetic tree revealing the relationships of these Ab alleles with respect to restriction site polymorphisms, but excluding the retroposon insertions, demonstrated that these lineages have distinctive genomic structures beyond the retroposon polymorphisms. In summary, these mouse Ab genes were produced from successive retroposon insertion events. Lineage 1 and 2 were detected in a variety of Mus species, including Mus caroli, indicating that these lineages diverged more than 2 million years ago. Lineage 3 alleles were found only in the Mus musculus subspecies, suggesting that it diverged from lineage 2 more recently. These results indicate that all three lineages of Ab have persisted through several speciation events in the genus Mus.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M60562     

Convergent evolution of major histocompatibility complex molecules in humans and New World monkeys

  In both Old World and New World monkeys Mhc-DRB sequences have been found which resemble human DRB1*03 and DRB3 genes in their second exon. The resemblance is shared sequence motifs and clustering of the genes or the encoded proteins in phylogenetic trees. This similarity could be due to common ancestry, convergence at the molecular level, or chance. To test which of these three explanations applies, we sequenced segments of New World monkey and macaque genes which encompass the entire second exon and large parts of both flanking introns. The test strongly supports the monophyly of New World monkey DRB intron sequences. The phylogenies of introns 1 and 2 from DRB1*03-like and DRB3-like genes are congruent, but both are incongruent with the exon 2-based phylogeny. The matching of intron 1- and intron 2-based phylogenies with each other suggests that reciprocal recombination has not played a major role in exon 2 evolution. Statistical comparisons of exon 2 from different DRB1*03 and DRB3 lineages indicate that it was neither gene conversion (descent), nor chance, but molecular convergence that has shaped their characteristic motifs. The demonstration of convergence in anthropoid Mhc-DRB genes has implications for the classification, age, and mechanism of generation of DRB allelic lineages. Received: 30 August 1999 / Revised: 19 October 1999     

Nucleotide sequence of the ribosomal RNA gene of Physarum polycephalum: intron 2 and its flanking regions of the 26S rRNA gene.

The 26S ribosomal RNA gene of Physarum polycephalum is interrupted by two introns, and we have previously determined the sequence of one of them (intron 1) (Nomiyama et al. Proc.Natl.Acad.Sci.USA 78, 1376-1380, 1981). In this study we sequenced the second intron (intron 2) of about 0.5 kb length and its flanking regions, and found that one nucleotide at each junction is identical in intron 1 and intron 2, though the junction regions share no other sequence homology. Comparison of the flanking exon sequences to E. coli 23S rRNA sequences shows that conserved sequences are interspersed with tracts having little homology. In particular, the region encompassing the intron 2 interruption site is highly conserved. The E. coli ribosomal protein L1 binding region is also conserved.     

Typing of artiodactyl MHC-DRB genes with the help of intronic simple repeated DNA sequences

An efficient oligonucleotide typing method for the highly polymorphic MHC-DRB genes is described for artiodactyls like cattle, sheep and goat. By means of the polymerase chain reaction, the second exon of MHC-DRB is amplified as well as part of the adjacent intron containing a mixed simple repeat sequence. Using this primer combination we were able to amplify the MHC-DRB exons 2 and adjacent introns from all of the investigated 10 species of the family of Bovidae and giraffes. Therefore, the DRB genes of novel artiodactyl species can also be readily studied. Oligonucleotide probes specific for the polymorphisms of ungulate DRB genes are used with which sequences differing in at least one single base can be distinguished. Exonic polymorphism was found to be correlated with the allele lengths and the patterns of the repeat structures. Hence oligonucleotide probes specific for different simple repeats and polymorphic positions serve also for typing across species barriers. The strict correlation of sequence length and exonic polymorphism permits a preselection of specific oligonucleotides for hybridization. Thus more than 20 alleles can already be differentiated from each of the three species.     

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.     

Distribution and evolution of introns in drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60–80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral. Received: 23 October 1995 / Accepted: 5 March 1996