Sample size for detecting differentially expressed genes in microarray experiments

Background  

Microarray experiments are often performed with a small number of biological replicates, resulting in low statistical power for detecting differentially expressed genes and concomitant high false positive rates. While increasing sample size can increase statistical power and decrease error rates, with too many samples, valuable resources are not used efficiently. The issue of how many replicates are required in a typical experimental system needs to be addressed. Of particular interest is the difference in required sample sizes for similar experiments in inbred vs. outbred populations (e.g. mouse and rat vs. human).     

DETECTING DESERTIFICATION PROCESSES USING TM AND ETM+ DATA, NORTH OF ISFAHAN, IRAN

 Aims Desertification results in ecological and biological diminution of the earth, and can happen naturally or cause by anthropogenic activities. This process especially affects arid and semi-arid regions, such as the Isfahan region, where the spread of desertification is reaching critical proportions. The aim of this study is to use remotely sensed data to review the trend of desertification in the northern of Isfahan, Iran. Methods Multi-temporal images were employed to evaluate the trend of desertification, specifically the TM and ETM+ data of September, 1990 and September, 2001. Geometric and radiometric corrections were applied to each image prior to image processing and supervised classification, and vegetation indices were applied to produce a land use map of each image in nine classes. The land use classification s in the two map images were compared and changes between land use classes were detected over the 11 year period using a fuzzy and post-classification technique. Important findings The maps and their comparison with false color composite images showed the differences efficiently. With the fuzzy and post-classification method the land use changes were sited on the map. Fuzzy confirmed 53% changed area and 47% unchanged areas in the study region. The results verify the desertification expansion in the study areas. Because of poor land management, agricultural lands converted to desert and abandoned areas, and some marginal pasture lands had to be changed to agricultural land which are desertification spreading according to United Nations Conference on Desertification (UNCOD). Also farmland and pastures have been converted to urban and industrial areas, and the rangelands have been spoiled due to opencast mine excavations. With the mine margins eroding as well as their debris accumulating on the pasture lands, desertification has become worse. Three areas of less-elevated mountains have remained unchanged. This study confirmed that the anthropogenic activities accelerated the desertification process and severely endangered the remaining areas.     

Bayesian mixture model based clustering of replicated microarray data

MOTIVATION: Identifying patterns of co-expression in microarray data by cluster analysis has been a productive approach to uncovering molecular mechanisms underlying biological processes under investigation. Using experimental replicates can generally improve the precision of the cluster analysis by reducing the experimental variability of measurements. In such situations, Bayesian mixtures allow for an efficient use of information by precisely modeling between-replicates variability. RESULTS: We developed different variants of Bayesian mixture based clustering procedures for clustering gene expression data with experimental replicates. In this approach, the statistical distribution of microarray data is described by a Bayesian mixture model. Clusters of co-expressed genes are created from the posterior distribution of clusterings, which is estimated by a Gibbs sampler. We define infinite and finite Bayesian mixture models with different between-replicates variance structures and investigate their utility by analyzing synthetic and the real-world datasets. Results of our analyses demonstrate that (1) improvements in precision achieved by performing only two experimental replicates can be dramatic when the between-replicates variability is high, (2) precise modeling of intra-gene variability is important for accurate identification of co-expressed genes and (3) the infinite mixture model with the 'elliptical' between-replicates variance structure performed overall better than any other method tested. We also introduce a heuristic modification to the Gibbs sampler based on the 'reverse annealing' principle. This modification effectively overcomes the tendency of the Gibbs sampler to converge to different modes of the posterior distribution when started from different initial positions. Finally, we demonstrate that the Bayesian infinite mixture model with 'elliptical' variance structure is capable of identifying the underlying structure of the data without knowing the 'correct' number of clusters. AVAILABILITY: The MS Windows based program named Gaussian Infinite Mixture Modeling (GIMM) implementing the Gibbs sampler and corresponding C++ code are available at http://homepages.uc.edu/~medvedm/GIMM.htm SUPPLEMENTAL INFORMATION: http://expression.microslu.washington.edu/expression/kayee/medvedovic2003/medvedovic_bioinf2003.html     

Studies on Calf Diarrhoea in Mozambique: Prevalence of Bacterial Pathogens

The prevalence of diarrhoea in calves was investigated in 8 dairy farms in Mozambique at 4 occasions during 2 consecutive years. A total of 1241 calves up to 6 months of age were reared in the farms, and 63 (5%) of them had signs of diarrhoea. Two farms had an overall higher prevalence (13% and 21%) of diarrhoea. Faecal samples were collected from all diarrhoeal calves (n = 63) and from 330 healthy calves and analysed for Salmonella species, Campylobacter jejuni and enterotoxigenic Escherichia coli (ETEC). Salmonella spp. was isolated in only 2% of all calves. Campylobacter was isolated in 11% of all calves, irrespective of health condition, and was more frequent (25%) in one of the 2 diarrhoeal farms (p = 0.001). 80% of the isolates were identified as C. jejuni. No ETEC strains were detected among the 55 tested strains from diarrhoeal calves, but 22/55 (40%) strains from diarrhoeal calves and 14/88 (16%) strains from healthy calves carried the K99 adhesin (p = 0.001). 6,757 E. coli isolates were typed with a biochemical fingerprinting method (the PhenePlate?) giving the same E. coli diversity in healthy and diarrhoeal calves. Thus it was concluded: i) the overall prevalence of diarrhoea was low, but 2 farms had a higher prevalence that could be due to an outbreak situation, ii) Salmonella did not seem to be associated with diarrhoea, iii) Campylobacter jejuni was common at one of the 2 diarrhoeal farms and iv) ETEC strains were not found, but K99 antigen was more prevalent in E. coli strains from diarrhoeal calves than from healthy, as well as more prevalent in one diarrhoeal farm.     

Advanced age promotes colonic dysfunction and gut‐derived lung infection after stroke

Bacterial infection a leading cause of death among patients with stroke, with elderly patients often presenting with more debilitating outcomes. The findings from our retrospective study, supported by previous clinical reports, showed that increasing age is an early predictor for developing fatal infectious complications after stroke. However, exactly how and why older individuals are more susceptible to infection after stroke remains unclear. Using a mouse model of transient ischaemic stroke, we demonstrate that older mice (>12 months) present with greater spontaneous bacterial lung infections compared to their younger counterparts (7–10 weeks) after stroke. Importantly, we provide evidence that older poststroke mice exhibited elevated intestinal inflammation and disruption in gut barriers critical in maintaining colonic integrity following stroke, including reduced expression of mucin and tight junction proteins. In addition, our data support the notion that the localized pro‐inflammatory microenvironment driven by increased tumour necrosis factor‐α production in the colon of older mice facilitates the translocation and dissemination of orally inoculated bacteria to the lung following stroke onset. Therefore, findings of this study demonstrate that exacerbated dysfunction of the intestinal barrier in advanced age promotes translocation of gut‐derived bacteria and contributes to the increased risk to poststroke bacterial infection.     

A preliminary investigation into the relationship between functional movement screen scores and athletic physical performance in female team sport athletes

There is little research investigating relationships between the Functional Movement Screen (FMS) and athletic performance in female athletes. This study analyzed the relationships between FMS (deep squat; hurdle step [HS]; in-line lunge [ILL]; shoulder mobility; active straight-leg raise [ASLR]; trunk stability push-up; rotary stability) scores, and performance tests (bilateral and unilateral sit-and-reach [flexibility]; 20-m sprint [linear speed]; 505 with turns from each leg; modified T-test with movement to left and right [change-of-direction speed]; bilateral and unilateral vertical and standing broad jumps; lateral jumps [leg power]). Nine healthy female recreational team sport athletes (age = 22.67 ± 5.12 years; height = 1.66 ± 0.05 m; body mass = 64.22 ± 4.44 kilograms) were screened in the FMS and completed the afore-mentioned tests. Percentage between-leg differences in unilateral sit-and-reach, 505 turns and the jumps, and difference between the T-test conditions, were also calculated. Spearman''s correlations (p ≤ 0.05) examined relationships between the FMS and performance tests. Stepwise multiple regressions (p ≤ 0.05) were conducted for the performance tests to determine FMS predictors. Unilateral sit-and-reach positive correlated with the left-leg ASLR (r = 0.704-0.725). However, higher-scoring HS, ILL, and ASLR related to poorer 505 and T-test performance (r = 0.722-0.829). A higher-scored left-leg ASLR related to a poorer unilateral vertical and standing broad jump, which were the only significant relationships for jump performance. Predictive data tended to confirm the correlations. The results suggest limitations in using the FMS to identify movement deficiencies that could negatively impact athletic performance in female team sport athletes.     

Comparison of Major and Minor Viral SNPs Identified through Single Template Sequencing and Pyrosequencing in Acute HIV-1 Infection

Massively parallel sequencing (MPS) technologies, such as 454-pyrosequencing, allow for the identification of variants in sequence populations at lower levels than consensus sequencing and most single-template Sanger sequencing experiments. We sought to determine if the greater depth of population sampling attainable using MPS technology would allow detection of minor variants in HIV founder virus populations very early in infection in instances where Sanger sequencing detects only a single variant. We compared single nucleotide polymorphisms (SNPs) during acute HIV-1 infection from 32 subjects using both single template Sanger and 454-pyrosequencing. Pyrosequences from a median of 2400 viral templates per subject and encompassing 40% of the HIV-1 genome, were compared to a median of five individually amplified near full-length viral genomes sequenced using Sanger technology. There was no difference in the consensus nucleotide sequences over the 3.6kb compared in 84% of the subjects infected with single founders and 33% of subjects infected with multiple founder variants: among the subjects with disagreements, mismatches were found in less than 1% of the sites evaluated (of a total of nearly 117,000 sites across all subjects). The majority of the SNPs observed only in pyrosequences were present at less than 2% of the subject’s viral sequence population. These results demonstrate the utility of the Sanger approach for study of early HIV infection and provide guidance regarding the design, utility and limitations of population sequencing from variable template sources, and emphasize parameters for improving the interpretation of massively parallel sequencing data to address important questions regarding target sequence evolution.     

Hierarchical gene expression profiles of HUVEC stimulated by different lipid A structures obtained from Porphyromonas gingivalis and Escherichia coli

The ability of lipid A structural variants to elicit unique endothelial cell gene expression was examined by measuring global gene expression profiles in human umbilical cord vein endothelial cells (HUVEC) using Affymetrix full genome chips. Two lipid A structural variants obtained from Porphyromonas gingivalis designated PgLPS(1435/1449) and PgLPS(1690) as well as LPS obtained from Escherichia coli wild type and an E. coli msbB mutant (missing myristic acid in the lipid A) were examined. Each of these lipid A structures has been shown to interact with TLR4; however, PgLPS(1435/1449) and E. coli msbB LPS have been shown to be TLR4 antagonists while PgLPS(1690) and wild-type E. coli LPS are TLR4 agonists. It was found that PgLPS(1435/1449) and PgLPS(1690) as well as E. coli msbB LPS activated a subset of those genes significantly transcribed in response to E. coli wild-type LPS. Furthermore, the subset of genes expressed in response to the different lipid A structural forms were those most significantly activated by wild-type E. coli LPS demonstrating a hierarchy in TLR4-dependent endothelial cell gene activation. A unique gene expression profile for the weak TLR4 agonist PgLPS(1690) was observed and represents a TLR4 hierarchy in endothelial cell gene activation.