High homologous gene conservation despite extreme autopolyploid redundancy in sugarcane

Modern sugarcane (Saccharum spp.) is the leading sugar crop and a primary energy crop. It has the highest level of 'vertical' redundancy (2n=12x=120) of all polyploid plants studied to date. It was produced about a century ago through hybridization between two autopolyploid species, namely S. officinarum and S. spontaneum. In order to investigate the genome dynamics in this highly polyploid context, we sequenced and compared seven hom(oe)ologous haplotypes (bacterial artificial chromosome clones). Our analysis revealed a high level of gene retention and colinearity, as well as high gene structure and sequence conservation, with an average sequence divergence of 4% for exons. Remarkably, all of the hom(oe)ologous genes were predicted as being functional (except for one gene fragment) and showed signs of evolving under purifying selection, with the exception of genes within segmental duplications. By contrast, transposable elements displayed a general absence of colinearity among hom(oe)ologous haplotypes and appeared to have undergone dynamic expansion in Saccharum, compared with sorghum, its close relative in the Andropogonea tribe. These results reinforce the general trend emerging from recent studies indicating the diverse and nuanced effect of polyploidy on genome dynamics.     

Construction of a composite sorghum genome map and comparison with sugarcane, a related complex polyploid

 A sorghum composite linkage map was constructed with two recombinant inbred line populations using heterologous probes already mapped on maize and sugarcane. This map includes 199 loci revealed by 188 probes and distributed on 13 linkage groups. A comparison based on 84 common probes was performed between the sorghum composite map and a map of a sugarcane (Saccharum spp.) cultivar being developed and presently comprising 10 tentative linkage groups. A straight synteny was observed for 2 pairs of linkage groups; in two cases, 1 sorghum linkage group corresponded to 2 or 3 sugarcane linkage groups, respectively; in two cases 1 sugarcane link- age group corresponded to 2 separate sorghum linkage groups; for 2 sorghum linkage groups, no complete correspondance was found in the sugarcane genome. In most cases loci appeared to be colinear between homoeologous chromosomal segments in sorghum and sugarcane. These results are discussed in relation to published data on sorghum genomic maps, with specific reference to the genetic organization of sugarcane cultivars, and they, illustrate how investigations on relatively simple diploid genomes as sorghum will facilitate the mapping of related polyploid species such as sugarcane. Received: 12 August 1996 / Accepted: 30 August 1996     

Diploid/polyploid syntenic shuttle mapping and haplotype-specific chromosome walking toward a rust resistance gene (Bru1) in highly polyploid sugarcane (2n approximately 12x approximately 115)

The genome of modern sugarcane cultivars is highly polyploid ( approximately 12x), aneuploid, of interspecific origin, and contains 10 Gb of DNA. Its size and complexity represent a major challenge for the isolation of agronomically important genes. Here we report on the first attempt to isolate a gene from sugarcane by map-based cloning, targeting a durable major rust resistance gene (Bru1). We describe the genomic strategies that we have developed to overcome constraints associated with high polyploidy in the successive steps of map-based cloning approaches, including diploid/polyploid syntenic shuttle mapping with two model diploid species (sorghum and rice) and haplotype-specific chromosome walking. Their applications allowed us (i) to develop a high-resolution map including markers at 0.28 and 0.14 cM on both sides and 13 markers cosegregating with Bru1 and (ii) to develop a physical map of the target haplotype that still includes two gaps at this stage due to the discovery of an insertion specific to this haplotype. These approaches will pave the way for the development of future map-based cloning approaches for sugarcane and other complex polyploid species.     

Spatio‐temporal patterns of genome evolution in allotetraploid species of the genus Oryza

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct ‘genome types’ (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1–Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio‐temporal dynamics of genome organization and evolution of this region in ‘natural’ polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage‐specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.     

HapTree: A Novel Bayesian Framework for Single Individual Polyplotyping Using NGS Data

As the more recent next-generation sequencing (NGS) technologies provide longer read sequences, the use of sequencing datasets for complete haplotype phasing is fast becoming a reality, allowing haplotype reconstruction of a single sequenced genome. Nearly all previous haplotype reconstruction studies have focused on diploid genomes and are rarely scalable to genomes with higher ploidy. Yet computational investigations into polyploid genomes carry great importance, impacting plant, yeast and fish genomics, as well as the studies of the evolution of modern-day eukaryotes and (epi)genetic interactions between copies of genes. In this paper, we describe a novel maximum-likelihood estimation framework, HapTree, for polyploid haplotype assembly of an individual genome using NGS read datasets. We evaluate the performance of HapTree on simulated polyploid sequencing read data modeled after Illumina sequencing technologies. For triploid and higher ploidy genomes, we demonstrate that HapTree substantially improves haplotype assembly accuracy and efficiency over the state-of-the-art; moreover, HapTree is the first scalable polyplotyping method for higher ploidy. As a proof of concept, we also test our method on real sequencing data from NA12878 (1000 Genomes Project) and evaluate the quality of assembled haplotypes with respect to trio-based diplotype annotation as the ground truth. The results indicate that HapTree significantly improves the switch accuracy within phased haplotype blocks as compared to existing haplotype assembly methods, while producing comparable minimum error correction (MEC) values. A summary of this paper appears in the proceedings of the RECOMB 2014 conference, April 2–5.     

The identification and characterisation of alleles of sucrose phosphate synthase gene family III in sugarcane

Little is known about the extent of allelic diversity of genes in the complex polyploid, sugarcane. Using sucrose phosphate synthase (SPS) Gene (SPS) Family III as an example, we have amplified and sequenced a 400 nt region from this gene from two sugarcane lines that are parents of a mapping population. Ten single nucleotide polymorphisms (SNPs) were identified within the 400 nt region of which seven were present in both lines. In the elite commercial cultivar Q165^A, 10 sequence haplotypes were identified, with four haplotypes recovered at 9% or greater frequency. Based on SNP presence, two clusters of haplotypes were observed. In IJ76-514, a Saccharum officinarum accession, 8 haplotypes were identified with 4 haplotypes recovered at 13% or greater frequency. Again, two clusters of haplotypes were observed. The results suggest that there may be two SPS Gene Family III genes per genome in sugarcane, each with different numbers of different alleles. This suggestion is supported by sequencing results in an elite parental sorghum line, 403463-2-1, in which 4 haplotypes, corresponding to two broad types, were also identified. Primers were designed to the sugarcane SNPs and screened over bulked DNA from high and low Sucrose-containing progeny from a cross between Q165^A and IJ76-514. The SNP frequency did not vary in the two bulked DNA samples, suggesting that these SNPs from this SPS gene family are not associated with variation in sucrose content. Using an ecotilling approach, two of the SPS Gene Family III haplotypes were mapped to two different linkage groups in homology group 1 in Q165^A. Both haplotypes mapped near QTLs for increased sucrose content but were not themselves associated with any sugar-related trait.     

ConPADE: Genome Assembly Ploidy Estimation from Next-Generation Sequencing Data

As a result of improvements in genome assembly algorithms and the ever decreasing costs of high-throughput sequencing technologies, new high quality draft genome sequences are published at a striking pace. With well-established methodologies, larger and more complex genomes are being tackled, including polyploid plant genomes. Given the similarity between multiple copies of a basic genome in polyploid individuals, assembly of such data usually results in collapsed contigs that represent a variable number of homoeologous genomic regions. Unfortunately, such collapse is often not ideal, as keeping contigs separate can lead both to improved assembly and also insights about how haplotypes influence phenotype. Here, we describe a first step in avoiding inappropriate collapse during assembly. In particular, we describe ConPADE (Contig Ploidy and Allele Dosage Estimation), a probabilistic method that estimates the ploidy of any given contig/scaffold based on its allele proportions. In the process, we report findings regarding errors in sequencing. The method can be used for whole genome shotgun (WGS) sequencing data. We also show applicability of the method for variant calling and allele dosage estimation. Results for simulated and real datasets are discussed and provide evidence that ConPADE performs well as long as enough sequencing coverage is available, or the true contig ploidy is low. We show that ConPADE may also be used for related applications, such as the identification of duplicated genes in fragmented assemblies, although refinements are needed.     

Development of microsatellite markers in autopolyploid sugarcane and comparative analysis of conserved microsatellites in sorghum and sugarcane

Sugarcane has become an increasingly important first-generation biofuel crop in tropical and subtropical regions. It has a large, complex, polyploid genome that has hindered the progress of genomic research and marker-assisted selection. Genetic mapping and ultimately genome sequence assembly require a large number of DNA markers. Simple sequence repeats (SSRs) are widely used in genetic mapping because of their abundance, high rates of polymorphism, and ease of use. The objectives of this study were to develop SSR markers for construction of a saturated genetic map and to characterize the frequency and distribution of SSRs in a polyploid genome. SSR markers were mined from expressed sequence tag (EST), reduced representation library genomic sequences, and bacterial artificial chromosome (BAC) sequences. A total of 5,675 SSR markers were surveyed in a segregating population. The overall successful amplification and polymorphic rates were 87.9 and 16.4%, respectively. The trinucleotide repeat motifs were most abundant, with tri- and hexanucleotide motifs being the most abundant for the ESTs. BAC and genomic SSRs were mostly AT-rich while the ESTs were relatively GC-rich due to codon bias. These markers were also aligned to the sorghum genome, resulting in 1,203 markers mapped in the sorghum genome. This set of SSRs conserved in sugarcane and sorghum would be the most informative for mapping quantitative trait loci in sugarcane and for comparative genomic analyses. This large collection of SSR markers is a valuable resource for sugarcane genomic research and crop improvement.     

A genome-wide BAC end-sequence survey of sugarcane elucidates genome composition, and identifies BACs covering much of the euchromatin

BAC-end sequences (BESs) of hybrid sugarcane cultivar R570 are presented. A total of 66,990 informative BESs were obtained from 43,874 BAC clones. Similarity search using a variety of public databases revealed that 13.5 and 42.8 % of BESs match known gene-coding and repeat regions, respectively. That 11.7 % of BESs are still unmatched to any nucleotide sequences in the current public databases despite the fact that a close relative, sorghum, is fully sequenced, indicates that there may be many sugarcane-specific or lineage-specific sequences. We found 1,742 simple sequence repeat motifs in 1,585 BESs, spanning 27,383 bp in length. As simple sequence repeat markers derived from BESs have some advantages over randomly generated markers, these may be particularly useful for comparing BAC-based physical maps with genetic maps. BES and overgo hybridization information was used for anchoring sugarcane BAC clones to the sorghum genome sequence. While sorghum and sugarcane have extensive similarity in terms of genomic structure, only 2,789 BACs (6.4 %) could be confidently anchored to the sorghum genome at the stringent threshold of having both-end information (BESs or overgos) within 300 Kb. This relatively low rate of anchoring may have been caused in part by small- or large-scale genomic rearrangements in the Saccharum genus after two rounds of whole genome duplication since its divergence from the sorghum lineage about 7.8 million years ago. Limiting consideration to only low-copy matches, 1,245 BACs were placed to 1,503 locations, covering ~198 Mb of the sorghum genome or about 78 % of the estimated 252 Mb of euchromatin. BESs and their analyses presented here may provide an early profile of the sugarcane genome as well as a basis for BAC-by-BAC sequencing of much of the basic gene set of sugarcane.     

Lr34 multi-pathogen resistance ABC transporter: molecular analysis of homoeologous and orthologous genes in hexaploid wheat and other grass species

The Triticum aestivum (bread wheat) disease resistance gene Lr34 confers durable, race non-specific protection against three fungal pathogens, and has been a highly relevant gene for wheat breeding since the green revolution. Lr34, located on chromosome 7D, encodes an ATP-binding cassette (ABC) transporter. Both wheat cultivars with and without Lr34-based resistance encode a putatively functional protein that differ by only two amino acid polymorphisms. In this study, we focused on the identification and characterization of homoeologous and orthologous Lr34 genes in hexaploid wheat and other grasses. In hexaploid wheat we found an expressed and putatively functional Lr34 homoeolog located on chromosome 4A, designated Lr34-B. Another homoeologous Lr34 copy, located on chromosome 7A, was disrupted by the insertion of repetitive elements. Protein sequences of LR34-B and LR34 were 97% identical. Orthologous Lr34 genes were detected in the genomes of Oryza sativa (rice) and Sorghum bicolor (sorghum). Zea mays (maize), Brachypodium distachyon and Hordeum vulgare (barley) lacked Lr34 orthologs, indicating independent deletion of this particular ABC transporter. Lr34 was part of a gene-rich island on the wheat D genome. We found gene colinearity on the homoeologous A and B genomes of hexaploid wheat, but little microcolinearity in other grasses. The homoeologous LR34-B protein and the orthologs from rice and sorghum have the susceptible haplotype for the two critical polymorphisms distinguishing the LR34 proteins from susceptible and resistant wheat cultivars. We conclude that the particular Lr34-haplotype found in resistant wheat cultivars is unique. It probably resulted from functional gene diversification that occurred after the polyploidization event that was at the origin of cultivated bread wheat.     

Changes in twelve homoeologous genomic regions in soybean following three rounds of polyploidy

With the advent of high-throughput sequencing, the availability of genomic sequence for comparative genomics is increasing exponentially. Numerous completed plant genome sequences enable characterization of patterns of the retention and evolution of genes within gene families due to multiple polyploidy events, gene loss and fractionation, and differential evolutionary pressures over time and across different gene families. In this report, we trace the changes that have occurred in 12 surviving homoeologous genomic regions from three rounds of polyploidy that contributed to the current Glycine max genome: a genome triplication before the origin of the rosids (~130 to 240 million years ago), a genome duplication early in the legumes (~58 million years ago), and a duplication in the Glycine lineage (~13 million years ago). Patterns of gene retention following the genome triplication event generally support predictions of the Gene Balance Hypothesis. Finally, we find that genes in networks with a high level of connectivity are more strongly conserved than those with low connectivity and that the enrichment of these highly connected genes in the 12 highly conserved homoeologous segments may in part explain their retention over more than 100 million years and repeated polyploidy events.     

Comparative genetic mapping between duplicated segments on maize chromosomes 3 and 8 and homoeologous regions in sorghum and sugarcane

Comparative mapping within maize, sorghum and sugarcane has previously revealed the existence of syntenic regions between the crops. In the present study, mapping on the sorghum genome of a set of probes previously located on the maize and sugarcane maps allow a detailed analysis of the relationship between maize chromosomes 3 and 8 and sorghum and sugarcane homoeologous regions. Of 49 loci revealed by 46 (4 sugarcane and 42 maize) polymorphic probes in sorghum, 42 were linked and were assigned to linkage groups G (28), E (10) and I (4). On the basis of common probes, a complete co-linearity is observed between sorghum linkage group G and the two sugarcane linkage groups II and III. The comparison between the consensus sorghum/sugarcane map (G/II/III) and the maps of maize chromosomes 3 and 8 reveals a series of linkage blocks within which gene orders are conserved. These blocks are interspersed with non-homoeologous regions corresponding to the central part of the two maize chromosomes and have been reshuffled, resulting in several inversions in maize compared to sorghum and sugarcane. The results emphasize the fact that duplication will considerably complicate precise comparative mapping at the whole genome scale between maize and other Poaceae.     

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class ( Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue ( Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.     

Ploidy in cyanobacteria

A recently developed real-time PCR method for the determination of genome copy numbers was optimized for the application to cyanobacteria. Three species were chosen to represent a fresh water species, a salt water species, and two strains of a widely used laboratory species. Synechococcus PCC 7942 and Synechococcus WH7803 were found to contain 3-4 genome copies per cell and are thus oligoploid, confirming earlier publications. In contrast, Synechocystis PCC 6803 is highly polyploid. The motile wild-type strain contains 218 genome copies in exponential phase and 58 genome copies in linear and in stationary growth phase. The GT wild-type strain contains 142 genome copies in exponential phase and 42 genome copies in linear and stationary growth phase. These are the highest numbers found for any cyanobacterial species. Notably these values are much higher than the value of 12 genome copies published for the 'Kazusa' strain more than 20 years ago. The results reveal that for Synechocystis PCC 6803 strain differences exist and that the ploidy level is highly growth phase-regulated. A compilation of the ploidy levels of all investigated cyanobacterial species gives an overview of the genome copy number distribution and shows that monoploid, oligoploid, and polyploid cyanobacteria exist.     

Genomic organization of the S core region and the S flanking regions of a class-II S haplotype in Brassica rapa

The nucleotide sequence of an 86.4-kb region that includes the SP11, SRK, and SLG genes of Brassica rapa S-60 (a class-II S haplotype) was determined. In the sequenced region, 13 putative genes were found besides SP11-60, SRK-60, and SLG-60. Five of these sequences were isolated as cDNAs, five were homologues of known genes, cDNAs, or ORFs, and three are hypothetical ORFs. Based on their nucleotide sequences, however, some of them are thought to be non-functional. Two regions of colinearity between the class-II S-60 and Brassica class-I S haplotypes were identified, i.e., S flanking region 1 which shows partial colinearity of non-genic sequences and S flanking region 2 which shows a high level of colinearity. The observed colinearity made it possible to compare the order of SP-11, SRK, and SLG genes in the S locus between the five sequenced S haplotypes. It emerged that the order of SRK and SLG in class-II S-60 is the reverse of that in the four class-I S haplotypes reported so far, and the order of SP11, SRK and SLG is the opposite of that in the class-I haplotype S-910. The possible gene designated as SAN1 (S locus Anther-expressed Non-coding RNA like-1), which is located in the region between SP11-60 and SRK-60, has features reminiscent of genes for non-coding RNAs (ncRNAs), but no homologous sequences were found in the databases. This sequence is transcribed in anthers but not in stigmas or leaves. These features of the genomic structure of S-60 are discussed with special reference to the characteristics of class-II S haplotypes.     

Large intraspecific haplotype variability at the Rph7 locus results from rapid and recent divergence in the barley genome

To study genome evolution and diversity in barley (Hordeum vulgare), we have sequenced and compared more than 300 kb of sequence spanning the Rph7 leaf rust disease resistance gene in two barley cultivars. Colinearity was restricted to five genic and two intergenic regions representing <35% of the two sequences. In each interval separating the seven conserved regions, the number and type of repetitive elements were completely different between the two homologous sequences, and a single gene was absent in one cultivar. In both cultivars, the nonconserved regions consisted of approximately 53% repetitive sequences mainly represented by long-terminal repeat retrotransposons that have inserted <1 million years ago. PCR-based analysis of intergenic regions at the Rph7 locus and at three other independent loci in 41 H. vulgare lines indicated large haplotype variability in the cultivated barley gene pool. Together, our data indicate rapid and recent divergence at homologous loci in the genome of H. vulgare, possibly providing the molecular mechanism for the generation of high diversity in the barley gene pool. Finally, comparative analysis of the gene composition in barley, wheat (Triticum aestivum), rice (Oryza sativa), and sorghum (Sorghum bicolor) suggested massive gene movements at the Rph7 locus in the Triticeae lineage.     

The colinearity of the Sh2/A1 orthologous region in rice, sorghum and maize is interrupted and accompanied by genome expansion in the triticeae

The Sh2/A1 orthologous region of maize, rice, and sorghum contains five genes in the order Sh2, X1, X2, and two A1 homologs in tandem duplication. The Sh2 and A1 homologs are separated by approximately 20 kb in rice and sorghum and by approximately 140 kb in maize. We analyzed the fate of the Sh2/A1 region in large-genome species of the Triticeae (wheat, barley, and rye). In the Triticeae, synteny in the Sh2/A1 region was interrupted by a break between the X1 and X2 genes. The A1 and X2 genes remained colinear in homeologous chromosomes as in other grasses. The Sh2 and X1 orthologs also remained colinear but were translocated to a nonhomeologous chromosome. Gene X1 was duplicated on two nonhomeologous chromosomes, and surprisingly, a paralog shared homology much higher than that of the orthologous copy to the X1 gene of other grasses. No tandem duplication of A1 homologs was detected but duplication of A1 on a nonhomeologous barley chromosome 6H was observed. Intergenic distances expanded greatly in wheat compared to rice. Wheat and barley diverged from each other 12 million years ago and both show similar changes in the Sh2/A1 region, suggesting that the break in colinearity as well as X1 duplications and genome expansion occurred in a common ancestor of the Triticeae species.     

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.     

Different patterns of evolution in the centromeric and telomeric regions of group A and B haplotypes of the human killer cell Ig-like receptor locus

The fast evolving human KIR gene family encodes variable lymphocyte receptors specific for polymorphic HLA class I determinants. Nucleotide sequences for 24 representative human KIR haplotypes were determined. With three previously defined haplotypes, this gave a set of 12 group A and 15 group B haplotypes for assessment of KIR variation. The seven gene-content haplotypes are all combinations of four centromeric and two telomeric motifs. 2DL5, 2DS5 and 2DS3 can be present in centromeric and telomeric locations. With one exception, haplotypes having identical gene content differed in their combinations of KIR alleles. Sequence diversity varied between haplotype groups and between centromeric and telomeric halves of the KIR locus. The most variable A haplotype genes are in the telomeric half, whereas the most variable genes characterizing B haplotypes are in the centromeric half. Of the highly polymorphic genes, only the 3DL3 framework gene exhibits a similar diversity when carried by A and B haplotypes. Phylogenetic analysis and divergence time estimates, point to the centromeric gene-content motifs that distinguish A and B haplotypes having emerged ~6 million years ago, contemporaneously with the separation of human and chimpanzee ancestors. In contrast, the telomeric motifs that distinguish A and B haplotypes emerged more recently, ~1.7 million years ago, before the emergence of Homo sapiens. Thus the centromeric and telomeric motifs that typify A and B haplotypes have likely been present throughout human evolution. The results suggest the common ancestor of A and B haplotypes combined a B-like centromeric region with an A-like telomeric region.