CNVassoc: Association analysis of CNV data using R

Background: Copy number variants (CNV) are a potentially important component of the genetic contribution to risk of common complex diseases. Analysis of the association between CNVs and disease requires that uncertainty in CNV copy-number calls, which can be substantial, be taken into account; failure to consider this uncertainty can lead to biased results. Therefore, there is a need to develop and use appropriate statistical tools. To address this issue, we have developed CNVassoc, an R package for carrying out association analysis of common copy number variants in population-based studies. This package includes functions for testing for association with different classes of response variables (e.g. class status, censored data, counts) under a series of study designs (case-control, cohort, etc) and inheritance models, adjusting for covariates. The package includes functions for inferring copy number (CNV genotype calling), but can also accept copy number data generated by other algorithms (e.g. CANARY, CGHcall, IMPUTE). Results: Here we present a new R package, CNVassoc, that can deal with different types of CNV arising from different platforms such as MLPA o aCGH. Through a real data example we illustrate that our method is able to incorporate uncertainty in the association process. We also show how our package can also be useful when analyzing imputed data when analyzing imputed SNPs. Through a simulation study we show that CNVassoc outperforms CNVtools in terms of computing time as well as in convergence failure rate. Conclusions: We provide a package that outperforms the existing ones in terms of modelling flexibility, power, convergence rate, ease of covariate adjustment, and requirements for sample size and signal quality. Therefore, we offer CNVassoc as a method for routine use in CNV association studiesThis work has been supported by the Spanish Ministry of Science and Innovation (MTM2008-02457 to JRG, BIO2009-12458 to RD-U and statistical genetics network MTM2010-09526-E (subprograma MTM) to JRG, IS, GL and RD-U). GL is supported by the Juan de la Cierva Program of the Spanish Ministry of Science and Innovation

Molecular Analysis of Precursor Lesions in Familial Pancreatic Cancer

PMCID: PMC3553106This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Measurement and Interpretation of Fermion-Pair Production at LEP energies above the Z Resonance

This paper presents DELPHI measurements and interpretations of cross-sections, forward-backward asymmetries, and angular distributions, for the e+e- -> ffbar process for centre-of-mass energies above the Z resonance, from sqrt(s) ~ 130 - 207 GeV at the LEP collider. The measurements are consistent with the predictions of the Standard Model and are used to study a variety of models including the S-Matrix ansatz for e+e- -> ffbar scattering and several models which include physics beyond the Standard Model: the exchange of Z' bosons, contact interactions between fermions, the exchange of gravitons in large extra dimensions and the exchange of sneutrino in R-parity violating supersymmetry.Comment: 79 pages, 16 figures, Accepted by Eur. Phys. J.

A Determination of the Centre-of-Mass Energy at LEP2 using Radiative 2-fermion Events

Using e+e- -> mu+mu-(gamma) and e+e- -> qqbar(gamma) events radiative to the Z pole, DELPHI has determined the centre-of-mass energy, sqrt{s}, using energy and momentum constraint methods. The results are expressed as deviations from the nominal LEP centre-of-mass energy, measured using other techniques. The results are found to be compatible with the LEP Energy Working Group estimates for a combination of the 1997 to 2000 data sets.Comment: 20 pages, 6 figures, Accepted by Eur. Phys. J.

A Measurement of the Tau Hadronic Branching Ratios

The exclusive and semi-exclusive branching ratios of the tau lepton hadronic decay modes (h- v_t, h- pi0 v_t, h- pi0 pi0 v_t, h- \geq 2pi0 v_t, h- \geq 3pi0 v_t, 2h- h+ v_t, 2h- h+ pi0 v_t, 2h- h+ \geq 2pi0 v_t, 3h- 2h+ v_t and 3h- 2h+ \geq 1pi0 v_t) were measured with data from the DELPHI detector at LEP.Comment: 53 pages, 18 figures, Accepted by Eur. Phys. J.

Search for supersymmetric particles in scenarios with a gravitino LSP and stau NLSP

Sleptons, neutralinos and charginos were searched for in the context of scenarios where the lightest supersymmetric particle is the gravitino. It was assumed that the stau is the next-to-lightest supersymmetric particle. Data collected with the DELPHI detector at a centre-of-mass energy near 189 GeV were analysed combining the methods developed in previous searches at lower energies. No evidence for the production of these supersymmetric particles was found. Hence, limits were derived at 95% confidence level.Comment: 31 pages, 14 figure

Haplotype block structure study of the CFTR gene. Most variants are associated with the M470 allele in several European populations

An average of about 1700 CFTR (cystic fibrosis transmembrane conductance regulator) alleles from normal individuals from different European populations were extensively screened for DNA sequence variation. A total of 80 variants were observed: 61 coding SNSs (results already published), 13 noncoding SNSs, three STRs, two short deletions, and one nucleotide insertion. Eight DNA variants were classified as non-CF causing due to their high frequency of occurrence. Through this survey the CFTR has become the most exhaustively studied gene for its coding sequence variability and, though to a lesser extent, for its noncoding sequence variability as well. Interestingly, most variation was associated with the M470 allele, while the V470 allele showed an 'extended haplotype homozygosity' (EHH). These findings make us suggest a role for selection acting either on the M470V itself or through an hitchhiking mechanism involving a second site. The possible ancient origin of the V allele in an 'out of Africa' time frame is discussed

Loss of Ribosomal Protein L11 Affects Zebrafish Embryonic Development through a p53-Dependent Apoptotic Response

Ribosome is responsible for protein synthesis in all organisms and ribosomal proteins (RPs) play important roles in the formation of a functional ribosome. L11 was recently shown to regulate p53 activity through a direct binding with MDM2 and abrogating the MDM2-induced p53 degradation in response to ribosomal stress. However, the studies were performed in cell lines and the significance of this tumor suppressor function of L11 has yet to be explored in animal models. To investigate the effects of the deletion of L11 and its physiological relevance to p53 activity, we knocked down the rpl11 gene in zebrafish and analyzed the p53 response. Contrary to the cell line-based results, our data indicate that an L11 deficiency in a model organism activates the p53 pathway. The L11-deficient embryos (morphants) displayed developmental abnormalities primarily in the brain, leading to embryonic lethality within 6–7 days post fertilization. Extensive apoptosis was observed in the head region of the morphants, thus correlating the morphological defects with apparent cell death. A decrease in total abundance of genes involved in neural patterning of the brain was observed in the morphants, suggesting a reduction in neural progenitor cells. Upregulation of the genes involved in the p53 pathway were observed in the morphants. Simultaneous knockdown of the p53 gene rescued the developmental defects and apoptosis in the morphants. These results suggest that ribosomal dysfunction due to the loss of L11 activates a p53-dependent checkpoint response to prevent improper embryonic development

Infrared Spectroscopic Studies of Cells and Tissues: Triple Helix Proteins as a Potential Biomarker for Tumors

In this work, the infrared (IR) spectra of living neural cells in suspension, native brain tissue, and native brain tumor tissue were investigated. Methods were developed to overcome the strong IR signal of liquid water so that the signal from the cellular biochemicals could be seen. Measurements could be performed during surgeries, within minutes after resection. Comparison between normal tissue, different cell lineages in suspension, and tumors allowed preliminary assignments of IR bands to be made. The most dramatic difference between tissues and cells was found to be in weaker IR absorbances usually assigned to the triple helix of collagens. Triple helix domains are common in larger structural proteins, and are typically found in the extracellular matrix (ECM) of tissues. An algorithm to correct offsets and calculate the band heights and positions of these bands was developed, so the variance between identical measurements could be assessed. The initial results indicate the triple helix signal is surprisingly consistent between different individuals, and is altered in tumor tissues. Taken together, these preliminary investigations indicate this triple helix signal may be a reliable biomarker for a tumor-like microenvironment. Thus, this signal has potential to aid in the intra-operational delineation of brain tumor borders. © 2013 Stelling et al