A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Rapid parallel expression in E-coli and insect cells: Analysis of five lef gene products of the Autographa californica multiple nuclear polyhedrosis virus (AcMNPV)

A number of strategies are emerging for the high throughput (HTP) expression of recombinant proteins to enable structural and functional study. Here we describe a workable HTP strategy based on parallel protein expression in E. coli and insect cells. Using this system we provide comparative expression data for five proteins derived from the Autographa californica polyhedrosis virus genome that vary in amino acid composition and in molecular weight. Although the proteins are part of a set of factors known to be required for viral late gene expression, the precise function of three of the five, late expression factors (lefs) 6, 7 and 10, is unknown. Rapid expression and characterisation has allowed the determination of their ability to bind DNA and shown a cellular location consistent with their properties. Our data point to the utility of a parallel expression strategy to rapidly obtain workable protein expression levels from many open reading frames (ORFs).

A Bayesian approach to disease gene location using allelic association

A Bayesian approach to analysing data from family-based association studies is developed. This permits direct assessment of the range of possible values of model parameters, such as the recombination frequency and allelic associations, in the light of the data. In addition, sophisticated comparisons of different models may be handled easily, even when such models are not nested. The methodology is developed in such a way as to allow separate inferences to be made about linkage and association by including theta, the recombination fraction between the marker and disease susceptibility locus under study, explicitly in the model. The method is illustrated by application to a previously published data set. The data analysis raises some interesting issues, notably with regard to the weight of evidence necessary to convince us of linkage between a candidate locus and disease.

Characterisation of an amphioxus Fringe gene and the evolution of the vertebrate segmentation clock

In mouse and chick embryos, cyclic expression of lunatic fringe has an important role in the regulation of mesoderm segmentation. We have isolated a Fringe gene from the protochordate amphioxus. Amphioxus is the closest living relative of the vertebrates, and has mesoderm that is definitively segmented in a manner that is similar to, and probably homologous with, that of vertebrates. AmphiFringe is placed basal to vertebrate Fringe genes in molecular phylogenetic analyses, indicating that the duplications that formed radical-, manic- and lunatic fringe are specific to the vertebrate lineage. AmphiFringe expression was detected in the anterior neural plate of early neurulae, where it resolved into a series of segmental patches by the mid-neurulae stage. No AmphiFringe transcripts were detected in the mesoderm. Based on these observations, we propose a model depicting a successive recruitment of Fringe in the maintenance then regulation of segmentation during vertebrate evolution.

Global iron-dependent gene regulation in Escherichia coli - A new mechanism for iron homeostasis

Organisms generally respond to iron deficiency by increasing their capacity to take up iron and by consuming intracellular iron stores. Escherichia coli, in which iron metabolism is particularly well understood, contains at least 7 iron-acquisition systems encoded by 35 iron-repressed genes. This Fe-dependent repression is mediated by a transcriptional repressor, Fur ( ferric uptake regulation), which also controls genes involved in other processes such as iron storage, the Tricarboxylic Acid Cycle, pathogenicity, and redox-stress resistance. Our macroarray-based global analysis of iron- and Fur-dependent gene expression in E. coli has revealed several novel Fur-repressed genes likely to specify at least three additional iron- transport pathways. Interestingly, a large group of energy metabolism genes was found to be iron and Fur induced. Many of these genes encode iron- rich respiratory complexes. This iron- and Fur-dependent regulation appears to represent a novel iron-homeostatic mechanism whereby the synthesis of many iron- containing proteins is repressed under iron- restricted conditions. This mechanism thus accounts for the low iron contents of fur mutants and explains how E. coli can modulate its iron requirements. Analysis of Fe-55-labeled E. coli proteins revealed a marked decrease in iron- protein composition for the fur mutant, and visible and EPR spectroscopy showed major reductions in cytochrome b and d levels, and in iron- sulfur cluster contents for the chelator-treated wild-type and/or fur mutant, correlating well with the array and quantitative RT-PCR data. In combination, the results provide compelling evidence for the regulation of intracellular iron consumption by the Fe2+-Fur complex.

Effects of chronic and acute consumption of fruit- and vegetable-puree-based drinks on vasodilation, risk factors for CVD and the response as a result of the eNOS G298T polymorphism

The average UK adult consumes less than three portions of fruit and vegetables daily, despite evidence to suggest that consuming five portions daily could help prevent chronic diseases. It is recommended that fruit juice should only count as one of these portions, as juicing removes fibre and releases sugars. However, fruit juices contain beneficial compounds such as vitamin C and flavonoids and could be a useful source of dietary phytochemicals. Two randomised controlled cross-over intervention studies investigating the effects of chronic and acute consumption of commercially-available fruit- and vegetable-puree-based drinks (FVPD) on bioavailability, antioxidant status and CVD risk factors are described. Blood and urine samples were collected during both studies and vascular tone was measured using laser Doppler imaging. In the chronic intervention study FVPD consumption was found to significantly increase dietary carotenoids (P = 0.001) and vitamin C (P = 0.003). Plasma carotenoids were increased (P = 0.001), but the increase in plasma vitamin C was not significant. There were no significant effects on oxidative stress, antioxidant status and other CVD risk factors. In the acute intervention study FVPD were found to increase total plasma nitrate and nitrite (P = 0.001) and plasma vitamin C (P = 0.002). There was no effect on plasma lipids or uric acid, but there was a lower glucose and insulin peak concentration after consumption of the FVPD compared with the sugar-matched control. There was a trend towards increased vasodilation following both chronic and acute FVPD consumption. All volunteers were retrospectively genotyped for the eNOS G298T polymorphism and the effect of genotype on the measurements is discussed. Overall, there was a non-significant trend towards increased endothelium-dependent vasodilation following both acute and chronic FVPD consumption. However, there was a significant time x treatment effect (P < 0.05) of acute FVPD consumption in individuals with the GG variant of the eNOS gene.

Assays for enhanced activity of low efficacy partial agonists at the D-2 dopamine receptor

Background and purpose: Low efficacy partial agonists at the D-2 dopamine receptor may be useful for treating schizophrenia. In this report we describe a method for assessing the efficacy of these compounds based on stimulation of [S-35]GTP gamma S binding. Experimental approach: Agonist efficacy was assessed from [S-35]GTP gamma S binding to membranes of CHO cells expressing D2 dopamine receptors in buffers with and without Na+. Effects of Na+ on receptor/G protein coupling were assessed using agonist/[H-3] spiperone competition binding assays. Key results: When [S-35]GTP gamma S binding assays were performed in buffers containing Na+, some agonists (aripiprazole, AJ-76, UH-232) exhibited very low efficacy whereas other agonists exhibited measurable efficacy. When Na+ was substituted by N-methyl D-glucamine, the efficacy of all agonists increased (relative to that of dopamine) but particularly for aripiprazole, aplindore, AJ-76, (-)-3-PPP and UH-232. In ligand binding assays, substitution of Na+ by N-methyl D-glucamine increased receptor/G protein coupling for some agonists -. aplindore, dopamine and (-)-3-PPP-but for aripiprazole, AJ-76 and UH-232 there was little effect on receptor/G protein coupling. Conclusions and implications: Substitution of Na+ by NMDG increases sensitivity in [S-35] GTPgS binding assays so that very low efficacy agonists were detected clearly. For some agonists the effect seems to be mediated via enhanced receptor/G protein coupling whereas for others the effect is mediated at another point in the G protein activation cycle. AJ-76, aripiprazole and UH-232 seem particularly sensitive to this change in assay conditions. This work provides a new method to discover these very low efficacy agonists.

Profiling dehydrin gene sequence and physiological parameters in drought tolerant and susceptible spring wheat cultivars.

Physiological and yield traits such as stomatal conductance (mmol m-2s-1), Leaf relative water content (RWC %) and grain yield per plant were studied in a separate experiment. Results revealed that five out of sixteen cultivars viz. Anmol, Moomal, Sarsabz, Bhitai and Pavan, appeared to be relatively more drought tolerant. Based on morphophysiological results, studies were continued to look at these cultivars for drought tolerance at molecular level. Initially, four well recognized primers for dehydrin genes (DHNs) responsible for drought induction in T. durum L., T. aestivum L. and O. sativa L. were used for profiling gene sequence of sixteen wheat cultivars. The primers amplified the DHN genes variably like Primer WDHN13 (T. aestivum L.) amplified the DHN gene in only seven cultivars whereas primer TdDHN15 (T. durum L.) amplified all the sixteen cultivars with even different DNA banding patterns some showing second weaker DNA bands. Third primer TdDHN16 (T. durum L.) has shown entirely different PCR amplification prototype, specially showing two strong DNA bands while fourth primer RAB16C (O. sativa L.) failed to amplify DHN gene in any of the cultivars. Examination of DNA sequences revealed several interesting features. First, it identified the two exon/one intron structure of this gene (complete sequences were not shown), a feature not previously described in the two database cDNA sequences available from T. aestivum L. (gi|21850). Secondly, the analysis identified several single nucleotide polymorphisms (SNPs), positions in gene sequence. Although complete gene sequence was not obtained for all the cultivars, yet there were a total of 38 variable positions in exonic (coding region) sequence, from a total gene length of 453 nucleotides. Matrix of SNP shows these 37 positions with individual sequence at positions given for each of the 14 cultivars (sequence of two cultivars was not obtained) included in this analysis. It demonstrated a considerable diversity for this gene with only three cultivars i.e. TJ-83, Marvi and TD-1 being similar to the consensus sequence. All other cultivars showed a unique combination of SNPs. In order to prove a functional link between these polymorphisms and drought tolerance in wheat, it would be necessary to conduct a more detailed study involving directed mutation of this gene and DHN gene expression.

Development and validation of a qPCR-based method for quantifying Shiga toxin-encoding and other lambdoid bacteriophages

P>To address whether seasonal variability exists among Shiga toxin-encoding bacteriophage (Stx phage) numbers on a cattle farm, conventional plaque assay was performed on water samples collected over a 17 month period. Distinct seasonal variation in bacteriophage numbers was evident, peaking between June and August. Removal of cattle from the pasture precipitated a reduction in bacteriophage numbers, and during the winter months, no bacteriophage infecting Escherichia coli were detected, a surprising occurrence considering that 1031 tailed-bacteriophages are estimated to populate the globe. To address this discrepancy a culture-independent method based on quantitative PCR was developed. Primers targeting the Q gene and stx genes were designed that accurately and discriminately quantified artificial mixed lambdoid bacteriophage populations. Application of these primer sets to water samples possessing no detectable phages by plaque assay, demonstrated that the number of lambdoid bacteriophage ranged from 4.7 x 104 to 6.5 x 106 ml-1, with one in 103 free lambdoid bacteriophages carrying a Shiga toxin operon (stx). Specific molecular biological tools and discriminatory gene targets have enabled virus populations in the natural environment to be enumerated and similar strategies could replace existing propagation-dependent techniques, which grossly underestimate the abundance of viral entities.

Development of a miniaturised microarray-based assay for the rapid identification of antimicrobial resistance genes in Gram-negative bacteria

We describe the development of a miniaturised microarray for the detection of antimicrobial resistance genes in Gram-negative bacteria. Included on the array are genes encoding resistance to aminoglycosides, trimethoprim, sulphonamides, tetracyclines and beta-lactams, including extended-spectrum beta-lactamases. Validation of the array with control strains demonstrated a 99% correlation between polymerase chain reaction and array results. There was also good correlation between phenotypic and genotypic results for a large panel of Escherichia coli and Salmonella isolates. Some differences were also seen in the number and type of resistance genes harboured by E. coli and Salmonella strains. The array provides an effective, fast and simple method for detection of resistance genes in clinical isolates suitable for use in diagnostic laboratories, which in future will help to understand the epidemiology of isolates and to detect gene linkage in bacterial populations. (C) 2008 Published by Elsevier B.V. and the International Society of Chemotherapy.

Evaluation and comparison of SYBR Green I Real-Time PCR and TaqMan Real-Time PCR methods for quantitative assay of Listeria monocytogenes in nutrient broth and milk

Specific traditional plate count method and real-time PCR systems based on SYBR Green I and TaqMan technologies using a specific primer pair and probe for amplification of iap-gene were used for quantitative assay of Listeria monocytogenes in seven decimal serial dilution series of nutrient broth and milk samples containing 1.58 to 1.58×107 cfu /ml and the real-time PCR methods were compared with the plate count method with respect to accuracy and sensitivity. In this study, the plate count method was performed using surface-plating of 0.1 ml of each sample on Palcam Agar. The lowest detectable level for this method was 1.58×10 cfu/ml for both nutrient broth and milk samples. Using purified DNA as a template for generation of standard curves, as few as four copies of the iap-gene could be detected per reaction with both real-time PCR assays, indicating that they were highly sensitive. When these real-time PCR assays were applied to quantification of L. monocytogenes in decimal serial dilution series of nutrient broth and milk samples, 3.16×10 to 3.16×105 copies per reaction (equals to 1.58×103 to 1.58×107 cfu/ml L. monocytogenes) were detectable. As logarithmic cycles, for Plate Count and both molecular assays, the quantitative results of the detectable steps were similar to the inoculation levels.

Partial identification of spirochaetes from two dairy cows with digital dermatitis by polymerase chain reaction analysis of the 16S ribosomal RNA gene

Specimens taken postmortem from typical lesions of digital dermatitis in two dairy cows were tested by the polymerase chain reaction (PCR) for the presence of a spirochaetal 16S rRNA gene. Seven different assays detected the gene in the samples from both cows. Two of the PCR products were sequenced and a comparison of the nucleotide sequences revealed that the spirochaete belonged to the genus Treponema and was closely related to Treponema denticola. A PCR specific for the detection of the digital dermatitis-associated treponeme was developed.

Acetone enhances the direct analysis of procyanidin- and prodelphinidin-based condensed tannins in lotus species by the butanol-HCl-iron assay

The butanol-HCl spectrophotometric assay is widely used for quantifying extractable and insoluble condensed tannins (CT, syn. proanthocyanidins) in foods, feeds, and foliage of herbaceous and woody plants, but the method underestimates total CT content when applied directly to plant material. To improve CT quantitation, we tested various cosolvents with butanol-HCl and found that acetone increased anthocyanidin yields from two forage Lotus species having contrasting procyanidin and prodelphinidin compositions. A butanol-HCl-iron assay run with 50% (v/v) acetone gave linear responses with Lotus CT standards and increased estimates of total CT in Lotus herbage and leaves by up to 3.2-fold over the conventional method run without acetone. The use of thiolysis to determine the purity of CT standards further improved quantitation. Gel-state 13C and 1H–13C HSQC NMR spectra of insoluble residues collected after butanol-HCl assays revealed that acetone increased anthocyanidin yields by facilitating complete solubilization of CT from tissue.

Construction of a Festuca pratensis BAC library for map-based cloning in Festulolium substitution lines

Introgression in Festulolium is a potentially powerful tool to isolate genes for a large number of traits which differ between Festuca pratensis Huds. and Lolium perenne L. Not only are hybrids between the two species fertile, but the two genomes can be distinguished by genomic in situ hybridisation and a high frequency of recombination occurs between homoeologous chromosomes and chromosome segments. By a programme of introgression and a series of backcrosses, L. perenne lines have been produced which contain small F. pratensis substitutions. This material is a rich source of polymorphic markers targeted towards any trait carried on the F. pratensis substitution not observed in the L. perenne background. We describe here the construction of an F. pratensis BAC library, which establishes the basis of a map-based cloning strategy in L. perenne. The library contains 49,152 clones, with an average insert size of 112 kbp, providing coverage of 2.5 haploid genome equivalents. We have screened the library for eight amplified fragment length polymorphism (AFLP) derived markers known to be linked to an F. pratensis gene introgressed into L. perenne and conferring a staygreen phenotype as a consequence of a mutation in primary chlorophyll catabolism. While for four of the markers it was possible to identify bacterial artificial chromosome (BAC) clones, the other four AFLPs were too repetitive to enable reliable identification of locus-specific BACs. Moreover, when the four BACs were partially sequenced, no obvious coding regions could be identified. This contrasted to BACs identified using cDNA sequences, when multiple genes were identified on the same BAC.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.