A confidence set inference procedure for gene mapping using markers with incomplete polymorphism

A recent approach for gene mapping based on confidence set inference (CSI) promises several advantages, including avoidance of corrections for multiple tests, availability of confidence intervals with known statistical properties, and sufficient localizations of disease genes. This paper proposes an extended CSI procedure that can handle markers with incomplete polymorphism, thereby increasing the applicability of the set of CSI methods in practical situations. Simulation studies show that the new procedure retains the main advantages of the original CSI. Although it generally requires more data to achieve a similar power, this increase is moderate for markers with 80% heterozygosity or higher. We also investigate the effects of relative risk estimates and disease models. Our analyses show that perturbation from actual relative risks or multilocus disease models generally leads to reduction in power or inflation in type I error, as expected. Nevertheless, for certain classes of two-locus disease models, CSI can still perform well, with reasonably high actual coverage probabilities for at least one of the disease loci. Application of CSI to the data provided by the Genetic Analysis Workshop 13 yields encouraging results, as they compare favorably to those obtained from GENEHUNTER using its NPL sib-pair method.     

Characterization and regional mapping of new anonymous chromosome 20-specific DNA markers isolated from a flow-sorted DNA library

Employing the flow-sorted chromosome 20-specific DNA library LL20NS01, we isolated seven novel unique poly- and monomorphic DNA markers specific to human chromosome 20. Initially, 201 phage clones were analyzed regarding insert size and repetitivity. By testing 14 single- and low-copy number clones for their ability to detect RFLPs, three polymorphisms were revealed by two probes, pFMS22-1.4 [D20S22] and pFMS76 [D20S23]. Seven of twenty probes (35%) were assigned to chromosome 20 using a somatic cell hybrid DNA panel. Five of them were regionally mapped by in situ hybridization. Three DNA markers, pFMS51 [D20S29], pFMS76 [D20S23], and pFMS106 [D20S30], were assigned to 20p11.2-p12, and two markers, pFMS22-1.4 [D20S22] and pFMS135 [D20S31], to 20q12-q13.3. Our new chromosome 20-specific DNA markers should be useful for the molecular characterization of this rather underpopulated human chromosome.     

A mapped set of DNA markers for human chromosome 17

We have developed and mapped by genetic linkage a primary set of markers for chromosome 17. The map consists of 21 loci derived from 27 probe/enzyme systems, including eight highly informative markers at loci containing a variable number of tandemly repeated DNA sequences (VNTRs). The map is continuous from the telomeric region of the short arm to the telomeric region of the long arm, covering estimated genetic distances of 218 cM in males and 279 cM in females. The average heterozygosity among all 21 loci in the population sample analyzed is 58%; 77% heterozygosity was observed among the eight VNTR markers that were highly informative. This map will make it possible to detect by linkage the location of genetic defects associated with chromosome 17 and will also provide anchor points for a high-resolution map of this chromosome.     

Additional microsatellite markers for mouse genome mapping

Mouse sequence information from the EMBL and GenBank databases, published sequences and genomic clones have been analyzed for simple repetitive elements or microsatellites. Each microsatellite has been amplified by the polymerase chain reaction (PCR) as a single locus marker. PCR primers were designed from unique sequence flanking each repeat. Size variation of PCR products less than 750 base pairs (bp) between mouse strains has been determined using ethidium bromide-stained acrylamide or agarose gels. A further 74 newly characterized microsatellites are presented in this paper, bringing to 185 the total we have analyzed. Of these, 157/185 (85%) have more than one allele, 143/178 (80%) vary in length between C57BL/6J and Mus spretus, and 82/168 (49%) vary between DBA/2J and C57BL/6J. Microsatellites provide informative single locus probes for linkage analysis in the construction of a genetic map of the mouse genome.     

Polymorphic molecular markers from anonymous nuclear DNA for genetic analysis of populations

A simple method is presented for developing polymorphic, anonymous DNA markers suitable for population genetic studies. Anonymous DNA fragments are screened for sequence variability using a common mutation detection technique (single strand conformation polymorphism analysis; SSCP) and locus-specific PCR primers are designed for polymorphic DNA fragments. Detection of the markers by SSCP analysis coupled with sequence analysis of SSCP variants allows rapid screening while retaining information about the genealogical relationship among alleles. Variability detected for six markers was assessed in rainbow trout Oncorhynchus mykiss and was compared with variability detected by similar analysis of intron loci. Between three and 12 distinct alleles were observed at each marker locus, and average within-population heterozygosity ranged from 0.12 to 0.44. Advantages and limitations of the methodology for population genetic analysis are discussed.     

Fine mapping of gene probes and anonymous DNA fragments to the long arm of chromosome 16

The fragile site, FRA16B, at 16q22.100 and four different translocations with breakpoints at 16q22.102, 16q22.105, 16q22.108, and 16q22.3 were used to locate and order DNA probes. This was achieved by Southern analysis of a somatic cell hybrid panel containing portions of chromosome 16 and by in situ hybridization. The anonymous DNA fragments D16S6, D16S10, and D16S11 were proximal to FRA16B and located at 16q13----q22.100. D16S4 and LCAT were located at 16q22.100----q22.102. TAT and HP were located at 16q22.105----q22.108. CTRB was located distal to 16q22.105 and therefore is in the distal half of 16q22. The order of markers in this region was determined as centromere-D16S6, D16S11, D16S10, MT-FRA16B-D16S4, LCAT-HP,TAT,CTRB-APRT- telomere. Linkage studies to determine map distances between the closest markers flanking the fragile site are now in progress.     

Complex probes for high-throughput parallel genetic mapping of genomic mouse BAC clones

We describe a novel approach for the identification and mapping of polymorphic markers. Amplicons are generated by ligation of double-stranded adaptor molecules to genomic DNA cleaved with a restriction enzyme. Using primers that extend beyond the restriction site, reduced-complexity subsets of fragments are generated by PCR. Differences in the composition of complex probes generated from DNA of different strains are revealed through hybridization against high-density filter grids of large-insert genomic clones. Genetic mapping of genomic clones is achieved by hybridizing complex probes derived from backcross animals against the polymorphic clones. The mouse was chosen as a model system to test the feasibility of this technique because of the general availability of backcross resources and genomic libraries. Nevertheless, we would expect the method to be of particular use to generate markers for species that have not yet been extensively studied, because a substantial number of easy-to-use markers can be recruited in a relatively short period of time. Received: 20 January 1998 / Accepted: 21 April 1998     

Isolation and regional localization of 35 unique anonymous DNA markers for human chromosome 22

Thirty-five new, unique, DNA probes have been isolated and each has been assigned to one of five regions on chromosome 22. The distribution of probes along the chromosome is what would be expected based on the estimated size of each region with the exception of the short arm (22p). RFLP analysis was performed using 13 different restriction enzymes and over 50% of the probes were found to have useful polymorphisms. Probes mapping to 22q11 were further characterized by pulsed-field gel analysis and it has been possible to link several on large restriction fragments. These 35 new probes will be useful reagents for producing genetic and physical maps of chromosome 22.     

Construction and mapping of safflower chloroplast DNA recombinants and location of selected gene markers

Summary DNA was isolated and purified from chloroplasts of safflower (Carthamus tinctorius L.), digested with HindIII restriction endonuclease, and ligated into the HindIII site of the plasmid pUC9. Recombinant DNAs were isolated from ampicillin resistant white colonies which grew in the presence of the appropriate indicator, digested with HindIII, and then identified by comparison of agarose gel electrophoretic mobilities. HindIII digests of chloroplast DNA were used as a standard. Such recombinants were radiolabeled and hybridized with Southern blots of PstI, SalI, KpnI, and HindIII single and double digests of safflower chloroplast DNA. A physical map was subsequently generated showing the location of each recombinant on the circular plastid genome. Recombinants containing heterologous chloroplast gene markers from spinach or Euglena were also radiolabeled and mapped. The relative mapping positions of these genes are in good agreement with those which have previously been published for spinach and several other higher plants.     

Interpreting anonymous DNA samples from mass disasters--probabilistic forensic inference using genetic markers

MOTIVATION: The problem of identifying victims in a mass disaster using DNA fingerprints involves a scale of computation that requires efficient and accurate algorithms. In a typical scenario there are hundreds of samples taken from remains that must be matched to the pedigrees of the alleged victim's surviving relatives. Moreover the samples are often degraded due to heat and exposure. To develop a competent method for this type of forensic inference problem, the complicated quality issues of DNA typing need to be handled appropriately, the matches between every sample and every family must be considered, and the confidence of matches need to be provided. RESULTS: We present a unified probabilistic framework that efficiently clusters samples, conservatively eliminates implausible sample-pedigree pairings, and handles both degraded samples (missing values) and experimental errors in producing and/or reading a genotype. We present a method that confidently exclude forensically unambiguous sample-family matches from the large hypothesis space of candidate matches, based on posterior probabilistic inference. Due to the high confidentiality of disaster DNA data, simulation experiments are commonly performed and used here for validation. Our framework is shown to be robust to these errors at levels typical in real applications. Furthermore, the flexibility in the probabilistic models makes it possible to extend this framework to include other biological factors such as interdependent markers, mitochondrial sequences, and blood type. AVAILABILITY: The software and data sets are available from the authors upon request.     

Isolation of anonymous DNA markers for human chromosome 22q11 from a flow-sorted library,and mapping using hybrids from patients with DiGeorge syndrome

Summary DiGeorge syndrome (DGS) is a human developmental defect of the structures derived from the third and fourth pharyngeal pouches. It apparently arises due to deletion of 22q11. We describe a strategy for the isolation of DNA probes for this region. A deleted chromosome 22, which includes 22q11, was flow-sorted from a lymphoblastoid cell line of a patient with cat eye syndrome and used as the source of DNA. A DNA library was constructed from this chromosome by cloning into the EcoR1 site of the vector Lambda gt10. Inserts were amplified by PCR and mapped using a somatic cell hybrid panel of this region. Out of 32 probes, 14 were mapped to 22q11. These probes were further sublocalised within the region by dosage analysis of DGS patients, and by the use of two new hybrid cell lines which we have produced from DGS patients. One of these lines (7939B662) contains the altered human chromosome segregated from its normal homologue. This chromosome 22 contains an interstitial deletion in 22q11, and will be useful for localising further probes to the DGS region.     

Genetic mapping of 40 cDNA clones on the mouse genome by PCR

We recently proposed a new PCR-based genetic marker assay for the mouse genome that exploits sequence differences in the 3-untranslated region (UTR) of cDNAs between different mouse strains, called biallelic polymorphic expressed sequence tags (bESTs). The specific use of 3-UTR has several advantages: (1) frequent sequence polymorphism between different mouse strains, (2) most commonly uninterrupted by introns, (3) usually unique sequence even among closely related gene family members. In this paper, we identify additional genetic loci defined by bEST and determine their location on the mouse genetic map by using interspecific backross mapping panels between C57BL/6J and Mus spretus. Of 136 markers tested, 86 produced unique PCR products from C57BL/6J and M. spretus genomic DNAs. We then sequenced these 86 PCR products from C57BL/6J and M. spretus and found that 59 markers have sequence polymorphisms. Of these, we mapped 36 by restriction fragment length polymorphism (RFLP) of the PCR products and 4 by length polymorphism (LP) of the PCR products. We discuss the possibility of a large-scale application of this method for cDNA mapping.     

Genetic mapping of Y-chromosomal DNA markers in Pacific salmon

Sex chromosomes in fish provide an intriguing view of how sex-determination mechanisms evolve in vertebrates. Many fish species with single-factor sex-determination systems do not have cytogenetically-distinguishable sex chromosomes, suggesting that few sex-specific sequences or chromosomal rearrangements are present and that sex-chromosome evolution is thus at an early stage. We describe experiments examining the linkage arrangement of a Y-chromosomal GH pseudogene (GH-Y) sequence in four species of salmon (chum, Oncorhynchus keta; pink, O. gorbuscha; coho, O. kisutch; chinook, O. tshawytscha). Phylogenetic analysis indicates that GH-Y arose early in Oncorhynchus evolution, after this genus had diverged from Salmo and Salvelinus. However, GH-Y has not been detected in some Oncorhynchus species (O. nerka, O. mykiss and O. clarki), consistent with this locus being deleted in some lineages. GH-Y is tightly linked genetically to the sex-determination locus on the Y chromosome and, in chinook salmon, to another Y-linked DNA marker OtY1. GH-Y is derived from an ancestral GH2 gene, but this latter functional GH locus is autosomal or pseudoautosomal. YY chinook salmon are viable and fertile, indicating the Y chromosome is not deficient of vital genetic functions present on the X chromosome, consistent with sex chromosomes that are in an early stage of divergence.     

A novel set of single-copy nuclear gene markers in white oak and implications for species delimitation

Oaks have often been the focus of research on plant evolution owing to their propensity to intercross and their important role in ecology and economy. Compared with traditional molecular markers, such as amplified fragment length polymorphisms (AFLPs) and simple sequence repeats (SSRs), multiple single-copy nuclear genes (SCNGs) are of greater utility in inferring evolutionary processes in oaks. Nineteen primer pairs were developed from expressed sequence tags (ESTs) of Quercus mongolica and Q. robur that could produce orthologous products in Chinese white oaks (section Quercus). These SCNG markers showed a moderate to high level of nucleotide polymorphism in 42 individuals of two closely related white oaks, Q. mongolica and Q. liaotungensis, and demonstrated high transferability across seven white oaks, four oaks from section Cerris, and one oak from section Lobatae. A phylogenetic tree based on these SCNGs provided resolution at deep nodes and robust support for delimiting populations of Q. mongolica and Q. liaotungensis; Bayesian analysis clustered individuals into their respective species with high probability and no admixture. When the same individuals were used, Bayesian clusters based on either 194 AFLPs or 19 SSRs gave comparable results, but one or several individuals respectively were identified as having admixed ancestry. This indicates that the hybridization rate between these two oaks may have been overestimated using SSR markers due to the occurrence of homoplasy. The SCNGs are powerful for species delimitation of white oaks, and these markers could be useful for future phylogenetics and phylogeography research in white oaks.     

Strategies for characterizing highly polymorphic markers in human gene mapping.

Before new markers are thoroughly characterized, they are usually screened for high polymorphism on the basis of a small panel of individuals. Four commonly used screening strategies are compared in terms of their power to correctly classify a marker as having heterozygosity of 70% or higher. A small number of typed individuals (10, say) are shown to provide good discrimination power between low- and high-heterozygosity markers when the markers have a small number of alleles. Characterizing markers in more detail requires larger sample sizes (e.g., at least 80-100 individuals) if there is to be a high probability of detecting most or all alleles. For linkage analyses involving highly polymorphic markers, the practice of arbitrarily assuming equal gene frequencies can cause serious trouble. In the presence of untyped individuals, when gene frequencies are unequal but are assumed to be equal in the analysis, recombination-fraction estimates tend to be badly biased, leading to strong false-positive evidence for linkage.     

A set of BAC clones spanning the human genome

Using the human bacterial artificial chromosome (BAC) fingerprint-based physical map, genome sequence assembly and BAC end sequences, we have generated a fingerprint-validated set of 32855 BAC clones spanning the human genome. The clone set provides coverage for at least 98% of the human fingerprint map, 99% of the current assembled sequence and has an effective resolving power of 79 kb. We have made the clone set publicly available, anticipating that it will generally facilitate FISH or array-CGH-based identification and characterization of chromosomal alterations relevant to disease.