The effects of inbreeding on DNA profile frequency estimates using PCR-based loci

Estimates of inbreeding were determined using Wright's F_ST for loci used for PCR-based forensic analyses. The populations analyzed were African Americans, Caucasians, Hispanics, and Orientals. In most cases the F_ST values at each locus were less than 0.01. The F_ST values over all loci for African Americans, Caucasians, and Orientals ranged from 0.0015 to 0.0048. No substantial differences were observed for DNA profile frequency estimates when calculated under the assumption of independence or with the incorporation of F_ST.Editor's commentsThe author provides an empirical study to accompany the papers by Balding and Nichols and Brookfield in this volume. As may have been expected, F_ST values for large heterogeneous populations such as Caucasian or African American are very small.     

The evidential value in the DNA database search controversy and the two-stain problem

Summary . Does the evidential strength of a DNA match depend on whether the suspect was identified through database search or through other evidence (“probable cause”)? In Balding and Donnelly (1995, Journal of the Royal Statistical Society, Series A 158, 21–53) and elsewhere, it has been argued that the evidential strength is slightly larger in a database search case than in a probable cause case, while Stockmarr (1999 , Biometrics 55, 671–677) reached the opposite conclusion. Both these approaches use likelihood ratios. By making an excursion to a similar problem, the two‐stain problem, we argue in this article that there are certain fundamental difficulties with the use of a likelihood ratio, which can be avoided by concentrating on the posterior odds. This approach helps resolving the above‐mentioned conflict.     

Inferring relatedness and heritability using molecular markers in radiata pine

A set of eight unlinked microsatellite markers was used to estimate relatedness among 355 individuals of a Pinus radiata breeding population. The average performance of open-pollinated progeny of each individual, for wood density, was considered to represent the phenotype of all 355 individuals. Marker-based estimates of relationship were compared with the pedigree-based coefficients of relationships. The phenotypic similarity among all pairs of individuals was regressed on marker-estimated relatedness to estimate the inheritance of wood density. The marker-based estimate of heritability was compared with that obtained using classical quantitative genetic methods. Overall, a low correlation (0.13) was observed between marker-based and pedigree-based estimates of relatedness. After discarding negative estimates of relatedness, the average coefficient of relationship among known groups of maternal half-sibs, full-sibs and unrelated individuals, increased from 0.24 to 0.29 (0.25 expected), from 0.43 to 0.48 (0.50 expected) and from –0.04 to 0.15 (0 expected), respectively. Marker-based and conventional estimates of heritability of wood density were 0.79 and 0.38, respectively. However, by using only marker loci with expected Hardy–Weinberg frequencies, marker-based estimate of heritability was 0.33, which is very similar to that obtained from conventional approaches. The use of molecular markers to understand quantitative genetic variation is discussed.     

Sporadic cases in Duchenne muscular dystrophy

Summary A new estimation of the proportion of sporadic cases in Duchenne muscular dystrophy was attempted by means of segregation analysis in a sample of 988 sibships collected on a world-wide scale by different authors. Maximum likelihood estimates of ascertainment probability (), segregation frequency (p), and frequency of sporadic cases (x) were calculated by Morton's equations under different hypotheses. The best fit was found for p=0.454±0.024 and x=0.235±0.034. The possibility that the proportion of sporadic cases might be lower than the expected 1/3 is suggested.     

Evaluation of biometrical methods for estimating the number of genes

Summary Computer simulation of several genetic models was used to assess the effect of type I and type II statistical errors on estimating the number of genes by the inbred-backcross and genotype assay procedures. Depending upon the actual number of genes, heritability, and the probability of type I errors (), substantial upward and downward biases were observed in estimates of the number of genes from both methods. The estimated number of genes increased as was increased from 0.01 to 0.30 and as heritability increased. With high and/or high heritability, the estimated number of genes often exceeded the actual number. Downward biases occurred with low and low heritability, and tended to become greater as the number of genes in the model was increased. Large type II errors were associated with downward biases. The choice of had a greater impact on biases in estimates from the genotype assay procedure than from the inbred-backcross procedure. Increasing the number of backcrosses in the inbred-backcross procedure or delaying the assay generation in genotype assay increased the probability of upward biases in the estimated number of genes. Unbiased estimates can be obtained only by choice of an optimum . There is no known way to choose the optimum in practice. This fact reduces the value of estimates of the number of genes by genotype assay or by the inbred-backcross methods.Part of this study is based on the Ph.D. thesis of the senior author     

Evolutionary relationship and secondary structure predictions in four transport proteins ofSaccharomyces cerevisiae

Summary The comparison of the amino acid sequences of four yeast transport proteins indicates that there is a questionable relatedness between the uracil permease (FUR4) and the purine-cytosine permease (FCY2), whereas the arginine permease (CAN1) and the histidine permease (HIP1) clearly originated from a common molecular ancestor. The analysis of the primary structure of these transport proteins by two methods of secondary structure predictions suggests the presence of 9–12 membrane-spanning -helices in each polypeptide chain. These results are concordant in that 90% of the -helices were determined by both methods to be at the same positions. In the aligned sequencesHIP1 andCAN1, the postulated membrane-spanning -helices often start at corresponding sites, even though the overall sequence similarity of the two proteins is only 30%. In the aligned DNA coding sequences ofCAN1 andHIP1, synonymous nucleotide substitutions occur with very similar frequencies in regions where the replacement substitution (changing the amino acids) frequencies are widely different. Moreover, our data suggest that the replacement substitutions can be considered as neutral in the N-terminal segment, whereas the other regions are subject to a conservative selective pressure because, if compared to a random drift, the replacement substitutions are underrepresented.     

The Carbohydrate-containing Cell-Wall Polymers of Certain Species from the Cluster “Streptomyces lavendulae”

The type strains of the species of the cluster Streptomyces lavendulae with a low level of DNA–DNA relatedness were found to contain different cell-wall carbohydrate polymers, whereas the species of this cluster with a level of DNA–DNA relatedness of about 60% contain similar or identical carbohydrate polymers. The type strains Streptomyces katraeVKM Ac-1220^Tand S. polychromogenesVKM Ac-1207^Tsynthesize mannan with different amounts of -1,2- and -1,3-substituted mannopyranose units and a small number of 1,3-poly(glycerolphosphate) chains. The cell walls of S. lavendulocolorVKM Ac-215^Tand Streptomycessp. VKM Ac-2117 were found to contain a hitherto unknown teichuronic acid, whose repeating unit is a disaccharide consisting of diaminomannuronic acid and N-acetylgalactosamine: 4)--D-ManpNAc3NAcA-(1 3)--D-GalpNAc-(1 . In addition, the cell walls of these two streptomycetes contain -glucosylated 1,5-poly(ribitol phosphate). The cell walls of S. virginiaeVKM Ac-1218^Tand S. flavotriciniVKM Ac-1277^Tcontain the same poly(glucosyl-glycerolphosphate). The results presented in this paper are in accordance with the DNA–DNA relatedness data and indicate a taxonomic significance of the structure of the cell-wall polysaccharides for the delineation of phenetically related Streptomycesspecies.     

Relatedness measured by oligonucleotide probe DNA fingerprints and an estimate of the mating system of Sea Lavender (Limonium carolinianum)

Using DNA fingerprint markers within species and populations of wild plants requires information on the relationship between fingerprint similarity and relatedness. We identified a hypervariable marker based on oliog(GATA)₄-hybridization of DpnII-cut genomic DNA from Sea Lavender (Limonium carolinianum). Banding patterns were somatically stable and highly variable among unrelated individuals. Band molecular-weight sizing errors (as a percent of band molecular weight) were estimated at 0.44%±0.003 within gels and 0.76%±0.964 between gels. Band sizing errors defined a 99% confidence bin of ±0.95% (1.90% total) of molecular weight. Band-sharing estimates were based on this bin size and on variance estimates that compensate for non-independent comparisons. Band-sharing among nine unrelated individuals () was 0.198±0.O11. Experimental pollinations designed to produce selfed, fulland half-sib progeny groups led to five selfed progeny groups and no outcrossed progeny (mean band-sharing, ovS=0.468±0.074). A linear regression between band-sharing (S) and relatedness (r) assuming 17% inbreeding was r=0.006+0.914*S (R²=0.973) and established the maximum amount of inbreeding. ovS(0.392±0.022) estimated from wild pollinated seeds from four maternal families was intermediate to unrelated individuals and experimental selfed progeny, giving evidence for mixed mating in wild plants. More extensive plant pedigrees with known levels of inbreeding will be needed to measure variation in the relationship between S and r among populations and families.     

Assessing probability of paternity and the product rule in DNA systems

The genetic resolution of paternity disputes begins with an intricate detection of inherited traits and finishes with a statistical inference (the probability of paternity,W). Notwithstanding some initial fanfare, statistical inference is a necessary component of DNA-based paternity tests because band patterns may be rare but not yet unique, and even rare events in a vacuum are meaningless. The genetic match must be combined with other evidence for relevancy, thus a Bayesian approach is preferred when computingW. This paper reviews the standard model used to computeW and discusses the model's various properties and assumptions. The standard model is extended to include DNA systems in which alleles are operationally continuous due to measurement error. This extension avoids problems associated with matched/non-matched binned decisions. After outlining the model assumptions for a single DNA system, particular attention is given to the product rule-the procedure of multiplying intermediate probabilities across genetic loci to form a combinedW. An empirical alternative to the product rule is also assessed and correlated with standard procedures.Editor's commentsThe authors present a carefully reasoned account of Bayesian inference in the paternity context. Their presentation of an alternative to the product rule for Paternity Index values offers a means of avoiding objections to the rule.     

Localization of the metJBLF gene cluster of Escherichia coli in lambda met transducing phage

Summary The position of the metJBLF gene cluster in the transducing phage dmet102 was determined by ligation of its leftmost EcoRI fragment (102-1) to the BCDEF (nin5) EcoRI fragment of gtl (BC) and characterization of the resultant recombinant phage. The new transducing phage carries about 6kb of bacterial DNA which contains the entire met gene cluster including the promoter of its rightmost member metF. Reasonable estimates of the coding capacity required for the four genes indicate that most of the bacterial DNA of the recombinant phage is occupied by the met gene cluster.     

Nucleotide sequence comparison of the Adh gene in three drosophilids

Summary The alcohol dehydrogenase (Adh) gene has been isolated fromDrosophila simulans andD. mauritiana by screening clone libraries of each with a previously cloned Adh gene fromD. melanogaster. The isolated clones were subcloned and partially sequenced to determine the relatedness of these species and to examine details of evolutionary change in the structure of the Adh gene. We report the sequence of the first 704 nucleotides of each gene as well as 127 bases in the 5 untranslated region. When these sequences are compared,D. melanogaster differs fromD. simulans andD. mauritiana by 2.8% and 3.1%, respectively.D. simulans andD. mauritiana differ by only 1.8%, implying that they are more closely related to each other than either is toD. melanogaster. This is consistent with phylogenetic relationships established by a variety of genetic, biochemical, and morphological means and illustrates that DNA sequencing of a single gene may be used to assess the evolutionary relationships of species.     

Minimal entropy probability paths between genome families

We develop a metric for probability distributions with applications to biological sequence analysis. Our distance metric is obtained by minimizing a functional defined on the class of paths over probability measures on N categories. The underlying mathematical theory is connected to a constrained problem in the calculus of variations. The solution presented is a numerical solution, which approximates the true solution in a set of cases called rich paths where none of the components of the path is zero. The functional to be minimized is motivated by entropy considerations, reflecting the idea that nature might efficiently carry out mutations of genome sequences in such a way that the increase in entropy involved in transformation is as small as possible. We characterize sequences by frequency profiles or probability vectors, in the case of DNA where N is 4 and the components of the probability vector are the frequency of occurrence of each of the bases A, C, G and T. Given two probability vectors a and b, we define a distance function based as the infimum of path integrals of the entropy function H(p) over all admissible paths p(t), 0 t1, with p(t) a probability vector such that p(0)=a and p(1)=b. If the probability paths p(t) are parameterized as y(s) in terms of arc length s and the optimal path is smooth with arc length L, then smooth and rich optimal probability paths may be numerically estimated by a hybrid method of iterating Newtons method on solutions of a two point boundary value problem, with unknown distance L between the abscissas, for the Euler–Lagrange equations resulting from a multiplier rule for the constrained optimization problem together with linear regression to improve the arc length estimate L. Matlab code for these numerical methods is provided which works only for rich optimal probability vectors. These methods motivate a definition of an elementary distance function which is easier and faster to calculate, works on non–rich vectors, does not involve variational theory and does not involve differential equations, but is a better approximation of the minimal entropy path distance than the distance ||b–a||₂. We compute minimal entropy distance matrices for examples of DNA myostatin genes and amino-acid sequences across several species. Output tree dendograms for our minimal entropy metric are compared with dendograms based on BLAST and BLAST identity scores.Mathematics Subject Classification (2000): 92B05, 92D20     

Gene isolation through genomic complementation using an indexed library of Chlamydomonas reinhardtii DNA

Hundreds of mutants with defects in a variety of physiologically important functions, such as photosynthesis, respiration, flagellar motility, phototaxis, circadian rhythms and the cell cycle, have been isolated from cultures of Chlamydomonas reinhardtii. In only a few cases have the genes responsible for these mutations been cloned and sequenced. The development of efficient methods for transformation with nuclear genes [7] has allowed the recent demonstration of gene isolation through genomic complementation with a pooled library of C. reinhardtii DNA [9]. To improve the efficiency with which genes complementing a particular mutation can be isolated, we have established an indexed (ordered) cosmid library of 11,280 individual clones contained in the separate wells of 120 microtiter plates. The average insert size is ca. 38 kb. PCR analysis of five sequenced nuclear genes present in the Chlamydomonas library revealed a range from two copies for the ₂ and ₂ tubulin genes to at least seven copies for the agininosuccinate lyase gene. Overall, these five clones were represented an average of >-3.4 times in the library. Thus, the probability that any one particular nuclear gene of < 1000 bp will be found in the library is >-97%, and the probability that a gene of ca. 10 000 bp will be found in the library is ca. 92%. Rapid screening methods with cosmid DNAs pooled from individual microtiter dishes have been applied successfully. Bacteria containing clones of the argininosuccinate lyase gene have been identified through genomic complementation of a Chlamydomonas mutant bearing an inactive arginnosuccinate lyase gene.We are using the nomenclature of indexed library versus ordered library to avoid confusion of this library with a library of ordered contigs.     

Renibacterium salmoninarum isolates from different sources possess two highly conserved copies of the rRNA operon

The nucleotide sequences of the rRNA genes and the 5 flanking region were determined for R. salmoninarum ATCC 33209^T from overlapping products generated by PCR amplification from the genomic DNA. Comparison of the sequences with rRNA genes from a variety of bacteria demonstrated the close relatedness between R. salmoninarum and the high G+C group of the actinobacteria, in particular, Arthrobacter species. A regulatory element within the 5 leader of the rRNA operon was identical to an element, CL2, described for mycobacteria. PCR, DNA sequence analysis, and DNA hybridisation were performed to examine variation between isolates from diverse sources which represented the four 16S–23S rRNA intergenic spacer sequevars previously described for R. salmoninarum. Two 23S–5S rRNA intergenic spacer sequevars of identical length were found. DNA hybridisation using probes complementary to 23S rDNA and 16S rDNA identified two rRNA operons which were identical or nearly identical amongst 40 isolates sourced from a variety of countries.     

Preliminary quantification of immunological relationships among the transaldolases of the genus Bifidobacterium

The immunological relatedness among the transaldolases (dihydroxyacetone transferase, E.C. 2.2.1.2) of twenty species of the genus Bifidobacterium has been tested by the microcomplement fixation method, using B. thermophilum (B. ruminale) RU326 (=ATCC 25866), B. cuniculi RA93(=ATCC 27916) and B. minimum (DNA homology group) F392 (=ATCC 27538) as references. Based on the serological relationships of the transaldolases, expressed either as indices of dissimilarity or as immunological distances, the twenty species of the genus Bifidobacterium were arranged into clusters. These clusters generally coincided with the immunological groups obtained previously by the immunodiffusion method (Sgorbati and Scardovi, 1979).     

Genomic characterization of thermophilic Geobacillus species isolated from the deepest sea mud of the Mariana Trench

The thermophilic strains HTA426 and HTA462 isolated from the Mariana Trench were identified as Geobacillus kaustophilus and G. stearothermophilus, respectively, based on physiologic and phylogenetic analyses using 16S rDNA sequences and DNA–DNA relatedness. The genome size of HTA426 and HTA462 was estimated at 3.23–3.49 Mb and 3.7–4.49 Mb, respectively. The nucleotide sequences of three independent -phage inserts of G. stearothermophilus HTA462 have been determined. The organization of protein coding sequences (CDSs) in the two -phage inserts was found to differ from that in the contigs corresponding to each insert assembled by the shotgun clones of the G. kaustophilus HTA426 genome, although the CDS organization in another insert is identical to that in the HTA426 genome.     

Regional replication of the bacterial chromosome induced by derepression of prophage lambda

Summary We have demonstrated previously by DNA-DNA hybridization that induction of phage with wild type O and P genes results in an increase of bacterial DNA in the chromosomal region adjacent to the left of the prophage, that is a segment between gal and att (gal DNA) (Imae and Fukasawa, 1970). Evidence is presented in this report that such an increase of bacterial DNA is also seen in the region to the right of the prophage; a segment between bio and att (bio DNA). We postulate therefore that the bidirectional replication of DNA extends beyond the prophage and copies the neighboring host DNA until the prophage is excised. The model is verified by making use of excision-defective phages. The synthesis of gal DNA (or bio DNA) slows down to a halt within 40 min after the induction in the normal lysogens. The results are attributed to the prophage excision: (1) In lysogens for int, synthesis of the bacterial DNA continues for longer times. (2) The synthesis of the bacterial DNA slows down to a halt in lysogens for xis or b2 as in the control. However DNA synthesis also slows down in parallel so that the amount of the bacterial DNA relative to that of DNA synthesized by a given time stays constant from 20 min to 80 min. During that time the relative amount of the bacterial DNA rapidly decreases in the normal lysogen.The first article of this series is in J. molec. Biol. 54, 585 (1970).     

The effect of linkage on sample size determination for multiple trait selection

Sufficient sample sizes are needed in breeding programs to be confident, with a specified probability , of obtaining a specified number of plants of a desired genotype in segregating populations. We develop a method of determining the minimum sample size needed to produce, with specified probability , at least m individuals of a desired genotype. This method takes into consideration factors affecting differential selection of gametes, segregation at a single locus, and linkage among the loci of interest. We first consider the effects in the gametophyte (haploid level) of fitness and linkage on the frequencies of alleles at two linked loci, then at three or more linked loci. The probability of obtaining at least m successes, or occurrences of the desired allele, among n gametes is given by a formula based on the binomial distribution. This probability is affected by fitness and linkage through their impact on the probability that a single randomly chosen gamete is of the desired type. Using an extension of this approach, we examine the effects of the altered allelic frequencies on the likelihood of obtaining the desired genotype from a randomly chosen pair of gametes in the sporophyte (diploid level). A table and a figure show the sample size required to produce, with probability 0.95, m individuals of the desired g enotype or phenotype, as a function of m and the probability that a randomly selected individual is of the desired type.BU-1031-MC in the Technical Report Series of the Biometrics Unit, Cornell University, Ithaca, New York 14853     

Detection of segregation distortions in an indica-japonica rice cross using a high-resolution molecular map

We have constructed a high-resolution rice genetic map containing 1383 DNA markers covering 1575 cM on the 12 linkage groups of rice using 186 F₂ progeny from a cross between a japonica variety, Nipponbare, and an indica variety, Kasalath. Using this high-resolution molecular linkage map, we detected segregation distortion in a single wide cross of rice. The frequencies of genotypes for 1181 markers with more than 176 genotype data were plotted along this map to detect segregation distortion. Several types of distorted segregation were observed on 6 of the chromosomes. We could detect 11 major segregation distortions at ten positions on chromosomes 1, 3, 6, 8, 9, and 10. The strongest segregation distortion was at 107.2 cM on chromosome 3 and may be the gametophyte gene 2 (ga-2). The Kasalath genotype at this position was transmitted to the progeny with about a 95% probability through the pollen gamete. At least 8 out of the 11 segregation distortions detected here are new. The use of the high-resolution molecular linkage map for improving our understanding of the genetic nature and cause of these segregation distortions is discussed.     

Characterization of an unusual <Emphasis Type="Italic">Ds</Emphasis> transposable element in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>: Insertion of an abortive circular transposition intermediate

The maize Ac/Dstransposable elements, which belong to the hAT transposon superfamily, are widely used as insertional mutagens in numerous plant species. Molecular studies suggest that Ac/Ds elements transpose in a conservative non-replicative fashion; however the molecular mechanism of transposition remains unclear. We describe here the identification of an unusual Ds element, Ds-mmd1, in a transgenic Arabidopsis line. Ds-mmd1 is rearranged relative to the original Ds element, such that the original 5 and 3 ends are internal and previously internal sequences are the new 5 and 3 termini of Ds-mmd1. Short duplications of plant genomic DNA and Ds sequences are present at the Ds-mmd1 junctions, suggesting that a circular Dsmolecule was part of the events that created the Ds-mmd1 element. In addition, a revertant analysis on mmd1 plants demonstrated that Ds-mmd1 can be eliminated from the genome in an Ac-dependent process.