A Fast Algorithm for Robust Regression with Penalised Trimmed Squares

The presence of groups containing high leverage outliers makes linear regression a difficult problem due to the masking effect. The available high breakdown estimators based on Least Trimmed Squares often do not succeed in detecting masked high leverage outliers in finite samples. An alternative to the LTS estimator, called Penalised Trimmed Squares (PTS) estimator, was introduced by the authors in \cite{ZiouAv:05,ZiAvPi:07} and it appears to be less sensitive to the masking problem. This estimator is defined by a Quadratic Mixed Integer Programming (QMIP) problem, where in the objective function a penalty cost for each observation is included which serves as an upper bound on the residual error for any feasible regression line. Since the PTS does not require presetting the number of outliers to delete from the data set, it has better efficiency with respect to other estimators. However, due to the high computational complexity of the resulting QMIP problem, exact solutions for moderately large regression problems is infeasible. In this paper we further establish the theoretical properties of the PTS estimator, such as high breakdown and efficiency, and propose an approximate algorithm called Fast-PTS to compute the PTS estimator for large data sets efficiently. Extensive computational experiments on sets of benchmark instances with varying degrees of outlier contamination, indicate that the proposed algorithm performs well in identifying groups of high leverage outliers in reasonable computational time.Comment: 27 page

Clonality and evolutionary history of rhabdomyosarcoma.

To infer the subclonality of rhabdomyosarcoma (RMS) and predict the temporal order of genetic events for the tumorigenic process, and to identify novel drivers, we applied a systematic method that takes into account germline and somatic alterations in 44 tumor-normal RMS pairs using deep whole-genome sequencing. Intriguingly, we find that loss of heterozygosity of 11p15.5 and mutations in RAS pathway genes occur early in the evolutionary history of the PAX-fusion-negative-RMS (PFN-RMS) subtype. We discover several early mutations in non-RAS mutated samples and predict them to be drivers in PFN-RMS including recurrent mutation of PKN1. In contrast, we find that PAX-fusion-positive (PFP) subtype tumors have undergone whole-genome duplication in the late stage of cancer evolutionary history and have acquired fewer mutations and subclones than PFN-RMS. Moreover we predict that the PAX3-FOXO1 fusion event occurs earlier than the whole genome duplication. Our findings provide information critical to the understanding of tumorigenesis of RMS

Offspring sex ratio and gonadal irradiation in the British Childhood Cancer Survivor Study

We investigated offspring sex ratio among 6232 offspring born to 3218 survivors of childhood cancer in relation to therapeutic irradiation, and pooled our data with those from two other large-scale studies giving a total of 9685 offspring. Exposure to high-dose gonadal irradiation was not associated with a significant alteration in offspring sex ratio compared to low doses (men: P=0.58, women: P=0.66). There was also no evidence that the ratio varied with time since cancer diagnosis when comparing survivors treated with radiotherapy vs those without (men: P=0.51; women: P=0.46). This, the largest study to date, finds no evidence that exposure to radiation affects the offspring sex ratio among survivors of childhood cancer

Atlantic overturning in decline?

Global ocean circulation is an important factor in climate variability and change. In particular, changes in the strength of the Atlantic meridional overturning circulation (AMOC) have been implicated in ancient climate events, as well as in recent climate anomalies such as the rapid warming of the North Atlantic Ocean in the mid-1990s. A series of moored current meters and temperature sensors deployed in the Atlantic at 26° N known as the RAPID-MOCHA array has been used to monitor the strength of meridional overturning since 2004. The data indicate a decline in this strength over the period 2004–20123. Here, using additional observations and climate model simulations we suggest that this measured decline is not merely a short-term fluctuation, but is part of a substantial reduction in meridional overturning occurring on a decadal timescale

Intragenic DNA methylation: implications of this epigenetic mechanism for cancer research

Epigenetics is the study of all mechanisms that regulate gene transcription and genome stability that are maintained throughout the cell division, but do not include the DNA sequence itself. The best-studied epigenetic mechanism to date is DNA methylation, where methyl groups are added to the cytosine base within cytosine–guanine dinucleotides (CpG sites). CpGs are frequently clustered in high density (CpG islands (CGIs)) at the promoter of over half of all genes. Current knowledge of transcriptional regulation by DNA methylation centres on its role at the promoter where unmethylated CGIs are present at most actively transcribed genes, whereas hypermethylation of the promoter results in gene repression. Over the last 5 years, research has gradually incorporated a broader understanding that methylation patterns across the gene (so-called intragenic or gene body methylation) may have a role in transcriptional regulation and efficiency. Numerous genome-wide DNA methylation profiling studies now support this notion, although whether DNA methylation patterns are a cause or consequence of other regulatory mechanisms is not yet clear. This review will examine the evidence for the function of intragenic methylation in gene transcription, and discuss the significance of this in carcinogenesis and for the future use of therapies targeted against DNA methylation

Enrichment analysis of Alu elements with different spatial chromatin proximity in the human genome

Transposable elements (TEs) have no longer been totally considered as “junk DNA” for quite a time since the continual discoveries of their multifunctional roles in eukaryote genomes. As one of the most important and abundant TEs that still active in human genome, Alu, a SINE family, has demonstrated its indispensable regulatory functions at sequence level, but its spatial roles are still unclear. Technologies based on 3C(chromosomeconformation capture) have revealed the mysterious three-dimensional structure of chromatin, and make it possible to study the distal chromatin interaction in the genome. To find the role TE playing in distal regulation in human genome, we compiled the new released Hi-C data, TE annotation, histone marker annotations, and the genome-wide methylation data to operate correlation analysis, and found that the density of Alu elements showed a strong positive correlation with the level of chromatin interactions (hESC: r=0.9, P<2.2×1016; IMR90 fibroblasts: r = 0.94, P < 2.2 × 1016) and also have a significant positive correlation withsomeremote functional DNA elements like enhancers and promoters (Enhancer: hESC: r=0.997, P=2.3×10−4; IMR90: r=0.934, P=2×10−2; Promoter: hESC: r = 0.995, P = 3.8 × 10−4; IMR90: r = 0.996, P = 3.2 × 10−4). Further investigation involving GC content and methylation status showed the GC content of Alu covered sequences shared a similar pattern with that of the overall sequence, suggesting that Alu elements also function as the GC nucleotide and CpG site provider. In all, our results suggest that the Alu elements may act as an alternative parameter to evaluate the Hi-C data, which is confirmed by the correlation analysis of Alu elements and histone markers. Moreover, the GC-rich Alu sequence can bring high GC content and methylation flexibility to the regions with more distal chromatin contact, regulating the transcription of tissue-specific genes

The epidemiology of injuries across the weight-training sports

Background: Weight-training sports, including weightlifting, powerlifting, bodybuilding, strongman, Highland Games, and CrossFit, are weight-training sports that have separate divisions for males and females of a variety of ages, competitive standards, and bodyweight classes. These sports may be considered dangerous because of the heavy loads commonly used in training and competition. Objectives: Our objective was to systematically review the injury epidemiology of these weight-training sports, and, where possible, gain some insight into whether this may be affected by age, sex, competitive standard, and bodyweight class. Methods: We performed an electronic search using PubMed, SPORTDiscus, CINAHL, and Embase for injury epidemiology studies involving competitive athletes in these weight-training sports. Eligible studies included peer-reviewed journal articles only, with no limit placed on date or language of publication. We assessed the risk of bias in all studies using an adaption of the musculoskeletal injury review method. Results: Only five of the 20 eligible studies had a risk of bias score ≥75 %, meaning the risk of bias in these five studies was considered low. While 14 of the studies had sample sizes >100 participants, only four studies utilized a prospective design. Bodybuilding had the lowest injury rates (0.12–0.7 injuries per lifter per year; 0.24–1 injury per 1000 h), with strongman (4.5–6.1 injuries per 1000 h) and Highland Games (7.5 injuries per 1000 h) reporting the highest rates. The shoulder, lower back, knee, elbow, and wrist/hand were generally the most commonly injured anatomical locations; strains, tendinitis, and sprains were the most common injury type. Very few significant differences in any of the injury outcomes were observed as a function of age, sex, competitive standard, or bodyweight class. Conclusion: While the majority of the research we reviewed utilized retrospective designs, the weight-training sports appear to have relatively low rates of injury compared with common team sports. Future weight-training sport injury epidemiology research needs to be improved, particularly in terms of the use of prospective designs, diagnosis of injury, and changes in risk exposure

Increasing incidence of childhood leukaemia: a controversy re-examined

We provide evidence of a gradual increase in the incidence of childhood leukaemia over the twentieth century from examination of trends in both incidence and mortality in England and Wales. We conclude that much of the recorded increase is likely to be real

Caspase Inhibitors of the P35 Family Are More Active When Purified from Yeast than Bacteria

Many insect viruses express caspase inhibitors of the P35 superfamily, which prevent defensive host apoptosis to enable viral propagation. The prototypical P35 family member, AcP35 from Autographa californica M nucleopolyhedrovirus, has been extensively studied. Bacterially purified AcP35 has been previously shown to inhibit caspases from insect, mammalian and nematode species. This inhibition occurs via a pseudosubstrate mechanism involving caspase-mediated cleavage of a “reactive site loop” within the P35 protein, which ultimately leaves cleaved P35 covalently bound to the caspase's active site. We observed that AcP35 purifed from Saccharomyces cerevisae inhibited caspase activity more efficiently than AcP35 purified from Escherichia coli. This differential potency was more dramatic for another P35 family member, MaviP35, which inhibited human caspase 3 almost 300-fold more potently when purified from yeast than bacteria. Biophysical assays revealed that MaviP35 proteins produced in bacteria and yeast had similar primary and secondary structures. However, bacterially produced MaviP35 possessed greater thermal stability and propensity to form higher order oligomers than its counterpart purified from yeast. Caspase 3 could process yeast-purified MaviP35, but failed to detectably cleave bacterially purified MaviP35. These data suggest that bacterially produced P35 proteins adopt subtly different conformations from their yeast-expressed counterparts, which hinder caspase access to the reactive site loop to reduce the potency of caspase inhibition, and promote aggregation. These data highlight the differential caspase inhibition by recombinant P35 proteins purified from different sources, and caution that analyses of bacterially produced P35 family members (and perhaps other types of proteins) may underestimate their activity