An efficient algorithm to perform multiple testing in epistasis screening

Background: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn’s disease. Results: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn’s disease (CD) data. Conclusions: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn’s disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

Transcription factor PROX1 induces colon cancer progression by promoting the transition from benign to highly dysplastic phenotype.

The Drosophila transcription factor Prospero functions as a tumor suppressor, and it has been suggested that the human counterpart of Prospero, PROX1, acts similarly in human cancers. However, we show here that PROX1 promotes dysplasia in colonic adenomas and colorectal cancer progression. PROX1 expression marks the transition from benign colon adenoma to carcinoma in situ, and its loss inhibits growth of human colorectal tumor xenografts and intestinal adenomas in Apc(min/+) mice, while its transgenic overexpression promotes colorectal tumorigenesis. Furthermore, in intestinal tumors PROX1 is a direct and dose-dependent target of the beta-catenin/TCF signaling pathway, responsible for the neoplastic transformation. Our data underscore the complexity of cancer pathogenesis and implicate PROX1 in malignant tumor progression through the regulation of cell polarity and adhesion.

Schizotypal perceptual aberrations of time: correlation between score, behavior and brain activity.

A fundamental trait of the human self is its continuum experience of space and time. Perceptual aberrations of this spatial and temporal continuity is a major characteristic of schizophrenia spectrum disturbances--including schizophrenia, schizotypal personality disorder and schizotypy. We have previously found the classical Perceptual Aberration Scale (PAS) scores, related to body and space, to be positively correlated with both behavior and temporo-parietal activation in healthy participants performing a task involving self-projection in space. However, not much is known about the relationship between temporal perceptual aberration, behavior and brain activity. To this aim, we composed a temporal Perceptual Aberration Scale (tPAS) similar to the traditional PAS. Testing on 170 participants suggested similar performance for PAS and tPAS. We then correlated tPAS and PAS scores to participants' performance and neural activity in a task of self-projection in time. tPAS scores correlated positively with reaction times across task conditions, as did PAS scores. Evoked potential mapping and electrical neuroimaging showed self-projection in time to recruit a network of brain regions at the left anterior temporal cortex, right temporo-parietal junction, and occipito-temporal cortex, and duration of activation in this network positively correlated with tPAS and PAS scores. These data demonstrate that schizotypal perceptual aberrations of both time and space, as reflected by tPAS and PAS scores, are positively correlated with performance and brain activation during self-projection in time in healthy individuals along the schizophrenia spectrum.

Mineralocorticoid effects in the kidney: correlation between alphaENaC, GILZ, and Sgk-1 mRNA expression and urinary excretion of Na+ and K+.

Aldosterone exerts its effects through interactions with two types of binding sites, the mineralocorticoid (MR) and the glucocorticoid (GR) receptors. Although both receptors are known to be involved in the anti-natriuretic response to aldosterone, the mechanisms of signal transduction leading to modulation of electrolyte transport are not yet fully understood. This study measured the Na(+) and K(+) urinary excretion and the mRNA levels of three known aldosterone-induced transcripts, the serum and glucocorticoid-induced kinase (Sgk-1), the alpha subunit of the epithelial Na(+) channel (alphaENaC), and the glucocorticoid-induced-leucine-zipper protein (GILZ) in the whole kidney and in isolated cortical collecting tubules of adrenalectomized rats treated with low doses of aldosterone and/or dexamethasone. The resulting plasma concentrations of both steroids were close to 1 nmol/L. Aldosterone, given with or without dexamethasone, induced anti-natriuresis and kaliuresis, whereas dexamethasone alone did not. GILZ and alphaENaC transcripts were higher after treatment with either or both hormones, whereas the mRNA abundance of Sgk-1 was increased in the cortical collecting tubule by aldosterone but not by dexamethasone. We conclude the increased expression of Sgk-1 in the cortical collecting tubules is a primary event in the early antinatriuretic and kaliuretic responses to physiologic concentrations of aldosterone. Induction of alphaENaC and/or GILZ mRNAs may play a permissive role in the enhancement of the early and/or late responses; these effects may be necessary for a full response but do not by themselves promote early changes in urinary Na(+) and K(+) excretion.

ABCB1 and cytochrome P450 genotypes and phenotypes: influence on methadone plasma levels and response to treatment.

BACKGROUND AND OBJECTIVE: The in vivo implication of various cytochrome P450 (CYP) isoforms and of P-glycoprotein on methadone kinetics is unclear. We aimed to thoroughly examine the genetic factors influencing methadone kinetics and response to treatment. METHODS: Genotyping for CYP1A2, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, ABCB1, and UGT2B7 polymorphisms was performed in 245 patients undergoing methadone maintenance treatment. To assess CYP3A activity, the patients were phenotyped with midazolam. RESULTS: The patients with lower CYP3A activity presented higher steady-state trough (R,S)-methadone plasma levels (4.3, 3.0, and 2.3 ng/mL x mg for low, medium, and high activity, respectively; P = .0002). As previously reported, CYP2B6*6/*6 carriers had significantly higher trough (S)-methadone plasma levels (P = .0001) and a trend toward higher (R)-methadone plasma levels (P = .07). CYP2D6 ultrarapid metabolizers presented lower trough (R,S)-methadone plasma levels compared with the extensive or intermediate metabolizers (2.4 and 3.3 ng/mL x mg, respectively; P = .04), whereas CYP2D6 poor metabolizer status showed no influence. ABCB1 3435TT carriers presented lower trough (R,S)-methadone plasma levels (2.7 and 3.4 ng/mL . mg for 3435TT and 3435CC carriers, respectively; P = .01). The CYP1A2, CYP2C9, CYP2C19, CYP3A5, and UGT2B7 genotypes did not influence methadone plasma levels. Only CYP2B6 displayed a stereoselectivity in its activity. CONCLUSION: In vivo, CYP3A4 and CYP2B6 are the major CYP isoforms involved in methadone metabolism, with CYP2D6 contributing to a minor extent. ABCB1 genetic polymorphisms also contribute slightly to the interindividual variability of methadone kinetics. The genetic polymorphisms of these 4 proteins had no influence on the response to treatment and only a small influence on the dose requirement of methadone.

Copy number variants, diseases and gene expression.

Copy number variation (CNV) has recently gained considerable interest as a source of genetic variation likely to play a role in phenotypic diversity and evolution. Much effort has been put into the identification and mapping of regions that vary in copy number among seemingly normal individuals in humans and a number of model organisms, using bioinformatics or hybridization-based methods. These have allowed uncovering associations between copy number changes and complex diseases in whole-genome association studies, as well as identify new genomic disorders. At the genome-wide scale, however, the functional impact of CNV remains poorly studied. Here we review the current catalogs of CNVs, their association with diseases and how they link genotype and phenotype. We describe initial evidence which revealed that genes in CNV regions are expressed at lower and more variable levels than genes mapping elsewhere, and also that CNV not only affects the expression of genes varying in copy number, but also have a global influence on the transcriptome. Further studies are warranted for complete cataloguing and fine mapping of CNVs, as well as to elucidate the different mechanisms by which they influence gene expression.

Genotype and environment interaction on soybean yield in Mato Grosso State, Brazil

The objective of this work was to investigate the genotype-environment interaction in Mato Grosso State, MT. The relative importance of locations, years, sowing dates and cultivars and their interactions was analyzed from data collected in regional yield trials performed in a randomized complete block design with four replications, from 1994-1995 through 1999-2000, in nine locations and two sowing dates. Individual and pooled analyses of variance over years and locations were performed. Complementary analyses of variances partitioned MT State in two main and five smaller regions, respectively: North and South of Cuiabá; and MT-South-A (Pedra Preta area), MT-South-B (Rondonópolis and Itiquira), MT-East (Primavera do Leste and Campo Verde), MT-Central (Nova Mutum, Lucas do Rio Verde and Sorriso) and MT-Parecis (Campo Novo dos Parecis and Sapezal). Locations are relatively more important than years for yield testing soybeans in the MT State, therefore, investment should be made in increasing locations rather than years to improve experimental precision. Partitioning the MT State into regions has little impact on soybean yield testing results and, consequently, on the efficiency of the soybean breeding program in the State. Breeding genotypes with broad adaptation for the MT State is an efficient strategy for cultivar development.

Ultrasound evaluation of synovitis in RA: Correlation with clinical disease activity and sensitivity to change in an observational cohort study.

OBJECTIVE: To evaluate the correlation between clinical measures of disease activity and a ultrasound (US) scoring system for synovitis applied by many different ultrasonographers in a daily routine care setting within the Swiss registry for RA (SCQM) and further to determine the sensitivity to change of this US Score. METHODS: One hundred and eight Swiss rheumatologists were trained in performing the Swiss Sonography in Arthritis and Rheumatism (SONAR) score. US B-mode and Power Doppler (PwD) scores were correlated with DAS28 and compared between the clinical categories in a cross-sectional cohort of patients. In patients with a second US (longitudinal cohort), we investigated if change in US score correlated with change in DAS and evaluated the responsiveness of both methods. RESULTS: In the cross-sectional cohort with 536 patients, correlation between the B-mode score and DAS28 was significant but modest (Pearson coefficient r=0.41, P<0.0001). The same was true for the PwD score (r=0.41, P<0.0001). In the longitudinal cohort with 183 patients we also found a significant correlation between change in B-mode and in PwD score with change in DAS28 (r=0.54, P<0.0001 and r=0.46, P<0.0001, respectively). Both methods of evaluation (DAS and US) showed similar responsiveness according to standardized response mean (SRM). CONCLUSIONS: The SONAR Score is practicable and was applied by many rheumatologists in daily routine care after initial training. It demonstrates significant correlations with the degree of as well as change in disease activity as measured by DAS. On the level of the individual, the US score shows many discrepancies and overlapping results exist.

Immunohistochemistry as a valuable tool to assess CD30 expression in peripheral T-cell lymphomas: high correlation with mRNA levels.

The extended use of brentuximab-vedotin was reported for CD30(+) nonanaplastic peripheral T-cell lymphomas (PTCLs) with promising efficacy. CD30 status assessment is thus a critical factor for therapeutic decision, but the reliability of immunohistochemistry (IHC) in evaluating its expression remains to be defined. This prompted us to investigate the correlation between semiquantitative CD30 protein assessment by IHC and messenger RNA (mRNA) assessment by microarrays in a cohort of 376 noncutaneous PTCLs representative of the main entities. By IHC, CD30 expression was heterogeneous across and within entities and significantly associated with large tumor cell size. In addition to 100% anaplastic large-cell lymphomas, 57% of other PTCL entities were CD30-positive at a 5% threshold. CD30 protein expression was highly correlated to mRNA levels. mRNA levels were bimodal, separating high from low CD30-expressing PTCL cases. We conclude that IHC is a valuable tool in clinical practice to assess CD30 expression in PTCLs.

Visual cortex in Alzheimer's disease: occurrence of neuronal death and glial proliferation, and correlation with pathological hallmarks.

Visual areas 17 and 18 were studied with morphometric methods for numbers of neurons, glia, senile plaques (SP), and neurofibrillary tangles (NFT) in 13 cases of Alzheimer's disease (AD) as compared to 11 controls. In AD cases, the mean neuronal density was significantly decreased by about 30% in both areas 17 and 18, while the glial density was increased significantly only in area 17. The volume of area 17 was unchanged in AD cases but its total number of neurons was decreased by 33% and its total number of glia increased by 45% compared to controls. In AD the number of SP was similar in areas 17 and 18, while that of NFT was significantly higher in area 18. The number of neurons with NFT was only 2% in area 17 and about 10% in area 18. The discrepancy between the loss of neurons and the amount of NFT suggests that neuronal loss can occur without passing through NFT degeneration. The deposition of SP was correlated with glial proliferation, but not with neuronal loss or neurofibrillary degeneration.

Haemophilia registry of the medical committee of the Swiss Haemophilia Society. Update and annual survey 2009.

The Swiss Haemophilia Registry of the Medical Committee of the Swiss Haemophilia Society started in 1996 but was set as an internet-based, double password-protected facility in the year 2000. With the inclusion of patients' data from two new centres in 2009, we assume a coverage rate of about 90% of all patients with inherited bleeding disorders in our country. Data concerning the phenotype and genotype of the disorder, its severity, its therapy, the prevalence of inhibitors are readily available to the registered users, allowing quality control of haemophilia therapy at a national level, but also rapid care of the patient visiting the emergency room of another treatment centre. Basing on the available data, about two thirds of the WFH global survey can be answered; the mortality statistics shows that bleeding remains a cause of death in haemophiliacs, also in the 21th century. The Registry allows for comparisons with international datasets, especially with respect to treatment (prophylaxis vs. on-demand therapy), factor consumption and costs.

Polymorphisms in Toll-like receptor 9 influence the clinical course of HIV-1 infection.

BACKGROUND: The clinical course of HIV-1 infection is highly variable among individuals, at least in part as a result of genetic polymorphisms in the host. Toll-like receptors (TLRs) have a key role in innate immunity and mutations in the genes encoding these receptors have been associated with increased or decreased susceptibility to infections. OBJECTIVES: To determine whether single-nucleotide polymorphisms (SNPs) in TLR2-4 and TLR7-9 influenced the natural course of HIV-1 infection. METHODS: Twenty-eight SNPs in TLRs were analysed in HAART-naive HIV-positive patients from the Swiss HIV Cohort Study. The SNPs were detected using Sequenom technology. Haplotypes were inferred using an expectation-maximization algorithm. The CD4 T cell decline was calculated using a least-squares regression. Patients with a rapid CD4 cell decline, less than the 15th percentile, were defined as rapid progressors. The risk of rapid progression associated with SNPs was estimated using a logistic regression model. Other candidate risk factors included age, sex and risk groups (heterosexual, homosexual and intravenous drug use). RESULTS: Two SNPs in TLR9 (1635A/G and +1174G/A) in linkage disequilibrium were associated with the rapid progressor phenotype: for 1635A/G, odds ratio (OR), 3.9 [95% confidence interval (CI),1.7-9.2] for GA versus AA and OR, 4.7 (95% CI,1.9-12.0) for GG versus AA (P = 0.0008). CONCLUSION: Rapid progression of HIV-1 infection was associated with TLR9 polymorphisms. Because of its potential implications for intervention strategies and vaccine developments, additional epidemiological and experimental studies are needed to confirm this association.

Correlation of the Road Rater and the Dynatest Falling Weight Deflectometer, MLR-91-04, 1991

The Falling Weight Deflectometer (FWD) has become the "standard" for deflection testing of pavements. Iowa has used a Road Rater since 1976 to obtain deflection information. A correlation between the Road Rater and the FWD was needed if Iowa was going to continue with the Road Rater. Comparative deflection testing was done using a Road Rater Model 400 and a Pynatest 8000 FWD on 26 pavement sections. The SHRP contractor, Braun Intertec Pavement, Inc., provided the FWD testing. The r^2 for the linear correlations ranged from 0.90 to 0.99 for the different pavement types and sensor locations.