ABCtoolbox: a versatile toolkit for approximate Bayesian computations.

BACKGROUND: The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. RESULTS: Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. CONCLUSION: ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

Antigenic heterogeneity of clones and subclones from human melanoma cell lines demonstrated by a panel of monoclonal antibodies and flow microfluorometry analysis.

Cells from two melanoma cell lines, Me43 and GLL-19, were cloned in methylcellulose cultures and 20 randomly selected colonies from each line were picked up by micromanipulation, expanded in liquid cultures, and considered as clones of the original cell lines. The antigenic cell surface phenotype of these clones defined by panel of 12 monoclonal antibodies (MAb) was analyzed by flow microfluorometry (FMF) using a fluorescence-activated cell sorter (FACS II) and compared with the known stable phenotype of the parent cell line. The antibody panel consisted of eight MAb against melanoma-associated antigens, two MAb against monomorphic determinants of HLA-DR (la) and HLA-ABC, respectively, one MAb against the common acute lymphoblastic leukemia antigen (CALLA) and one MAb against carcinoembryonic antigen used as control. A remarkable heterogeneity in terms of qualitative and quantitative expression of the cell surface antigens studied was observed among and within the different clones. The single-cell origin of the clones was assessed by comparing the clonogenic cell frequency, determined by limiting dilutions in microculture plates, with the cloning efficiency observed in Petri dishes. Both techniques using methylcellulose medium gave the same percentages of growing colonies. Cells from four Me43 clones were recloned in methylcellulose and the phenotype of five randomly selected subclones from each clone was analysed using the same panel of monoclonal antibodies. Each subclone also displayed heterogeneity with individual phenotypes different from that of the original clone and from the parental Me43 cell line. The antigen expression by individual cells in situ within clones was analyzed on frozen sections from colonies using the same panel of MAb and a biotin-avidin immunoperoxidase method. The results confirmed the marked heterogeneity of antigen expression within and among colonies, as indicated by the FMF analysis.

Overcoming the heterologous bias: an in vivo functional analysis of multidrug efflux transporter, CgCdr1p in matched pair clinical isolates of Candida glabrata.

We have taken advantage of the natural milieu of matched pair of azole sensitive (AS) and azole resistant (AR) clinical isolates of Candida glabrata for expressing its major ABC multidrug transporter, CgCdr1p for structure and functional analysis. This was accomplished by tagging a green fluorescent protein (GFP) downstream of ORF of CgCDR1 and integrating the resultant fusion protein at its native chromosomal locus in AS and AR backgrounds. The characterization confirmed that in comparison to AS isolate, CgCdr1p-GFP was over-expressed in AR isolates due to its hyperactive native promoter and the GFP tag did not affect its functionality in either construct. We observed that in addition to Rhodamine 6 G (R6G) and Fluconazole (FLC), a recently identified fluorescent substrate of multidrug transporters Nile Red (NR) could also be expelled by CgCdr1p. Competition assays with these substrates revealed the presence of overlapping multiple drug binding sites in CgCdr1p. Point mutations employing site directed mutagenesis confirmed that the role played by unique amino acid residues critical to ATP catalysis and localization of ABC drug transporter proteins are well conserved in C. glabrata as in other yeasts. This study demonstrates a first in vivo novel system where over-expression of GFP tagged MDR transporter protein can be driven by its own hyperactive promoter of AR isolates. Taken together, this in vivo system can be exploited for the structure and functional analysis of CgCdr1p and similar proteins wherein the artefactual concerns encountered in using heterologous systems are totally excluded.

Similarities and dissimilarities between the Movement ABC-2 and the Zurich Neuromotor Assessment in children with suspected developmental coordination disorder.

An established tool for the assessment of motor performance in children with developmental coordination disorder (DCD) is the Movement-ABC-2 (M-ABC-2). The Zurich Neuromotor Assessment (ZNA) is also widely used for the evaluation of children's motor performance, but has not been compared with the M-ABC-2. Fifty-one children (39 males) between 5 and 7 years of age with suspected DCD were assessed using the M-ABC-2 and the ZNA. Rank correlations between scores of different test components were calculated. The structure of the tests was explored using canonical-correlation analysis. The correlation between total scores of the two motor tests was reasonable (0.66; p<0.001). However, ZNA scores were generally lower than those of M-ABC-2, due to poor performance in the fine motor adaptive component and increased contralateral associated movements (CAM). The canonical-correlation analysis revealed that ZNA measures components like pure motor skills and CAM that are not represented in the M-ABC-2. Furthermore, there was also no equivalent for the aiming and catching items of the M-ABC-2 in ZNA. The two tests measure different motor characteristics in children with suspected DCD and, thus, can be used complementary for the diagnosis of the disorder.

Regional variations in ABC transporter expression along the mouse intestinal tract.

The ATP-binding cassette (ABC) family of proteins comprise a group of membrane transporters involved in the transport of a wide variety of compounds, such as xenobiotics, vitamins, lipids, amino acids, and carbohydrates. Determining their regional expression patterns along the intestinal tract will further characterize their transport functions in the gut. The mRNA expression levels of murine ABC transporters in the duodenum, jejunum, ileum, and colon were examined using the Affymetrix MuU74v2 GeneChip set. Eight ABC transporters (Abcb2, Abcb3, Abcb9, Abcc3, Abcc6, Abcd1, Abcg5, and Abcg8) displayed significant differential gene expression along the intestinal tract, as determined by two statistical models (a global error assessment model and a classic ANOVA, both with a P < 0.01). Concordance with semiquantitative real-time PCR was high. Analyzing the promoters of the differentially expressed ABC transporters did not identify common transcriptional motifs between family members or with other genes; however, the expression profile for Abcb9 was highly correlated with fibulin-1, and both genes share a common complex promoter model involving the NFkappaB, zinc binding protein factor (ZBPF), GC-box factors SP1/GC (SP1F), and early growth response factor (EGRF) transcription binding motifs. The cellular location of another of the differentially expressed ABC transporters, Abcc3, was examined by immunohistochemistry. Staining revealed that the protein is consistently expressed in the basolateral compartment of enterocytes along the anterior-posterior axis of the intestine. Furthermore, the intensity of the staining pattern is concordant with the expression profile. This agrees with previous findings in which the mRNA, protein, and transport function of Abcc3 were increased in the rat distal intestine. These data reveal regional differences in gene expression profiles along the intestinal tract and demonstrate that a complete understanding of intestinal ABC transporter function can only be achieved by examining the physiologically distinct regions of the gut.

The half-size ABC transporters STR1 and STR2 are indispensable for mycorrhizal arbuscule formation in rice.

The central structure of the symbiotic association between plants and arbuscular mycorrhizal (AM) fungi is the fungal arbuscule that delivers minerals to the plant. Our earlier transcriptome analyses identified two half-size ABCG transporters that displayed enhanced mRNA levels in mycorrhizal roots. We now show specific transcript accumulation in arbusculated cells of both genes during symbiosis. Presently, arbuscule-relevant factors from monocotyledons have not been reported. Mutation of either of the Oryza sativa (rice) ABCG transporters blocked arbuscule growth of different AM fungi at a small and stunted stage, recapitulating the phenotype of Medicago truncatula stunted arbuscule 1 and 2 (str1 and str2) mutants that are deficient in homologous ABCG genes. This phenotypic resemblance and phylogenetic analysis suggest functional conservation of STR1 and STR2 across the angiosperms. Malnutrition of the fungus underlying limited arbuscular growth was excluded by the absence of complementation of the str1 phenotype by wild-type nurse plants. Furthermore, plant AM signaling was found to be intact, as arbuscule-induced marker transcript accumulation was not affected in str1 mutants. Strigolactones have previously been hypothesized to operate as intracellular hyphal branching signals and possible substrates of STR1 and STR2. However, full arbuscule development in the strigolactone biosynthesis mutants d10 and d17 suggested strigolactones to be unlikely substrates of STR1/STR2. Interestingly, rice STR1 is associated with a cis-natural antisense transcript (antiSTR1). Analogous to STR1 and STR2, at the root cortex level, the antiSTR1 transcript is specifically detected in arbusculated cells, suggesting unexpected modes of STR1 regulation in rice.