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

The secret world of shrimps: polarisation vision at its best

Animal vision spans a great range of complexity, with systems evolving to detect variations in optical intensity, distribution, colour, and polarisation. Polarisation vision systems studied to date detect one to four channels of linear polarisation, combining them in opponent pairs to provide intensity-independent operation. Circular polarisation vision has never been seen, and is widely believed to play no part in animal vision. Polarisation is fully measured via Stokes' parameters--obtained by combined linear and circular polarisation measurements. Optimal polarisation vision is the ability to see Stokes' parameters: here we show that the crustacean \emph{Gonodactylus smithii} measures the exact components required. This vision provides optimal contrast-enhancement, and precise determination of polarisation with no confusion-states or neutral-points--significant advantages. We emphasise that linear and circular polarisation vision are not different modalities--both are necessary for optimal polarisation vision, regardless of the presence of strongly linear or circularly polarised features in the animal's environment.Comment: 10 pages, 6 figures, 2 table

Chromatin establishes an immature version of neuronal protocadherin selection during the naive-to-primed conversion of pluripotent stem cells.

In the mammalian genome, the clustered protocadherin (cPCDH) locus provides a paradigm for stochastic gene expression with the potential to generate a unique cPCDH combination in every neuron. Here we report a chromatin-based mechanism that emerges during the transition from the naive to the primed states of cell pluripotency and reduces, by orders of magnitude, the combinatorial potential in the human cPCDH locus. This mechanism selectively increases the frequency of stochastic selection of a small subset of cPCDH genes after neuronal differentiation in monolayers, 10-month-old cortical organoids and engrafted cells in the spinal cords of rats. Signs of these frequent selections can be observed in the brain throughout fetal development and disappear after birth, except in conditions of delayed maturation such as Down's syndrome. We therefore propose that a pattern of limited cPCDH-gene expression diversity is maintained while human neurons still retain fetal-like levels of maturation

Recruitment of ethnic minority patients to a cardiac rehabilitation trial: The Birmingham Rehabilitation Uptake Maximisation (BRUM) study [ISRCTN72884263]

Background: Concerns have been raised about low participation rates of people from minority ethnic groups in clinical trials. However, the evidence is unclear as many studies do not report the ethnicity of participants and there is insufficient information about the reasons for ineligibility by ethnic group. Where there are data, there remains the key question as to whether ethnic minorities more likely to be ineligible (e.g. due to language) or decline to participate. We have addressed these questions in relation to the Birmingham Rehabilitation Uptake Maximisation (BRUM) study, a randomized controlled trial (RCT) comparing a home-based with a hospital-based cardiac rehabilitation programme in a multi-ethnic population in the UK. Methods: Analysis of the ethnicity, age and sex of presenting and recruited subjects for a trial of cardiac rehabilitation in the West-Midlands, UK. Participants: 1997 patients presenting post-myocardial infarction, percutaneous transluminal coronary angioplasty or coronary artery bypass graft surgery. Data collected: exclusion rates, reasons for exclusion and reasons for declining to participate in the trial by ethnic group. Results: Significantly more patients of South Asian ethnicity were excluded (52% of 'South Asian' v 36% 'White European' and 36% 'Other', p < 0.001). This difference in eligibility was primarily due to exclusion on the basis of language (i.e. the inability to speak English or Punjabi). Of those eligible, similar proportions were recruited from the different ethnic groups (white, South Asian and other). There was a marked difference in eligibility between people of Indian, Pakistani or Bangladeshi origin

Automatically extracting functionally equivalent proteins from SwissProt

In summary, FOSTA provides an automated analysis of annotations in UniProtKB/Swiss-Prot to enable groups of proteins already annotated as functionally equivalent, to be extracted. Our results demonstrate that the vast majority of UniProtKB/Swiss-Prot functional annotations are of high quality, and that FOSTA can interpret annotations successfully. Where FOSTA is not successful, we are able to highlight inconsistencies in UniProtKB/Swiss-Prot annotation. Most of these would have presented equal difficulties for manual interpretation of annotations. We discuss limitations and possible future extensions to FOSTA, and recommend changes to the UniProtKB/Swiss-Prot format, which would facilitate text-mining of UniProtKB/Swiss-Prot

Eigen-transitions in cantilever cylindrical shells subjected to vertical edge loads

A thin cantilever cylindrical shell subjected to a transverse shear force at the free end can experience two distinct modes of buckling, depending on its relative thickness and length. If the former parameter is fixed then a short cylinder buckles in a diffuse manner, while the eigenmodal deformation of a moderately long shell is localised, both axially and circumferentially, near its fixed end. Donnelltype buckling equations for cylindrical shells are here coupled with a non-symmetric membrane basic state to produce a linear boundary-value problem that is shown to capture the transition between the aforementioned instability modes. The main interest lies in exploring the approximate asymptotic separation of the independent variables in the corresponding stability equations, when the eigen-deformation is doubly localised. Comparisons with direct numerical simulations of the full buckling problem provide further insight into the accuracy and limitations of our approximations

Dietary garlic and hip osteoarthritis: evidence of a protective effect and putative mechanism of action

Background Patterns of food intake and prevalent osteoarthritis of the hand, hip, and knee were studied using the twin design to limit the effect of confounding factors. Compounds found in associated food groups were further studied in vitro. Methods Cross-sectional study conducted in a large population-based volunteer cohort of twins. Food intake was evaluated using the Food Frequency Questionnaire; OA was determined using plain radiographs. Analyses were adjusted for age, BMI and physical activity. Subsequent in vitro studies examined the effects of allium-derived compounds on the expression of matrix-degrading proteases in SW1353 chondrosarcoma cells. Results Data were available, depending on phenotype, for 654-1082 of 1086 female twins (median age 58.9 years; range 46-77). Trends in dietary analysis revealed a specific pattern of dietary intake, that high in fruit and vegetables, showed an inverse association with hip OA (p = 0.022). Consumption of 'non-citrus fruit' (p = 0.015) and 'alliums' (p = 0.029) had the strongest protective effect. Alliums contain diallyl disulphide which was shown to abrogate cytokine-induced matrix metalloproteinase expression. Conclusions Studies of diet are notorious for their confounding by lifestyle effects. While taking account of BMI, the data show an independent effect of a diet high in fruit and vegetables, suggesting it to be protective against radiographic hip OA. Furthermore, diallyl disulphide, a compound found in garlic and other alliums, represses the expression of matrix-degrading proteases in chondrocyte-like cells, providing a potential mechanism of action

Systematic and detailed analysis of behavioural tests in the rat Middle Cerebral Artery Occlusion (MCAO) model of stroke: tests for long-term assessment

In order to test therapeutics, functional assessments are required. In pre-clinical stroke research, there is little consensus regarding the most appropriate behavioural tasks to assess deficits; especially when testing over extended times in milder models with short occlusion times and small lesion volumes. In this study we comprehensively assessed 16 different behavioural tests, with the aim of identifying those that show robust, reliable and stable deficits for up to 2 months. These tasks are regularly used in stroke research, as well as being useful for examining striatal dysfunction in models of Huntington’s and Parkinson’s disease. Two cohorts of male Wistar rats underwent the intraluminal filament model of MCAO (30min) and were imaged 24hrs later. This resulted in primarily subcortical infarcts, with a small amount of cortical damage. Animals were tested, along with sham and naïve groups at 24hrs, 7 days, and 1 and 2 months. Following behavioural testing, brains were processed and striatal neuronal counts were performed alongside measurements of total brain and white matter atrophy. The staircase, adjusting steps, rotarod and apomorphine induced rotations were the most reliable for assessing long-term deficits in the 30 min transient MCAO model of stroke

Middleborns disadvantaged? testing birth-order effects on fitness in pre-industrial finns

Parental investment is a limited resource for which offspring compete in order to increase their own survival and reproductive success. However, parents might be selected to influence the outcome of sibling competition through differential investment. While evidence for this is widespread in egg-laying species, whether or not this may also be the case in viviparous species is more difficult to determine. We use pre-industrial Finns as our model system and an equal investment model as our null hypothesis, which predicts that (all else being equal) middleborns should be disadvantaged through competition. We found no overall evidence to suggest that middleborns in a family are disadvantaged in terms of their survival, age at first reproduction or lifetime reproductive success. However, when considering birth-order only among same-sexed siblings, first-, middle-and lastborn sons significantly differed in the number of offspring they were able to rear to adulthood, although there was no similar effect among females. Middleborn sons appeared to produce significantly less offspring than first-or lastborn sons, but they did not significantly differ from lastborn sons in the number of offspring reared to adulthood. Our results thus show that taking sex differences into account is important when modelling birth-order effects. We found clear evidence of firstborn sons being advantaged over other sons in the family, and over firstborn daughters. Therefore, our results suggest that parents invest differentially in their offspring in order to both preferentially favour particular offspring or reduce offspring inequalities arising from sibling competition