Quantifying single nucleotide variant detection sensitivity in exome sequencing

BACKGROUND: The targeted capture and sequencing of genomic regions has rapidly demonstrated its utility in genetic studies. Inherent in this technology is considerable heterogeneity of target coverage and this is expected to systematically impact our sensitivity to detect genuine polymorphisms. To fully interpret the polymorphisms identified in a genetic study it is often essential to both detect polymorphisms and to understand where and with what probability real polymorphisms may have been missed. RESULTS: Using down-sampling of 30 deeply sequenced exomes and a set of gold-standard single nucleotide variant (SNV) genotype calls for each sample, we developed an empirical model relating the read depth at a polymorphic site to the probability of calling the correct genotype at that site. We find that measured sensitivity in SNV detection is substantially worse than that predicted from the naive expectation of sampling from a binomial. This calibrated model allows us to produce single nucleotide resolution SNV sensitivity estimates which can be merged to give summary sensitivity measures for any arbitrary partition of the target sequences (nucleotide, exon, gene, pathway, exome). These metrics are directly comparable between platforms and can be combined between samples to give “power estimates” for an entire study. We estimate a local read depth of 13X is required to detect the alleles and genotype of a heterozygous SNV 95% of the time, but only 3X for a homozygous SNV. At a mean on-target read depth of 20X, commonly used for rare disease exome sequencing studies, we predict 5–15% of heterozygous and 1–4% of homozygous SNVs in the targeted regions will be missed. CONCLUSIONS: Non-reference alleles in the heterozygote state have a high chance of being missed when commonly applied read coverage thresholds are used despite the widely held assumption that there is good polymorphism detection at these coverage levels. Such alleles are likely to be of functional importance in population based studies of rare diseases, somatic mutations in cancer and explaining the “missing heritability” of quantitative traits

Measurement of D+- and D0 production in deep inelastic scattering using a lifetime tag at HERA

The production of D-+/-- and D-0-mesons has been measured with the ZEUS detector at HERA using an integrated luminosity of 133.6 pb(-1). The measurements cover the kinematic range 5 < Q(2) < 1000 GeV2, 0.02 < y < 0.7, 1.5 < p(T)(D) < 15 GeV and |eta(D)| < 1.6. Combinatorial background to the D-meson signals is reduced by using the ZEUS microvertex detector to reconstruct displaced secondary vertices. Production cross sections are compared with the predictions of next-to-leading-order QCD, which is found to describe the data well. Measurements are extrapolated to the full kinematic phase space in order to obtain the open-charm contribution, F-2(c (c) over bar), to the proton structure function, F-2

Cytogenetic characterization and genome size of the medicinal plant Catharanthus roseus (L.) G. Don

The genome size and organization of the important medicinal plant Catharanthus roseus is shown to correspond to 1C = 0.76 pg (~738 Mbps) and 2n = 16 chromosomes. The data provide a sound basis for future studies including cytogenetic mapping, genomics and breeding

Elliptic flow of charged particles in Pb-Pb collisions at 2.76 TeV

We report the first measurement of charged particle elliptic flow in Pb-Pb collisions at 2.76 TeV with the ALICE detector at the CERN Large Hadron Collider. The measurement is performed in the central pseudorapidity region (|$\eta$|<0.8) and transverse momentum range 0.2< $p_{\rm T}$< 5.0 GeV/$c$. The elliptic flow signal v$_2$, measured using the 4-particle correlation method, averaged over transverse momentum and pseudorapidity is 0.087 $\pm$ 0.002 (stat) $\pm$ 0.004 (syst) in the 40-50% centrality class. The differential elliptic flow v$_2(p_{\rm T})$ reaches a maximum of 0.2 near $p_{\rm T}$ = 3 GeV/$c$. Compared to RHIC Au-Au collisions at 200 GeV, the elliptic flow increases by about 30%. Some hydrodynamic model predictions which include viscous corrections are in agreement with the observed increase.Comment: 10 pages, 4 captioned figures, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/389