Genetic analysis of completely sequenced disease-associated MHC haplotypes identifies shuffling of segments in recent human history.

The major histocompatibility complex (MHC) is recognised as one of the most important genetic regions in relation to common human disease. Advancement in identification of MHC genes that confer susceptibility to disease requires greater knowledge of sequence variation across the complex. Highly duplicated and polymorphic regions of the human genome such as the MHC are, however, somewhat refractory to some whole-genome analysis methods. To address this issue, we are employing a bacterial artificial chromosome (BAC) cloning strategy to sequence entire MHC haplotypes from consanguineous cell lines as part of the MHC Haplotype Project. Here we present 4.25 Mb of the human haplotype QBL (HLA-A26-B18-Cw5-DR3-DQ2) and compare it with the MHC reference haplotype and with a second haplotype, COX (HLA-A1-B8-Cw7-DR3-DQ2), that shares the same HLA-DRB1, -DQA1, and -DQB1 alleles. We have defined the complete gene, splice variant, and sequence variation contents of all three haplotypes, comprising over 259 annotated loci and over 20,000 single nucleotide polymorphisms (SNPs). Certain coding sequences vary significantly between different haplotypes, making them candidates for functional and disease-association studies. Analysis of the two DR3 haplotypes allowed delineation of the shared sequence between two HLA class II-related haplotypes differing in disease associations and the identification of at least one of the sites that mediated the original recombination event. The levels of variation across the MHC were similar to those seen for other HLA-disparate haplotypes, except for a 158-kb segment that contained the HLA-DRB1, -DQA1, and -DQB1 genes and showed very limited polymorphism compatible with identity-by-descent and relatively recent common ancestry (<3,400 generations). These results indicate that the differential disease associations of these two DR3 haplotypes are due to sequence variation outside this central 158-kb segment, and that shuffling of ancestral blocks via recombination is a potential mechanism whereby certain DR-DQ allelic combinations, which presumably have favoured immunological functions, can spread across haplotypes and populations.

The DNA sequence of the human X chromosome.

The human X chromosome has a unique biology that was shaped by its evolution as the sex chromosome shared by males and females. We have determined 99.3% of the euchromatic sequence of the X chromosome. Our analysis illustrates the autosomal origin of the mammalian sex chromosomes, the stepwise process that led to the progressive loss of recombination between X and Y, and the extent of subsequent degradation of the Y chromosome. LINE1 repeat elements cover one-third of the X chromosome, with a distribution that is consistent with their proposed role as way stations in the process of X-chromosome inactivation. We found 1,098 genes in the sequence, of which 99 encode proteins expressed in testis and in various tumour types. A disproportionately high number of mendelian diseases are documented for the X chromosome. Of this number, 168 have been explained by mutations in 113 X-linked genes, which in many cases were characterized with the aid of the DNA sequence.

Genomic sequence of the class II region of the canine MHC: comparison with the MHC of other mammalian species.

The domestic dog, Canis familiaris, is an excellent model species in which to study complex inherited diseases, having over 200 recognized breeds, each of which represents a closed gene pool. Overlapping canine genomic BAC clones were sequenced to obtain 711,521 bp of the canine classical and extended MHC class II regions. Analysis and annotation of this sequence reveals that it contains 45 loci, of which 29 are predicted to be functionally expressed. Comparison of the DLA class II sequence with those of the cat, human, and mouse highlights regions of syntenic conservation and species-specific gene rearrangement and duplication and gives an insight into the evolution of the DR region in the order Carnivora. Elucidation of functionally important dog class II genes and the identification of 23 microsatellite markers spanning this region will contribute significantly to the study of canine diseases that have an immune component.

Identification of mammalian microRNA host genes and transcription units.

To derive a global perspective on the transcription of microRNAs (miRNAs) in mammals, we annotated the genomic position and context of this class of noncoding RNAs (ncRNAs) in the human and mouse genomes. Of the 232 known mammalian miRNAs, we found that 161 overlap with 123 defined transcription units (TUs). We identified miRNAs within introns of 90 protein-coding genes with a broad spectrum of molecular functions, and in both introns and exons of 66 mRNA-like noncoding RNAs (mlncRNAs). In addition, novel families of miRNAs based on host gene identity were identified. The transcription patterns of all miRNA host genes were curated from a variety of sources illustrating spatial, temporal, and physiological regulation of miRNA expression. These findings strongly suggest that miRNAs are transcribed in parallel with their host transcripts, and that the two different transcription classes of miRNAs ('exonic' and 'intronic') identified here may require slightly different mechanisms of biogenesis.

Complete MHC haplotype sequencing for common disease gene mapping.

The future systematic mapping of variants that confer susceptibility to common diseases requires the construction of a fully informative polymorphism map. Ideally, every base pair of the genome would be sequenced in many individuals. Here, we report 4.75 Mb of contiguous sequence for each of two common haplotypes of the major histocompatibility complex (MHC), to which susceptibility to >100 diseases has been mapped. The autoimmune disease-associated-haplotypes HLA-A3-B7-Cw7-DR15 and HLA-A1-B8-Cw7-DR3 were sequenced in their entirety through a bacterial artificial chromosome (BAC) cloning strategy using the consanguineous cell lines PGF and COX, respectively. The two sequences were annotated to encompass all described splice variants of expressed genes. We defined the complete variation content of the two haplotypes, revealing >18,000 variations between them. Average SNP densities ranged from less than one SNP per kilobase to >60. Acquisition of complete and accurate sequence data over polymorphic regions such as the MHC from large-insert cloned DNA provides a definitive resource for the construction of informative genetic maps, and avoids the limitation of chromosome regions that are refractory to PCR amplification.

Organisation of the pantothenate (vitamin B5) biosynthesis pathway in higher plants.

Pantothenate (vitamin B5) is the precursor for the biosynthesis of the phosphopantetheine moiety of coenzyme A and acyl carrier protein, and is synthesised in Escherichia coli by four enzymic reactions. Ketopantoate hydroxymethyltransferase (KPHMT) and pantothenate synthetase (PtS) catalyse the first and last steps, respectively. Two genes encoding KPHMT and one for PtS were identified in the Arabidopsis thaliana genome, and cDNAs for all three genes were amplified by PCR. The cDNAs were able to complement their respective E. coli auxotrophs, demonstrating that they encoded functional enzymes. Subcellular localisation of the proteins was investigated using green fluorescent protein (GFP) fusions and confocal microscopy. The two KPHMT-GFP fusion proteins were targeted exclusively to mitochondria, whereas PtS-GFP was found in the cytosol. This implies that there must be transporters for pathway intermediates. KPHMT enzyme activity could be measured in purified mitochondria from both pea leaves and Arabidopsis suspension cultures. We investigated whether Arabidopsis encoded homologues of the remaining two pantothenate biosynthesis enzymes from E. coli, l-aspartate-alpha-decarboxylase (ADC) and ketopantoate reductase (KPR). No homologue of ADC could be identified using either conventional blast or searches with the program fugue in which the structure of the E. coli ADC was compared to all the annotated proteins in Arabidopsis. ADC also appears to be absent from the genome of the yeast, Saccharomyces cerevisiae, by the same criteria. In contrast, a putative Arabidopsis oxidoreductase with some similarity to KPR was identified with fugue.

Gene annotation: prediction and testing.

Fifty years after the publication of DNA structure, the whole human genome sequence will be officially finished. This achievement marks the beginning of the task to catalogue every human gene and identify each of their function expression patterns. Currently, researchers estimate that there are about 30,000 human genes and approximately 70% of these can be automatically predicted using a combination of ab initio and similarity-based programs. However, to experimentally investigate every gene's function, the research community requires a high-quality annotation of alternative splicing, pseudogenes, and promoter regions that can only be provided by manual intervention. Manual curation of the human genome will be a long-term project as experimental data are continually produced to confirm or refine the predictions, and new features such as noncoding RNAs and enhancers have not been fully identified. Such a highly curated human gene-set made publicly available will be a great asset for the experimental community and for future comparative genome projects.