The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Phylogenetic approach reveals that virus genotype largely determines HIV set-point viral load.

HIV virulence, i.e. the time of progression to AIDS, varies greatly among patients. As for other rapidly evolving pathogens of humans, it is difficult to know if this variance is controlled by the genotype of the host or that of the virus because the transmission chain is usually unknown. We apply the phylogenetic comparative approach (PCA) to estimate the heritability of a trait from one infection to the next, which indicates the control of the virus genotype over this trait. The idea is to use viral RNA sequences obtained from patients infected by HIV-1 subtype B to build a phylogeny, which approximately reflects the transmission chain. Heritability is measured statistically as the propensity for patients close in the phylogeny to exhibit similar infection trait values. The approach reveals that up to half of the variance in set-point viral load, a trait associated with virulence, can be heritable. Our estimate is significant and robust to noise in the phylogeny. We also check for the consistency of our approach by showing that a trait related to drug resistance is almost entirely heritable. Finally, we show the importance of taking into account the transmission chain when estimating correlations between infection traits. The fact that HIV virulence is, at least partially, heritable from one infection to the next has clinical and epidemiological implications. The difference between earlier studies and ours comes from the quality of our dataset and from the power of the PCA, which can be applied to large datasets and accounts for within-host evolution. The PCA opens new perspectives for approaches linking clinical data and evolutionary biology because it can be extended to study other traits or other infectious diseases.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Fluorescent amplified fragment length polymorphism (fAFLP) analyses and genetic diversity in Litopenaeus vannamei (Penaeidae)

The Pacific white shrimp, Litopenaeus vannamei (Penaeidae), represents about 95% of all Brazilian shrimp production. The Brazilian L. vannamei foundation broodstock was made up of specimens collected from different American Pacific sites, but little information was collected on the genetic structure of the broodstock. We used the fluorescence amplified fragment length polymorphism (fAFLP) method to study the genetic diversity of L. vannamei broodstock lines 03CMF1 and 03CBF1 originally produced by breeder-shrimps imported mainly from Panama and Ecuador, although wild individuals from other localities may also have been used in producing these two lines. Our results showed a total of 93 polymorphic bands ranging from 50 to 500 bp, the mean Nei's genetic diversity calculated for the total sample was 13.4% and identity and genetic distance analyses indicated high genetic homogeneity within and between both the broodstock lineages studied which suggests that they had similar genetic structure. These results may represent an important tool for the appropriate management of L. vannamei broodstocks. Copyright by the Brazilian Society of Genetics.

Estimation of the incidence of a rare genetic disease through a two-tier mutation survey.

Recent attempts to detect mutations involving single base changes or small deletions that are specific to genetic diseases provide an opportunity to develop a two-tier mutation-screening program through which incidence of rare genetic disorders and gene carriers may be precisely estimated. A two-tier survey consists of mutation screening in a sample of patients with specific genetic disorders and in a second sample of newborns from the same population in which mutation frequency is evaluated. We provide the statistical basis for evaluating the incidence of affected and gene carriers in such two-tier mutation-screening surveys, from which the precision of the estimates is derived. Sample-size requirements of such two-tier mutation-screening surveys are evaluated. Considering examples of cystic fibrosis (CF) and medium-chain acyl-CoA dehydrogenase deficiency (MCAD), the two most frequent autosomal recessive disease in Caucasian populations and the two most frequent mutations (delta F508 and G985) that occur on these disease allele-bearing chromosomes, we show that, with 50-100 patients and a 20-fold larger sample of newborns screened for these mutations, the incidence of such diseases and their gene carriers in a population may be quite reliably estimated. The theory developed here is also applicable to rare autosomal dominant diseases for which disease-specific mutations are found.

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Complement facilitates early prion pathogenesis.

New-variant Creutzfeldt-Jakob disease and scrapie are typically initiated by extracerebral exposure to the causative agent, and exhibit early prion replication in lymphoid organs. In mouse scrapie, depletion of B-lymphocytes prevents neuropathogenesis after intraperitoneal inoculation, probably due to impaired lymphotoxin-dependent maturation of follicular dendritic cells (FDCs), which are a major extracerebral prion reservoir. FDCs trap immune complexes with Fc-gamma receptors and C3d/C4b-opsonized antigens with CD21/CD35 complement receptors. We examined whether these mechanisms participate in peripheral prion pathogenesis. Depletion of circulating immunoglobulins or of individual Fc-gamma receptors had no effect on scrapie pathogenesis if B-cell maturation was unaffected. However, mice deficient in C3, C1q, Bf/C2, combinations thereof or complement receptors were partially or fully protected against spongiform encephalopathy upon intraperitoneal exposure to limiting amounts of prions. Splenic accumulation of prion infectivity and PrPSc was delayed, indicating that activation of specific complement components is involved in the initial trapping of prions in lymphoreticular organs early after infection.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values &lt; or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values &lt; or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

Remodeling of DNA replication structures by the branch point translocase FANCM

Fanconi anemia (FA) is a genetically heterogeneous chromosome instability syndrome associated with congenital abnormalities, bone marrow failure, and cancer predisposition. Eight FA proteins form a nuclear core complex, which promotes tolerance of DNA lesions in S phase, but the underlying mechanisms are still elusive. We reported recently that the FA core complex protein FANCM can translocate Holliday junctions. Here we show that FANCM promotes reversal of model replication forks via concerted displacement and annealing of the nascent and parental DNA strands. Fork reversal by FANCM also occurs when the lagging strand template is partially single-stranded and bound by RPA. The combined fork reversal and branch migration activities of FANCM lead to extensive regression of model replication forks. These observations provide evidence that FANCM can remodel replication fork structures and suggest a mechanism by which FANCM could promote DNA damage tolerance in S phase

Neurological complications of infective endocarditis: risk factors, outcome, and impact of cardiac surgery: a multicenter observational study.

BACKGROUND The purpose of this study was to assess the incidence of neurological complications in patients with infective endocarditis, the risk factors for their development, their influence on the clinical outcome, and the impact of cardiac surgery. METHODS AND RESULTS This was a retrospective analysis of prospectively collected data on a multicenter cohort of 1345 consecutive episodes of left-sided infective endocarditis from 8 centers in Spain. Cox regression models were developed to analyze variables predictive of neurological complications and associated mortality. Three hundred forty patients (25%) experienced such complications: 192 patients (14%) had ischemic events, 86 (6%) had encephalopathy/meningitis, 60 (4%) had hemorrhages, and 2 (1%) had brain abscesses. Independent risk factors associated with all neurological complications were vegetation size ≥3 cm (hazard ratio [HR] 1.91), Staphylococcus aureus as a cause (HR 2.47), mitral valve involvement (HR 1.29), and anticoagulant therapy (HR 1.31). This last variable was particularly related to a greater incidence of hemorrhagic events (HR 2.71). Overall mortality was 30%, and neurological complications had a negative impact on outcome (45% of deaths versus 24% in patients without these complications; P<0.01), although only moderate to severe ischemic stroke (HR 1.63) and brain hemorrhage (HR 1.73) were significantly associated with a poorer prognosis. Antimicrobial treatment reduced (by 33% to 75%) the risk of neurological complications. In patients with hemorrhage, mortality was higher when surgery was performed within 4 weeks of the hemorrhagic event (75% versus 40% in later surgery). CONCLUSIONS Moderate to severe ischemic stroke and brain hemorrhage were found to have a significant negative impact on the outcome of infective endocarditis. Early appropriate antimicrobial treatment is critical, and transitory discontinuation of anticoagulant therapy should be considered.

Parallel changes in genetic diversity and species diversity following a natural disturbance.

We examined the spatial and temporal variation of species diversity and genetic diversity in a metacommunity comprising 16 species of freshwater gastropods. We monitored species abundance at five localities of the Ain river floodplain in southeastern France, over a period of four years. Using 190 AFLP loci, we monitored the genetic diversity of Radix balthica, one of the most abundant gastropod species of the metacommunity, twice during that period. An exceptionally intense drought occurred during the last two years and differentially affected the study sites. This allowed us to test the effect of natural disturbances on changes in both genetic and species diversity. Overall, local (alpha) diversity declined as reflected by lower values of gene diversity H(S) and evenness. In parallel, the among-sites (beta) diversity increased at both the genetic (F(ST)) and species (F(STC)) levels. These results suggest that disturbances can lead to similar changes in genetic and community structure through the combined effects of selective and neutral processes.

Effect of statin therapy in the outcome of bloodstream infections due to Staphylococcus aureus: a prospective cohort study.

INTRODUCTION Statins have pleiotropic effects that could influence the prevention and outcome of some infectious diseases. There is no information about their specific effect on Staphylococcus aureus bacteremia (SAB). METHODS A prospective cohort study including all SAB diagnosed in patients aged ≥18 years admitted to a 950-bed tertiary hospital from March 2008 to January 2011 was performed. The main outcome variable was 14-day mortality, and the secondary outcome variables were 30-day mortality, persistent bacteremia (PB) and presence of severe sepsis or septic shock at diagnosis of SAB. The effect of statin therapy at the onset of SAB was studied by multivariate logistic regression and Cox regression analysis, including a propensity score for statin therapy. RESULTS We included 160 episodes. Thirty-three patients (21.3%) were receiving statins at the onset of SAB. 14-day mortality was 21.3%. After adjustment for age, Charlson index, Pitt score, adequate management, and high risk source, statin therapy had a protective effect on 14-day mortality (adjusted OR = 0.08; 95% CI: 0.01-0.66; p = 0.02), and PB (OR = 0.89; 95% CI: 0.27-1.00; p = 0.05) although the effect was not significant on 30-day mortality (OR = 0.35; 95% CI: 0.10-1.23; p = 0.10) or presentation with severe sepsis or septic shock (adjusted OR = 0.89; CI 95%: 0.27-2.94; p = 0.8). An effect on 30-day mortality could neither be demonstrated on Cox analysis (adjusted HR = 0.5; 95% CI: 0.19-1.29; p = 0.15). CONCLUSIONS Statin treatment in patients with SAB was associated with lower early mortality and PB. Randomized studies are necessary to identify the role of statins in the treatment of patients with SAB.