Adjusting for regression effect in uncontrolled studies

When clinical studies require enrolled patients to have abnormal assays, the natural tendency of repeat measurements to regress toward the mean can lead to a false assessment of effectiveness of therapy. We propose a method to more accurately estimate the true effect of therapy by adjusting for a component of improvement that can be attributed to regression effect. The model we use allows for a combination of additive and/or multiplicative effects of the therapy.     

Why are translationally sub-optimal synonymous codons used in Escherichia coli?

Natural selection favors certain synonymous codons which aid translation in Escherichia coli, yet codons not favored by translational selection persist. We use the frequency distributions of synonymous polymorphisms to test three hypotheses for the existence of translationally sub-optimal codons: (1) selection is a relatively weak force, so there is a balance between mutation, selection, and drift; (2) at some sites there is no selection on codon usage, so some synonymous sites are unaffected by translational selection; and (3) translationally sub-optimal codons are favored by alternative selection pressures at certain synonymous sites. We find that when all the data is considered, model 1 is supported and both models 2 and 3 are rejected as sole explanations for the existence of translationally sub-optimal codons. However, we find evidence in favor of both models 2 and 3 when the data is partitioned between groups of amino acids and between regions of the genes. Thus, all three mechanisms appear to contribute to the existence of translationally sub-optimal codons in E. coli. Received: 18 July 2000 / Accepted: 17 April 2001     

Mutability of microsatellites developed for the ant Camponotus consobrinus

Five highly polymorphic (GA)n microsatellite loci are reported for the formicine ant Camponotus consobrinus. The occurrence of many nests with a simple family structure enabled a search for new mutations, 11 of which were found from 3055 informative typings. These mutations were not randomly distributed across loci, 10 of them occurring at the locus Ccon70. The spectrum of mutations across alleles at Ccon70 was also nonrandom, with all of them occurring in alleles in the upper half of the allele size distribution. Six of the Ccon70 mutations decreased allele size. The mutations observed fit the stepwise mutation model well, i.e. mutations could always be assigned to an allele which differed in size from them by one repeat unit. The parental origins of the Ccon70 mutations were established and appear more female biased than vertebrate mutations, significantly so compared with human haemophilia A and primate intron mutations. This result may indicate that the lack of meiosis in males (which are haploid in ants) reduces the mutation rate in that sex relative to species in which both sexes are diploid.     

Generalized estimating equations for ordinal categorical data: arbitrary patterns of missing responses and missingness in a key covariate

We propose methods for regression analysis of repeatedly measured ordinal categorical data when there is nonmonotone missingness in these responses and when a key covariate is missing depending on observables. The methods use ordinal regression models in conjunction with generalized estimating equations (GEEs). We extend the GEE methodology to accommodate arbitrary patterns of missingness in the responses when this missingness is independent of the unobserved responses. We further extend the methodology to provide correction for possible bias when missingness in knowledge of a key covariate may depend on observables. The approach is illustrated with the analysis of data from a study in diagnostic oncology in which multiple correlated receiver operating characteristic curves are estimated and corrected for possible verification bias when the true disease status is missing depending on observables.     

Codon Bias and Mutability in HIV Sequences

A survey of the patterns of synonymous codon preference in the HIV env gene reveals a correlation between the codon bias and the mutability requirements of different regions of the protein. At hypervariable regions in gp120 one finds a greater proportion of codons that tend to mutate nonsynonymously, but to a target that is similar in hydrophobicity and volume. We argue that this strategy results from a compromise between the selective pressure placed on the virus by the induced immune response, which favors amino acid substitutions in the complementarity determining regions, and the negative selection against missense mutations that violate structural constraints of the env protein. Received: 9 June 1997 / Accepted: 25 May 1998     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998     

Compositional Bias May Affect Both DNA-Based and Protein-Based Phylogenetic Reconstructions

It is now well-established that compositional bias in DNA sequences can adversely affect phylogenetic analysis based on those sequences. Phylogenetic analyses based on protein sequences are generally considered to be more reliable than those derived from the corresponding DNA sequences because it is believed that the use of encoded protein sequences circumvents the problems caused by nucleotide compositional biases in the DNA sequences. There exists, however, a correlation between AT/GC bias at the nucleotide level and content of AT- and GC-rich codons and their corresponding amino acids. Consequently, protein sequences can also be affected secondarily by nucleotide compositional bias. Here, we report that DNA bias not only may affect phylogenetic analysis based on DNA sequences, but also drives a protein bias which may affect analyses based on protein sequences. We present a striking example where common phylogenetic tools fail to recover the correct tree from complete animal mitochondrial protein-coding sequences. The data set is very extensive, containing several thousand sites per sequence, and the incorrect phylogenetic trees are statistically very well supported. Additionally, neither the use of the LogDet/paralinear transform nor removal of positions in the protein alignment with AT- or GC-rich codons allowed recovery of the correct tree. Two taxa with a large compositional bias continually group together in these analyses, despite a lack of close biological relatedness. We conclude that even protein-based phylogenetic trees may be misleading, and we advise caution in phylogenetic reconstruction using protein sequences, especially those that are compositionally biased. Received: 19 February 1998 / Accepted: 28 August 1998     

Testing for age-at-onset anticipation with affected parent-child pairs

The tendency for the onset of a genetic disease to occur at progressively earlier ages or with progressively greater severity in successive generations is known as anticipation. Following the discovery of trinucleotide repeat expansion as a plausible genetic mechanism for anticipation, interest in testing for anticipation has increased. Studies of anticipation can be biased when parents with late onset or children with early onset are preferentially ascertained. This paper presents a nonparametric approach to testing for age-at-onset anticipation that adjusts for such preferential ascertainment. The approach is illustrated through application to data on panic disorder.     

Deletion bias in avian introns over evolutionary timescales

The role that introns play in the function and evolution of nuclear genomes is not fully understood. Recent models of intron evolution suggest that selection and drift may interact to maintain introns in multicellular organisms. In addition, deletion mutations are more likely to become fixed than insertion mutations. Examination of indel substitutions over macroevolutionary timescales in pigeons and doves (Aves: Columbiformes) reveals that deletion substitutions outnumber insertion substitutions by over six times. The length of indel events is variable.     

Evolutionary pressures on apicoplast transit peptides

Malaria parasites (species of the genus Plasmodium) harbor a relict chloroplast (the apicoplast) that is the target of novel antimalarials. Numerous nuclear-encoded proteins are translocated into the apicoplast courtesy of a bipartite N-terminal extension. The first component of the bipartite leader resembles a standard signal peptide present at the N-terminus of secreted proteins that enter the endomembrane system. Analysis of the second portion of the bipartite leaders of P. falciparum, the so-called transit peptide, indicates similarities to plant transit peptides, although the amino acid composition of P. falciparum transit peptides shows a strong bias, which we rationalize by the extraordinarily high AT content of P. falciparum DNA. 786 plastid transit peptides were also examined from several other apicomplexan parasites, as well as from angiosperm plants. In each case, amino acid biases were correlated with nucleotide AT content. A comparison of a spectrum of organisms containing primary and secondary plastids also revealed features unique to secondary plastid transit peptides. These unusual features are explained in the context of secondary plastid trafficking via the endomembrane system.