Evaluating the statistical significance of single band profiles in VNTR analyses.

VNTR profiles may present either a single band or two bands. If two bands are present then the individual is a heterozygote for these two bands. However, if only one band is present there is ambiguity as to the true genotype of the individual. This person may be a homozygote in that he has two copies of the same allele, or he may be a heterozygote for two very close bands that cannot be separated on the gel. The second NRC report proposed the use of the '2p' rule, or Formula 4.10a in the sub-structure case, as a conservative upper bound in the statistical interpretation. However, further examination suggests that these formulae are not necessarily conservative. In this paper we examine this phenomenon by deriving a formula that contains both the corrections for null alleles and for subpopulation effects.

Intensive care unit syndrome: a dangerous misnomer.

The terms intensive care unit (ICU) syndrome and ICU psychosis have been used interchangeably to describe a cluster of psychiatric symptoms that are unique to the ICU environment. It is often postulated that aspects of the ICU, such as sleep deprivation and sensory overload or monotony, are causes of the syndrome. This article reviews the empirical support for these propositions. We conclude that ICU syndrome does not differ from delirium and that ICU syndrome is caused exclusively by organic stressors on the central nervous system. We argue further that the term ICU syndrome is dangerous because it impedes standardized communication and research and may reduce the vigilance necessary to promptly investigate and reverse the medical cause of the delirium. Directions for future research are suggested.

Use of combined frequencies for RFLP and PCR based loci in determining match probability.

Statistical analysis was performed on a subset of the Pennsylvania State Police Caucasian, African American and Hispanic database for the purpose of determining Hardy-Weinberg equilibrium and associations across the RFLP loci D1S7, D2S44, D4S139, D5S110, D10S28 and D17S79 and the PCR-based loci HLA-DQA1, LDLR, GYPA, HBGG, D7S8 and Gc. Overall, the statistical results are consistent with a population in equilibrium both within and between loci. The assumption for independence is valid.

Estimating the genetic architecture of quantitative traits.

Understanding and estimating the structure and parameters associated with the genetic architecture of quantitative traits is a major research focus in quantitative genetics. With the availability of a well-saturated genetic map of molecular markers, it is possible to identify a major part of the structure of the genetic architecture of quantitative traits and to estimate the associated parameters. Multiple interval mapping, which was recently proposed for simultaneously mapping multiple quantitative trait loci (QTL), is well suited to the identification and estimation of the genetic architecture parameters, including the number, genomic positions, effects and interactions of significant QTL and their contribution to the genetic variance. With multiple traits and multiple environments involved in a QTL mapping experiment, pleiotropic effects and QTL by environment interactions can also be estimated. We review the method and discuss issues associated with multiple interval mapping, such as likelihood analysis, model selection, stopping rules and parameter estimation. The potential power and advantages of the method for mapping multiple QTL and estimating the genetic architecture are discussed. We also point out potential problems and difficulties in resolving the details of the genetic architecture as well as other areas that require further investigation. One application of the analysis is to improve genome-wide marker-assisted selection, particularly when the information about epistasis is used for selection with mating.

The exact test for cytonuclear disequilibria.

We extend the analysis of the statistical properties of cytonuclear disequilibria in two major ways. First, we develop the asymptotic sampling theory for the nonrandom associations between the alleles at a haploid cytoplasmic locus and the alleles and genotypes at a diploid nuclear locus, when there are an arbitrary number of alleles at each marker. This includes the derivation of the maximum likelihood estimators and their sampling variances for each disequilibrium measure, together with simple tests of the null hypothesis of no disequilibrium. In addition to these new asymptotic tests, we provide the first implementation of Fisher's exact test for the genotypic cytonuclear disequilibria and some approximations of the exact test. We also outline an exact test for allelic cytonuclear disequilibria in multiallelic systems. An exact test should be used for data sets when either the marginal frequencies are extreme or the sample size is small. The utility of this new sampling theory is illustrated through applications to recent nuclear-mtDNA and nuclear-cpDNA data sets. The results also apply to population surveys of nuclear loci in conjunction with markers in cytoplasmically inherited microorganisms.

Mitochondrial gene divergence of Colombian Drosophila pseudoobscura.

Isolated populations of drosophila pseudoobscura, separated from North American populations by about 2,400 km, were found in Colombia in 1960. We compared for sequences of the small ribosomal RNA (srRNA) gene on the mitochondria between North American and Colombian D. pseudoobscura in order to clarify the age of the Colombian isolates. The North American populations were not genetically different from each other but were genetically different from the Colombian populations. The Mexican strains represent the area from which the Colombian founders might have come. The estimated net nucleotide divergence between Mexican and Colombian D. pseudoobscura indicates that the Colombian population is not an ancient lineage. Phylogenies using both distance and parsimony methodologies reinforced this conclusion. The Colombian samples group together with both methods but, according to the bootstrap analysis, not significantly. It appears that the populations have not been separated long enough for their DNA sequences to show much divergence.

Alcohol consumption by undergraduate students.

This study examines a prediction of alcohol consumption in a sample of 212 undergraduate students at a university in Sydney. Data were collected using self-administered questionnaires, given to students in Semesters one and two. During both semesters the college residents reported drinking significantly more alcohol than students living elsewhere, but during vacation the intake of the two groups was approximately equal. Higher alcohol consumption in Semester two was best predicted by higher Semester one intake, followed by more consumption by friends and higher parental occupation status. Female students residing on campus were those most likely to be drinking at hazardous levels. Results also indicated that a significant proportion of residential students were drinking at hazardous levels, and that prevention interventions may need to focus on changing the attitudes and behavior of the social group.

Cytonuclear genetic structure of a hybrid zone in lizards of the Sceloporus grammicus complex (Sauria, Phrynosomatidae).

Lizards of the Sceloporus grammicus complex are comprised of multiple chromosome races that form several zones of parapatric hybridization in central Mexico. We scored diagnostic mitochondrial DNA (mtDNA) haplotypes and autosomal chromosome markers in a sample of 342 lizards from one well-defined zone between 2n = 34 and 2n = 46 races. A two-part analysis was performed on this data set in an attempt to infer the predominant evolutionary forces shaping the cytonuclear structure of this zone. The complications posed by its spatial structure were addressed by analysing a hierarchical series of smaller subsamples chosen to approximate single mating units. Two critical conclusions were drawn from this first-stage analysis. First and foremost, the three chromosomes have largely concordant cytonuclear disequilibrium patterns within each subsample with adequate numbers of individuals for detecting nonrandom cytonuclear associations. This suggests that the cytonuclear structure of this zone is predominantly a result of deterministic genome-wide forces rather than genetic drift of deterministic forces specific to individual chromosomes or loci. Second, the fit of a series of migration models to the data shows that the cytonuclear structure of the subsamples is well accounted for by continued gene flow from the two parental races alone, with random mating with respect to cytonuclear genotype and no other evolutionary forces. These results motivate several further empirical and theoretical investigations to refine our understanding of the relative roles of migration and other potentially important forces such as natural selection and genetic drift, in this and other hybrid zones.

Constraints and normalized measures for cytonuclear disequilibria.

The full bounds are derived for cytonuclear disequilibria in two-locus systems with an arbitrary number of alleles at the cytoplasmic and nuclear markers. The associated marginal frequencies constrain the nonrandom associations between cytoplasmic alleles and nuclear genotypes in the same way that the allele frequencies constrain the linkage disequilibrium between two nuclear loci. Additional constraints are imposed on the nonrandom associations between cytoplasmic and nuclear alleles, however, by the marginal frequencies of nuclear genotypes carrying either two or no copies of the associated nuclear allele. These bounds are analysed and used to define normalized measures of cytonuclear disequilibria, whose practical utility is illustrated through applications to two sets of recent nuclear-mitochondrial data.

Sampling theory for cytonuclear disequilibria.

We examine the statistical properties of cytonuclear disequilibria within a system including one diploid nuclear locus and one haploid cytoplasmic locus, each with two alleles. The results provide practical guidelines for the design and interpretation of cytonuclear surveys seeking to utilize the novel evolutionary information recorded in the observed pattern of cytonuclear associations. Important applications include population studies of nuclear allozymes in conjunction with genes from mitochondria, chloroplasts, or cytoplasmically inherited microorganisms. Our attention focuses on the allelic and genotypic disequilibria, which respectively measure the nonrandom associations between the cytotypes and the nuclear alleles and genotypes. We first derive the maximum likelihood estimators and their approximate large sample variances for each disequilibrium measure. These are each in turn used to set up an asymptotic test of the null hypothesis of no disequilibrium. We then calculate the minimum sample sizes required to detect the disequilibria under specified alternate hypotheses. The work also incorporates the deviation from Hardy-Weinberg equilibrium at the nuclear locus, which can significantly affect the results. The practical utility of this new sampling theory is illustrated through applications to two nuclear-mitochondrial data sets.

Estimation of mitotic stability in conidial fungi: a theoretical framework.

Mitotic stability refers to the probability that genetic elements are transmitted to both daughters during mitosis. This is of practical importance in molecular genetics because autonomous cloning vectors should be transmitted at high frequency during mitosis. In filamentous coencytic fungi it is difficult to quantify mitotic stability because a fluctuation test is not feasible. We show how to get around this problem by formulating a general model of the transmission of nuclear genetic elements through the course of conidiogenesis. We derive formulas by two different methods for the expected proportion of conidiospores that retain the element as a function of its mitotic stability and the number of generations of spore production. An important by-product yields the exact probability distributions for the number of conidiospores retaining elements at each stage of conidiophore development. We outline, and illustrate through specific numerical examples, how to use these formulas to estimate mitotic stability. Although we use Aspergillus nidulans as our biological paradigm, the same general framework can be extended to other fungal species, and possibly to less closely related systems as well.

Simulation study of a multigene family, with special reference to the evolution of compensatory advantageous mutations.

We investigate the evolution of a multigene family incorporating the forces of drift, mutation, gene conversion, unequal crossing over and selection. The use of simulation studies is required due to the complexity of the model. Selection is modeled in two modes: positive selection as a function of the number of different beneficial alleles and negative selection against deleterious alleles. We assume that gene conversion is unbiased, and that all mutations are initially deleterious. Compensation between mutants creates beneficial and neutral alleles, and allowances are made for compensatory mutations either within or between the members of a multigene family. We find that gene conversion can enhance the rate of acquisition of compensatory advantageous mutations when genes are redundant.

Gene conversion generates hypervariability at the variable regions of kallikreins and their inhibitors.

The two mechanisms for generating hypervariability at the reactive center of serine proteases and their inhibitors are gene conversion followed by natural selection and natural selection for point mutation. One way to clarify the effects of these two mechanisms is to calculate separately the number of nonsynonymous substitutions and that of synonymous substitutions at the variable regions and at the conserved regions. Our data analysis shows that not only the number of nonsynonymous substitutions but also the number of synonymous substitutions at the variable regions exceed the corresponding numbers at the conserved regions. Thus gene conversion has provided needed variability at the variable regions of serine proteases and their inhibitors. Natural selection has helped perpetuate such variability.

A branching-process model for the evolution of transposable elements incorporating selection.

We have formulated a very general mathematical model to analyze the evolution of transposable genetic elements in prokaryotic populations. Transposable genetic elements are DNA sequences able to replicate and insert copies of themselves at new locations in the genome. This work characterizes the equilibrium distribution of copy number under the influence of copy number-dependent selection, transposition and deletion. Our principal results concern the equilibrium distribution of copy number in response to various selective regimes. For particular transposition patterns (e.g., unregulated transposition or copy number-dependent transposition), equilibrium distributions are calculated numerically for a variety of specific selection patterns. Selection is quantified through specification of the expected number of offspring for individuals of each type, which is generally a non-increasing function of copy number, in accord with the usual evolutionary speculations.

Effect of gene conversion on variances of digenic identity measures.

The variances and covariances of digenic descent measures are studied for a two-locus model incorporating mutation, gene conversion, recombination, drift, and finite sampling. Gene conversion can occur between allelic pairs of genes or between non-allelic pairs on the same or different gametes within individuals. Most interest therefore centers on pairs of genes, and five digenic identity measures are required. The behavior over time of these measures is studied, with an emphasis on the effects of gene conversion. Because of the stochastic nature of the forces of drift, recombination, mutation, and conversion, the actual identity status of gene pairs can vary from expectation among replicate populations. To study this variation we compute the expected variances and covariances of the measures, and show that this requires the introduction of trigenic and quadrigenic measures. Allowing for conversion between genes on different gametes requires a large number of these higher-order measures.

Sampling strategies for distances between DNA sequences.

An international effort is now underway to obtain the DNA sequence for the entire human genome (Watson and Jordan, 1989, Genomics 5, 654-656; Barnhart, 1989, Genomics 5, 657-660). This Human Genome Initiative will generate sequence data from several species other than humans, and will result in several copies per species of at least some regions of the genome. Although the project has generated much interest, it is but one aspect of the widespread effort to generate DNA sequence data. Published sequences are collected in common databases, and release 63 of GenBank in March 1990 contained 40,127,752 bases from 33,337 reported sequences (News from GenBank 3; Mountain View, California: Intelligenetics, Inc., 1990). Large though this database is, it is only about 1% of the number of bases in the human genome. Interpretations of data of such magnitude are going to require the collaborative efforts of biometricians and molecular biologists, and an aim of this paper is to show that there is also a role for readers of this journal in the design of surveys of DNA sequences. Discussion here will center on the use of sequence data in evolutionary studies, where some region of DNA is sequenced in several different species. The object is to infer the evolutionary history of that particular region, or of the species themselves. Statistical issues in the very important studies on sequences to locate and characterize regions responsible for human diseases will not be addressed here. We will discuss appropriate ways of measuring distances between DNA sequences and of predicting the sampling properties of the distances. There are procedures for inferring evolutionary histories for a set of elements that depend on a matrix of distances between each pair of elements, and the precision of resulting trees must be influenced by the precision of the distances. We will show that account needs to be taken of two sampling processes--the sampling of sequences by the investigator ("statistical sampling"), and the sampling of genetic material involved in the formation of offspring from a parental population ("genetic sampling").

The evolution of latent genes in subdivided populations.

We define latent genes as phenotypically silent DNA sequences which may be reactivated by various genetic mechanisms. Of interest is how they and their functional counterparts can be maintained at high frequency in the face of mutation and selection pressure. We propose a two-deme, three-allele model incorporating viability selection, mutation and migration in haploid populations. It is shown that polymorphism for the three alleles can be easily maintained for a wide range of biologically meaningful parameter values. Computer simulations were employed to gain qualitative insight into the global dynamics of the system. It was found that the dynamics of the latent allele is closely correlated with that of the functional allele. In addition, bias in the migration rates can strengthen or weaken selective conditions for preservation of the functional and latent alleles.