Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

Mynodbcsv: Lightweight Zero-Config Database Solution for Handling Very Large CSV Files.

Volumes of data used in science and industry are growing rapidly. When researchers face the challenge of analyzing them, their format is often the first obstacle. Lack of standardized ways of exploring different data layouts requires an effort each time to solve the problem from scratch. Possibility to access data in a rich, uniform manner, e.g. using Structured Query Language (SQL) would offer expressiveness and user-friendliness. Comma-separated values (CSV) are one of the most common data storage formats. Despite its simplicity, with growing file size handling it becomes non-trivial. Importing CSVs into existing databases is time-consuming and troublesome, or even impossible if its horizontal dimension reaches thousands of columns. Most databases are optimized for handling large number of rows rather than columns, therefore, performance for datasets with non-typical layouts is often unacceptable. Other challenges include schema creation, updates and repeated data imports. To address the above-mentioned problems, I present a system for accessing very large CSV-based datasets by means of SQL. It's characterized by: "no copy" approach - data stay mostly in the CSV files; "zero configuration" - no need to specify database schema; written in C++, with boost [1], SQLite [2] and Qt [3], doesn't require installation and has very small size; query rewriting, dynamic creation of indices for appropriate columns and static data retrieval directly from CSV files ensure efficient plan execution; effortless support for millions of columns; due to per-value typing, using mixed text/numbers data is easy; very simple network protocol provides efficient interface for MATLAB and reduces implementation time for other languages. The software is available as freeware along with educational videos on its website [4]. It doesn't need any prerequisites to run, as all of the libraries are included in the distribution package. I test it against existing database solutions using a battery of benchmarks and discuss the results.

Addressing trypsin bias in large scale (phospho)proteome analysis by size exclusion chromatography and secondary digestion of large post-trypsin peptides.

In the vast majority of bottom-up proteomics studies, protein digestion is performed using only mammalian trypsin. Although it is clearly the best enzyme available, the sole use of trypsin rarely leads to complete sequence coverage, even for abundant proteins. It is commonly assumed that this is because many tryptic peptides are either too short or too long to be identified by RPLC-MS/MS. We show through in silico analysis that 20-30% of the total sequence of three proteomes (Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Homo sapiens) is expected to be covered by Large post-Trypsin Peptides (LpTPs) with M(r) above 3000 Da. We then established size exclusion chromatography to fractionate complex yeast tryptic digests into pools of peptides based on size. We found that secondary digestion of LpTPs followed by LC-MS/MS analysis leads to a significant increase in identified proteins and a 32-50% relative increase in average sequence coverage compared to trypsin digestion alone. Application of the developed strategy to analyze the phosphoproteomes of S. pombe and of a human cell line identified a significant fraction of novel phosphosites. Overall our data indicate that specific targeting of LpTPs can complement standard bottom-up workflows to reveal a largely neglected portion of the proteome.

Expression of MYC, IgM, As Well As Non-Germinal Centre B-Cell Like Immunophenotype and Positive Immunofish Index Predict a Worse Progression Free Survival and Overall Survival in a Series of 670 De Novo Diffuse Large B-Cell Lymphomas Included in Clinical Trials: A GELA Study of the 2003 Program

Introduction: Diffuse large B-cell lymphomas (DLBCL) represent a heterogeneous disease with variable clinical outcome. Identifying phenotypic biomarkers of tumor cells on paraffin sections that predict different clinical outcome remain an important goal that may also help to better understand the biology of this lymphoma. Differentiating non-germinal centre B-cell-like (non-GCB) from Germinal Centre B-cell-like (GCB) DLBCL according to Hans algorithm has been considered as an important immunohistochemical biomarker with prognostic value among patients treated with R-CHOP although not reproducibly found by all groups. Gene expression studies have also shown that IgM expression might be used as a surrogate for the GCB and ABC subtypes with a strong preferential expression of IgM in ABC DLBCL subtype. ImmunoFISH index based on the differential expression of MUM-1, FOXP1 by immunohistochemistry and on the BCL6 rearrangement by FISH has been previously reported (C Copie-Bergman, J Clin Oncol. 2009;27:5573-9) as prognostic in an homogeneous series of DLBCL treated with R-CHOP. In addition, oncogenic MYC protein overexpression by immunohistochemistry may represent an easy tool to identify the consequences of MYC deregulation in DLBCL. Our aim was to analyse by immunohistochemistry the prognostic relevance of MYC, IgM, GCB/nonGCB subtype and ImmunoFISH index in a large series of de novo DLBCL treated with Rituximab (R)-chemotherapy (anthracyclin based) included in the 2003 program of the Groupe d'Etude des Lymphomes de l'Adulte (GELA) trials. Methods: The 2003 program included patients with de novo CD20+ DLBCL enrolled in 6 different LNH-03 GELA trials (LNH-03-1B, -B, -3B, 39B, -6B, 7B) stratifying patients according to age and age-adjusted IPI. Tumor samples were analyzed by immunohistochemistry using CD10, BCL6, MUM1, FOXP1 (according to Barrans threshold), MYC, IgM antibodies on tissue microarrays and by FISH using BCL6 split signal DNA probes. Considering evaluable Hans score, 670 patients were included in the study with 237 (35.4%) receiving intensive R-ACVBP regimen and 433 (64.6%) R-CHOP/R-mini-CHOP. Results: 304 (45.4%) DLBCL were classified as GCB and 366 (54.6%) as non-GCB according to Hans algorithm. 337/567 cases (59.4%) were positive for the ImmunoFISH index (i.e. two out of the three markers positive: MUM1 protein positive, FOXP1 protein Variable or Strong, BCL6 rearrangement). Immunofish index was preferentially positive in the non-GCB subtype (81.3%) compared to the GCB subtype (31.2%), (p<0.001). IgM was recorded as positive in tumor cells in 351/637 (52.4%) DLBCL cases with a preferential expression in non-GCB 195 (53.3%) vs GCB subtype 100(32.9%), p<0.001). MYC was positive in 170/577 (29.5%) cases with a 40% cut-off and in 44/577 (14.2%) cases with a cut-off of 70%. There was no preferential expression of MYC among GCB or non-GCB subtype (p>0.4) for both cut-offs. Progression-free Survival (PFS) was significantly worse among patients with high IPI score (p<0.0001), IgM positive tumor (p<0.0001), MYC positive tumor with a 40% threshold (p<0.001), ImmunoFISH positive index (p<0.002), non-GCB DLBCL subtype (p<0.0001). Overall Survival (OS) was also significantly worse among patients with high IPI score (p<0.0001), IgM positive tumor (p=0.02), MYC positive tumor with a 40% threshold (p<0.01), ImmunoFISH positive index (p=0.02), non-GCB DLBCL subtype (p<0.0001). All significant parameters were included in a multivariate analysis using Cox Model and in addition to IPI, only the GCB/non-GCB subtype according to Hans algorithm predicted significantly a worse PFS among non-GCB subgroup (HR 1.9 [1.3-2.8] p=0.002) as well as a worse OS (HR 2.0 [1.3-3.2], p=0.003). This strong prognostic value of non-GCB subtyping was confirmed considering only patients treated with R- CHOP for PFS (HR 2.1 [1.4-3.3], p=0.001) and for OS (HR 2.3 [1.3-3.8], p=0.002). Conclusion: Our study on a large series of patients included in trials confirmed the relevance of immunohistochemistry as a useful tool to identify significant prognostic biomarkers for clinical use. We show here that IgM and MYC might be useful prognostic biomarkers. In addition, we confirmed in this series the prognostic value of the ImmunoFISH index. Above all, we fully validated the strong and independent prognostic value of the Hans algorithm, daily used by the pathologists to subtype DLBCL.

Iron Supplementation Therapy in a Large Swiss Cohort of Inflammatory Bowel Disease Patients: Shift from Oral to Intravenous Therapy over Time

Background a nd A ims: T he 2 007 ECCO g uidelines o nanemia in inflammatory bowel disease (IBD) favour intravenous(iv) over oral (po) i ron supplementation due to bettereffectiveness and tolerance. Application of guidelines in clinicalpractice m ay r equire time. We a imed to determine thepercentage of IBD patients under iron supplementation therapyand its application mode over time in a large IBD cohort.Methods: Helsana, a leading Swiss health insurance companyprovides c overage f or approximately 18% of t he Swisspopulation, corresponding to about 1.2 million enrollees.Patients with Crohn's disease (CD) and ulcerative colitis (UC)were identified b y keyword search from t he a nonymisedHelsana database.Results: I n total, 6 29 CD ( 61% female) a nd 4 03 UC ( 56%female) patients w ere identified, mean retrospectiveobservation time w as 2 0.4 m onths f or CD and 13 m onths f orUC patients. Of t he entire study population, 29.3% wereprescribed iron. O ccurrence of iron prescription was 21.3% inmales a nd 31.2% in f emales ( odds r atio [OR] 1 .69, 95%-confidence interval [CI] 1.26-2.28). The prescription of iv i ronincreased from 2006/2007 ( 48.8% w ith iv i ron) to 2 008/2009(65.2% with iv iron) by a factor of 1.89.Conclusions: One third of the IBD population was treated withiron supplementation. A gradual s hift from oral t o iv iron wasobserved over time in a large Swiss IBD cohort. This switch inprescription habits g oes a long with the implementation of theECCO consensus guidelines on anemia in IBD.

Variation in thermal performance and reaction norms among populations of Drosophila melanogaster.

The major goal of evolutionary thermal biology is to understand how variation in temperature shapes phenotypic evolution. Comparing thermal reaction norms among populations from different thermal environments allows us to gain insights into the evolutionary mechanisms underlying thermal adaptation. Here, we have examined thermal adaptation in six wild populations of the fruit fly (Drosophila melanogaster) from markedly different natural environments by analyzing thermal reaction norms for fecundity, thorax length, wing area, and ovariole number under ecologically realistic fluctuating temperature regimes in the laboratory. Contrary to expectation, we found only minor differences in the thermal optima for fecundity among populations. Differentiation among populations was mainly due to differences in absolute (and partly also relative) thermal fecundity performance. Despite significant variation among populations in the absolute values of morphological traits, we observed only minor differentiation in their reaction norms. Overall, the thermal reaction norms for all traits examined were remarkably similar among different populations. Our results therefore suggest that thermal adaptation in D. melanogaster predominantly involves evolutionary changes in absolute trait values rather than in aspects of thermal reaction norms.

Reactivity of antiglioma monoclonal antibodies for a large panel of cultured gliomas and other neuroectoderm derived tumors.

We produced three monoclonal antibodies, BF7, GE2 and CG12, against cultured human glioma cells. Their specificity was tested by an indirect antibody-binding radioimmunoassay on a panel of glial and non-glial tumor cell lines. BF7 and GE2 react preferentially with glioma cells and, except for one colon carcinoma line, they do not bind to the control non-neuroectodermal cells; they appear to be directed against common malignant glioma associated antigens. CG12, the third monoclonal antibody, binds to the great majority of tumor cell lines of neuroectodermal origin and does not bind to any other cell lines tested.

Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways.

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have increased the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin concentration showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional analysis of these newly discovered loci will further improve our understanding of glycemic control.

Detection of HIV drug resistance during antiretroviral treatment and clinical progression in a large European cohort study.

OBJECTIVE(S): To investigate the relationship between detection of HIV drug resistance by 2 years from starting antiretroviral therapy and the subsequent risk of progression to AIDS and death. DESIGN: Virological failure was defined as experiencing two consecutive viral loads of more than 400 copies/ml in the time window between 0.5 and 2 years from starting antiretroviral therapy (baseline). Patients were grouped according to evidence of virological failure and whether there was detection of the International AIDS Society resistance mutations to one, two or three drug classes in the time window. METHODS: Standard survival analysis using Kaplan-Meier curves and Cox proportional hazards regression model with time-fixed covariates defined at baseline was employed. RESULTS: We studied 8229 patients in EuroSIDA who started antiretroviral therapy and who had at least 2 years of clinical follow-up. We observed 829 AIDS events and 571 deaths during 38,814 person-years of follow-up resulting in an overall incidence of new AIDS and death of 3.6 per 100 person-years of follow-up [95% confidence interval (CI):3.4-3.8]. By 96 months from baseline, the proportion of patients with a new AIDS diagnosis or death was 20.3% (95% CI:17.7-22.9) in patients with no evidence of virological failure and 53% (39.3-66.7) in those with virological failure and mutations to three drug classes (P = 0.0001). An almost two-fold difference in risk was confirmed in the multivariable analysis (adjusted relative hazard = 1.8, 95% CI:1.2-2.7, P = 0.005). CONCLUSION: Although this study shows an association between the detection of resistance at failure and risk of clinical progression, further research is needed to clarify whether resistance reflects poor adherence or directly increases the risk of clinical events via exhaustion of drug options.

Iron supplementation in a large Swiss Cohort of Inflammatory Bowel Disease patients demonstrates a shift form oral to intravenous therapy over time

Background and Aims: The 2007 European Crohn's and Colitis Organization guidelines on anemia in inflammatory bowel disease (IBD) favour intravenous (iv) over oral (po) iron supplementation due to better effectiveness and tolerance. We aimed to determine the percentage of IBD patients under iron supplementation therapy and the dynamics of prescription habits (iv versus po) over time. Methods: Helsana, a leading Swiss health insurance company provides coverage for approximately 18% of the Swiss population, corresponding to about 1.2 million enrollees. Patients with Crohn's disease (CD) and ulcerative colitis (UC) were analyzed from the anonymised Helsana database. Results: In total, 629 CD (61% female) and 398 UC (57% female) patients were identified, mean observation time was 31.8 months for CD and 31.0 months for UC patients. Of the entire study population, 27.1% were prescribed iron (21.1% in males and 31.1% in females). Patients treated with IBD-specific drugs (steroids, immunomodulators, anti-TNF agents) were more frequently treated with iron compared to patients without any medication (35.0% vs. 20.9%, OR 1.91, 95%-CI 1.41-2.61). The prescription of iv iron increased from 2006/2007 (48.8% of all patients receiving any iron priscription) to 65.2% in 2008/2009 by a factor of 1.89. Conclusions: One third of the IBD population was treated with iron supplementation. A gradual shift from oral to iv iron was observed over time. This switch in prescription habits goes along with the implementation of the ECCO consensus guidelines on anemia in IBD.

Geochemistry of the thermal springs and fumaroles of Basse-Terre Island, Guadeloupe, Lesser Antilles

The purpose of this work was to study jointly the volcanic-hydrothermal system of the high-risk volcano La Soufriere, in the southern part of Basse-Terre, and the geothermal area of Bouillante, on its western coast, to derive an all-embracing and coherent conceptual geochemical model that provides the necessary basis for adequate volcanic surveillance and further geothermal exploration. The active andesitic dome of La Soufriere has erupted eight times since 1660, most recently in 1976-1977. All these historic eruptions have been phreatic. High-salinity, Na-CI geothermal liquids circulate in the Bouillante geothermal reservoir, at temperatures close to 250 degrees C. These Na-CI solutions rise toward the surface, undergo boiling and mixing with groundwater and/or seawater, and feed most Na-CI thermal springs in the central Bouillante area. The Na-Cl thermal springs are surrounded by Na-HCO3 thermal springs and by the Na-Cl thermal spring of Anse a la Barque (a groundwater slightly mixed with seawater), which are all heated through conductive transfer. The two main fumarolic fields of La Soufriere area discharge vapors formed through boiling of hydrothermal aqueous solutions at temperatures of 190-215 degrees C below the ``Ty'' fault area and close to 260 degrees C below the dome summit. The boiling liquid producing the vapors of the Ty fault area has SD and delta(18)O values relatively similar to those of the Na-CI liquids of the Bouillante geothermal reservoir, whereas the liquid originating the vapors of the summit fumaroles is strongly enriched in O-18, due to input of magmatic fluids from below. This process is also responsible for the paucity of CH;I in the fumaroles. The thermal features around La Soufriere dome include: (a) Ca-SO4 springs, produced through absorption of hydrothermal vapors in shallow groundwaters; (b) conductively heated, Ca-Na-HCO3 springs; and (c) two Ca-Na-Cl springs produced through mixing of shallow Ca-SO4 waters and deep Na-Cl hydrothermal liquids. The geographical distribution of the different thermal features of La Soufriere area indicates the presence of: (a) a central zone dominated by the ascent of steam, which either discharges at the surface in the fumarolic fields or is absorbed in shallow groundwaters; and (b) an outer zone, where the shallow groundwaters are heated through conduction or addition of Na-Cl liquids coming from hydrothermal aquifer(s).

Recent human-induced trophic change in the large and deep perialpine Lake Lucerne (Switzerland) compared to historical geochemical variations

This study investigates the sedimentological and geochemical changes that occurred during the last 2200 years in the meromictic Lake Lucerne (Switzerland), one of the largest freshwater lakes of Central Europe. The stable isotope composition (delta C-13 and delta O-18 values) of bulk carbonates is compared to changes in grain-size distribution (clay and silt fraction), natural trace element input (titanium and thorium concentrations), and organic material abundance (C-org, nitrogen and phosphorus) and composition (C/N ratios and hydrogen and oxygen indexes). A drop in carbonate accumulation and in the delta O-18 values of sediments between ca. AD 500 and 700 followed a large and consistent rise in chemical weathering, marked by increases in the silicate-clay fraction and in crustal element concentrations. During the following millennium, there was a long-term decreasing trend in the lithogenic trace element input and in the phosphorus loading, suggesting decreasing terrigeneous input from runoff water. The major sedimentological change over the studied period occurred after ca. AD 1800 with a significant increase in the erosion-driven silt-fraction and in the sedimentation rate. During the last century, human-induced increase in nutrient input to the lake highly enhanced the accumulation of organic matter in sediment. Changes in nutrients and oxygen conditions in the hypolimnion of Lake Lucerne during the eutrophication period (i.e., the last 40 years) highly modified the geochemical fluxes compared to the relatively stable oligotrophic conditions that prevailed during the previous 2000 years. Before the 19th century, climate driven meromixis had a limited impact on the organic matter flux to the sediments, but the accumulation of carbonate considerably decreased during periods of lower mechanical erosion rates and high chemical weathering rates. (C) 2012 Elsevier B.V. All rights reserved.