Description of three Pristionchus species (Nematoda: Diplogastridae) from Japan that form a cryptic species complex with the model organism P. pacificus

Three new species of Pristionchus (P. exspectatus, P. arcanus, and P. japonicus) are described from Japan. They are morphologically similar, with P. exspectatus and P. arcanus being almost indistinguishable from the model organism P. pacificus. Reproductive isolation, namely the inability to produce interfertile F1 hybrids, separates all species pairs in the species complex. Additionally, all three new species are distinguished from P. pacificus Sommer, Carta, Kim, and Sternberg, 1996 by having a gonochoristic instead of hermaphroditic mode of reproduction. In addition to its reproductive isolation, P. japonicus is distinct from other Pristionchus species by its arrangement of genital papillae. All species in the complex are separated from each other by molecular sequence divergence, as indicated by analysis of 27 nuclear protein-coding genes and unique sequences of the small subunit ribosomal RNA gene. The identification of a species complex that includes P. pacificus is invaluable for studies of population genetics, speciation, and macroevolution, particularly the evolution of hermaphroditism in the genus.

Parapristionchus giblindavisi n. gen., n. sp. (Rhabditida: Diplogastridae) isolated from stag beetles (Coleoptera: Lucanidae) in Japan

A new species of diplogastrid nematode, isolated in a previous survey of nematodes associated with stag beetles in Japan, is described as Parapristionchus giblindavisi n. gen., n. sp. Parapristionchus n. gen. differs from other diplogastrid genera chiefly by its stomatal morphology. Distinguishing the genus are the presence of a claw-like dorsal tooth in both the eurystomatous and stenostomatous forms and the division of the cheilostom into 12 plates lacking apical flaps. According to phylogenetic analysis of nine ribosomal protein gene sequences, Parapristionchus n. gen. shows deep divergence from other known genera. Molecular evidence strongly supports P. giblindavisi n. gen., n. sp. + Pristionchus spp. as monophyletic with respect to all other diplogastrids examined. Congruent with a clade of P. giblindavisi n. gen., n. sp. + Pristionchus spp. is the shared presence of a bifurcate P7 genital papilla. Discovery and description of a close sister group to Pristionchus, a model biological system, enables character polarisation in macroevolutionary studies of Pristionchus nematodes.

Horizontal gene transfer of microbial cellulases into nematode genomes is associated with functional assimilation and gene turnover

BACKGROUND

Natural acquisition of novel genes from other organisms by horizontal or lateral gene transfer is well established for microorganisms. There is now growing evidence that horizontal gene transfer also plays important roles in the evolution of eukaryotes. Genome-sequencing and EST projects of plant and animal associated nematodes such as Brugia, Meloidogyne, Bursaphelenchus and Pristionchus indicate horizontal gene transfer as a key adaptation towards parasitism and pathogenicity. However, little is known about the functional activity and evolutionary longevity of genes acquired by horizontal gene transfer and the mechanisms favoring such processes.

RESULTS

We examine the transfer of cellulase genes to the free-living and beetle-associated nematode Pristionchus pacificus, for which detailed phylogenetic knowledge is available, to address predictions by evolutionary theory for successful gene transfer. We used transcriptomics in seven Pristionchus species and three other related diplogastrid nematodes with a well-defined phylogenetic framework to study the evolution of ancestral cellulase genes acquired by horizontal gene transfer. We performed intra-species, inter-species and inter-genic analysis by comparing the transcriptomes of these ten species and tested for cellulase activity in each species. Species with cellulase genes in their transcriptome always exhibited cellulase activity indicating functional integration into the host's genome and biology. The phylogenetic profile of cellulase genes was congruent with the species phylogeny demonstrating gene longevity. Cellulase genes show notable turnover with elevated birth and death rates. Comparison by sequencing of three selected cellulase genes in 24 natural isolates of Pristionchus pacificus suggests these high evolutionary dynamics to be associated with copy number variations and positive selection.

CONCLUSION

We could demonstrate functional integration of acquired cellulase genes into the nematode's biology as predicted by theory. Thus, functional assimilation, remarkable gene turnover and selection might represent key features of horizontal gene transfer events in nematodes.

Quantitative assessment of the nematode fauna present on Geotrupes dung beetles reveals species-rich communities with a heterogeneous distribution

Pristionchus spp. nematodes exhibit several traits that might serve as pre-adaptations to parasitism. Under harsh environmental conditions, these nematodes can arrest development and form dauer larvae. In addition, they have been shown to live in necromenic association with a range of beetles, including dung beetles ( Geotrupes stercorosus ) on which, for example, Pristionchus entomophagus is commonly found. It has been argued that the formation of dauer larvae and the association with invertebrates represent intermediate steps towards parasitism. To better understand necromenic associations, and to gain information on Pristionchus spp. abundance and the general species composition on dung beetles, we extracted all the nematode fauna present on 114 individuals of G. stercorosus. By direct sequencing using the 18S SSU, we provide a barcode for all nematodes isolated from the beetle samples. In total, 5,002 dauer-stage nematodes were sequenced, which included Pristionchus spp., Koerneria spp. (Diplogastridae), Pelodera spp. (Rhabditidae), and Strongyloidea as well as Spirurida. Intensities of infection varied from over 1,000 nematodes isolated from a single G. stercorosus to none, with Pelodera spp. being the most abundant group isolated. This study presents the first quantitative data on the Pristionchus spp. infection of beetles.

The Strongyloides (Nematoda) of sheep and the predominant Strongyloides of cattle form at least two different, genetically isolated populations

Strongyloides sp. (Nematoda) are very wide spread small intestinal parasites of vertebrates that can form a facultative free-living generation. Most authors considered all Strongyloides of farm ruminants to belong to the same species, namely Strongyloides papillosus (Wedl, 1856). Here we show that, at least in southern Germany, the predominant Strongyloides found in cattle and the Strongyloides found in sheep belong to separate, genetically isolated populations. While we did find mixed infections in cattle, one form clearly dominated. This variety, in turn, was never found in sheep, indicating that the two forms have different host preferences. We also present molecular tools for distinguishing the two varieties, and an analysis of their phylogenetic relationship with the human parasite Strongyloides stercoralis and the major laboratory model species Strongyloides ratti. Based on our findings we propose that Strongyloides from sheep and the predominant Strongyloides from cattle should be considered separate species as it had already been proposed by [Brumpt, E., 1921. Recherches sur le determinisme des sexes et de l'evolution des Anguillules parasites (Strongyloides). Comptes rendu hebdomadaires des séances et mémoires de la Société de Biologie et de ses filiales 85, 149-152], but was largely ignored by later authors. For nomenclature, we follow [Brumpt, E., 1921. Recherches sur le determinisme des sexes et de l'evolution des Anguillules parasites (Strongyloides). Comptes rendu hebdomadaires des séances et mémoires de la Société de Biologie et de ses filiales 85, 149-152] and use the name S. papillosus for the Strongyloides of sheep and the name Strongyloides vituli for the predominant Strongyloides of cattle.

The nematode Pristionchus pacificus (Nematoda: Diplogastridae) is associated with the oriental beetle Exomala orientalis (Coleoptera: Scarabaeidae) in Japan

Pristionchus pacificus has been developed as a nematode satellite organism in evolutionary developmental biology. Detailed studies of vulva development revealed multiple differences in genetic and molecular control in P. pacificus compared to the model organism Caenorhabditis elegans. To place evolutionary developmental biology in a comprehensive evolutionary context, such studies have to be complemented with ecology. In recent field studies in western Europe and eastern North America we found 11 Pristionchus species that are closely associated with scarab beetles and the Colorado potato beetle. However, P. pacificus was not commonly found in association with scarab beetles in these studies. Here, we describe the results of a similar survey of scarab beetles in Japan. Pristionchus pacificus was the most common Pristionchus species on scarab beetles in Japan, with 40 out of 43 (93%) isolates. The other Pristionchus isolates represent three novel species, which we refer to as Pristionchus sp. 11, Pristionchus sp. 14, and Pristionchus sp. 15. Thirty-seven of the established P. pacificus strains were found on the oriental beetle Exomala orientalis. Laboratory studies with the sex pheromone (Z)-7-tetradecen-2-one of the oriental beetle revealed that P. pacificus shows strong olfactory attraction to the beetle's sex pheromone, which provides a potential mechanism for the recognition and interaction of P. pacificus and E. orientalis. Together, this study identifies P. pacificus as the most common Pristionchus nematode in field studies in Japan, identifies E. orientalis as an important host species, and provides the basis for the ecology of P. pacificus.

Distinct patterns of genetic variation in Pristionchus pacificus and Caenorhabditis elegans, two partially selfing nematodes with cosmopolitan distribution

Hermaphroditism has evolved several times independently in nematodes. The model organism Caenorhabditis elegans and Pristionchus pacificus are self-fertile hermaphrodites with rare facultative males. Both species are members of different families: C. elegans belongs to the Rhabditidae and P. pacificus to the Diplogastridae. Also, both species differ in their ecology: C. elegans is a soil-dwelling nematode that is often found in compost heaps. In contrast, field studies in Europe and North America indicate that Pristionchus nematodes are closely associated with scarab beetles. In C. elegans, several recent studies have found low genetic diversity and rare out-crossing events. Little is known about diversity levels and population structure in free-living hermaphroditic nematodes outside the genus Caenorhabditis. Taking a comparative approach, we analyse patterns of molecular diversity and linkage disequilibrium in 18 strains of P. pacificus from eight countries and four continents. Mitochondrial sequence data of P. pacificus isolates reveal a substantially higher genetic diversity on a global scale when compared to C. elegans. A mitochondrial-derived hermaphrodite phylogeny shows little geographic structuring, indicating several worldwide dispersal events. Amplified fragment length polymorphism and single strand conformation polymorphism analyses demonstrate a high degree of genome-wide linkage disequilibrium, which also extends to the mitochondrial genome. Together, these findings indicate distinct patterns of genetic variation of the two species. The low level of genetic diversity observed in C. elegans might reflect a recent human-associated dispersal, whereas the P. pacificus diversity might reflect a long-lasting and ongoing insect association. Thus, despite similar lifestyle characteristics in the laboratory, the reproductive mode of hermaphroditism with rare facultative males can result in distinct genetic variability patterns in different ecological settings.

The free-living generation of the nematode Strongyloides papillosus undergoes sexual reproduction

The nematode genus Strongyloides consists of parasites that live as parthenogenetic females in the small intestines of their hosts. They can also form a facultative free-living generation with males and females. Recently, research on Strongyloides cellular and molecular biology has concentrated on Strongyloides ratti and Strongyloides stercoralis. We propose that the related nematode Strongyloides papillosus, a common parasite of ruminants, is well suited for comparative and evolutionary studies and we show that it is phylogentically basal to S. ratti and S. stercoralis. Based on cytological observations several reports have proposed that Strongyloides males do not contribute genetically to the next generation, leaving open the question of why males still exist. In contrast, the only study employing molecular markers showed that S. ratti males do pass on genetic material. Here, we demonstrate that in S. papillosus males also contribute molecular genetic markers to the next generation. This is interesting for two reasons. First, it shows that S. papillosus is amenable to genetic analysis and second, it indicates that sexual reproduction is more common in Strongyloides than previously assumed.

Sex, bugs and Haldane's rule: the nematode genus Pristionchus in the United States

Background

The nematode Pristionchus pacificus has been developed as a satellite organism in evolutionary developmental biology for comparison to Caenorhabditis elegans. Comparative studies have revealed major differences in the regulation of developmental processes between P. pacificus and C. elegans. To place evolutionary developmental biology and the observed developmental differences between species in a comprehensive evolutionary context, such studies have to be complemented with ecological aspects. Knowledge about the ecology of the organism in question might indicate specific environmental conditions that can result in developmental adaptations and could account for species differences in development. To this end, we have started to investigate the ecology of Pristionchus nematodes. In recent field studies in Western Europe we found six Pristionchus species that are closely associated with scarab beetles and the Colorado potato beetle. This Pristionchus – beetle association provides the unique opportunity to combine research in evolutionary developmental biology with ecology. However, it remains unknown how general these findings from Europe are on a global scale.

Results

Here, we describe the Pristionchus species associated with scarab and Colorado potato beetles in the Eastern United States and show striking transatlantic differences and unexpected evolutionary and ecological patterns. Twohundredeighty of 285 (98%) isolates from American scarab beetles belong to five Pristionchus species, all of which are different from the European species. We describe four of them as novel Pristionchus species. The five American Pristionchus species fall into a single phylogenetic clade and have a male-female (gonochoristic) mode of reproduction, whereas the majority of European isolates are hermaphroditic. Crosses between the two most closely related species, P. aerivorus and P. pseudaerivorus n. sp., follow Haldane's rule in that heterogametic F1 males are inviable. We observed P. aerivorus and P. pseudaerivorus n. sp. coexisting on the same scarab beetle and obtained two cases of F1 hybrids from wild beetles. Finally, the Colorado potato beetle is associated with the same nematode, P. uniformis in the United States and Europe. Given the introduction of the Colorado potato beetle to Europe in 1877, our results suggest that P. uniformis was introduced together with its beetle vector.

Conclusion

Taken together, the Pristionchus – beetle association provides a powerful tool for studying biodiversity, biogeography, speciation and species invasion on a global scale.

Nematodes of the genus Pristionchus are closely associated with scarab beetles and the Colorado potato beetle in Western Europe

Evolutionary developmental biology examines how changes in developmental programmes give rise to developmental and, ultimately, morphological novelty. To this end, comparisons of related but distinct organisms have to be performed. The diplogastrid nematode Pristionchus pacificus has been developed as a satellite system for a detailed comparison of various developmental processes to the model organism Caenorhabditis elegans, a rhabditid nematode. In addition to developmental and genetic studies, a genomic platform has been established to analyse the biology of this organism. However, only little is known about where and how Pristionchus pacificus and its relatives live in the wild. Here we show that nematodes of the genus Pristionchus live in close association with scarabaeoid beetles and the Colorado potato beetle. In total, we generated 371 isogenic female lines from 4242 beetles collected at 25 sampling sites all over Europe. Isogenic female lines were subjected to sequence analysis and mating experiments for species determination. The 371 isolates fell into six species. Two hermaphroditic species account for about 60% of the collected nematodes. We found Pristionchus maupasi almost exclusively on cockchafers and Pristionchus entomophagus predominantly on dung beetles. Colorado potato beetles carried the gonochoristic species Pristionchus uniformis, which was only rarely observed on scarabaeoid beetles. We describe the initial evidence for the association of Pristionchus nematodes with beetles and provide a phylogeny based on sequence analysis of the small subunit ribosomal RNA gene.

Organization, alternative splicing, polymorphism, and phylogenetic position of lamprey CD45 gene

CD45 of jawed vertebrates is a receptor-type protein tyrosine phosphatase regulating lymphocyte development and activation. To shed light on the evolution of the CD45 gene, the organization of its orthologue in the lamprey, a jawless vertebrate, was determined. Compared to its mammalian and fugu counterparts, the lamprey gene was found to be lacking several exons in the segment encoding the extracellular part of the protein. In consequence, this part contains only one instead of the two or three fibronectin type III domains typical of the mammalian molecules. The lamprey transcripts of the CD45 gene occur in several variants originating by alternative splicing, including some not observed previously in other vertebrates. Most remarkable of these are splice variants generated by the use of intra-exonic splicing signals and thus lacking one half, one third, or two thirds of an exon and yet apparently translated in the correct reading frame. The lamprey gene contains polymorphic sites not only in the segment encoding the extracellular portion but also in the segment specifying the cytoplasmic part of the molecule. Polymorphism is generated by both mutations and recombination. Some of the alleles may have persisted long enough to represent transspecies polymorphism presumably maintained by positive selection. Phylogenetic analysis suggests that ancestors of the CD45 gene may have existed before the divergence of coelomate from pseudocoelomate metazoans.

Prototypic T cell receptor and CD4-like coreceptor are expressed by lymphocytes in the agnathan sea lamprey

All jawed vertebrates have highly diverse lymphocyte receptors, which allow discrimination between self and nonself antigens as well as the recognition of potential pathogens. Key elements of the anticipatory recombinatorial immune system in jawed vertebrates are the TCR, Ig, and MHC genes, but their ancestral genes have not been found in more basal vertebrates. In this study, we extended our analysis of the transcriptome of lymphocyte-like cells in the lamprey to identify the TCR-like and CD4-like genes. The structural features of these genes and their preferential expression in lymphocytes make them attractive candidates for ancestral TCR and CD4 genes. The TCR-like gene contains both V (variable) and J (joining) sequences in its first exon and exists as a single-copy gene that is invariant. Thus, the TCR-like gene cannot account for the receptor diversity that is required for the immune responses reported for lamprey, but it could have been easily modified to serve as an evolutionary precursor of modern TCR and Ig genes.

The Ig-Like Domain of Tapasin Influences Intermolecular Interactions

Presentation of antigenic peptides to T lymphocytes by MHC class I molecules is regulated by events involving multiple endoplasmic reticulum proteins, including tapasin. By studying the effects of substitutions in the tapasin Ig-like domain, we demonstrated that H-2L(d)/tapasin association can be segregated from reconstitution of folded L(d) surface expression. This finding suggests that peptide acquisition by L(d) is influenced by tapasin functions that are independent of L(d) binding. We also found that the presence of a nine-amino acid region in the Ig-like domain of mouse or human tapasin is required for association with L(d), and certain point substitutions in this sequence abrogate human, but not mouse, tapasin association with L(d). These data are consistent with a higher overall affinity between L(d) and mouse tapasin compared with human tapasin. In addition, we found that other point mutations in the same region of the tapasin Ig-like domain affect MHC class I surface expression and Ag presentation. Finally, we showed that the cysteine residues in the Ig-like domain of tapasin influence tapasin's stability, its interaction with the MHC class I H chain, and its stabilization of TAP. Mutagenesis of these cysteines decreases tapasin's electrophoretic mobility, suggesting that these residues form an intramolecular disulfide bond. Taken together, these results reveal a critical role for the tapasin Ig-like domain in tapasin function.

Phylogenetic relationships among East African haplochromine fish as revealed by short interspersed elements (SINEs)

Genomic DNA libraries were prepared from two endemic species of Lake Victoria haplochromine (cichlid) fish and used to isolate and characterize a set of short interspersed elements (SINEs). The distribution and sequences of the SINEs were used to infer phylogenetic relationships among East African haplochromines. The SINE-based classification divides the fish into four groups, which, in order of their divergence from a stem lineage, are the endemic Lake Tanganyika flock (group 1); fish of the nonendemic, monotypic, widely distributed genus Astatoreochromis (group 2); the endemic Lake Malawi flock (group 3); and group 4, which contains fish from widely dispersed East African localities including Lakes Victoria, Edward, George, Albert, and Rukwa, as well as many rivers. The group 4 haplochromines are characterized by a subset of polymorphic SINEs, each of which is present in some individuals and absent in others of the same population at a given locality, the same morphologically defined species, and the same mtDNA-defined haplogroup. SINE-defined group 4 contains six of the seven previously described mtDNA haplogroups. One of the polymorphic SINEs appears to be fixed in the endemic Lake Victoria flock; four others display the presence-or-absence polymorphism within the species of this flock. These findings have implications for the origin of Lake Victoria cichlids and for their founding population sizes.

Clustering of C-Type Lectin Natural Killer Receptor-Like Loci in the Bony Fish Oreochromis niloticus

The genome of the cichlid (teleost) fish Oreochromis niloticus contains a set of genes which encode group V C-type lectin proteins homologous to the mammalian NKG2/CD94 family of natural killer (NK) cell receptors. To determine the genomic organization of these killer cell-like receptor (KLR) genes, an O. niloticus BAC library was screened with a cDNA probe derived previously from an expressed sequence tag of the related cichlid species Paralabidochromis chilotes. Four distinct KLR-bearing BAC clones were analysed, three of which could be assembled into a contig. One of the clones was sequenced in its entirety, whereas the others were partially sequenced to identify the KLR loci borne by them. Altogether, 28 distinct KLR loci were identified, of which at least 26 occupy a single chromosomal region, the KLR complex. One half of the loci appear to be occupied by pseudogenes. Compared to the human NK cell receptor complex, the Oreochromis KLR complex is more compact and, apart from transposons, appears to contain only KLR loci. The gene density of the complex is one KLR locus per 18 kb of sequence. All the KLR loci constituting the complex are derived from a single most recent common ancestor, which is estimated to have existed 7.7 million years ago. The 180 kb of the determined sequence is a mosaic of blocks of similar segments reflecting a complex history of duplications, deletions and rearrangements. The transposons found in the sequenced part belong to the TC1, Xena, CR1 and TX1 families.

Phylogeny of Antigen-Processing Enzymes: Cathepsins of a Cephalochordate, an Agnathan and a Bony Fish

Cathepsins are enzymes that have been cleaving peptide bonds of lysosomal proteins probably since lysosomes appeared in early eucaryotes. When the adaptive system emerged in gnathostomes, cathepsins were recruited to produce peptides for loading onto the major histocompatibility complex class II molecules and for degrading the class II-associated invariant chain just before the loading. The circumstances under which this recruitment took place are unclear because the knowledge about vertebrate cathepsins is limited largely to mammals. To shed light on the recruitment, 10 amphioxus, one lamprey and one cichlid fish cathepsin cDNA clone were characterized and analysed phylogenetically. Disregarding cathepsin O, whose phylogenetic position is uncertain, the analysis confirms the existence of two old lines of descent, the B and the L lineages of cathepsins, which diverged from each other early in the evolution of eucaryotes. The B lineage encompasses cathepsins B, C and Z (X). The L lineage splits off sublineages encompassing cathepsins F and W before the plant-animal separation and cathepsin H early in the evolution of the metazoa. The remaining cathepsins belonging to the L lineage diverged from one another during the evolution of vertebrates: S, K and L before the emergence of bony fishes, and the group of rodent placentally expressed cathepsins [J (P), M, Q, R, 3, 6, 7 and 8] as well as the testis/ova-expressed cathepsins (testins) probably after the divergence of rodents from primates. The part possibly played by the adaptive immune system in some of these divergences is discussed.

Origin and speciation of haplochromine fishes in East African crater lakes investigated by the analysis of their mtDNA, Mhc genes, and SINEs

The Western Branch of the East African Great Rift Valley is pocketed with craters of extinct or dormant volcanoes. Many of the craters are filled with water, and the lakes are inhabited by fishes. The objective of the present study was to determine the amount and nature of genetic variation in haplochromine fishes inhabiting two of these crater lakes, Lake Lutoto and Lake Nshere, and to use this information to infer the origin and history of the two populations. To this end, sequences of mitochondrial (mt) DNA control region, exon 2 of major histocompatibility complex (Mhc) class II B genes, and short interspersed elements (SINEs) were analyzed. The results indicate that the Lake Nshere and Lake Lutoto fishes originated from different but related large founding populations derived from the Kazinga Channel, which connects Lake Edward and Lake George. Some of the genetic polymorphism that existed in the ancestral populations was lost in the populations of the two lakes. The polymorphism that has been retained has persisted for some 50000 generations (years). During this time, new mutations arose and became fixed in each of the two populations in the mtDNA, giving rise to sets of diagnostic substitutions. Each population evolved in isolation after the colonization of the lakes less than 50000 years ago. There appears to be no population structure within the crater lake fishes, and their present effective population sizes are in the order of 104 to 105 individuals. Comparisons with the endemic haplochromine species of Lake Victoria reveal interesting parallels, as well as differences, which may help to understand the nature of the speciation process.

A molecule bearing an immunoglobulin-like V region of the CTX subfamily in amphioxus

An expressed sequence tag with significant similarity to a vertebrate T-cell receptor (Tcr) sequence was found in a cDNA library prepared from the posterior part of the adult amphioxus, Branchiostoma lanceolatum ( Brla). Characterization of the corresponding cDNA clone revealed the presence of an open reading frame encoding a 351 amino acid residue-long polypeptide. The putative protein, tentatively designated Brla-VDB for "V-domain bearing", appears to consist of two domains. The N-terminal domain begins with a putative leader peptide followed by a sequence resembling the V domain of the CTX protein originally found to be expressed in the cortical thymocytes of the clawed frog Xenopus. The C-terminal domain of the VDB protein does not show significant sequence similarity to any entry in the databases and contains five hydrophobic segments separated by short intervening hydrophilic stretches. It may therefore belong to a protein that crosses the plasma membrane five times. These findings support the notion that V domains resembling those found in Tcrs evolved in nonvertebrates before the emergence of the adaptive immune system and may have participated in functions not involved directly in immunity.

Genes encoding putative natural killer cell C-type lectin receptors in teleostean fishes

Mammalian natural killer (NK) cells are cytotoxic lymphocytes that express receptors specific for MHC class I molecules. The NK cell receptors belong to two structurally unrelated families, the killer cell Ig-like receptors and the killer cell C-type lectin receptors. We describe a cDNA clone derived from the bony (cichlid) fish Paralabidochromis chilotes and show that it encodes a protein related to the CD94/NK cell group 2 (NKG2) subfamily of the killer cell C-type lectin receptors. The gene encoding this receptor in a related species, Oreochromis niloticus, has a similar structure to the human CD94/NKG2 genes and is a member of a multigene cluster that resembles the mammalian NK cell gene complex. Thus, the CD94/NKG2 subfamily of NK cell receptors must have arisen before the divergence of fish and tetrapods and may have retained its function (possibly monitoring the expression of MHC class I molecules) for >400 million years.

Macrophage migration inhibitory factor (MIF) of jawed and jawless fishes: implications for its evolutionary origin

The macrophage migration inhibitory factor (MIF) is a cytokine produced by T lymphocytes and macrophages in response to inflammatory stimuli. We sequenced MIF cDNA clones of two jawless fishes, the sea lamprey (Petromyzon marinus) and the North Atlantic hagfish (Myxine glutinosa), as well as of the jawed (cichlid) fish Paralabidochromis chilotes. The fish MIF-encoding genes have the same exon-intron organization as the mammalian MIF genes and are present in one copy per haploid genome. Secondary and tertiary structure predictions suggest that the fish MIF proteins have a topology characteristic of the entire MIF-family of proteins. Phylogenetic analysis separates the known nematode members of the family into two groups, one having a sister group relationship with the mammalian D-dopachrome tautomerase (DDT) proteins and the other being related to vertebrate MIFs. It also reveals a high degree of convergent evolution among the members of the family. Finally, it suggests that the divergence of MIF and DDT occurred before the emergence of nematodes in metazoan evolution.

Identification and characterization of a TAP-family gene in the lamprey

An expressed sequence tag obtained from a sea lamprey ( Petromyzon marinus) cDNA library was used to obtain a full-length coding sequence showing significant similarity to ABCB transporter proteins. The sequence is closely related to the mammalian ABCB9 protein and the TAP1 and TAP2 proteins that transport peptides for loading onto nascent Mhc class I molecules. The Pema-ABCB9 gene has an exon-intron organization similar to that of the mammalian TAP genes, with the exception of exon 2, which in the lamprey is split into two by a 949-bp long intron. The gene probably occurs in a single copy in the haploid lamprey genome. The ABCB9 genes appear to be evolving four-to-ten times slower than the TAP1 and TAP2 genes. Six putative transmembrane helices and the nucleotide-binding domain of the lamprey ABCB9 protein show high sequence similarity with the TAP1 and TAP2 molecules. The lamprey protein also contains sequence stretches that resemble the putative peptide interacting parts of the TAP1 and TAP2 molecules, but are peppered with ABCB9-specific residues.

Identification of chemokines and a chemokine receptor in cichlid fish, shark, and lamprey

Chemokines are small, inducible, structurally related proteins that guide cells expressing the right chemokine receptors to sites of immune response. They have been identified and studied extensively in mammals, but little is known about their presence in other vertebrate groups. Here we describe seven new chemokines in bony fish and one in a cartilaginous fish, as well as one chemokine receptor in a jawless vertebrate. All eight chemokines belong to the SCYA (CC) subfamily characterized by four conserved cysteine residues of which the first two are adjacent. The chemokine receptor is of the CXCR4 type. Phylogenetic analysis does not reveal any clear evidence of orthology of fish and human chemokines. Although the divergence of the subfamilies began before the fish-tetrapod split, much of the divergence within the subfamilies took place separately in the two vertebrate groups. The existence of a chemokine receptor in the lamprey indicates that chemokines are apparently also present in the Agnatha.

Persistence of Mhc heterozygosity in homozygous clonal killifish, Rivulus marmoratus: implications for the origin of hermaphroditism

The mangrove killifish Rivulus marmoratus, a neotropical fish in the order Cyprinodontiformes, is the only known obligatorily selfing, synchronous hermaphroditic vertebrate. To shed light on its population structure and the origin of hermaphroditism, major histocompatibility complex (Mhc) class I genes of the killifish from seven different localities in Florida, Belize, and the Bahamas were cloned and sequenced. Thirteen loci and their alleles were identified and classified into eight groups. The loci apparently arose approximately 20 million years ago (MYA) by gene duplications from a single common progenitor in the ancestors of R. marmoratus and its closest relatives. Distinct loci were found to be restricted to different populations and different individuals in the same population. Up to 44% of the fish were heterozygotes at Mhc loci, as compared to near homozygosity at non-Mhc loci. Large genetic distances between some of the Mhc alleles revealed the presence of ancestral allelic lineages. Computer simulation designed to explain these findings indicated that selfing is incomplete in R. marmoratus populations, that Mhc allelic lineages must have diverged before the onset of selfing, and that the hermaphroditism arose in a population containing multiple ancestral Mhc lineages. A model is proposed in which hermaphroditism arose stage-wise by mutations, each of which spread through the entire population and was fixed independently in the emerging clones.

The effect of selection on a long wavelength-sensitive (LWS) opsin gene of Lake Victoria cichlid fishes

In East African Lake Victoria >200 endemic species of haplochromine fishes have been described on the basis of morphological and behavioral differences. Yet molecular analysis has failed to reveal any species-specific differences among these fishes in either mitochondrial or nuclear genes. Although the genes could be shown to vary, the variations represent trans-species polymorphisms not yet assorted along species lines. Nevertheless, fixed genetic differences must exist between the species at loci responsible for the adaptive characters distinguishing the various forms from one another. Here we describe variation and fixation at the long wavelength-sensitive (LWS) opsin locus, which is selection-driven, adaptive, and if not species- then at least population-specific. Because color is one of the characters distinguishing species of haplochromine fishes and color perception plays an important part in food acquisition and mate choice, we suggest that the observed variation and fixation at the LWS opsin locus may have been involved in the process that has led to the spectacular species divergence of haplochromine fishes in Lake Victoria.

Isolation and characterization of lymphocyte-like cells from a lamprey

Lymphocyte-like cells in the intestine of the sea lamprey, Petromyzon marinus, were isolated by flow cytometry under light-scatter conditions used for the purification of mouse intestinal lymphocytes. The purified lamprey cells were morphologically indistinguishable from mammalian lymphocytes. A cDNA library was prepared from the lamprey lymphocyte-like cells, and more than 8,000 randomly selected clones were sequenced. Homology searches comparing these ESTs with sequences deposited in the databases led to the identification of numerous genes homologous to those predominantly or characteristically expressed in mammalian lymphocytes, which included genes controlling lymphopoiesis, intracellular signaling, proliferation, migration, and involvement of lymphocytes in innate immune responses. Genes closely related to those that in gnathostomes control antigen processing and transport of antigenic peptides could be ascertained, although no sequences with significant similarity to MHC, T cell receptor, or Ig genes were found. The data suggest that the evolution of lymphocytes in the lamprey has reached a stage poised for the emergence of adaptive immunity.

Lamprey lymphocyte-like cells express homologs of genes involved in immunologically relevant activities of mammalian lymphocytes

To shed light on the origin of adaptive immunity, a cDNA library was prepared from purified lymphocyte-like cells of a jawless vertebrate, the sea lamprey (Petromyzon marinus). Randomly selected cDNA clones were sequenced, and their homologies to proteins in the databases were determined. Of the sequences homologous to proteins involved in immune responses, five were selected for further characterization. Their encoding genes corresponded to loci that in jawed vertebrates are essential for activities of lymphocytes. These activities include regulation of T and B cell stimulation and proliferation (CD45); stabilization of molecular complexes involved in lymphocyte activation, adhesion, migration, and differentiation (CD9/CD81); adaptor functions in signaling leading to the activation of B lymphocytes (BCAP) and T lymphocytes (CAST); and amino acid transport associated with cell activation (CD98). The presence of these genes in the lamprey genome and their expression in lymphocyte-like cells support the notion that these cells perform many of the functions of gnathostome lymphocytes. It reopens the question of the stage jawless fishes reached in the evolution of their immune system.

Linkage relationships of genes coding for alpha2-macroglobulin, C3 and C4 in the zebrafish: implications for the evolution of the complement and Mhc systems

The alpha2-macroglobulin (A2M) and the complement components C3 and C4 are related proteins derived from a common ancestor. Theoretically, this derivation could have occurred either by tandem duplications of their encoding genes or by polyploidization involving chromosomal segments, a chromosome or the whole genome. In tetrapods the A2M-, C3- and C4-encoding genes are generally each located on a different chromosome. This observation has been interpreted as supporting their origin by polyploidization. We identified and mapped (with the help of a radiation hybrid panel of cell lines) the A2M, C3 and C4 loci in the zebrafish, Danio rerio. Each of the three types of loci is present in the zebrafish in multiple copies, but all of the identified copies of a given type map to the same region in linkage groups 1 (C3) and 15 (A2M, C4). The A2M and C4 loci are mapped in the same region not linked to any of the class I or class II major histocompatibility complex (Mhc) loci. These observations are interpreted as supporting the origin of the A2M family of genes by tandem duplications, followed by the dispersal of the copies to different chromosomes. It is also argued that the association of C4 with the class I/II loci in tetrapods is accidental and without functional significance.

Identification of two Ikaros-like transcription factors in lamprey

The jawless Agnatha (lampreys and hagfishes) represent the phylogenetically oldest order of vertebrates that are believed to lack the adaptive immune system of jawed vertebrates. In order to search for molecular markers specific for cellular components of the adaptive immune system in lampreys, we used the polymerase chain reaction (PCR) to identify genes for transcription factors of the Ikaros family in genomic DNA and cDNA libraries from two species of lampreys, Petromyzon marinus and Lampetra fluviatilis. The mammalian Ikaros-like family of transcription factors consists of five members, Ikaros, Helios, Aiolos, Eos and Pegasus, of which the first three appear to be essential for lymphocyte development. Two different Ikaros-like genes, named IKLF1 and IKLF2, were identified in lamprey. They both have the conserved exon-intron structure of seven exons and show alternative splicing like their counterparts in jawed vertebrates. The genes code for predicted proteins of 589 and 513 amino acid residues, respectively. The proteins contain six highly conserved zinc finger motifs that are 83-91% identical to the mammalian members of the Ikaros-like family. The remaining parts of the sequences are, however, mostly unalignable. Phylogenetic analysis based on the alignable segments of the sequences does not identify the orthologous gene in jawed vertebrates but rather shows equidistance of the lamprey Ikaros-like factors to each other and to Ikaros, Helios, Aiolos and Eos. Expression studies by reverse transcription (RT)-PCR and in situ hybridization (ISH), however, provide evidence for moderate expression in presumed lymphoid tissues like the gut epithelium and for high levels of expression in the gonads, especially in the ovary.

Evolution of Mhc class II B genes in Darwin's finches and their closest relatives: birth of a new gene

The 15 extant species of Darwin's finches on the Galápagos and Cocos Islands are the products of an unfinished adaptive radiation from a founder flock of birds related to the South American species Tiaris obscura. Molecular characterization of their major histocompatibility complex ( Mhc) class II B genes has revealed the existence of several related groups of sequences (presumably encoded in distinct loci) from which one (group 5) stands out because of its low divergence over extended time periods. Analysis of group 5 exon 2 and intron 2 sequences has revealed that the encoding locus apparently arose 2-3 million years ago in the Tiaris group of South and Central American Thraupini. The locus shows no evidence of inactivation, but displays a very low degree of polymorphism, both in terms of number of alleles and genetic distances between alleles. Some of the polymorphism, however, appears to be trans-specific. All the observed intergenic differences can be explained by point mutations and most of the exon 2 changes represent non-synonymous substitutions, although the rate of non-synonymous and synonymous substitutions appears to be the same. The origin of the new locus is explained by the birth-and-death model of Mhc evolution with two important extensions. First, the ancestor of the group 5 genes may have arisen without new gene duplication and second, the birth of the new group may have been brought about by a switch from balancing to directional selection. The ancestor of the group 5 genes may have been a classical class II B allele (one of many) which directional selection fixed in the ancestral population and drove into the category of nonclassical genes.

Mhc class I genes of swordtail fishes, Xiphophorus: variation in the number of loci and existence of ancient gene families

Swordtail fishes and platies in the genus Xiphophorus (order Cyprinodontiformes, Teleostei) encompass 22 closely related species which are the products of a recent adaptive radiation in the streams of Central America. To investigate the evolution of the major histocompatibility complex (Mhc) genes in the period immediately following speciation, the class I genes from 20 of the 22 species were cloned and characterized by sequencing. The analysis revealed the existence of multiple loci (at least seven in some individuals) whose numbers vary among the different species and probably also among individuals of the same species. The variation does not seem to bear any relationship to the taxonomy of the genus. Genes at the different loci are distinguished by their intron sequences and by the presence of characteristic motifs in exons 2 and 3. The variation in copy number of loci may have been effected in part by unequal crossing over occurring between introns of misaligned closely related genes. The sequences of the genes fall into two groups, A and B, which represent ancient lineages. The groups define two families of loci, which diverged from each other an estimated 85 million years ago, before the separation of the Acanthopterygii from the Paracanthopterygii of the advanced bony fishes. Evolution of the genes within each family can be explained by the birth-and-death process driven by gene duplications and mutational differentiation.

Classification and phylogenetic relationships of African tilapiine fishes inferred from mitochondrial DNA sequences

African cichlid fishes are composed of two major lineages, the haplochromines and the tilapiines. Whereas the phylogenetic relationships of the haplochromines have been studied extensively, primarily because of their spectacular adaptive radiations in the Great Lakes of East Africa, little is known about the relationships among the tilapiine species, despite the fact that they have become an important component of African, indeed world, aquaculture. To remedy this situation, molecular phylogenetic analysis of tilapiine fishes was undertaken. A segment of mitochondrial DNA encompassing the terminal part of the tRNA(Pro) gene and the most variable part of the control region was amplified by the polymerase chain reaction with DNA samples isolated from 42 tilapiine species, and the amplification products were subjected to heteroduplex analysis and sequencing. Phylogenetic trees based on 68 sequences revealed the existence of 11 sequence groups and 11 single-sequence branches. The groups, designated Ti1 through Ti11, were distinguished by specific combinations of diagnostic substitutions, formation of monophyletic clusters, and separation by genetic distances in excess of 0.04. Although the relationships among the groups could not be resolved, the sequences separated Oreochromis and Sarotherodon from Tilapia, as defined by Trewavas. The Oreochromis sequences clustered with the Sarotherodon sequences and thus supported the hypothesis that the mouthbrooding behavior of the tilapiine fishes evolved only once from the substrate-spawning behavior. Since on phylogenetic trees the O. alcalicus (sub)species were always separated from O. amphimelas by other Oreochromis species, it was concluded that the adaptation to life in water with a high salt concentration and high pH values evolved independently at least twice in the tilapiine fishes. The tilapiines diverged from the haplochromines more than 8 million years ago; most of the intragroup divergences among the tilapiines took place an estimated 1.1 to 6 million years ago.

Isolation and sequencing of cDNA clones coding for the catalytic unit of glucose-6-phosphatase from two haplochromine cichlid fishes

Complementary DNA clones coding for the catalytic unit of the enzyme glucose-6-phosphatase (G6Pase) were obtained from Haplochromis nubilus and Haplochromis xenognathus, two cichlid fish species from Lake Victoria. The translated sequence of these two cDNAs identifies a polypeptide consisting of 352 amino acid residues and showing a 54.4% similarity to the human form of G6Pase. The amino acid sequences of the two fish species are identical. The comparison of the fish amino acid sequence with the corresponding sequences of rat, mouse, and human G6Pase revealed that the amino acid residues, which are involved in G6Pase catalysis in humans, are also conserved in fish G6Pase. Northern blot analysis showed that G6Pase is expressed at the same level in 6- and 10-day-old fish. A three base pair insertion/deletion polymorphism was found in the 3'-untranslated region of the fish G6Pase gene. The polymorphism will be a useful marker in a phylogenetic study of Lake Victoria cichlids.

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 6,000 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.

The origin and age of haplochromine fishes in Lake Victoria, east Africa

According to a widely held view, the more than 300 species of haplochromine cichlid fishes in Lake Victoria (LV), East Africa, originated from a single founder species in less than 12,000 years. This view, however, does not follow from the published geological and molecular evidence. The former does indeed suggest that the LV basin dried out less than 15,000 years ago, but it does not provide any information about the species that re-colonized the new lake or that remained in the rivers draining the area. The molecular evidence is inconclusive with respect to the origin of the LV haplochromines because cichlids from critical regions around LV were not adequately sampled; and as far as the age of the LV haplochromines is concerned, it in fact led to an estimate of 250,000-750,000 years old. In the present study, mitochondrial DNA (control region) variation was determined by heteroduplex and sequencing analyses of more than 670 specimens collected at widely distributed East African riverine and lacustrine localities. The analyses revealed the existence of seven haplogroups (I-VII) distinguishable by characteristic substitutions. All endemic LV samples tested fell into one of these haplogroups (V) which, however, was also found to be present at various other localities, both riverine and lacustrine, outside LV. Within this haplogroup, four subgroups (VA through VD) could be distinguished, two of which (VB and VC) were represented in LV and at other localities. The great majority of the LV haplochromine species could be classified as belonging to the VC subgroup, which was found only in LV and in the rivers draining into it. Hence, while the endemic haplochromine species of LV could not have originated from a single founding population, the lake does harbour a large species flock which probably arose in situ.

Molecular phylogeny of Philippine freshwater sardines based on mitochondrial DNA analysis

The commercially important Sardinella species (family Clupeidae or herrings) usually thrive in marine environments. An exception is Sardinella tawilis of Taal Lake, Batangas, Philippines, the only known freshwater sardine. This species is believed to have immigrated from Balayan Bay to the lake when it was formed in the course of volcanic eruptions some 240 years ago. To determine the relationship of S. tawilis to the marine species S. albella, S. fimbriata, and S. longiceps from the Balayan Bay we sequenced 358 bp of the cytochrome b gene and the mitochondrial control region. The cytochrome b gene was highly conserved and contained little phylogenetic information. The control region sequences, however, demonstrated two highly diversified main haplotypes grouping S. tawilis with S. albella, as shown by maximum parsimony and neighbor-joining analysis. The haplotypes are characterized by the presence of an 81 bp indel and up to eight 35 bp tandem repeat elements. The repeat copy number varied within individuals of S. tawilis and S. albella, thus showing heteroplasmy in these two species only. The analysis of two subpopulations of S. tawilis revealed restricted substitutions that may indicate the beginning of genetic differentiation of the two subpopulations.

Ancient allelism at the cytosolic chaperonin-alpha-encoding gene of the zebrafish

The T-complex protein 1, TCP1, gene codes for the CCT-alpha subunit of the group II chaperonins. The gene was first described in the house mouse, in which it is closely linked to the T locus at a distance of approximately 11 cM from the Mhc. In the zebrafish, Danio rerio, in which the T homolog is linked to the class I Mhc loci, the TCP1 locus segregates independently of both the T and the Mhc loci. Despite its conservation between species, the zebrafish TCP1 locus is highly polymorphic. In a sample of 15 individuals and the screening of a cDNA library, 12 different alleles were found, and some of the allelic pairs were found to differ by up to nine nucleotides in a 275-bp-long stretch of sequence. The substitutions occur in both translated and untranslated regions, but in the former they occur predominantly at synonymous codon sites. Phylogenetically, the alleles fall into two groups distinguished also by the presence or absence of a 10-bp insertion/deletion in the 3' untranslated region. The two groups may have diverged as long as 3.5 mya, and the polymorphic differences may have accumulated by genetic drift in geographically isolated populations.

Fish trypanosomes: their position in kinetoplastid phylogeny and variability as determined from 12S rRNA kinetoplast sequences

Fish trypanosomes have traditionally been classified according to the host species from which they were isolated, each isolate being regarded as a distinct species. To test the soundness of this practice, the genetic variabilities of the kinetoplast 12S rRNA-encoding genes of different fish trypanosomes isolates were compared. The DNAs were extracted from trypanosomes cloned from blood samples of 15 donors representing ten different fish species in four orders from waters of three major river systems of Central and Northern Europe. Comparison with other trypanosomatid sequences revealed that the fish trypanosomes form a monophyletic group with Trypanosoma brucei as a sister group. Pairwise comparisons of genetic distances yielded a wide range of continuous variation with no indication of any discontinuities attributable to barriers to gene flow. The genetic distances did not correlate with either the identity of the host species or geography. The host specificity of fish trypanosomes appears to be limited.

cDNA sequence coding for the alpha'-chain of the third complement component in the African lungfish

cDNA clones coding for almost the entire C3 alpha-chain of the African lungfish (Protopterus aethiopicus), a representative of the Sarcopterygii (lobe-finned fishes), were sequenced and characterized. From the sequence it is deduced that the lungfish C3 molecule is probably a disulphide-bonded alpha:beta dimer similar to that of the C3 components of other jawed vertebrates. The deduced sequence contains conserved sites presumably recognized by proteolytic enzymes (e.g. factor I) involved in the activation and inactivation of the component. It also contains the conserved thioester region and the putative site for binding properdin. However, the site for the interaction with complement receptor 2 and factor H are poorly conserved. Either complement receptor 2 and factor H are not present in the lungfish or they bind to different residues at the same or a different site than mammalian complement receptor 2 and factor H. The C3 alpha-chain sequences faithfully reflect the phylogenetic relationships among vertebrate classes and can therefore be used to help to resolve the long-standing controversy concerning the origin of the tetrapods.

Persistence of neutral polymorphisms in Lake Victoria cichlid fish

Phylogenetic trees for groups of closely related species often have different topologies, depending on the genes used. One explanation for the discordant topologies is the persistence of polymorphisms through the speciation phase, followed by differential fixation of alleles in the resulting species. The existence of transspecies polymorphisms has been documented for alleles maintained by balancing selection but not for neutral alleles. In the present study, transspecific persistence of neutral polymorphisms was tested in the endemic haplochromine species flock of Lake Victoria cichlid fish. Putative noncoding region polymorphisms were identified at four randomly selected nuclear loci and tested on a collection of 12 Lake Victoria species and their putative riverine ancestors. At all loci, the same polymorphism was found to be present in nearly all the tested species, both lacustrine and riverine. Different polymorphisms at these loci were found in cichlids of other East African lakes (Malawi and Tanganyika). The Lake Victoria polymorphisms must have therefore arisen after the flocks now inhabiting the three great lakes diverged from one another, but before the riverine ancestors of the Lake Victoria flock colonized the Lake. Calculations based on the mtDNA clock suggest that the polymorphisms have persisted for about 1.4 million years. To maintain neutral polymorphisms for such a long time, the population size must have remained large throughout the entire period.

Phylogeny of African cichlid fishes as revealed by molecular markers

The species flocks of cichlid fish in the three great East African Lakes, Victoria, Malawi, and Tanganyika, have arisen in each lake by explosive adaptive radiation. Various questions concerning their phylogeny have not yet been answered. In particular, the identity of the ancestral founder species and the monophyletic origin of the haplochromine cichlids from the East African lakes have not been established conclusively. In the present study, we used the anonymous nuclear DNA marker DXTU1 as a step towards answering these questions. A 280 bp-fragment of the DXTU1 locus was amplified by the polymerase chain reaction from East African lacustrine species, the East African riverine cichlid species Haplochromis bloyeti, H. burtoni and H. sparsidens, and other African cichlids. Sequencing revealed several indels and substitutions that were used as cladistically informative markers to support a phylogenetic tree constructed by the neighbor-joining method. The topology, although not supported by high bootstrap values, corresponds well to the geographical distribution and previous classification of the cichlids. Markers could be defined that: (i) differentiate East African from West African cichlids; (ii) distinguish the riverine and Lake Victoria/Malawi haplochromines from Lake Tanganyika cichlids; and (iii) indicate the existence of a monophyletic Lake Victoria cichlid superflock which includes haplochromines from satellite lakes and East African rivers. In order to resolve further the relationship of East African riverine and lacustrine species, mtDNA cytochrome b and control region segments were sequenced. The mtDNA-based trees support the notion of the monophyly of the Lake Victoria superflock but are ambiguous with respect to the phylogenetic position of the Lake Malawi flock.

Isolation of candidate hybrid sterility 1 genes by cDNA selection in a 1.1 megabase pair region on mouse chromosome 17

The Hybrid sterility 1 (Hst1) gene causes male infertility in crosses between certain inbred strains of the laboratory and wild mouse, Mus musculus. To identify the causative gene, we have searched YAC clones encompassing the Hst1 region for testis-expressed sequences, using the cDNA selection method. We isolated 12 non-overlapping cDNA clones, sequenced them, and placed them on a physical map based on the analysis of YAC clones and total genomic DNA. The cDNA clones map to ten loci. Three cDNA sequences correspond to the proteasome subunit C5 (locus Psmb1), ornithine decarboxylase (Odc-rs15), and penta-zinc finger (Zfp91-rs1) transcripts. Three of the ten testis-expressed loci described in this report (D17Ph4e, Psmb1, and Zfp91-rs1) co-segregate with all Hst1 recombinants and, together with the Tbp gene, are therefore potential candidates for the Hst1 gene. The presented physical and genetic mapping data indicate there are no gross rearrangements distinguishing the Hst1(f) and Hst1(s) alleles.

HLA-DRB intron 1 sequences: implications for the evolution of HLA-DRB genes and haplotypes

Human DRB genes encode beta chains of the major histocompatibility complex (MHC) class II molecules. Although nine DRB loci have been mapped to the short arm of chromosome 6, an individual chromosome contains only one to five loci and is classified into one of five major haplotypes. To elucidate the origin of human DRB loci and haplotypes, intron 1 sequences approximately 5000 bp in length were determined for three DRB1 alleles (DRB1*03, DRB1*04, and DRB1*15) and five DRB genes (DRB2, DRB3, DRB4, DRB5, and DRB7). The sequences were subjected to phylogenetic analyses together with previously determined intron 4 and 5 sequences. The sequences provided two sources of information: Nucleotide substitutions that could be used to construct phylogenetic trees and to estimate divergence times and a set of insertions (mostly Alu elements) that reveal the order of splitting of duplicated genes. The combined data indicate that the ancestor of the human DRB genes was HLA-DRB1*04-like and that the DRB2, DRB7, DRB5, and DRB3 genes arose from this ancestor by four rounds of duplication 58, 56, 53, and 36 million years (MY) ago, respectively. The DRB4 gene may have arisen 46 MY ago by a deletion from the DRB1 and DRB2 genes and the DRB6 gene is probably an allele at the DRB2 locus. During the course of its evolution, the DRB1*04 gene acquired an intron 1 segment (including two Alu elements) from a gene that became the ancestor of DRB1*03. The present-day HLA-DR haplotypes were derived from three principal ancestral haplotypes: DRB1-DRB2, DRB1-DRB5, and DRB1-DRB7.

Evolutionary relationship of HLA-DRB genes inferred from intron sequences

The major histocompatibility complex (Mhc) consists of class I and class II genes. In the human Mhc (HLA) class II genes, nine DRB loci have been identified. To elucidate the origin of these duplicated loci and allelic divergences at the most polymorphic DRB1 locus, introns 4 and 5 as well as the 3' untranslated region (altogether approximately 1,000 base pairs) of seven HLA-DRB loci, three HLA-DRB1 alleles, and nine nonhuman primate DRB genes were examined. It is shown that there were two major diversification events in HLA-DRB genes, each involving gene duplications and allelic divergences. Approximately 50 million years (my) ago, DRB1*04 and an ancestor of the DRB1*03 cluster (DRB1*03, DRB1*15, and DRB3) diverged from each other and DRB5, DRB7, DRB8, and an ancestor of the DRB2 cluster (DRB2, DRB4, and DRB6) arose by gene duplication. Later, about 25 my ago, DRB1*15 diverged from DRB1*03, and DRB3 was duplicated from DRB1*03. Then, some 20 my ago, the lineage leading to the DRB2 cluster produced two new loci, DRB4 and DRB6. The DRB1*03 and DRB1*04 allelic lineages are extraordinarily old and have persisted longer than some duplicated genes. The orthologous relationships of DRB genes between human and Old World monkeys are apparent, but those between Catarrhini and New World monkeys are equivocal because of a rather rapid expansion and contraction of primate DRB genes by duplication and deletion.

A complement factor B-like cDNA clone from the zebrafish (Brachydanio rerio)

An important molecule in the activation of the complement system in vertebrates is factor B, a serine protease with a molecular mass of 95,000. Factor B and the complement component C2 are thought to have arisen by gene duplication. In mammals and in Xenopus the factor B gene is linked to the major histocompatibility complex (MHC), whereas in domestic fowl it segregates independently of the MHC. Here we describe the isolation of a cDNA clone coding for factor B in the zebrafish, Brachydanio rerio. The deduced protein sequence exhibits a characteristic mosaic structure consisting of the short consensus repeat (SCR), the von Willebrand factor, and the serine protease domains. The estimated time of factor B and C2 divergence (approximately 350 million years ago), combined with the fact that C2 has thus far been found only in mammals, suggest that the factor B-C2 gene duplication occurred after the divergence of mammal-like reptiles from other reptiles and hence also birds. After the duplication, the C2 component evolved significantly faster than factor B.

Mhc-DRB genes and the origin of New World monkeys

The major histocompatibility complex (Mhc) is a family of loci characterized by its relatively rapid evolutionary turnover, large genetic distances between genes, and long persistence of allelic lineages effected by balancing selection. These features render the Mhc highly suitable for answering questions concerning speciation and adaptive radiation. The aim of the present study was to use Mhc-DRB genes to make inferences about the founding population of the Platyrrhini. Three segments, each approximately 300 base pairs in length, of the platyrrhine DRB genes were amplified by the polymerase chain reaction and sequenced. The segments were derived from intron 2, exon 3, and exon 6 of DRB genes from different species of New World monkeys. The results of the study have revealed that on a phylogenetic tree, all of the tested platyrrhine genes appear to form a single cluster, while all catarrhine DRB genes form a distinct cluster, although the bootstrap values fail to provide statistically significant support for the separation of these two clades. This observation suggests that the multiple platyrrhine genes originated from a single ancestral gene after the divergence of the Platyrrhini and Catarrhini and thus contradicts the results of an earlier study in which some exon 2 DRB sequences appeared to predate the split of the two primate groups. The inconsistency in the DRB gene phylogeny can be explained by postulating convergent evolution for the peptide-binding region of the DRB exon 2 sequences. The phylogeny of the platyrrhine DRB genes (except for exon 2) is relatively "shallow"; the distances between genes are relatively short (in comparison to the catarrhine DRB genes), and there is a tendency for sequences of individual species to cluster together. The phylogeny of the platyrrhine DRB genes is consistent with the postulate that a small population founded the group and that there is an ongoing adaptive radiation from small, relatively isolated founding populations.

Phylogenetic analysis of cichlid fishes using nuclear DNA markers

The recent explosive adaptive radiation of cichlids in the great lakes of Africa has attracted the attention of both morphologists and molecular biologists. To decipher the phylogenetic relationships among the various taxa within the family Cichlidae is a prerequisite for answering some fundamental questions about the nature of the speciation process. In the present study, we used the random amplification of polymorphic DNA (RAPD) technique to obtain sequence differences between selected cichlid species. We then designed specific primers based on these sequences and used them to amplify template DNA from a large number of species by the polymerase chain reaction (PCR). We sequenced the amplified products and searched the sequences for indels and shared substitutions. We identified a number of such characters at three loci--DXTU1, DXTU2, and DXTU3--and used them for phylogenetic and cladistic analysis of the relationships among the various cichlid groups. Our studies assign an outgroup position to Neotropical cichlids in relation to African cichlids, provide evidence for a sister-group relationship of tilapiines to the haplochromines, group Cyphotilapia frontosa with the lamprologines of Lake Tanganyika, place Astatoreochromis alluaudi to an outgroup position with respect to other haplochromines of Lakes Victoria and Malawi, and provide additional support for the monophyly of the remaining Lake Victoria haplochromines and the Lake Malawi haplochromines. The described approach holds great promise for further resolution of cichlid phylogeny.

A cDNA clone from the sea lamprey Petromyzon marinus coding for a scavenger receptor Cys-rich (SRCR) domain protein

Our knowledge of the immune system in the early vertebrates, the Agnatha, and the molecules involved in their immune reactions is fragmentary. By serendipity we discovered a cDNA clone in a library made from gut poly(A)+RNA of the sea lamprey, Petromyzon marinus (Pema), that translates into the SREG (SRCR-EGF, see below) protein which resembles cell-membrane proteins of mammalian immune cells. The putative translated product is a type-I integral membrane glycoprotein which contains two scavenger receptor Cys-rich (SRCR) domains flanking five epidermal growth factor (EGF)-like repeats. The two SRCR domains are closely related to CD6 (expressed on human lymphocytes), WC1 (expressed on mammalian CD4-CD8(-)-gamma delta T cells) and M130 (expressed on human macrophages). The Pema-SREG may therefore be involved in intercellular contacts and cell activation or differentiation in the immune system. It is thus a potential marker that can be used to investigate the lamprey immune system.

Inactivation of E alpha and E beta expression in inbred and wild mice by multiple distinct mutations, some of which predate speciation within Mus species

The H-2 MHC of mice encodes two functional class II heterodimeric proteins: A alpha A beta (A) and E alpha E beta (E). While failure to express the A protein has not been reported, a significant proportion of of H-2 haplotypes in both inbred and wild mice do not express E proteins. We and others have previously characterized the molecular basis for defective E expression in haplotypes from Mus domesticus (b, f, q, s, from inbred strains) and M. castaneus (w17, wild-derived) species, identifying six distinct defects in the genes for E alpha or E beta. In this report we have extended these studies to other E- haplotypes, including several from t-haplotype-bearing M. domesticus mice (w29, w57, w302) and one derived from the Asian species M. bactrianus (w301). Analyses at the protein, RNA and DNA levels were employed to identify the defects in the genes for Ea and Eb. At least one new defect was identified that prevents E beta expression in a t-associated H-2 haplotypes (w57), bringing the number of distinct mutations causing the E- phenotype to seven. Another t-associated haplotype, w302, was found to share the same E beta defect with mice of the inbred q haplotype and of the w17 haplotype from M. castaneus, while its Ea gene contains the deletion carried also by the inbred b and s haplotypes and by a number of wild haplotypes. The mutations in the Ea and Eb genes of the w301 haplotype from M. bactrianus were found to be identical to those of the inbred f haplotype. This indicates that the origin of the mutations in the Eb genes of the q, w17 and w302 haplotypes and in the Ea and Eb genes of the f and w301 haplotypes, predated speciation within Mus, thought to have occurred approximately 0.35-1 million years ago. Their maintenance in mouse populations suggests that in certain conditions the failure to express E alpha E beta proteins may be advantageous and selected for.