Evolutionary relationships among the primate Mhc-DQA1 and DQA2 alleles

The variation of the Mhc-DQA1 and DQA2 loci of ten different primate species (hominoids and Old World monkeys) was studied in order to obtain an insight in the processes that generate polymorphism of major histocompatibility complex (Mhc) class II genes and to establish the evolutionary relationships of their alleles. To that end nucleotide sequences of 36 Mhc class II DQA1 and seven DQA2 second exons were determined and phylogenetic trees that illustrate their evolutionary relationships were constructed. We demonstrate the existence of four primate Mhc-DQA1 allele lineages, two of which probably existed before the separation of the ancestors of the hominoids and Old World monkeys (approximately 22–28 million years ago). Mhc-DQA2 sequences were found only in the hominoid species and showed little diversity. We found no evidence for a major contribution of recombinational events to the generation of allelic diversity of the primate Mhc-DQA1 locus. Instead, our data suggest that the primate Mhc-DQA1 and DQA2 loci are relatively stable entities that mutated primarily as a result of point mutations.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M76186-M76229.     

猕猴MHC-DPB1基因外显子2多态性研究

猕猴(Macaca mulatta)是最理想的医学实验灵长类动物, 且为国家二级保护动物。为了解中国猕猴主要组织相容复合体(Major histocompatibility complex, MHC)基因的遗传多态性背景, 为它们在生物医学研究中的应用及其遗传资源的保护提供一定的科学依据, 文章采用变性梯度凝胶电泳(Denaturing gradient gel electrophoresis, DGGE)和克隆测序技术分析了106个四川野生猕猴MHC-DPB1基因的exon 2, 共检测到21个Mamu-DPB1等位基因, 其中有15个为本研究中首次发现的新等位基因; 从整个大的猕猴群体(106个个体)来看, 等位基因频率最高的是Mamu-DPB1*30(0.1120); 单独从不同地理群体来看, 最高等位基因频率分别为: 小金-DPB1*30 (0.1120), 黑水-DPB1*04 (0.1702), 巴中-DPB1*32 (0.1613), 汉源-DPB1*30(0.1120), 九龙-DPB1*04(0.1139); 氨基酸序列比对发现, 猕猴Mamu-DPB1等位基因编码的氨基酸序列中, 有12个氨基酸残基变异位点表现出物种特异性, 其中有9个位于新发现的15个Mamu-DPB1等位基因氨基酸序列中; 不同物种来源的DPB1等位基因系统发生树表明, 猕猴与其近缘物种食蟹猴(Macaca fascicularis)的DPB1等位基因间存在着跨种多态(Trans-species polymorphism)现象。研究还表明, MHC-DPB1等位基因在中国猕猴群体和先前为主要研究对象的印度猕猴群体间具有较大的差异。     

Molecular study of Mhc-DRB in wild chacma baboons reveals high variability and evidence for trans-species inheritance

The MHC class II genes of many primate species were investigated extensively in recent years. However, while Mhc-DRB genes were studied in Old World monkeys such as rhesus macaques, the Mhc-DRB of baboons was only studied in a limited way. Because of their close anatomical and physiological relationship to humans, baboons are often used as models for reproduction and transplantation research. Baboons are also studied as a model species in behavioural ecology. Thus, identification of MHC genes would provide a foundation for studies of Mhc, biology and behaviour. Here, we describe the use of PCR, cloning, denaturing gradient gel electrophoresis (DGGE) and sequencing to identify Mhc-DRB sequences in wild chacma baboons (Papio ursinus). We amplified the highly variable second exon of baboon Mhc-DRB sequences using generic DRB primers. To validate and optimize the DGGE protocol, four DNA samples were initially studied using cloning and sequencing. Clones were screened using a novel RFLP approach to increase the number of clones identified for each individual. Results from cloning and sequencing were used to optimise DGGE conditions for Mhc-DRB genotyping of the remaining study subjects. Using these techniques, we identified 16 Paur-DRB sequences from 30 chacma baboons. On the basis of phylogenetic tree analyses, representatives of the Mhc-DRB1 and Mhc-DRB5 loci, and 13 different DRB lineages were identified. Evidence for trans-species inheritance of some Mhc-DRB sequences comes from high identity between the new Paur-DRB sequences and sequences from Papio cynocephalus, Macaca mulatta and possibly Galago moholi.Nucleotide sequence data reported are available in the GenBank/EMBL/DDBJ databases under the accession numbers DQ339722–DQ339737.     

High levels of diversity characterize mandrill (Mandrillus sphinx) Mhc-DRB sequences

The major histocompatibility complex (MHC) is highly polymorphic in most primate species studied thus far. The rhesus macaque (Macaca mulatta) has been studied extensively and the Mhc-DRB region demonstrates variability similar to humans. The extent of MHC diversity is relatively unknown for other Old World monkeys (OWM), especially among genera other than Macaca. A molecular survey of the Mhc-DRB region in mandrills (Mandrillus sphinx) revealed extensive variability, suggesting that other OWMs may also possess high levels of Mhc-DRB polymorphism. In the present study, 33 Mhc-DRB loci were identified from only 13 animals. Eleven were wild-born and presumed to be unrelated and two were captive-born twins. Two to seven different sequences were identified for each individual, suggesting that some mandrills may have as many as four Mhc-DRB loci on a single haplotype. From these sequences, representatives of at least six Mhc-DRB loci or lineages were identified. As observed in other primates, some new lineages may have arisen through the process of gene conversion. These findings indicate that mandrills have Mhc-DRB diversity not unlike rhesus macaques and humans.     

Characterization of the rhesus macaque (Macaca mulatta) equivalent of HLA-F

Nucleotide sequence analysis of rhesus macaque major histocompatibility complex class I cDNAs allowed the identification of the orthologue of HLA-F, designated Mamu-F. Comparison of Mamu-F with earlier published human and chimpanzee orthologues demonstrated that these sequences share a high degree of similarity, both at the nucleotide and amino acid level, whereas a New World monkey (cotton-top tamarin) equivalent is more distantly related. Exon 7, encoding one of the cytoplasmatic domains, is absent for all primate Mhc-F cDNA sequences analyzed so far. In contrast to the human, chimpanzee, and rhesus macaque equivalents, the cotton-top tamarin Saoe-F gene seems to have accumulated far more nonsynomynous than synonymous differences.The nucleotide sequence data reported in this paper have been submitted to the Genbank nucleotide sequence database and have been assigned the accession number Z 21819.     

Comprehensive identification of MHC class II alleles in a cohort of Chinese rhesus macaques

Rhesus macaque is a very important animal model for various human diseases, especially for AIDS and vaccine research. The susceptibility and/or resistance to some of these diseases are related to the major histocompatibility complex (MHC). To gain insight into the MHC background and to facilitate the experimental use of Chinese rhesus macaques, Mamu-DPB1, Mamu-DQB1, and Mamu-DRB alleles were investigated in 30 Chinese rhesus macaques through gene cloning and sequencing. A total of 66 alleles were identified in this study, including 14 Mamu-DPB1, 20 Mamu-DQB1, and 30 Mamu-DRB alleles as well as 2 high-frequency Mamu-DPB1 alleles. Interestingly, one of the high-frequency Mamu-DPB1 alleles had been undocumented in earlier studies. Eleven of the other alleles, including four Mamu-DPB1, three Mamu-DQB1, and four Mamu-DRB alleles were also novel. Importantly, like MHC-DRB, more than two Mamu-DPB1 sequences per animal were detected in 13 monkeys, which suggested that they might represent gene duplication. Our data also indicated quite a few differences in the distribution of MHC class II alleles between the Chinese rhesus macaques and the previously reported Indian rhesus macaques. To our knowledge, our results revealed comprehensively the combination of MHC II alleles. This information will not only promote the understanding of Chinese rhesus macaque MHC polymorphism but will also facilitate the use of Chinese rhesus macaques in studies of human disease.     

The sheep CD1 gene family contains at least four CD1B homologues

We identified four cDNA sequences encoding sheep homologues of the CD1 molecule. The sheep sequences were selected from λgt11 thymocyte cDNA libraries by hybridization with a humanCD1C probe and a homologous sheep probe. TheSCD1B-42 andSCD1A25 sequences encode complete CD1 molecules. The third sequence,SCD1B-52, which is closely related toSCD1B-42 and may be an allele, has the sequence encoding the α3 region precisely deleted. The fourth sequence,SCD1T10, is truncated at the 5′ end. All four sequences are related to the humanCD1B and domestic rabbitCD1B-like sequences at both nucleotide and amino acid level. Comparison of the derivedCD1 amino acid sequences with the sequence of major histocompatibility complex class I molecules showed that the sheep CD1 molecules, like human CD1 molecules, lack most of the conserved class I residues known to be involved in interaction with 132-microglobulin and the CD8 molecule. They do not contain the peptide docking residues involved in anchoring peptides in the peptide binding groove of class I molecules. Southern hybridization of sheep DNA with a sheepCD1 exon 4/ga3 probe showed that the sheep genome encodes at least sevenCD1 genes. The implications of these analyses for CD1 function are discussed. The nucleotide sequence data reported in this paper have been submitted to the EMBL/GenBank nucleotide sequence databases and have been assigned the accession numbers Z36890 (SCD1A25), Z36891 (SCD1B-42), Z36892 (SCD1B-52), and X90567 (SCDIT10)     

Identification of major histocompatibility complex genes in the guppy,Poecilia reticulata

The guppy, Poecilia reticulata, a teleostean fish of the order Cyprinodontiformes, has been used extensively in studies of host-parasite interactions, courtship behavior, and mating preference, as well as in ecological and evolutionary genetics. A related species was among the first poikilotherm vertebrates to be used in the study of histocompatibility genes. All these studies could benefit from the identification and characterization of the guppy major histocompatibility complex (Mhc) genes. Here, both class I and class II genes of the guppy are described. The number of expressed loci, as determined by representation of clones in a cDNA library, sequencing, and Southern blot analysis, may be low in both Mhc classes: combined evidence suggests that there may be one expressed class II locus only and one or two expressed class I loci. The variability of aquaristic guppy stocks is very low: only three and two genes have been detected at the class I and class II loci, respectively, in the stocks examined. This genetic paucity is most likely the consequence of breeding practices employed by aquarists and commercial establishments. Limited sampling of wild guppy populations revealed extensive Mhc polymorphism at loci of both classes in nature. Comparison of guppy Mhc sequences with those of other vertebrates has revealed the existence of a set of insertions/deletions which can be used as characters in cladistic analysis to infer phylogenetic relationships among vertebrate taxa and the Mhc genes themselves. These indels are particularly frequent in the regions coding for the loops of 1 and 2 domains of class I proteins.The nucleotide sequence data reported in this publication have been submitted to the EMBL nucleotide sequence database and have been assigned the accession numbers Z54076-Z54095     

Analysis of the sequence variations in the Mhc DRB1-like gene of the endangered Humboldt penguin (Spheniscus humboldti)

The Major Histocompatibility Complex (Mhc) genomic region of many vertebrates is known to contain at least one highly polymorphic class II gene that is homologous in sequence to one or other of the human Mhc DRB1 class II genes. The diversity of the avian Mhc class II gene sequences have been extensively studied in chickens, quails, and some songbirds, but have been largely ignored in the oceanic birds, including the flightless penguins. We have previously reported that several penguin species have a high degree of polymorphism on exon 2 of the Mhc class II DRB1-like gene. In this study, we present for the first time the complete nucleotide sequences of exon 2, intron 2, and exon 3 of the DRB1-like gene of 20 Humboldt penguins, a species that is presently vulnerable to the dangers of extinction. The Humboldt DRB1-like nucleotide and amino acid sequences reveal at least eight unique alleles. Phylogenetic analysis of all the available avian DRB-like sequences showed that, of five penguin species and nine other bird species, the sequences of the Humboldt penguins grouped most closely to the Little penguin and the mallard, respectively. The present analysis confirms that the sequence variations of the Mhc class II gene, DRB1, are useful for discriminating among individuals within the same penguin population as well those within different penguin population groups and species.The nucleotide sequence and amino acid sequence data reported in this paper have been submitted to the DDBJ database and have been assigned the accession numbers AB088371–AB088374, AB089199, AB154393–AB154399, and AB162144.     

Mapping and nucleotide sequence of a newHLA class II light chain gene,DQB3

A genomic clone specifying a new HLA class II antigen β chain,DQB3, was isolated from a human genomic phage library using aDQB1 cDNA probe under low stringency conditions. Southern hybridization and nucleotide sequence analyses identified the β2 domain exon (exon 3) with several deleterious mutations and the CP-TM-CY exon [connecting peptide, transmembrane, and cytoplasmic regions, (exon 4)], but the first, second, and fifth exons encoding the 5′ UT-leader, the β1 domain, and the 3′ UT domain of normal β chains, respectively, were entirely missing. The nucleotide sequences of these two exons were distinct from those of other class II β chain genes, but slightly more related to theDQB1 andDQB2 genes than to other class II genes. TheDQB3 sequence mapped betweenDQA2 andDQB1, 15 kb upstream fromDQA2, by analysis of overlapping cosmid clones. This mapping was supported by the fact thatTaq I,Msp I, andBam HIDQB3 polymorphisms were perfectly correlated with theDQA2 polymorphism and not with any polymorphisms in theDR orDQ subregion, suggesting the presence of a hot spot for recombination betweenDQB3 andDQB1. The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M26577.     

Detection of novel sequence heterogeneity and haplotypic diversity of HLA class II genes

Nucleic acid sequences of the second exons of HLA-DRB1, –DRB3/4/5, –DQB1, and –DQA1 genes were determined from 43 homozygous cell lines, representing each of the known class II haplotypes, and from 30 unrelated Caucasian subjects, comprising 60 haplotypes. This systematic sequence analysis was undertaken in order to a) determine the existence of sequence microheterogeneity among cell lines which type as identical by methods other than sequencing; b) determine whether direct sequencing of class II genes will identify the presence of more extensive sequence polymorphism at the population level than that identified with other typing methods; c) accurately determine the molecular composition of the known class II haplotypes; and d) study their evolutionary relatedness by maximum parsimony analysis. The identification of seven previously unidentified haplotypes carrying five new allelic amino acid sequences suggests that sequence microheterogeneity at the population level may be more frequent than previously thought. Maximum parsimony analysis of these haplotypes allowed their evolutionary classification and indicates that the higher mutation rate at DRB1 compared to DQB1 loci in most haplotypic groups is inversed in specific haplotype lineages. Furthermore, the extent and localization of gene conversions and point mutations at class II loci in the evolution of these haplotypes is significantly different at each locus. Identification of additional HLA class II molecular microheterogeneity suggests that direct sequence analysis of class II HLA genes can uncover new allelic sequences in the population and may represent a useful alternative to current typing methodologies to study the effects of sequence allelism in organ transplantation.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M35890 through M35953.     

Extensive sharing of MHC class II alleles between rhesus and cynomolgus macaques

In contrast to rhesus monkeys, substantial knowledge on cynomolgus monkey major histocompatibility complex (MHC) class II haplotypes is lacking. Therefore, 17 animals, including one pedigreed family, were thoroughly characterized for polymorphic Mhc class II region genes as well as their mitochondrial DNA (mtDNA) sequences. Different cynomolgus macaque populations appear to exhibit unique mtDNA profiles reflecting their geographic origin. Within the present panel, 10 Mafa-DPB1, 14 Mafa-DQA1, 12 Mafa-DQB1, and 35 Mafa-DRB exon 2 sequences were identified. All of these alleles cluster into lineages that were previously described for rhesus macaques. Moreover, about half of the Mafa-DPB1, Mafa-DQA1, and Mafa-DQB1 alleles and one third of the Mafa-DRB exon 2 sequences are identical to rhesus macaque orthologues. Such a high level of Mhc class II allele sharing has not been reported for primate species. Pedigree analysis allowed the characterization of nine distinct Mafa class II haplotypes, and seven additional ones could be deduced. Two of these haplotypes harbor a duplication of the Mafa-DQB1 locus. Despite extensive allele sharing, rhesus and cynomolgus monkeys do not appear to possess identical Mhc class II haplotypes, thus illustrating that new haplotypes were generated after speciation by recombination-like processes.     

Characterization of the major histocompatibility complex class II <Emphasis Type="Italic">DOB,DPB1</Emphasis>, and <Emphasis Type="Italic">DQB1</Emphasis> alleles in cynomolgus macaques of Vietnamese origin

Major histocompatibility complex (MHC) molecules play an important role in the susceptibility and/or resistance to many diseases. To gain an insight into the MHC background and to facilitate the experimental use of cynomolgus macaques, the second exon of the MhcMafa-DOB, -DPB1, and -DQB1 genes from 143 cynomolgus macaques were characterized by cloning to sequencing. A total of 16 Mafa-DOB, 16 Mafa-DPB1, and 34 Mafa-DQB1 alleles were identified, which revealed limited, moderate, and marked allelic polymorphism at DOB, DPB1, and DQB1, respectively, in a cohort of cynomolgus macaques of Vietnamese origin. In addition, 16 Mafa-DOB, 5 Mafa-DPB1, and 8 Mafa-DQB1 alleles represented novel sequences that had not been reported in earlier studies. Almost of the sequences detected at the DOB and DQB1 locus in the present study belonged to DOB*01 (100%) and DQB1*06 (62%) lineages, respectively. Interestingly, four, three, and one high-frequency alleles were detected at Mafa-DOB, -DPB1, and -DQB1, respectively, in this monkeys. The alleles with the highest frequency among these monkeys were Mafa-DOB*010102, Mafa-DPB1*13, and Mafa-DQB1*0616, and these were found in 33 (25.6%) of 129 monkeys, 32 (31.37%) of 102 monkeys, and 30 (31%) of 143 monkeys, respectively. The high-frequency alleles may represent high priority targets for additional characterization of immune function. We also carried out evolutionary and population analyses using these sequences to reveal population-specific alleles. This information will not only promote the understanding of MHC diversity and polymorphism in the cynomolgus macaque but will also increase the value of this species as a model for biomedical research.     

Major histocompatibility complex class II A gene polymorphism in the striped bass

Adaptions of the polymerase chain reaction were used to isolate cDNA sequences encoding the Major histocompatibility complex(Mhc) class II A gene(s) of the striped bass (Morone saxatilis). Four complete Mhc class II A genes were cloned and sequenced from a specimen originating in the Roanoke River, North Carolina, and another three A genes from a specimen originating from the Santee-Cooper Reservoir, South Carolina, identifying a total of seven unique sequences. The sequence suggests the presence of at least two Mhc class II A loci. The extensive sequence variability observed between the seven different Mhc class II clones was concentrated in the 1 encoding domain. The encoded 2, transmembrane, and cytoplasmic regions of all seven striped genes correlated well with those of known vertebrate Mhc class II proteins. Overall, the striped bass sequences showed greatest similarity to the Mhc class II A genes of the zebrafish. Southern blot analysis demonstrated extensive polymorphism in the Mhc class II A genes in members of a Roanoke river-caught population of striped bass versus a lesser degree of polymorphism in an aquacultured Santee-Cooper population of striped bass.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers (Mosa-A-S5) L35062, (Mosa-A-S8) L35066, (Mosa-A-R7) L35067, and (Mosa-A-S7) L35072 L35066, (Mossa-A-R7) L35067, and (Mosa-A-S7) L35072     

C4 gene polymorphism in primates: evolution,generation, and Chido and Rodgers antigenicity

Eleven new C4d genomic primate sequences of the fourth complement factor (C4) have been obtained. Seven of them belong to five species not yet explored for this gene: Pan paniscus (pygmy chimpanzee), Cercopithecus aethiops (green monkey), Macaca mulatta (rhesus monkey), Macaca fascicularis (cynomolgus), and Saguinus oedipus (cotton top tamarin). The New World monkeys (tamarins, four individuals) sequenced for C4 have a single C4d sequence only, which shows a B isotypic specificity and a Rodgers 3 (Rg3), Chido 1 (Ch1) antigenicity. Rg3 and Ch1 could thus be the oldest Rg/Ch specificity (at least 50 million years old) and Rg1, Rg2, Ch3, and Ch6 could be more recent human-specific antigens. Mechanisms of C4d polymorphism generation were analyzed by compiling all the presently available sequences. Examples of both point mutations and crossing-over events among C4d primate sequences could be detected. The problem of a possible trans-species inheritance of C4d polymorphism was addressed and two apparently contradicting dendrograms were obtained. One of them, constructed by using both exon and intron sequences, does not support trans-species evolution, but supports the proposed theory of extensive homogenization of the C4 genes occurring within each species, because alleles from each primate species cluster together. Another completely different dendrogram, obtained by using exon sequences only, suggests the existence of trans-species evolution for C4d polymorphism, because alleles belonging to different species cluster together in a way similar to that found for HLA class I or II alleles. However, orangutan sequences group together in both kinds of C4d sequence dendrograms and seem to have arisen from an ancestor different from that of chimpanzee, gorilla and man C4d sequences. Finally, further data have been obtained that support trans-species conservation of A-ness and B-ness and the existence of trans-specifically conserved allelic motifs, both in intronic and exonic sequences.     

Ancestral Population Sizes and Species Divergence Times in the Primate Lineage on the Basis of Intron and BAC End Sequences

The effective sizes of ancestral populations and species divergence times of six primate species (humans, chimpanzees, gorillas, orangutans, and representatives of Old World monkeys and New World monkeys) are estimated by applying the two-species maximum likelihood (ML) method to intron sequences of 20 different loci. Examination of rate heterogeneity of nucleotide substitutions and intragenic recombination identifies five outrageous loci (ODC1, GHR, HBE, INS, and HBG). The estimated ancestral polymorphism ranges from 0.21 to 0.96% at major divergences in primate evolution. One exceptionally low polymorphism occurs when African and Asian apes diverged. However, taking into consideration the possible short generation times in primate ancestors, it is concluded that the ancestral population size in the primate lineage was no smaller than that of extant humans. Furthermore, under the assumption of 6 million years (myr) divergence between humans and chimpanzees, the divergence time of humans from gorillas, orangutans, Old World monkeys, and New World monkeys is estimated as 7.2, 18, 34, and 65 myr ago, respectively, which are generally older than traditional estimates. Beside the intron sequences, three other data sets of orthologous sequences are used between the human and the chimpanzee comparison. The ML application to these data sets including 58,156 random BAC end sequences (BES) shows that the nucleotide substitution rate is as low as 0.6–0.8 × 10^–9 per site per year and the extent of ancestral polymorphism is 0.33–0.51%. With such a low substitution rate and short generation time, the relatively high extent of polymorphism suggests a fairly large effective population size in the ancestral lineage common to humans and chimpanzees.[Reviewing Editor: Dr. Magnus Nordborg]     

Maintenance of DQB1 polymorphisms in primates.

To understand the evolution of the class II major histocompatibility complex (MHC) DQB1 locus in primates, the second exons of seven DQB1 alleles from five non-human primate species were amplified by polymerase chain reaction. Comparisons of these and other primate sequences show that no between-species diversity is greater than within-species diversity, suggesting maintenance of DQB1 alleles through the history of Old-World primates. There is a preponderance of nonsynonymous nucleotide substitutions at antigen-binding-site codons; this pattern is in marked contrast to what is seen at the closely related, presumably nonfunctional DQB2 gene. The results support the hypothesis that DQB1 polymorphism is maintained by overdominant selection relating to antigen presentation.     

Molecular diversity at the self-incompatibility locus is a salient feature in natural populations of wild tomato (Lycopersicon peruvianum)

A cDNA encoding a stylar protein was cloned from flowers of self-incompatible wild tomato (Lycopersicon peruvianum). The corresponding gene was mapped to the S locus, which is responsible for self-incompatibility. The nucleotide sequence was determined for this allele, and compared to other S-related sequences in the Solanaceae. The S allele was used to probe DNA from 92 plants comprising 10 natural populations of Lycopersicon peruvianum. Hybridization was conducted under moderate and permissive stringencies in order to detect homologous sequences. Few alleles were detected, even under permissive conditions, underscoring the great sequence diversity at this locus. Those alleles that were detected are highly homologous. Sequences could not be detected in self-incompatible Nicotiana alata, self-compatible L. esculentum (cultivated tomato) or self-compatible L. hirsutum. However, hybridization to an individual of self-incompatible L. hirsutum revealed a closely related sequence that maps to the S locus in this reproductively isolated species. This supports the finding that S locus polymorphism predates speciation. The extraordinarily high degree of sequence diversity present in the gametophytic self-incompatibility system is discussed in the context of other highly divergent systems representing several kingdoms.