Mutations in TNFRSF13B encoding TACI are associated with common variable immunodeficiency in humans.

The functional interaction of BAFF and APRIL with TNF receptor superfamily members BAFFR, TACI and BCMA is crucial for development and maintenance of humoral immunity in mice and humans. Using a candidate gene approach, we identified homozygous and heterozygous mutations in TNFRSF13B, encoding TACI, in 13 individuals with common variable immunodeficiency. Homozygosity with respect to mutations causing the amino acid substitutions S144X and C104R abrogated APRIL binding and resulted in loss of TACI function, as evidenced by impaired proliferative response to IgM-APRIL costimulation and defective class switch recombination induced by IL-10 and APRIL or BAFF. Family members heterozygous with respect to the C104R mutation and individuals with sporadic common variable immunodeficiency who were heterozygous with respect to the amino acid substitutions A181E, S194X and R202H had humoral immunodeficiency. Although signs of autoimmunity and lymphoproliferation are evident, the human phenotype differs from that of the Tnfrsf13b-/- mouse model.

Co-mapping studies of QTLs for fruit acidity and candidate genes of organic acid metabolism and proton transport in sweet melon (Cucumis melo L.).

Sweet melon cultivars contain a low level of organic acids and, therefore, the quality and flavor of sweet melon fruit is determined almost exclusively by fruit sugar content. However, genetic variability for fruit acid levels in the Cucumis melo species exists and sour fruit accessions are characterized by acidic fruit pH of <5, compared to the sweet cultivars that are generally characterized by mature fruit pH values of >6. In this paper, we report results from a mapping population based on recombinant inbred lines (RILs) derived from the cross between the non-sour 'Dulce' variety and the sour PI 414323 accession. Results show that a single major QTL for pH co-localizes with major QTLs for the two predominant organic acids in melon fruit, citric and malic, together with an additional metabolite which we identified as uridine. While the acidic recombinants were characterized by higher citric and malic acid levels, the non-acidic recombinants had a higher uridine content than did the acidic recombinants. Additional minor QTLs for pH, citric acid and malic acid were also identified and for these the increased acidity was unexpectedly contributed by the non-sour parent. To test for co-localization of these QTLs with genes encoding organic acid metabolism and transport, we mapped the genes encoding structural enzymes and proteins involved in organic acid metabolism, transport and vacuolar H+ pumps. None of these genes co-localized with the major pH QTL, indicating that the gene determining melon fruit pH is not one of the candidate genes encoding this primary metabolic pathway. Linked markers were tested in two additional inter-varietal populations and shown to be linked to the pH trait. The presence of the same QTL in such diverse segregating populations suggests that the trait is determined throughout the species by variability in the same gene and is indicative of a major role of the evolution of this gene in determining the important domestication trait of fruit acidity within the species.

A genetic map of melon highly enriched with fruit quality QTLs and EST markers, including sugar and carotenoid metabolism genes.

A genetic map of melon enriched for fruit traits was constructed, using a recombinant inbred (RI) population developed from a cross between representatives of the two subspecies of Cucumis melo L.: PI 414723 (subspecies agrestis) and 'Dulce' (subspecies melo). Phenotyping of 99 RI lines was conducted over three seasons in two locations in Israel and the US. The map includes 668 DNA markers (386 SSRs, 76 SNPs, six INDELs and 200 AFLPs), of which 160 were newly developed from fruit ESTs. These ESTs include candidate genes encoding for enzymes of sugar and carotenoid metabolic pathways that were cloned from melon cDNA or identified through mining of the International Cucurbit Genomics Initiative database (http://www.icugi.org/). The map covers 1,222 cM with an average of 2.672 cM between markers. In addition, a skeleton physical map was initiated and 29 melon BACs harboring fruit ESTs were localized to the 12 linkage groups of the map. Altogether, 44 fruit QTLs were identified: 25 confirming QTLs described using other populations and 19 newly described QTLs. The map includes QTLs for fruit sugar content, particularly sucrose, the major sugar affecting sweetness in melon fruit. Six QTLs interacting in an additive manner account for nearly all the difference in sugar content between the two genotypes. Three QTLs for fruit flesh color and carotenoid content were identified. Interestingly, no clear colocalization of QTLs for either sugar or carotenoid content was observed with over 40 genes encoding for enzymes involved in their metabolism. The RI population described here provides a useful resource for further genomics and metabolomics studies in melon, as well as useful markers for breeding for fruit quality.

Mapping of a gene for Parkinson's disease to chromosome 4q21-q23.

Parkinson's disease (PD) is the second most common neurodegenerative disorder after Alzheimer's disease, affecting approximately 1 percent of the population over age 50. Recent studies have confirmed significant familial aggregation of PD and a large number of large multicase families have been documented. Genetic markers on chromosome 4q21-q23 were found to be linked to the PD phenotype in a large kindred with autosomal dominant PD, with a Zmax = 6.00 for marker D4S2380. This finding will facilitate identification of the gene and research on the pathogenesis of PD.

A 4,103 marker integrated physical and comparative map of the horse genome.

A comprehensive second-generation whole genome radiation hybrid (RH II), cytogenetic and comparative map of the horse genome (2n = 64) has been developed using the 5000rad horse x hamster radiation hybrid panel and fluorescence in situ hybridization (FISH). The map contains 4,103 markers (3,816 RH; 1,144 FISH) assigned to all 31 pairs of autosomes and the X chromosome. The RH maps of individual chromosomes are anchored and oriented using 857 cytogenetic markers. The overall resolution of the map is one marker per 775 kilobase pairs (kb), which represents a more than five-fold improvement over the first-generation map. The RH II incorporates 920 markers shared jointly with the two recently reported meiotic maps. Consequently the two maps were aligned with the RH II maps of individual autosomes and the X chromosome. Additionally, a comparative map of the horse genome was generated by connecting 1,904 loci on the horse map with genome sequences available for eight diverse vertebrates to highlight regions of evolutionarily conserved syntenies, linkages, and chromosomal breakpoints. The integrated map thus obtained presents the most comprehensive information on the physical and comparative organization of the equine genome and will assist future assemblies of whole genome BAC fingerprint maps and the genome sequence. It will also serve as a tool to identify genes governing health, disease and performance traits in horses and assist us in understanding the evolution of the equine genome in relation to other species.

First radiation hybrid map of the river buffalo X chromosome (BBUX) and comparison with BTAX.

We report the first radiation hybrid map of the river buffalo X chromosome generated from a recently constructed river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel (BBURH(5000)). This map contains a total of 33 cattle-derived markers, including 10 genes, four ESTs and 19 microsatellites. The markers are distributed in two linkage groups: LG1 contains eight markers spanning 125.6 cR, and LG2 contains 25 markers spanning 366.3 cR. LG1 contains six markers in common with bovine sequence assembly build 3.1. With the exception of BMS2152, the order of these markers on our BBUX map is shuffled when compared to the cow X chromosome (Bos taurus; BTAX). From LG2, two markers (AMELX and BL22) map to a more distal portion of BTAX compared to BBUX. In addition, two pairs of LG2 markers exhibit inversions compared to BTAX (ILSTS017 and ATRX; XBM38 and PPEF1). Alternatively, when compared to the most recent bovine RH map (Bov-Gen 3000rads), BL1098 and BMS2227 from LG1 as well as PLS3 and BMS1820 from LG2 showed inverted positions on the BBUX map. These discrepancies in buffalo and cattle maps may reflect evolutionary divergence of the chromosomes or mapping errors in one of the two species. Although the set of mapped markers does not cover the entire X chromosome, this map is a starting point for the construction of a high-resolution map, which is necessary for characterization of small rearrangements that might have occurred between the Bubalus bubalis and Bos taurus X chromosomes.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 1 (BBU1).

The largest chromosome in the river buffalo karyotype, BBU1, is a submetacentric chromosome with reported homology between BBU1q and bovine chromosome 1 and between BBU1p and BTA27. We present the first radiation hybrid map of this chromosome containing 69 cattle derived markers including 48 coding genes, 17 microsatellites and four ESTs distributed in two linkage groups spanning a total length of 1330.1 cR(5000). The RH map was constructed based on analysis of a recently developed river buffalo-hamster whole genome radiation hybrid (BBURH(5000)) panel. The retention frequency of individual markers across the panel ranged from 17.8 to 52.2%. With few exceptions, the order of markers within linkage groups is identical to the order established for corresponding cattle RH maps. The BBU1 map provides a starting point for comparison of gene order rearrangements between river buffalo chromosome 1 and its bovine homologs.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 7 and comparative mapping to the cattle and human genomes.

A preliminary radiation hybrid (RH) map containing 50 loci on chromosome 7 of the domestic river buffalo Bubalus bubalis (BBU; 2n = 50) was constructed based on a comparative mapping approach. The RH map of BBU7 includes thirty-seven gene markers and thirteen microsatellites. All loci have been previously assigned to Bos taurus (BTA) chromosome BTA6, which is known for its association with several economically important milk production traits in cattle. The map consists of two linkage groups spanning a total length of 627.9 cR(5,000). Comparative analysis of the BBU7 RH(5,000) map with BTA6 in cattle gave new evidence for strong similarity between the two chromosomes over their entire length and exposed minor differences in locus order. Comparison of the BBU7 RH(5,000) map with the Homo sapiens (HSA) genome revealed similarity with a large chromosome segment of HSA4. Comparative analysis of loci in both species revealed more variability than previously known in gene order and several chromosome rearrangements including centromere relocation. The data obtained in our study define the evolutionarily conserved segment on BBU7 and HSA4 to be between 3.5 megabases (Mb) and 115.8 Mb in the HSA4 (genome build 36) DNA sequence.

A 1.3-Mb interval map of equine homologs of HSA2.

A comparative approach that utilizes information from more densely mapped or sequenced genomes is a proven and efficient means to increase our knowledge of the structure of the horse genome. Human chromosome 2 (HSA2), the second largest human chromosome, comprising 243 Mb, and containing 1246 known genes, corresponds to all or parts of three equine chromosomes. This report describes the assignment of 140 new markers (78 genes and 62 microsatellites) to the equine radiation hybrid (RH) map, and the anchoring of 24 of these markers to horse chromosomes by FISH. The updated equine RH maps for ECA6p, ECA15, and ECA18 resulting from this work have one, two, and three RH linkage groups, respectively, per chromosome/chromosome-arm. These maps have a three-fold increase in the number of mapped markers compared to previous maps of these chromosomes, and an increase in the average marker density to one marker per 1.3 Mb. Comparative maps of ECA6p, ECA15, and ECA18 with human, chimpanzee, dog, mouse, rat, and chicken genomes reveal blocks of conserved synteny across mammals and vertebrates.

Improving the accuracy of PSI-BLAST protein database searches with composition-based statistics and other refinements.

PSI-BLAST is an iterative program to search a database for proteins with distant similarity to a query sequence. We investigated over a dozen modifications to the methods used in PSI-BLAST, with the goal of improving accuracy in finding true positive matches. To evaluate performance we used a set of 103 queries for which the true positives in yeast had been annotated by human experts, and a popular measure of retrieval accuracy (ROC) that can be normalized to take on values between 0 (worst) and 1 (best). The modifications we consider novel improve the ROC score from 0.758 +/- 0.005 to 0.895 +/- 0.003. This does not include the benefits from four modifications we included in the 'baseline' version, even though they were not implemented in PSI-BLAST version 2.0. The improvement in accuracy was confirmed on a small second test set. This test involved analyzing three protein families with curated lists of true positives from the non-redundant protein database. The modification that accounts for the majority of the improvement is the use, for each database sequence, of a position-specific scoring system tuned to that sequence's amino acid composition. The use of composition-based statistics is particularly beneficial for large-scale automated applications of PSI-BLAST.

Construction of evolutionary tree models for renal cell carcinoma from comparative genomic hybridization data.

Renal cell carcinoma is characterized by an accumulation of complex chromosomal alterations during tumor progression. Chromosome 3p deletions are known to occur early in the carcinogenesis, but the nature of subsequent events, their interrelationships, and their sequence is poorly understood, as one usually only obtains a single "view" of the dynamic process of tumor development in a particular cancer patient. To address this limitation, we used comparative genomic hybridization analysis in combination with a distance-based and a branching-tree method to search for tree models of the oncogenesis process of 116 conventional (clear cell) renal carcinomas. This provides a means to analyze and model cancer development processes based on a more dynamic model, including the presence of multiple pathways, as compared with the fixed linear model first proposed by Vogelstein et al. (N. Engl. J. Med., 319: 525-532, 1988) for colorectal cancer. The most common DNA losses involved 3p (61%), 4q (50%), 6q (40%), 9p (35%), 13q (37%), and Xq (21%). The most common gains were seen at chromosome 17p and 17q (20%). The tree model derived from the distance-based method is consistent with the established theory that -3p is an important early event in conventional (clear cell) renal cancer and supports the prediction made from the branching tree that -4q is another important early event. Both tree models suggest that there may be two groups of clear cell renal cancers: one characterized by -6q, +17q, and + 17p, and another by -9p, -13q, and -18q. Putative prognostic parameters were -9p and -13q. The distance-based tree clarifies that -8p (present in 12% of tumors) is a late event, largely independent of other events. In summary, tree modeling of comparative genomic hybridization data provided new information on the interrelationships of genetic changes in renal cancer and their possible order, as well as a clustering of these events. Using tree analysis, one can derive a more in-depth understanding of the renal cancer development process than is possible by simply focusing on the frequencies of genetic events in a given cancer type.

Genetic modifiers of the insulin resistance phenotype in mice.

Insulin resistance can result from genetic interactions among susceptibility alleles. To identify genetic loci predisposing to insulin resistance, we used crosses between different strains of mice with a targeted null allele of the insulin receptor gene. On the genetic background of B6 mice, the insulin receptor gene mutation causes mild hyperinsulinemia. In contrast, on the genetic background of 129/Sv mice, the same mutation causes severe hyperinsulinemia, suggesting that the 129/Sv strain harbors alleles that interact with the insulin receptor mutation and predispose to insulin resistance. As a first step to identify these alleles, we generated an F2 intercross between insulin receptor heterozygous mutant mice on B6 and 129/Sv backgrounds (B6IR x 129IR) and performed a genome-wide scan with polymorphic markers at a 20-cM resolution. We report the identification of loci on chromosomes 2 (logarithm of odds [LOD] 5.58) and 10 (LOD 5.58) that show significant evidence for linkage to plasma insulin levels as a quantitative trait. These findings indicate that targeted mutations in knockout mice can be used to unravel the complex genetic interactions underlying insulin resistance.

Sucrose-to-Starch Metabolism in Tomato Fruit Undergoing Transient Starch Accumulation.

Immature green tomato (Lycopersicon esculentum) fruits undergo a period of transient starch accumulation characterized by developmental changes in the activities of key enzymes in the sucrose (Suc)-to-starch metabolic pathway. Activities of Suc synthase, fructokinase, ADP-glucose (Glc) pyrophosphorylase, and soluble and insoluble starch synthases decline dramatically in parallel to the decrease in starch levels in the developing fruit. Comparison of "maximal" in vitro activities of the enzymes in the Suc-to-starch pathway suggests that these same enzymes are limiting to the rate of starch accumulation. In contrast, activities of invertase, UDP-Glc pyrophosphorylase, nucleoside diphosphate kinase, phosphoglucoisomerase, and phosphoglucomutase do not exhibit dramatic decreases in activity and appear to be in excess of starch accumulation rates. Starch accumulation is spatially localized in the inner and radial pericarp and columella, whereas the outer pericarp and seed locule contain little starch. The seed locule is characterized by lower activities of Suc synthase, UDP-Glc pyrophosphorylase, phosphoglucomutase, ADP-Glc pyrophosphorylase, and soluble and insoluble starch synthases. The outer pericarp exhibits comparatively lower activities of ADP-Glc pyrophosphorylase and insoluble starch synthase only. These data are discussed in terms of the developmental and tissue-specific coordinated control of Suc-to-starch metabolism.

Improving the practical space and time efficiency of the shortest-paths approach to sum-of-pairs multiple sequence alignment.

The MSA program, written and distributed in 1989, is one of the few existing programs that attempts to find optimal alignments of multiple protein or DNA sequences. The MSA program implements a branch-and-bound technique together with a variant of Dijkstra's shortest paths algorithm to prune the basic dynamic programming graph. We have made substantial improvements in the time and space usage of MSA. The improvements make feasible a variety of problem instances that were not feasible previously. On some runs we achieve an order of magnitude reduction in space usage and a significant multiplicative factor speedup in running time. To explain how these improvements work, we give a much more detailed description of MSA than has been previously available. In practice, MSA rarely produces a provably optimal alignment and we explain why.

Approximation algorithms for a genetic diagnostics problem.

We define and study a combinatorial problem called WEIGHTED DIAGNOSTIC COVER (WDC) that models the use of a laboratory technique called genotyping in the diagnosis of an important class of chromosomal aberrations. An optimal solution to WDC would enable us to define a genetic assay that maximizes the diagnostic power for a specified cost of laboratory work. We develop approximation algorithms for WDC by making use of the well-known problem SET COVER for which the greedy heuristic has been extensively studied. We prove worst-case performance bounds on the greedy heuristic for WDC and for another heuristic we call directional greedy. We implemented both heuristics. We also implemented a local search heuristic that takes the solutions obtained by greedy and dir-greedy and applies swaps until they are locally optimal. We report their performance on a real data set that is representative of the options that a clinical geneticist faces for the real diagnostic problem. Many open problems related to WDC remain, both of theoretical interest and practical importance.

Distance-based reconstruction of tree models for oncogenesis.

Comparative genomic hybridization (CGH) is a laboratory method to measure gains and losses in the copy number of chromosomal regions in tumor cells. It is hypothesized that certain DNA gains and losses are related to cancer progression and that the patterns of these changes are relevant to the clinical consequences of the cancer. It is therefore of interest to develop models which predict the occurrence of these events, as well as techniques for learning such models from CGH data. We continue our study of the mathematical foundations for inferring a model of tumor progression from a CGH data set that we started in Desper et al. (1999). In that paper, we proposed a class of probabilistic tree models and showed that an algorithm based on maximum-weight branching in a graph correctly infers the topology of the tree, under plausible assumptions. In this paper, we extend that work in the direction of the so-called distance-based trees, in which events are leaves of the tree, in the style of models common in phylogenetics. Then we show how to reconstruct the distance-based trees using tree-fitting algorithms developed by researchers in phylogenetics. The main advantages of the distance-based models are that 1) they represent information about co-occurrences of all pairs of events, instead of just some pairs, 2) they allow quantitative predictions about which events occur early in tumor progression, and 3) they bring into play the extensive methodology and software developed in the context of phylogenetics. We illustrate the distance-based tree method and how it complements the branching tree method, with a CGH data set for renal cancer.

Inferring tree models for oncogenesis from comparative genome hybridization data.

Comparative genome hybridization (CGH) is a laboratory method to measure gains and losses of chromosomal regions in tumor cells. It is believed that DNA gains and losses in tumor cells do not occur entirely at random, but partly through some flow of causality. Models that relate tumor progression to the occurrence of DNA gains and losses could be very useful in hunting cancer genes and in cancer diagnosis. We lay some mathematical foundations for inferring a model of tumor progression from a CGH data set. We consider a class of tree models that are more general than a path model that has been developed for colorectal cancer. We derive a tree model inference algorithm based on the idea of a maximum-weight branching in a graph, and we show that under plausible assumptions our algorithm infers the correct tree. We have implemented our methods in software, and we illustrate with a CGH data set for renal cancer.

Radiation hybrid mapping of 304 novel microsatellites in the domestic cat genome.

Effective utilization of the domestic cat as an animal model for hereditary and infectious disease requires the development and implementation of high quality gene maps incorporating microsatellites and conserved coding gene markers. Previous feline linkage and radiation hybrid maps have lacked sufficient microsatellite coverage on all chromosomes to make effective use of full genome scans. Here we report the isolation and genomic mapping of 304 novel polymorphic repeat loci in the feline genome. The new loci were mapped in the domestic cat radiation hybrid panel using an automated fluorescent TAQ-Man based assay. The addition of these 304 microsatellites brings the total number of microsatellites mapped in the feline genome to 580, and the total number of loci placed onto the RH map to 1,126. Microsatellites now span every autosome with an average spacing of roughly one polymorphic STR every five centimorgans, and full genome coverage of one marker every 2.7 megabases. These loci now provide a useful tool for undertaking full-genome scans to identify genes associated with phenotypes of interest, such as those relating to hereditary disease, coat color, patterning and morphology. These resources can also be extended to the remaining 36 species of the cat family for population genetic and evolutionary genomic analyses.

Towards a complete North American Anabaptist genealogy: A systematic approach to merging partially overlapping genealogy resources.

Genetic studies on closed populations benefit from comprehensive, searchable genealogy resources. To support studies of Anabaptists, a computerized database has been developed that merges two large genealogy books, the "Fisher Family History" and the "Amish and Amish Mennonite Genealogies." The former is more current but the latter is more comprehensive for the 18th and 19th centuries. Therefore, the merger of the two books is significantly more useful than either book alone. We demonstrate the utility of the merged database with two results: an increase in inbreeding coefficients and the identification of parent-child relationships that do not exist in either book. Am. J. Med. Genet. 86:156-161, 1999. Published 1999 Wiley-Liss, Inc.