Classes of Multiple Decision Functions Strongly Controlling FWER and FDR

This paper provides two general classes of multiple decision functions where each member of the first class strongly controls the family-wise error rate (FWER), while each member of the second class strongly controls the false discovery rate (FDR). These classes offer the possibility that an optimal multiple decision function with respect to a pre-specified criterion, such as the missed discovery rate (MDR), could be found within these classes. Such multiple decision functions can be utilized in multiple testing, specifically, but not limited to, the analysis of high-dimensional microarray data sets.Comment: 19 page

Caterpillars and fungal pathogens: two co-occurring parasites of an ant-plant mutualism

In mutualisms, each interacting species obtains resources from its partner that it would obtain less efficiently if alone, and so derives a net fitness benefit. In exchange for shelter (domatia) and food, mutualistic plant-ants protect their host myrmecophytes from herbivores, encroaching vines and fungal pathogens. Although selective filters enable myrmecophytes to host those ant species most favorable to their fitness, some insects can by-pass these filters, exploiting the rewards supplied whilst providing nothing in return. This is the case in French Guiana for Cecropia obtusa (Cecropiaceae) as Pseudocabima guianalis caterpillars (Lepidoptera, Pyralidae) can colonize saplings before the installation of their mutualistic Azteca ants. The caterpillars shelter in the domatia and feed on food bodies (FBs) whose production increases as a result. They delay colonization by ants by weaving a silk shield above the youngest trichilium, where the FBs are produced, blocking access to them. This probable temporal priority effect also allows female moths to lay new eggs on trees that already shelter caterpillars, and so to occupy the niche longer and exploit Cecropia resources before colonization by ants. However, once incipient ant colonies are able to develop, they prevent further colonization by the caterpillars. Although no higher herbivory rates were noted, these caterpillars are ineffective in protecting their host trees from a pathogenic fungus, Fusarium moniliforme (Deuteromycetes), that develops on the trichilium in the absence of mutualistic ants. Therefore, the Cecropia treelets can be parasitized by two often overlooked species: the caterpillars that shelter in the domatia and feed on FBs, delaying colonization by mutualistic ants, and the fungal pathogen that develops on old trichilia. The cost of greater FB production plus the presence of the pathogenic fungus likely affect tree growth

Large-scale associations between the leukocyte transcriptome and BOLD responses to speech differ in autism early language outcome subtypes.

Heterogeneity in early language development in autism spectrum disorder (ASD) is clinically important and may reflect neurobiologically distinct subtypes. Here, we identified a large-scale association between multiple coordinated blood leukocyte gene coexpression modules and the multivariate functional neuroimaging (fMRI) response to speech. Gene coexpression modules associated with the multivariate fMRI response to speech were different for all pairwise comparisons between typically developing toddlers and toddlers with ASD and poor versus good early language outcome. Associated coexpression modules were enriched in genes that are broadly expressed in the brain and many other tissues. These coexpression modules were also enriched in ASD-associated, prenatal, human-specific, and language-relevant genes. This work highlights distinctive neurobiology in ASD subtypes with different early language outcomes that is present well before such outcomes are known. Associations between neuroimaging measures and gene expression levels in blood leukocytes may offer a unique in vivo window into identifying brain-relevant molecular mechanisms in ASD

Biological-effective versus conventional dose volume histograms correlated with late genitourinary and gastrointestinal toxicity after external beam radiotherapy for prostate cancer: a matched pair analysis

BACKGROUND: To determine whether the dose-volume histograms (DVH's) for the rectum and bladder constructed using biological-effective dose (BED-DVH's) better correlate with late gastrointestinal (GI) and genitourinary (GU) toxicity after treatment with external beam radiotherapy for prostate cancer than conventional DVH's (C-DVH's). METHODS: The charts of 190 patients treated with external beam radiotherapy with a minimum follow-up of 2 years were reviewed. Six patients (3.2%) were found to have RTOG grade 3 GI toxicity, and similarly 6 patients (3.2%) were found to have RTOG grade 3 GU toxicity. Average late C-DVH's and BED-DVH's of the bladder and rectum were computed for these patients as well as for matched-pair control patients. For each matched pair the following measures of normalized difference in the DVH's were computed: (a) δ(AUC )= (Area Under Curve [AUC] in grade 3 patient – AUC in grade 0 patient)/(AUC in grade 0 patient) and (b) δ(V60 )= (Percent volume receiving = 60 Gy [V60] in grade 3 patient – V60 in grade 0 patient)/(V60 in grade 0 patient). RESULTS: As expected, the grade 3 curve is to the right of and above the grade 0 curve for all four sets of average DVH's – suggesting that both the C-DVH and the BED-DVH can be used for predicting late toxicity. δ(AUC )was higher for the BED-DVH's than for the C-DVH's – 0.27 vs 0.23 (p = 0.036) for the rectum and 0.24 vs 0.20 (p = 0.065) for the bladder. δ(V60 )was also higher for the BED-DVH's than for the C-DVH's – 2.73 vs 1.49 for the rectum (p = 0.021) and 1.64 vs 0.71 (p = 0.021) for the bladder. CONCLUSIONS: When considering well-established dosimetric endpoints used in evaluating treatment plans, BED-DVH's for the rectum and bladder correlate better with late toxicity than C-DVH's and should be considered when attempting to minimize late GI and GU toxicity after external beam radiotherapy for prostate cancer

The cost of large numbers of hypothesis tests on power, effect size and sample size

Advances in high-throughput biology and computer science are driving an exponential increase in the number of hypothesis tests in genomics and other scientific disciplines. Studies using current genotyping platforms frequently include a million or more tests. In addition to the monetary cost, this increase imposes a statistical cost owing to the multiple testing corrections needed to avoid large numbers of false-positive results. To safeguard against the resulting loss of power, some have suggested sample sizes on the order of tens of thousands that can be impractical for many diseases or may lower the quality of phenotypic measurements. This study examines the relationship between the number of tests on the one hand and power, detectable effect size or required sample size on the other. We show that once the number of tests is large, power can be maintained at a constant level, with comparatively small increases in the effect size or sample size. For example at the 0.05 significance level, a 13% increase in sample size is needed to maintain 80% power for ten million tests compared with one million tests, whereas a 70% increase in sample size is needed for 10 tests compared with a single test. Relative costs are less when measured by increases in the detectable effect size. We provide an interactive Excel calculator to compute power, effect size or sample size when comparing study designs or genome platforms involving different numbers of hypothesis tests. The results are reassuring in an era of extreme multiple testing

The Evolution of Gene Expression QTL in Saccharomyces cerevisiae

Understanding the evolutionary forces that influence patterns of gene expression variation will provide insights into the mechanisms of evolutionary change and the molecular basis of phenotypic diversity. To date, studies of gene expression evolution have primarily been made by analyzing how gene expression levels vary within and between species. However, the fundamental unit of heritable variation in transcript abundance is the underlying regulatory allele, and as a result it is necessary to understand gene expression evolution at the level of DNA sequence variation. Here we describe the evolutionary forces shaping patterns of genetic variation for 1206 cis-regulatory QTL identified in a cross between two divergent strains of Saccharomyces cerevisiae. We demonstrate that purifying selection against mildly deleterious alleles is the dominant force governing cis-regulatory evolution in S. cerevisiae and estimate the strength of selection. We also find that essential genes and genes with larger codon bias are subject to slightly stronger cis-regulatory constraint and that positive selection has played a role in the evolution of major trans-acting QTL

Sample size calculation for microarray experiments with blocked one-way design

<p>Abstract</p> <p>Background</p> <p>One of the main objectives of microarray analysis is to identify differentially expressed genes for different types of cells or treatments. Many statistical methods have been proposed to assess the treatment effects in microarray experiments.</p> <p>Results</p> <p>In this paper, we consider discovery of the genes that are differentially expressed among <it>K </it>(> 2) treatments when each set of <it>K </it>arrays consists of a block. In this case, the array data among <it>K </it>treatments tend to be correlated because of block effect. We propose to use the blocked one-way ANOVA <it>F</it>-statistic to test if each gene is differentially expressed among <it>K </it>treatments. The marginal p-values are calculated using a permutation method accounting for the block effect, adjusting for the multiplicity of the testing procedure by controlling the false discovery rate (FDR). We propose a sample size calculation method for microarray experiments with a blocked one-way design. With FDR level and effect sizes of genes specified, our formula provides a sample size for a given number of true discoveries.</p> <p>Conclusion</p> <p>The calculated sample size is shown via simulations to provide an accurate number of true discoveries while controlling the FDR at the desired level.</p

Complexity in the genetic architecture of leukoaraiosis in hypertensive sibships from the GENOA Study

<p>Abstract</p> <p>Background</p> <p>Subcortical white matter hyperintensity on magnetic resonance imaging (MRI) of the brain, referred to as leukoaraiosis, is associated with increased risk of stroke and dementia. Hypertension may contribute to leukoaraiosis by accelerating the process of arteriosclerosis involving penetrating small arteries and arterioles in the brain. Leukoaraiosis volume is highly heritable but shows significant inter-individual variability that is not predicted well by any clinical covariates (except for age) or by single SNPs.</p> <p>Methods</p> <p>As part of the Genetics of Microangiopathic Brain Injury (GMBI) Study, 777 individuals (74% hypertensive) underwent brain MRI and were genotyped for 1649 SNPs from genes known or hypothesized to be involved in arteriosclerosis and related pathways. We examined SNP main effects, epistatic (gene-gene) interactions, and context-dependent (gene-environment) interactions between these SNPs and covariates (including conventional and novel risk factors for arteriosclerosis) for association with leukoaraiosis volume. Three methods were used to reduce the chance of false positive associations: 1) false discovery rate (FDR) adjustment for multiple testing, 2) an internal replication design, and 3) a ten-iteration four-fold cross-validation scheme.</p> <p>Results</p> <p>Four SNP main effects (in <it>F3</it>, <it>KITLG</it>, <it>CAPN10</it>, and <it>MMP2</it>), 12 SNP-covariate interactions (including interactions between <it>KITLG </it>and homocysteine, and between <it>TGFB3 </it>and both physical activity and C-reactive protein), and 173 SNP-SNP interactions were significant, replicated, and cross-validated. While a model containing the top single SNPs with main effects predicted only 3.72% of variation in leukoaraiosis in independent test samples, a multiple variable model that included the four most highly predictive SNP-SNP and SNP-covariate interactions predicted 11.83%.</p> <p>Conclusion</p> <p>These results indicate that the genetic architecture of leukoaraiosis is complex, yet predictive, when the contributions of SNP main effects are considered in combination with effects of SNP interactions with other genes and covariates.</p

Height and timing of growth spurt during puberty in young people living with vertically acquired HIV in Europe and Thailand.

OBJECTIVE: The aim of this study was to describe growth during puberty in young people with vertically acquired HIV. DESIGN: Pooled data from 12 paediatric HIV cohorts in Europe and Thailand. METHODS: One thousand and ninety-four children initiating a nonnucleoside reverse transcriptase inhibitor or boosted protease inhibitor based regimen aged 1-10 years were included. Super Imposition by Translation And Rotation (SITAR) models described growth from age 8 years using three parameters (average height, timing and shape of the growth spurt), dependent on age and height-for-age z-score (HAZ) (WHO references) at antiretroviral therapy (ART) initiation. Multivariate regression explored characteristics associated with these three parameters. RESULTS: At ART initiation, median age and HAZ was 6.4 [interquartile range (IQR): 2.8, 9.0] years and -1.2 (IQR: -2.3 to -0.2), respectively. Median follow-up was 9.1 (IQR: 6.9, 11.4) years. In girls, older age and lower HAZ at ART initiation were independently associated with a growth spurt which occurred 0.41 (95% confidence interval 0.20-0.62) years later in children starting ART age 6 to 10 years compared with 1 to 2 years and 1.50 (1.21-1.78) years later in those starting with HAZ less than -3 compared with HAZ at least -1. Later growth spurts in girls resulted in continued height growth into later adolescence. In boys starting ART with HAZ less than -1, growth spurts were later in children starting ART in the oldest age group, but for HAZ at least -1, there was no association with age. Girls and boys who initiated ART with HAZ at least -1 maintained a similar height to the WHO reference mean. CONCLUSION: Stunting at ART initiation was associated with later growth spurts in girls. Children with HAZ at least -1 at ART initiation grew in height at the level expected in HIV negative children of a comparable age