Haplotype variation and linkage disequilibrium in 313 human genes

Variation within genes has important implications for all biological traits. We identified 3899 single nucleotide polymorphisms (SNPs) that were present within 313 genes from 82 unrelated individuals of diverse ancestry, and we organized the SNPs into 4304 different haplotypes. Each gene had several variable SNPs and haplotypes that were present in all populations, as well as a number that were population-specific. Pairs of SNPs exhibited variability in the degree of linkage disequilibrium that was a function of their location within a gene, distance from each other, population distribution, and population frequency. Haplotypes generally had more information content (heterozygosity) than did individual SNPs. Our analysis of the pattern of variation strongly supports the recent expansion of the human population.

SPROUT: recent developments in the de novo design of molecules

SPROUT is a computer program for constrained structure generation. It is designed to generate molecules for a range of applications in molecular recognition. The program uses a number of approximations that enable a wide variety of diverse structures to be generated. Practical use of the program is demonstrated in two examples. The first demonstrates the ability of the program to generate candidate inhibitors for a receptor site of known 3D structure, specifically the GDP binding site of p21. In the second example, structures are generated to fit a pharmacophore hypothesis that models morphine agonists.

Protein secondary structure of the isolated photosystem II reaction center and conformational changes studied by Fourier transform infrared spectroscopy

The secondary structure of the photosystem II (PSII) reaction center isolated from pea chloroplasts has been characterized by Fourier transform infrared (FTIR) spectroscopy. Spectra were recorded in aqueous buffers containing H2O or D2O; the detergent present for most measurements was dodecyl maltoside. The broad amide I and amide II bands were analyzed by using second-derivative and deconvolution procedures. Absorption bands were assigned to the presence of alpha-helices, beta-sheets, turns, or random structure. Quantitative analysis revealed that this complex contained a high proportion of alpha-helices (67%) and some antiparallel beta-sheets (9%) and turns (11%). An irreversible decrease in the intensity of the band associated with the alpha-helices occurs upon exposure of the isolated PSII reaction center to bright illumination. This loss of alpha-helical content gave rise to an increase in other secondary structures, particularly beta-sheets. After similar pretreatment with light, sodium dodecyl sulfate polyacrylamide gel electrophoresis reveals lower mobility and solubility of constituent D1 and D2 polypeptides of the PSII reaction center. Some degradation of these polypeptides also occurs. In contrast, there is no change in the mobility of the two subunits of cytochrome b559. In the absence of illumination, the PSII reaction center exchanged into dodecyl maltoside shows good thermal stability as compared with samples in Triton X-100. Only at a temperature of about 60 degrees C do spectral changes take place that are indicative of denaturation.

Evolutionary dynamics of non-coding sequences within the class II region of the human MHC

About 40% (350 kb) of the human MHC class II region has been sequenced and a coordinated effort to sequence the entire MHC is underway. In addition to the coding information (22 genes/pseudogenes), the non-coding sequences reveal novel information on the organisation and evolution of the MHC as demonstrated here by the example of a 200 kb contig that has been analysed for local and global features. In conjunction with cross-species comparisons, our results present new evidence on the structure of isochores, the evolutionary dynamics of repeat-mediated recombination and its effect on certain MHC encoded genes, and a higher than average degree of natural polymorphism that has implications for sequencing the human genome. We also report the finding of a class I-related pseudogene (HLA-ZI) in the middle of the class II region, which provides the first direct evidence for DNA exchange between these two related regions in man.

Cross sections and oscillator strengths for electron-impact excitation of the ÃB11 electronic state of water

The authors report absolute differential and integral cross section measurements for electron-impact excitation of the A (1)B(1) electronic state of water. This is an important channel for the production of the OH (X (2)Pi) radical, as well as for understanding the origin of the atmospheric Meinel [Astrophys. J. 111, 555 (1950)] bands. The incident energy range of our measurements is 20-200 eV, while the angular range of the differential cross section data is 3.5 degrees -90 degrees . This is the first time such data are reported in the literature and, where possible, comparison to existing theoretical work, and new scaled Born cross sections calculated as a part of the current study, is made. The scaled Born cross sections are in good agreement with the integral cross sections deduced from the experimental differential cross sections. In addition they report (experimental) generalized oscillator strength data at the incident energies of 100 and 200 eV. These data are used to derive a value for the optical oscillator strength which is found to be in excellent agreement with that from an earlier dipole (e,e) experiment and an earlier photoabsorption experiment.

Norbornane: An investigation into its valence electronic structure using electron momentum spectroscopy, and density functional and Green’s function theories

We report on the results of an exhaustive study of the valence electronic structure of norbornane (C(7)H(12)), up to binding energies of 29 eV. Experimental electron momentum spectroscopy and theoretical Green's function and density functional theory approaches were all utilized in this investigation. A stringent comparison between the electron momentum spectroscopy and theoretical orbital momentum distributions found that, among all the tested models, the combination of the Becke-Perdew functional and a polarized valence basis set of triple-zeta quality provides the best representation of the electron momentum distributions for all of the 20 valence orbitals of norbornane. This experimentally validated quantum chemistry model was then used to extract some chemically important properties of norbornane. When these calculated properties are compared to corresponding results from other independent measurements, generally good agreement is found. Green's function calculations with the aid of the third-order algebraic diagrammatic construction scheme indicate that the orbital picture of ionization breaks down at binding energies larger than 22.5 eV. Despite this complication, they enable insights within 0.2 eV accuracy into the available ultraviolet photoemission and newly presented (e,2e) ionization spectra, except for the band associated with the 1a(2) (-1) one-hole state, which is probably subject to rather significant vibronic coupling effects, and a band at approximately 25 eV characterized by a momentum distribution of "s-type" symmetry, which Green's function calculations fail to reproduce. We note the vicinity of the vertical double ionization threshold at approximately 26 eV.

Identification and analysis of multigene families by comparison of exon fingerprints

Gene families are often recognised by sequence homology using similarity searching to find relationships, however, genomic sequence data provides gene architectural information not used by conventional search methods. In particular, intron positions and phases are expected to be relatively conserved features, because mis-splicing and reading frame shifts should be selected against. A fast search technique capable of detecting possible weak sequence homologies apparent at the intron/exon level of gene organization is presented for comparing spliceosomal genes and gene fragments. FINEX compares strings of exons delimited by intron/exon boundary positions and intron phases (exon fingerprint) using a global dynamic programming algorithm with a combined intron phase identity and exon size dissimilarity score. Exon fingerprints are typically two orders of magnitude smaller than their nucleic acid sequence counterparts giving rise to fast search times: a ranked search against a library of 6755 fingerprints for a typical three exon fingerprint completes in under 30 seconds on an ordinary workstation, while a worst case largest fingerprint of 52 exons completes in just over one minute. The short "sequence" length of exon fingerprints in comparisons is compensated for by the large exon alphabet compounded of intron phase types and a wide range of exon sizes, the latter contributing the most information to alignments. FINEX performs better in some searches than conventional methods, finding matches with similar exon organization, but low sequence homology. A search using a human serum albumin finds all members of the multigene family in the FINEX database at the top of the search ranking, despite very low amino acid percentage identities between family members. The method should complement conventional sequence searching and alignment techniques, offering a means of identifying otherwise hard to detect homologies where genomic data are available.

Efficient high-resolution genetic mapping of mouse interspersed repetitive sequence PCR products, toward integrated genetic and physical mapping of the mouse genome

The ability to carry out high-resolution genetic mapping at high throughput in the mouse is a critical rate-limiting step in the generation of genetically anchored contigs in physical mapping projects and the mapping of genetic loci for complex traits. To address this need, we have developed an efficient, high-resolution, large-scale genome mapping system. This system is based on the identification of polymorphic DNA sites between mouse strains by using interspersed repetitive sequence (IRS) PCR. Individual cloned IRS PCR products are hybridized to a DNA array of IRS PCR products derived from the DNA of individual mice segregating DNA sequences from the two parent strains. Since gel electrophoresis is not required, large numbers of samples can be genotyped in parallel. By using this approach, we have mapped > 450 polymorphic probes with filters containing the DNA of up to 517 backcross mice, potentially allowing resolution of 0.14 centimorgan. This approach also carries the potential for a high degree of efficiency in the integration of physical and genetic maps, since pooled DNAs representing libraries of yeast artificial chromosomes or other physical representations of the mouse genome can be addressed by hybridization of filter representations of the IRS PCR products of such libraries.

Construction of genetic maps using distance geometry

The techniques of distance geometry, which generate coordinates from observed interpoint distances, have been applied to the problem of determining the relative positions of linked genetic loci from observed interlocus distances. Only the most precise data needed to join the loci are used, with missing distances substituted by sums of precise intermediate distances. Good initial positions (and therefore the order) of loci on a linear map are obtained in an operation of complexity O(N3). The method can therefore be used to generate good initial framework maps for the large numbers of markers encountered in current mapping projects. The locus positions can be subsequently refined to maximize the agreement with the originally observed distances, taking account of the weights of individual interlocus distances. By choosing only small distances from which to construct the map, the method reduces any error due to an incorrect choice of mapping function. It also prevents undue expansion of the map due to error-prone markers, since such markers are accommodated in higher dimensions. The method estimates the error in the positions of individual markers on the final map and identifies well- and ill-defined regions of the map.

MHCDB: database of the human MHC (release 2)

The second release of the human major histocompatibility complex (MHC) database is now publicly available. It contains an updated physical map and considerably more genomic sequence. cDNA sequences of all current alleles are accessible as individual sequence entries. The variability of different genes is displayed graphically as static and dynamic images accessible from the database. Known disease-serotype associations have also been incorporated, together with data from the MHCPEP database of eluted peptides.

Estimation of distances and map construction using radiation hybrids

A method of estimating distances between pairs of genetic markers is described that directly uses their observed joint frequency distribution in a panel of radiation hybrids (RHs). The distance measure is based on the strength of association between marker pairs, which is high for close markers and decays with distance. These distances are then submitted to a previous method that generates linear coordinates for the markers directly from the intermarker distance matrix. This method of map building from RH data is simpler than others, because it uses only the observed joint frequency distributions of markers in the panel, and does not attempt to model unobserved quantities such as the retention of different sized fragments that contain the markers. It also incorporates directly the observed variation in retention of different markers, without needing a model for differential fragment retention dependent on chromosomal location, which is generally not known. Only small, precise distances are used in map construction, thereby reducing any effects of different fragment retention frequencies and local variations in X-ray sensitivity. The method is tested by simulation, and known marker distances and locations are successfully recovered from RH raw data. The method is also applied to publicly available data sets related to the recent transcript map of the human genome.

MHCDB--database of the human MHC

Genetic and physical data relating to the human major histocompatibility complex (MHC) were compiled and analyzed, using a genome analysis program. The current contents of the database include: 1) location of over 100 genes and other markers; 2) location of over 250 YAC and cosmid clones; 3) 150 kilobases of genomic DNA sequence including full annotation (exon/intron boundaries, repeats, promoters, etc.); 4) cDNA sequences of currently-known class I and class II alleles; and 5) accompanying descriptive data--references, comments, laboratory addresses, and so on. The database allows rapid access, retrieval, and display of all these data, which should make it a useful tool for the study of the human MHC. MHCDB is publicly available and will be updated as new data become available.