Improved statistical tests for differential gene expression by shrinking variance components estimates

Combining information across genes in the statistical analysis of microarray data is desirable because of the relatively small number of data points obtained for each individual gene. Here we develop an estimator of the error variance that can borrow information across genes using the James-Stein shrinkage concept. A new test statistic (FS) is constructed using this estimator. The new statistic is compared with other statistics used to test for differential expression: the gene-specific F test (F1), the pooled-variance F statistic (F3), a hybrid statistic (F2) that uses the average of the individual and pooled variances, the regularized t-statistic, the posterior odds statistic B, and the SAM t-test. The FS-test shows best or nearly best power for detecting differentially expressed genes over a wide range of simulated data in which the variance components associated with individual genes are either homogeneous or heterogeneous. Thus FS provides a powerful and robust approach to test differential expression of genes that utilizes information not available in individual gene testing approaches and does not suffer from biases of the pooled variance approach.     

DNA binding of a non-sequence-specific HMG-D protein is entropy driven with a substantial non-electrostatic contribution

The thermal properties of two forms of the Drosophila melanogaster HMG-D protein, with and without its highly basic 26 residue C-terminal tail (D100 and D74) and the thermodynamics of their non-sequence-specific interaction with linear DNA duplexes were studied using scanning and titration microcalorimetry, spectropolarimetry, fluorescence anisotropy and FRET techniques at different temperatures and salt concentrations. It was shown that the C-terminal tail of D100 is unfolded at all temperatures, whilst the state of the globular part depends on temperature in a rather complex way, being completely folded only at temperatures close to 0 degrees C and unfolding with significant heat absorption at temperatures below those of the gross denaturational changes. The association constant and thus Gibbs energy of binding for D100 is much greater than for D74 but the enthalpies of their association are similar and are large and positive, i.e. DNA binding is a completely entropy-driven process. The positive entropy of association is due to release of counterions and dehydration upon forming the protein/DNA complex. Ionic strength variation showed that electrostatic interactions play an important but not exclusive role in the DNA binding of the globular part of this non-sequence-specific protein, whilst binding of the positively charged C-terminal tail of D100 is almost completely electrostatic in origin. This interaction with the negative charges of the DNA phosphate groups significantly enhances the DNA bending. An important feature of the non-sequence-specific association of these HMG boxes with DNA is that the binding enthalpy is significantly more positive than for the sequence-specific association of the HMG box from Sox-5, despite the fact that these proteins bend the DNA duplex to a similar extent. This difference shows that the enthalpy of dehydration of apolar groups at the HMG-D/DNA interface is not fully compensated by the energy of van der Waals interactions between these groups, i.e. the packing density at the interface must be lower than for the sequence-specific Sox-5 HMG box.     

Genetics of colitis susceptibility in IL-10-deficient mice: backcross versus F2 results contrasted by principal component analysis

Interleukin-10-deficient (Il10(-/-)) mice on a C3H/HeJBir genetic background develop more severe colitis than those on a C57BL/6J background. We performed genome screens for quantitative trait loci (QTLs) regulating colitis susceptibility in this model system using two first backcross populations derived from these two strains. To reduce the complexity of this analysis, the information from numerous histologic phenotypes was summarized by principal component analysis. A similar approach was applied to previously published data from an F2 intercross (involving the same progenitor strains), which allowed us to ascertain all six previously reported cytokine-deficiency-induced colitis susceptibility loci (Cdcs1-6) with main and/or interacting effects on chromosomes 3, 1, 2, 8, 17, and 18. The colitogenic effect of Cdcs1 was confirmed in the backcross to C3H/HeJBir-Il10(-/-). Its effect was epistatically modified by another locus on chromosome 12. In addition, three main effect QTLs on chromosomes 4, 5, and 12 were identified in the backcross to C57BL/6J-Il10(-/-). Analyses of the modes of inheritance in these crosses revealed colitogenic contributions by both parental genomes. These findings show the complexity of inheritance underlying susceptibility to colitis and illustrate why detection of human inflammatory bowel disease loci has proven to be so difficult.     

The influence of short-term fasting on the quality of small bowel graft preservation

INTRODUCTION: Donor nutritional status may be a determinant of small bowel (SB) quality following storage. In this study, we investigated the effect of donor nutritional status and a proven nutrient-rich preservation solution on graft quality following cold storage. METHODS: Rats were fasted (12-14 h) or non-fasted. SB (n=6) was flushed vascularly with modified University of Wisconsin (UW) solution and flushed luminally with UW or an amino acid-rich (AA) solution as follows: Fasted. UWV, none; UWL, UW solution; AAL, AA solution. Non-fasted. UWV, none; UWL, UW solution; AAL, AA solution. Energetics, peroxidation (malondialdehyde; MDA), glutathione and histology were assessed over 24 h at 4 degrees C. RESULTS: Energetics (ATP, ATP/ADP, and energy charge) were significantly higher in AAL (fasted and non-fasted) groups than other groups. However, there were no differences in energetics parameters between fasted and non-fasted animals in all groups. MDA was higher in fasted groups than non-fasted tissues; interestingly, AAL values were up to 10-fold lower than other groups. Higher glutathione levels were detected in non-fasted AAL tissues. Mucosal integrity was markedly superior in luminally treated tissues (UWL and AAL) in fasted and non-fasted states. Most noteably, AAL tissues from fasted animals exhibited grade 2 injury (villus clefting), whereas normal mucosa was observed in non-fasted tissues (grade 0). CONCLUSION: Luminal flushing and a nutrient-rich preservation solution improve energetics, oxidative stress, and mucosal integrity during storage. Poorer donor nutritional status does not affect energetics throughout storage, but causes mucosal injury as a result of increased oxidative stress, even after a brief period of donor fasting.     

Combining data from multiple inbred line crosses improves the power and resolution of quantitative trait loci mapping

Rodent inbred line crosses are widely used to map genetic loci associated with complex traits. This approach has proven to be powerful for detecting quantitative trait loci (QTL); however, the resolution of QTL locations, typically approximately 20 cM, means that hundreds of genes are implicated as potential candidates. We describe analytical methods based on linear models to combine information available in two or more inbred line crosses. Our strategy is motivated by the hypothesis that common inbred strains of the laboratory mouse are derived from a limited ancestral gene pool and thus QTL detected in multiple crosses are likely to represent shared ancestral polymorphisms. We demonstrate that the combined-cross analysis can improve the power to detect weak QTL, can narrow support intervals for QTL regions, and can be used to separate multiple QTL that colocalize by chance. Moreover, combined-cross analysis can establish the allelic states of a QTL among a set of parental lines, thus providing critical information for narrowing QTL regions by haplotype analysis.     

Selective phenotyping for increased efficiency in genetic mapping studies

The power of a genetic mapping study depends on the heritability of the trait, the number of individuals included in the analysis, and the genetic dissimilarity among them. In experiments that involve microarrays or other complex physiological assays, phenotyping can be expensive and time-consuming and may impose limits on the sample size. A random selection of individuals may not provide sufficient power to detect linkage until a large sample size is reached. We present an algorithm for selecting a subset of individuals solely on the basis of genotype data that can achieve substantial improvements in sensitivity compared to a random sample of the same size. The selective phenotyping method involves preferentially selecting individuals to maximize their genotypic dissimilarity. Selective phenotyping is most effective when prior knowledge of genetic architecture allows us to focus on specific genetic regions. However, it can also provide modest improvements in efficiency when applied on a whole-genome basis. Importantly, selective phenotyping does not reduce the efficiency of mapping as compared to a random sample in regions that are not considered in the selection process. In contrast to selective genotyping, inferences based solely on a selectively phenotyped population of individuals are representative of the whole population. The substantial improvement introduced by selective phenotyping is particularly useful when phenotyping is difficult or costly and thus limits the sample size in a genetic mapping study.     

Calcium signalling by nicotinic acid adenine dinucleotide phosphate (NAADP)

Nicotinic acid adenine dinucleotide phosphate (NAADP) is a recently described Ca2+ mobilizing messenger, and probably the most potent. We briefly review its unique properties as a Ca2+ mobilizing agent. We present arguments for its action in targeting acidic calcium stores rather than the endoplasmic reticulum. Finally, we discuss possible biosynthetic pathways for NAADP in cells and candidates for its target Ca2+ release channel, which has eluded identification so far.     

DNA binding and bending by HMG boxes: energetic determinants of specificity

To clarify the physical basis of DNA binding specificity, the thermodynamic properties and DNA binding and bending abilities of the DNA binding domains (DBDs) of sequence-specific (SS) and non-sequence-specific (NSS) HMG box proteins were studied with various DNA recognition sequences using micro-calorimetric and optical methods. Temperature-induced unfolding of the free DBDs showed that their structure does not represent a single cooperative unit but is subdivided into two (in the case of NSS DBDs) or three (in the case of SS DBDs) sub-domains, which differ in stability. Both types of HMG box, most particularly SS, are partially unfolded even at room temperature but association with DNA results in stabilization and cooperation of all the sub-domains. Binding and bending measurements using fluorescence spectroscopy over a range of ionic strengths, combined with calorimetric data, allowed separation of the electrostatic and non-electrostatic components of the Gibbs energies of DNA binding, yielding their enthalpic and entropic terms and an estimate of their contributions to DNA binding and bending. In all cases electrostatic interactions dominate non-electrostatic in the association of a DBD with DNA. The main difference between SS and NSS complexes is that SS are formed with an enthalpy close to zero and a negative heat capacity effect, while NSS are formed with a very positive enthalpy and a positive heat capacity effect. This indicates that formation of SS HMG box-DNA complexes is specified by extensive van der Waals contacts between apolar groups, i.e. a more tightly packed interface forms than in NSS complexes. The other principal difference is that DNA bending by the NSS DBDs is driven almost entirely by the electrostatic component of the binding energy, while DNA bending by SS DBDs is driven mainly by the non-electrostatic component. The basic extensions of both categories of HMG box play a similar role in DNA binding and bending, making solely electrostatic interactions with the DNA.     

Structure of the Pseudomonas aeruginosa acyl-homoserinelactone synthase LasI

The LasI/LasR quorum-sensing system plays a pivotal role in virulence gene regulation of the opportunistic human pathogen, Pseudomonas aeruginosa. Here we report the crystal structure of the acyl-homoserine lactone (AHL) synthase LasI that produces 3-oxo-C12-AHL from the substrates 3-oxo-C12-acyl-carrier protein (acyl-ACP) and S-adenosyl-L-methionine. The LasI six-stranded beta sheet platform, buttressed by three alpha helices, forms a V-shaped substrate-binding cleft that leads to a tunnel passing through the enzyme that can accommodate the acyl-chain of acyl-ACP. This tunnel places no apparent restriction on acyl-chain length, in contrast to a restrictive hydrophobic pocket seen in the AHL-synthase EsaI. Interactions of essential conserved N-terminal residues, Arg23, Phe27 and Trp33, suggest that the N-terminus forms an enclosed substrate-binding pocket for S-adenosyl-L-methionine. Analysis of AHL-synthase surface residues identified a binding site for acyl-ACP, a role that was supported by in vivo reporter assay analysis of the mutated residues, including Arg154 and Lys150. This structure and the novel explanation of AHL-synthase acyl-chain-length selectivity promise to guide the design of Pseudomonas aeruginosa-specific quorum-sensing inhibitors as antibacterial agents.