A locus on chromosome 7 determines myocardial cell necrosis and calcification (dystrophic cardiac calcinosis) in mice.

Dystrophic cardiac calcinosis, an age-related cardiomyopathy that occurs among certain inbred strains of mice, involves myocardial injury, necrosis, and calcification. Using a complete linkage map approach and quantitative trait locus analysis, we sought to identify genetic loci determining dystrophic cardiac calcinosis in an F2 intercross of resistant C57BL/6J and susceptible C3H/HeJ inbred strains. We identified a single major locus, designated Dyscalc, located on proximal chromosome 7 in a region syntenic with human chromosomes 19q13 and 11p15. The statistical significance of Dyscalc (logarithm of odds score 14.6) was tested by analysis of permuted trait data. Analysis of BxH recombinant inbred strains confirmed the mapping position. The inheritance pattern indicated that this locus influences susceptibility of cells both to enter necrosis and to subsequently undergo calcification.

A complete genome screen for genes predisposing to severe bipolar disorder in two Costa Rican pedigrees

Bipolar mood disorder (BP) is a debilitating syndrome characterized by episodes of mania and depression. We designed a multistage study to detect all major loci predisposing to severe BP (termed BP-I) in two pedigrees drawn from the Central Valley of Costa Rica, where the population is largely descended from a few founders in the 16th–18th centuries. We considered only individuals with BP-I as affected and screened the genome for linkage with 473 microsatellite markers. We used a model for linkage analysis that incorporated a high phenocopy rate and a conservative estimate of penetrance. Our goal in this study was not to establish definitive linkage but rather to detect all regions possibly harboring major genes for BP-I in these pedigrees. To facilitate this aim, we evaluated the degree to which markers that were informative in our data set provided coverage of each genome region; we estimate that at least 94% of the genome has been covered, at a predesignated threshold determined through prior linkage simulation analyses. We report here the results of our genome screen for BP-I loci and indicate several regions that merit further study, including segments in 18q, 18p, and 11p, in which suggestive lod scores were observed for two or more contiguous markers. Isolated lod scores that exceeded our thresholds in one or both families also occurred on chromosomes 1, 2, 3, 4, 5, 7, 13, 15, 16, and 17. Interesting regions highlighted in this genome screen will be followed up using linkage disequilibrium (LD) methods.

Linkage disequilibrium mapping in isolated populations: The example of Finland revisited

Linkage disequilibrium analysis can provide high resolution in the mapping of disease genes because it incorporates information on recombinations that have occurred during the entire period from the mutational event to the present. A circumstance particularly favorable for high-resolution mapping is when a single founding mutation segregates in an isolated population. We review here the population structure of Finland in which a small founder population some 100 generations ago has expanded into 5.1 million people today. Among the 30-odd autosomal recessive disorders that are more prevalent in Finland than elsewhere, several appear to have segregated for this entire period in the “panmictic” southern Finnish population. Linkage disequilibrium analysis has allowed precise mapping and determination of genetic distances at the 0.1-cM level in several of these disorders. Estimates of genetic distance have proven accurate, but previous calculations of the confidence intervals were too small because sampling variation was ignored. In the north and east of Finland the population can be viewed as having been “founded” only after 1500. Disease mutations that have undergone such a founding bottleneck only 20 or so generations ago exhibit linkage disequilibrium and haplotype sharing over long genetic distances (5–15 cM). These features have been successfully exploited in the mapping and cloning of many genes. We review the statistical issues of fine mapping by linkage disequilibrium and suggest that improved methodologies may be necessary to map diseases of complex etiology that may have arisen from multiple founding mutations.

Complete resolution of the solid-state NMR spectrum of a uniformly 15N-labeled membrane protein in phospholipid bilayers

Complete resolution of the amide resonances in a three-dimensional solid-state NMR correlation spectrum of a uniformly 15N-labeled membrane protein in oriented phospholipid bilayers is demonstrated. The three orientationally dependent frequencies, 1H chemical shift, 1H–15N dipolar coupling, and 15N chemical shift, associated with each amide resonance are responsible for resolution among resonances and provide sufficient angular restrictions for protein structure determination. Because the protein is completely immobilized by the phospholipids on the relevant NMR time scales (10 kHz), the linewidths will not degrade in the spectra of larger proteins. Therefore, these results demonstrate that solid-state NMR experiments can overcome the correlation time problem and extend the range of proteins that can have their structures determined by NMR spectroscopy to include uniformly 15N-labeled membrane proteins in phospholipid bilayers.

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Plasmodium falciparum antigenic diversity: Evidence of clonal population structure

Plasmodium falciparum, the agent of malignant malaria, is one of mankind’s most severe scourges. Efforts to develop preventive vaccines or remedial drugs are handicapped by the parasite’s rapid evolution of drug resistance and protective antigens. We examine 25 DNA sequences of the gene coding for the highly polymorphic antigenic circumsporozoite protein. We observe total absence of silent nucleotide variation in the two nonrepeated regions of the gene. We propose that this absence reflects a recent origin (within several thousand years) of the world populations of P. falciparum from a single individual; the amino acid polymorphisms observed in these nonrepeat regions would result from strong natural selection. Analysis of these polymorphisms indicates that: (i) the incidence of recombination events does not increase with nucleotide distance; (ii) the strength of linkage disequilibrium between nucleotides is also independent of distance; and (iii) haplotypes in the two nonrepeat regions are correlated with one another, but not with the central repeat region they span. We propose two hypotheses: (i) variation in the highly polymorphic central repeat region arises by mitotic intragenic recombination, and (ii) the population structure of P. falciparum is clonal—a state of affairs that persists in spite of the necessary stage of physiological sexuality that the parasite must sustain in the mosquito vector to complete its life cycle.

Complete mitochondrial genome suggests diapsid affinities of turtles

Despite more than a century of debate, the evolutionary position of turtles (Testudines) relative to other amniotes (reptiles, birds, and mammals) remains uncertain. One of the major impediments to resolving this important evolutionary problem is the highly distinctive and enigmatic morphology of turtles that led to their traditional placement apart from diapsid reptiles as sole descendants of presumably primitive anapsid reptiles. To address this question, the complete (16,787-bp) mitochondrial genome sequence of the African side-necked turtle (Pelomedusa subrufa) was determined. This molecule contains several unusual features: a (TA)n microsatellite in the control region, the absence of an origin of replication for the light strand in the WANCY region of five tRNA genes, an unusually long noncoding region separating the ND5 and ND6 genes, an overlap between ATPase 6 and COIII genes, and the existence of extra nucleotides in ND3 and ND4L putative ORFs. Phylogenetic analyses of the complete mitochondrial genome sequences supported the placement of turtles as the sister group of an alligator and chicken (Archosauria) clade. This result clearly rejects the Haematothermia hypothesis (a sister-group relationship between mammals and birds), as well as rejecting the placement of turtles as the most basal living amniotes. Moreover, evidence from both complete mitochondrial rRNA genes supports a sister-group relationship of turtles to Archosauria to the exclusion of Lepidosauria (tuatara, snakes, and lizards). These results challenge the classic view of turtles as the only survivors of primary anapsid reptiles and imply that turtles might have secondarily lost their skull fenestration.

Genome scan of human systemic lupus erythematosus: Evidence for linkage on chromosome 1q in African-American pedigrees

Systemic lupus erythematosus (SLE) is an autoimmune disorder characterized by production of autoantibodies against intracellular antigens including DNA, ribosomal P, Ro (SS-A), La (SS-B), and the spliceosome. Etiology is suspected to involve genetic and environmental factors. Evidence of genetic involvement includes: associations with HLA-DR3, HLA-DR2, Fcγ receptors (FcγR) IIA and IIIA, and hereditary complement component deficiencies, as well as familial aggregation, monozygotic twin concordance >20%, λs > 10, purported linkage at 1q41–42, and inbred mouse strains that consistently develop lupus. We have completed a genome scan in 94 extended multiplex pedigrees by using model-based linkage analysis. Potential [log10 of the odds for linkage (lod) > 2.0] SLE loci have been identified at chromosomes 1q41, 1q23, and 11q14–23 in African-Americans; 14q11, 4p15, 11q25, 2q32, 19q13, 6q26–27, and 12p12–11 in European-Americans; and 1q23, 13q32, 20q13, and 1q31 in all pedigrees combined. An effect for the FcγRIIA candidate polymorphism) at 1q23 (lod = 3.37 in African-Americans) is syntenic with linkage in a murine model of lupus. Sib-pair and multipoint nonparametric analyses also support linkage (P < 0.05) at nine loci detected by using two-point lod score analysis (lod > 2.0). Our results are consistent with the presumed complexity of genetic susceptibility to SLE and illustrate racial origin is likely to influence the specific nature of these genetic effects.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

Comparison of parametric and nonparametric methods to map oligogenes by linkage

A sample of 95 sib pairs affected with insulin-dependent diabetes and typed with their normal parents for 28 markers on chromosome 6 has been analyzed by several methods. When appropriate parameters are efficiently estimated, a parametric model is equivalent to the β model, which is superior to nonparametric alternatives both in single point tests (as found previously) and in multipoint tests. Theory is given for meta-analysis combined with allelic association, and problems that may be associated with errors of map location and/or marker typing are identified. Reducing by multipoint analysis the number of association tests in a dense map can give a 3-fold reduction in the critical lod, and therefore in the cost of positional cloning.

Glucose-induced Autophagy of Peroxisomes in Pichia pastoris Requires a Unique E1-like Protein

Cytosolic and peroxisomal enzymes necessary for methanol assimilation are synthesized when Pichia pastoris is grown in methanol. Upon adaptation from methanol to a glucose environment, these enzymes are rapidly and selectively sequestered and degraded within the yeast vacuole. Sequestration begins when the vacuole changes shape and surrounds the peroxisomes. The opposing membranes then fuse, engulfing the peroxisome. In this study, we have characterized a mutant cell line (glucose-induced selective autophagy), gsa7, which is defective in glucose-induced selective autophagy of peroxisomes, and have identified the GSA7 gene. Upon glucose adaptation, gsa7 cells were unable to degrade peroxisomal alcohol oxidase. We observed that the peroxisomes were surrounded by the vacuole, but complete uptake into the vacuole did not occur. Therefore, we propose that GSA7 is not required for initiation of autophagy but is required for bringing the opposing vacuolar membranes together for homotypic fusion, thereby completing peroxisome sequestration. By sequencing the genomic DNA fragment that complemented the gsa7 phenotype, we have found that GSA7 encodes a protein of 71 kDa (Gsa7p) with limited sequence homology to a family of ubiquitin-activating enzymes, E1. The knockout mutant gsa7Δ had an identical phenotype to gsa7, and both mutants were rescued by an epitope-tagged Gsa7p (Gsa7-hemagglutinin [HA]). In addition, a GSA7 homolog, APG7, a protein required for autophagy in Saccharomyces cerevisiae, was capable of rescuing gsa7. We have sequenced the human homolog of GSA7 and have shown many regions of identity between the yeast and human proteins. Two of these regions align to the putative ATP-binding domain and catalytic site of the family of ubiquitin activating enzymes, E1 (UBA1, UBA2, and UBA3). When either of these sites was mutated, the resulting mutants [Gsa7(ΔATP)-HA and Gsa7(C518S)-HA] were unable to rescue gsa7 cells. We provide evidence to suggest that Gsa7-HA formed a thio-ester linkage with a 25–30 kDa protein. This conjugate was not observed in cells expressing Gsa7(ΔATP)-HA or in cells expressing Gsa7(C518S)-HA. Our results suggest that this unique E1-like enzyme is required for homotypic membrane fusion, a late event in the sequestration of peroxisomes by the vacuole.

Whole genome amplification using a degenerate oligonucleotide primer allows hundreds of genotypes to be performed on less than one nanogram of genomic DNA

Genetic analysis of limiting quantities of genomic DNA play an important role in DNA forensics, paleoarcheology, genetic disease diagnosis, genetic linkage analysis, and genetic diversity studies. We have tested the ability of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to amplify picogram quantities of human genomic DNA for the purpose of increasing the amount of template for genotyping with microsatellite repeat markers. DNA was uniformly amplified at a large number of typable loci throughout the human genome with starting template DNAs from as little as 15 pg to as much as 400 ng. A much greater-fold enrichment was seen for the smaller genomic DOP-PCRs. All markers tested were amplified from starting genomic DNAs in the range of 0.6–40 ng with amplifications of 200- to 600-fold. The DOP-PCR-amplified genomic DNA was an excellent and reliable template for genotyping with microsatellites, which give distinct bands with no increase in stutter artifact on di-, tri-, and tetranucleotide repeats. There appears to be equal amplification of genomic DNA from 55 of 55 tested discrete microsatellites implying near complete coverage of the human genome. Thus, DOP-PCR appears to allow unbiased, hundreds-fold whole genome amplification of human genomic DNA for genotypic analysis.

Native display of complete foreign protein domains on the surface of hepatitis B virus capsids

The nucleocapsid of hepatitis B virus (HBV), or HBcAg, is a highly symmetric structure formed by multiple dimers of a single core protein that contains potent T helper epitopes in its 183-aa sequence. Both factors make HBcAg an unusually strong immunogen and an attractive candidate as a carrier for foreign epitopes. The immunodominant c/e1 epitope on the capsid has been suggested as a superior location to convey high immunogenicity to a heterologous sequence. Because of its central position, however, any c/e1 insert disrupts the core protein’s primary sequence; hence, only peptides, or rather small protein fragments seemed to be compatible with particle formation. According to recent structural data, the epitope is located at the tips of prominent surface spikes formed by the very stable dimer interfaces. We therefore reasoned that much larger inserts might be tolerated, provided the individual parts of a corresponding fusion protein could fold independently. Using the green fluorescent protein (GFP) as a model insert, we show that the chimeric protein efficiently forms fluorescent particles; hence, all of its structurally important parts must be properly folded. We also demonstrate that the GFP domains are surface-exposed and that the chimeric particles elicit a potent humoral response against native GFP. Hence, proteins of at least up to 238 aa can be natively displayed on the surface of HBV core particles. Such chimeras may not only be useful as vaccines but may also open the way for high resolution structural analyses of nonassembling proteins by electron microscopy.

An Escherichia coli strain with all chromosomal rRNA operons inactivated: Complete exchange of rRNA genes between bacteria

Current global phylogenies are built predominantly on rRNA sequences. However, an experimental system for studying the evolution of rRNA is not readily available, mainly because the rRNA genes are highly repeated in most experimental organisms. We have constructed an Escherichia coli strain in which all seven chromosomal rRNA operons are inactivated by deletions spanning the 16S and 23S coding regions. A single E. coli rRNA operon carried by a multicopy plasmid supplies 16S and 23S rRNA to the cell. By using this strain we have succeeded in creating microorganisms that contain only a foreign rRNA operon derived from either Salmonella typhimurium or Proteus vulgaris, microorganisms that have diverged from E. coli about 120–350 million years ago. We also were able to replace the E. coli rRNA operon with an E. coli/yeast hybrid one in which the GTPase center of E. coli 23S rRNA had been substituted by the corresponding domain from Saccharomyces cerevisiae. These results suggest that, contrary to common belief, coevolution of rRNA with many other components in the translational machinery may not completely preclude the horizontal transfer of rRNA genes.

Type 2 diabetes: Evidence for linkage on chromosome 20 in 716 Finnish affected sib pairs

We are conducting a genome scan at an average resolution of 10 centimorgans (cM) for type 2 diabetes susceptibility genes in 716 affected sib pairs from 477 Finnish families. To date, our best evidence for linkage is on chromosome 20 with potentially separable peaks located on both the long and short arms. The unweighted multipoint maximum logarithm of odds score (MLS) was 3.08 on 20p (location, x̂ = 19.5 cM) under an additive model, whereas the weighted MLS was 2.06 on 20q (x̂ = 57 cM, recurrence risk, λ̂s = 1.25, P = 0.009). Weighted logarithm of odds scores of 2.00 (x̂ = 69.5 cM, P = 0.010) and 1.92 (x̂ = 18.5 cM, P = 0.013) were also observed. Ordered subset analyses based on sibships with extreme mean values of diabetes-related quantitative traits yielded sets of families who contributed disproportionately to the peaks. Two-hour glucose levels in offspring of diabetic individuals gave a MLS of 2.12 (P = 0.0018) at 9.5 cM. Evidence from this and other studies suggests at least two diabetes-susceptibility genes on chromosome 20. We have also screened the gene for maturity-onset diabetes of the young 1, hepatic nuclear factor 4-a (HNF-4α) in 64 affected sibships with evidence for high chromosomal sharing at its location on chromosome 20q. We found no evidence that sequence changes in this gene accounted for the linkage results we observed.