Bias-corrected maximum likelihood estimator of the negative binomial dispersion parameter

We derive a first-order bias-corrected maximum likelihood estimator for the negative binomial dispersion parameter. This estimator is compared, in terms of bias and efficiency, with the maximum likelihood estimator investigated by Piegorsch (1990, Biometrics46, 863-867), the moment and the maximum extended quasi-likelihood estimators investigated by Clark and Perry (1989, Biometrics45, 309-316), and a double-extended quasi-likelihood estimator. The bias-corrected maximum likelihood estimator has superior bias and efficiency properties in most instances. For ease of comparison we give results for the two-parameter negative binomial model. However, an example involving negative binomial regression is given.     

fog-2 and the Evolution of Self-Fertile Hermaphroditism in Caenorhabditis

Somatic and germline sex determination pathways have diverged significantly in animals, making comparisons between taxa difficult. To overcome this difficulty, we compared the genes in the germline sex determination pathways of Caenorhabditis elegans and C. briggsae, two Caenorhabditis species with similar reproductive systems and sequenced genomes. We demonstrate that C. briggsae has orthologs of all known C. elegans sex determination genes with one exception: fog-2. Hermaphroditic nematodes are essentially females that produce sperm early in life, which they use for self fertilization. In C. elegans, this brief period of spermatogenesis requires FOG-2 and the RNA-binding protein GLD-1, which together repress translation of the tra-2 mRNA. FOG-2 is part of a large C. elegans FOG-2-related protein family defined by the presence of an F-box and Duf38/FOG-2 homogy domain. A fog-2-related gene family is also present in C. briggsae, however, the branch containing fog-2 appears to have arisen relatively recently in C. elegans, post-speciation. The C-terminus of FOG-2 is rapidly evolving, is required for GLD-1 interaction, and is likely critical for the role of FOG-2 in sex determination. In addition, C. briggsae gld-1 appears to play the opposite role in sex determination (promoting the female fate) while maintaining conserved roles in meiotic progression during oogenesis. Our data indicate that the regulation of the hermaphrodite germline sex determination pathway at the level of FOG-2/GLD-1/tra-2 mRNA is fundamentally different between C. elegans and C. briggsae, providing functional evidence in support of the independent evolution of self-fertile hermaphroditism. We speculate on the convergent evolution of hermaphroditism in Caenorhabditis based on the plasticity of the C. elegans germline sex determination cascade, in which multiple mutant paths yield self fertility.     

Life and death of lymphocytes: a role in immunesenescence

Human aging is associated with progressive decline in immune functions, increased frequency of infections. Among immune functions, a decline in T cell functions during aging predominates. In this review, we will discuss the molecular signaling in two major pathways of apoptosis, namely death receptor pathway and mitochondrial pathway, and their alterations in both T and B lymphocytes in human aging with a special emphasis on naïve and different memory subsets of CD8+ T cells. We will also discuss a possible role of lymphocyte apoptosis in immune senescence.     

Restrained molecular dynamics studies on complex of adriamycin with DNA hexamer sequence d-CGATCG

Adriamycin is an anthracycline anticancer drug used widely for solid tumors in spite of its adverse side effects. The solution structure of 2:1 adriamycin-d-(CGATCG)(2) complex has been studied by restrained molecular dynamics simulations. The restraint data set consists of several intramolecular and intermolecular nuclear Overhauser enhancement cross-peaks obtained from two-dimensional nuclear magnetic resonance spectroscopy data. The drug is found to intercalate between CG and GC base pairs at two d-CpG sites. The drug-DNA complex is stabilized via specific hydrogen bonding and van der Waal's interactions involving 4OCH(3), O5, 6OH, and NH(3)(+) moiety of daunosamine sugar, and rings A protons. The O-glycosidic bond C7-O7-C1'-C2' lies in the range 138 degrees -160 degrees during the course of simulations. The O6-H6...O5 hydrogen bond is stable while O11-H11...O12 hydrogen bond is not favored. The intercalating base pairs are buckled and minor groove is wider in the complex. The phosphate on one strand at intercalation site C1pG2 is in B(I) conformation and the phosphates directly lying on opposite strand is in B(II) conformation. The phosphorus on adjacent site G2pA3 is in B(II) conformation and hence a distinct pattern of B(I) and B(II) conformations is induced and stabilized. The role of various functional groups by which the molecular action is mediated has been discussed and correlated to the available biochemical evidence.     

Protein stabilization by introduction of cross-strand disulfides

Disulfides cross-link residues in a protein that are separated in primary sequence and stabilize the protein through entropic destabilization of the unfolded state. While the removal of naturally occurring disulfides leads to protein destabilization, introduction of engineered disulfides does not always lead to significant stabilization of a protein. We have analyzed naturally occurring disulfides that span adjacent antiparallel strands of beta sheets (cross-strand disulfides). Cross-strand disulfides have recently been implicated as redox-based conformational switches in proteins such as gp120 and CD4. The propensity of these disulfides to act as conformational switches was postulated on the basis of the hypothesis that this class of disulfide is conformationally strained. In the present analysis, there was no evidence to suggest that cross-strand disulfides are more strained compared to other disulfides as assessed by their torsional energy. It was also observed that these disulfides occur solely at non-hydrogen-bonded (NHB) registered pairs of adjacent antiparallel strands and not at hydrogen-bonded (HB) positions as suggested previously. One of the half-cystines involved in cross-strand disulfide formation often occurs at an edge strand. Experimental confirmation of the stabilizing effects of such disulfides was carried out in Escherichia coli thioredoxin. Four pairs of cross-strand cysteines were introduced, two at HB and two at NHB pairs. Disulfides were formed in all four cases. However, as predicted from our analysis, disulfides at NHB positions resulted in an increase in melting temperature of 7-10 degrees C, while at HB positions there was a corresponding decrease of -7 degrees C. The reduced state of all proteins had similar stability.     

A sequence-based genetic map of Medicago truncatula and comparison of marker colinearity with M. sativa

A core genetic map of the legume Medicago truncatula has been established by analyzing the segregation of 288 sequence-characterized genetic markers in an F(2) population composed of 93 individuals. These molecular markers correspond to 141 ESTs, 80 BAC end sequence tags, and 67 resistance gene analogs, covering 513 cM. In the case of EST-based markers we used an intron-targeted marker strategy with primers designed to anneal in conserved exon regions and to amplify across intron regions. Polymorphisms were significantly more frequent in intron vs. exon regions, thus providing an efficient mechanism to map transcribed genes. Genetic and cytogenetic analysis produced eight well-resolved linkage groups, which have been previously correlated with eight chromosomes by means of FISH with mapped BAC clones. We anticipated that mapping of conserved coding regions would have utility for comparative mapping among legumes; thus 60 of the EST-based primer pairs were designed to amplify orthologous sequences across a range of legume species. As an initial test of this strategy, we used primers designed against M. truncatula exon sequences to rapidly map genes in M. sativa. The resulting comparative map, which includes 68 bridging markers, indicates that the two Medicago genomes are highly similar and establishes the basis for a Medicago composite map.     

Gene expression intensity shapes evolutionary rates of the proteins encoded by the vertebrate genome

Natural selection leaves its footprints on protein-coding sequences by modulating their silent and replacement evolutionary rates. In highly expressed genes in invertebrates, these footprints are seen in the higher codon usage bias and lower synonymous divergence. In mammals, the highly expressed genes have a shorter gene length in the genome and the breadth of expression is known to constrain the rate of protein evolution. Here we have examined how the rates of evolution of proteins encoded by the vertebrate genomes are modulated by the amount (intensity) of gene expression. To understand how natural selection operates on proteins that appear to have arisen in earlier and later phases of animal evolution, we have contrasted patterns of mouse proteins that have homologs in invertebrate and protist genomes (Precambrian genes) with those that do not have such detectable homologs (vertebrate-specific genes). We find that the intensity of gene expression relates inversely to the rate of protein sequence evolution on a genomic scale. The most highly expressed genes actually show the lowest total number of substitutions per polypeptide, consistent with cumulative effects of purifying selection on individual amino acid replacements. Precambrian genes exhibit a more pronounced difference in protein evolutionary rates (up to three times) between the genes with high and low expression levels as compared to the vertebrate-specific genes, which appears to be due to the narrower breadth of expression of the vertebrate-specific genes. These results provide insights into the differential relationship and effect of the increasing complexity of animal body form on evolutionary rates of proteins.     

FtsZ collaborates with penicillin binding proteins to generate bacterial cell shape in Escherichia coli

The mechanisms by which bacteria adopt and maintain individual shapes remain enigmatic. Outstanding questions include why cells are a certain size, length, and width; why they are uniform or irregular; and why some branch while others do not. Previously, we showed that Escherichia coli mutants lacking multiple penicillin binding proteins (PBPs) display extensive morphological diversity. Because defective sites in these cells exhibit the structural and functional characteristics of improperly localized poles, we investigated the connection between cell division and shape. Here we show that under semipermissive conditions the temperature-sensitive FtsZ84 protein produces branched and aberrant cells at a high frequency in mutants lacking PBP 5, and this phenotype is exacerbated by the loss of additional peptidoglycan endopeptidases. Surprisingly, certain ftsZ84 strains lyse at the nonpermissive temperature instead of filamenting, and inhibition of wild-type FtsZ forces some mutants into tightly wound spirillum-like morphologies. The results demonstrate that significant aspects of bacterial shape are dictated by a previously unrecognized relationship between the septation machinery and ostensibly minor peptidoglycan-modifying enzymes and that under certain circumstances improper FtsZ function can destroy the structural integrity of the cell.