RNA-based phylogenetic methods: application to mammalian mitochondrial RNA sequences

The PHASE software package allows phylogenetic tree construction with a number of evolutionary models designed specifically for use with RNA sequences that have conserved secondary structure. Evolution in the paired regions of RNAs occurs via compensatory substitutions, hence changes on either side of a pair are correlated. Accounting for this correlation is important for phylogenetic inference because it affects the likelihood calculation. In the present study we use the complete set of tRNA and rRNA sequences from 69 complete mammalian mitochondrial genomes. The likelihood calculation uses two evolutionary models simultaneously for different parts of the sequence: a paired-site model for the paired sites and a single-site model for the unpaired sites. We use Bayesian phylogenetic methods and a Markov chain Monte Carlo algorithm is used to obtain the most probable trees and posterior probabilities of clades. The results are well resolved for almost all the important branches on the mammalian tree. They support the arrangement of mammalian orders within the four supra-ordinal clades that have been identified by studies of much larger data sets mainly comprising nuclear genes. Groups such as the hedgehogs and the murid rodents, which have been problematic in previous studies with mitochondrial proteins, appear in their expected position with the other members of their order. Our choice of genes and evolutionary model appears to be more reliable and less subject to biases caused by variation in base composition than previous studies with mitochondrial genomes.     

The extent of linkage disequilibrium and haplotype sharing around a polymorphic site

Various expressions related to the length of a conserved haplotype around a polymorphism of known frequency are derived. We obtain exact expressions for the probability that no recombination has occurred in a sample or subsample. We obtain an approximation for the probability that no recombination that could give rise to a detectable recombination event (through the four-gamete test) has occurred. The probabilities can be used to obtain approximate distributions for the length of variously defined haplotypes around a polymorphic site. The implications of our results for data analysis, and in particular for detecting selection, are discussed.     

Tryparedoxins from Crithidia fasciculata and Trypanosoma brucei: photoreduction of the redox disulfide using synchrotron radiation and evidence for a conformational switch implicated in function

Tryparedoxin (TryX) is a member of the thioredoxin (TrX) fold family involved in the regulation of oxidative stress in parasitic trypanosomatids. Like TrX, TryX carries a characteristic Trp-Cys-Xaa-Xaa-Cys motif, which positions a redox-active disulfide underneath a tryptophan lid. We report the structure of a Crithidia fasciculata tryparedoxin isoform (CfTryX2) in two crystal forms and compare them with structures determined previously. Efforts to chemically generate crystals of reduced TryX1 were unsuccessful, and we carried out a novel experiment to break the redox-active disulfide, formed between Cys-40 and Cys-43, utilizing the intense x-radiation from a third generation synchrotron undulator beamline. A time course study of the S-S bond cleavage is reported with the structure of a TryX1 C43A mutant as the control. When freed from the constraints of a disulfide link to Cys-43, Cys-40 pivots to become slightly more solvent-accessible. In addition, we have determined the structure of Trypanosoma brucei TryX, which, influenced by the molecular packing in the crystal lattice, displays a significantly different orientation of the active site tryptophan lid. This structural change may be of functional significance when TryX interacts with tryparedoxin peroxidase, the final protein in the trypanothione-dependent peroxidase pathway. Comparisons with chloroplast TrX and its substrate fructose 1,6-bisphosphate phosphatase suggest that this movement may represent a general feature of redox regulation in the trypanothione and thioredoxin peroxidase pathways.     

Decorin binds fibrinogen in a Zn2+-dependent interaction

We have previously shown that decorin, a member of the small leucine-rich proteoglycan family of extracellular matrix proteoglycans/glycoproteins is a Zn(2+) metalloprotein at physiological Zn(2+) concentrations (Yang, V. W-C., LaBrenz, S. R., Rosenberg, L. C., McQuillan, D., and H??k, M. (1999) J. Biol. Chem. 274, 12454-12460). We now report that the decorin proteoglycan binds fibrinogen in the presence of Zn(2+). The fibrinogen-binding site is located in the N-terminal domain of the decorin core protein and a 45-amino acid peptide representing this domain binds to the fibrinogen D fragment with an apparent K(D) of 1.7 x 10(-6) m, as determined from fluorescence polarization data. Furthermore, we show that Zn(2+) promotes the self-association of decorin. The N-terminal domain of the core protein also mediates this activity. The results of solid-phase binding assays and gel filtration chromatography suggest that the N-terminal domain of decorin, when present at low micromolar concentrations, forms an oligomer in a Zn(2+)-dependent manner. Thus, Zn(2+) appears to play a pivotal role in the interactions and biological function of decorin.     

Human bronchoalveolar lavage: biofluid analysis with special emphasis on sample preparation

Respiratory diseases are an important health problem throughout the world. Whether caused by industrial pollutants, infections, smoking, cancer or metabolic diseases, damage to the lungs and airways often lead to morbidity or death. Bronchoalveolar lavage (BAL) obtained by fiber-optic bronchoscopy is a biofluid mirroring the expression of normally secreted pulmonary proteins and the products of activated cells and destructive processes. The characterization of the proteome within this compartment provides an opportunity to establish temporal and prognostic indicators of airway disease. The objective of this study was to develop methods of analysis of BAL samples, which achieved the highest level of annotation of the expression map of this proteome. We have optimized the process of sample preparation after investigating a variety of techniques including dialysis, ultramembrane filtration, precipitation and gel filtration. We have further studied methods to remove albumin from BAL in order to unmask proteins hidden on two-dimensional gels. In a pilot application of the method, BAL protein profiles obtained from healthy nonsmokers and smokers at risk for developing chronic obstructive pulmonary disease showed distinct differences.     

Haplotype block structure and its applications to association studies: power and study designs

Recent studies have shown that the human genome has a haplotype block structure, such that it can be divided into discrete blocks of limited haplotype diversity. In each block, a small fraction of single-nucleotide polymorphisms (SNPs), referred to as "tag SNPs," can be used to distinguish a large fraction of the haplotypes. These tag SNPs can potentially be extremely useful for association studies, in that it may not be necessary to genotype all SNPs; however, this depends on how much power is lost. Here we develop a simulation study to quantitatively assess the power loss for a variety of study designs, including case-control designs and case-parental control designs. First, a number of data sets containing case-parental or case-control samples are generated on the basis of a disease model. Second, a small fraction of case and control individuals in each data set are genotyped at all the loci, and a dynamic programming algorithm is used to determine the haplotype blocks and the tag SNPs based on the genotypes of the sampled individuals. Third, the statistical power of tests was evaluated on the basis of three kinds of data: (1) all of the SNPs and the corresponding haplotypes, (2) the tag SNPs and the corresponding haplotypes, and (3) the same number of randomly chosen SNPs as the number of tag SNPs and the corresponding haplotypes. We study the power of different association tests with a variety of disease models and block-partitioning criteria. Our study indicates that the genotyping efforts can be significantly reduced by the tag SNPs, without much loss of power. Depending on the specific haplotype block-partitioning algorithm and the disease model, when the identified tag SNPs are only 25% of all the SNPs, the power is reduced by only 4%, on average, compared with a power loss of approximately 12% when the same number of randomly chosen SNPs is used in a two-locus haplotype analysis. When the identified tag SNPs are approximately 14% of all the SNPs, the power is reduced by approximately 9%, compared with a power loss of approximately 21% when the same number of randomly chosen SNPs is used in a two-locus haplotype analysis. Our study also indicates that haplotype-based analysis can be much more powerful than marker-by-marker analysis.     

Statistical analysis of the influence of conspecifics on the dispersal of a soil collembola

The evidence for dispersal activity among soil-living invertebrates comes mainly from observations of their movement on artificial substrates or of colonisation of defaunated soils in the field. In an attempt to elucidate the dispersal pattern of soil collembolans in the presence of conspecifics, statistical analyses were undertaken to describe and simulate the movement of groups of Onychiurus armatus released in trays of homogeneous soil. A chi(2) test was used to reject the null hypothesis that individuals moved independently of each other and uniformly in all directions. The mean radial distance moved (1-2 cm day(-1)) and the radial standard deviation varied temporally and with the density of conspecifics. To capture the interaction between the moving individuals, four dispersal models (pure diffusion, diffusion with drift interaction, drift interaction and synchronised diffusion, and drift interaction and behavioural mood), were formulated as stochastic differential equations. The parameters of the models were estimated by minimising the deviance between the observed replicates and replicates that were simulated using the models. The dynamics of movement were best described by modelling the drift interaction as dependent on whether individuals were in a social or an asocial mood.     

Creating arrays by centrifugation

We describe afast, low-cost, and reliable way of creating arrays from sample molecules of interest present within microformatted sample vessels (such as 1536-well microplates). The principle involves simple centrifugal transfer of molecules of interest onto a solid planar or membrane surfaces placed over the initial sample vessel. Tools and procedures are presented that validate the robustness and precision of this facile solution to an otherwise difficult problem in modern molecular genetics. The availability of transferred DNA molecules for hybridization is also demonstrated. In conclusion, this "centrifugal-array" concept should help research studies to be applied on ever-greater scales with very simple machinery.