GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

CD28-negative cytolytic effector T cells frequently express NK receptors and are present at variable proportions in circulating lymphocytes from healthy donors and melanoma patients.

In humans, NK receptors are expressed by natural killer cells and some T cells, the latter of which are preferentially alphabetaTCR+ CD8+ cytolytic T lymphocytes (CTL). In this study we analyzed the expression of nine NK receptors (p58.1, p58.2, p70, p140, ILT2, NKRP1A, ZIN176, CD94 and CD94/NKG2A) in PBL from both healthy donors and melanoma patients. The percentages of NK receptor-positive T cells (NKT cells) varied strongly, and this variation was more important between individual patients than between individual healthy donors. In all the individuals, the NKT cells were preferentially CD28-, and a significant correlation was found between the percentage of CD28- T cells and the percentage of NK receptor+ T cells. Based on these data and the known activated phenotype of CD28- T cells, we propose that the CD28- CD8+ T cell pool represents or contains the currently active CTL population, and that the frequent expression of NK receptors reflects regulatory mechanisms modulating the extent of CTL effector function. Preliminary results indicate that some tumor antigen-specific T cells may indeed be CD28- and express NK receptors in vivo.

The value of intraoperative ultrasound in oblique corpectomy for cervical spondylotic myelopathy and ossified posterior longitudinal ligament.

Intraoperative ultrasound (IOUS) has been described to be useful during central corpectomy for compressive cervical myelopathy. This study aimed at documenting the utility of IOUS in oblique cervical corpectomy (OCC). Prospective data from 24 patients undergoing OCC for cervical spondylotic myelopathy and ossified posterior longitudinal ligament (OPLL) were collected. Patients had a preoperative cervical spine magnetic resonance (MR) image, IOUS and a postoperative cervical CT scan. Retrospective data from 16 historical controls that underwent OCC without IOUS were analysed to compare the incidence of residual compression between the two groups. IOUS identified the vertebral artery in all cases, detected residual cord compression in six (27%) and missed compression in two cases (9%). In another two cases with OPLL, IOUS was sub-optimal due to shadowing. IOUS measurement of the corpectomy width correlated well with these measurements on the postoperative CT. The extent of cord expansion noted on IOUS after decompression showed no correlation with immediate or 6-month postoperative neurological recovery. No significant difference in residual compression was noted in the retrospective and prospective groups of the study. Craniocaudal spinal cord motion was noted after the completion of the corpectomy. IOUS is an inexpensive and simple real-time imaging modality that may be used during OCC for cervical spondylotic myelopathy. It is helpful in identifying the vertebral artery and determining the trajectory of approach, however, it has limited utility in patients with OPLL due to artifacts from residual ossification.

The effect of obesity, age, puberty and gender on resting metabolic rate in children and adolescents.

During puberty fat-free mass (FFM) and fat mass (FM) change quickly and these changes are influenced by sex and obesity. Since it is not completely known how these changes affect resting metabolic rate (RMR), the aim of the present study was to investigate the effect of body composition, age, sex and pubertal development of postabsorptive RMR in 9.5- to 16.5- year-old obese and non-obese children. Postabsorptive RMR was measured in a sample of 371 pre- and postpubertal children comprising 193 males (116 non-obese and 77 obese) and 178 females (119 non-obese and 59 obese). RMR was assessed by indirect calorimetry using a ventilated hood system for 45 min after an overnight fast. Body composition (FFM and FM) was estimated from skinfold measurements. The mean (+/- SD) RMR was significantly (P < 0.001) lower in non-obese (males: 5600 +/- 972 kJ/24 h; females: 5112 +/- 632 kJ/24 h) than in obese (males: 7223 +/- 1220 kJ/24 h; females: 6665 +/- 1106 kJ/24 h) children. This difference became non-significant when RMR was adjusted for body composition (FFM+FM). However, the difference between the genders still remained significant (control male: 6118 +/- 507, control female: 5652 +/- 507, P < 0.001; obese male: 6256 +/- 507, obese female: 5818 +/- 507 kJ/24 h, P < 0.001). The main determinant of RMR was FFM. In the whole cohort. FFM explained 79.8% of the variation in RMR, followed by age, gender and FM adding further 3.8%, 1.1% and 0.8% to the predictability of RMR, respectively. No significant contribution for study group (obese, non-obese), pubertal stage, or fat distribution was found in the regression for RMR. The adjusted value of RMR (for FFM and FM) slightly, but significantly (P < 0.01) decreased between the age of 10-16 years, demonstrating the important effect of age on RMR. CONCLUSIONS: The resting metabolic rate of obese and control children is not different when adjusted for body composition. The main determinant of RMR is the fat-free mass, however, age, gender and fat mass are also significant factors. Pubertal development and fat distribution do not influence RMR independently from the changes in body composition.

Genome-wide association study identifies six new loci influencing pulse pressure and mean arterial pressure.

Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans. We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N = 74,064) and follow-up studies (N = 48,607), we identified at genome-wide significance (P = 2.7 × 10(-8) to P = 2.3 × 10(-13)) four new PP loci (at 4q12 near CHIC2, 7q22.3 near PIK3CG, 8q24.12 in NOV and 11q24.3 near ADAMTS8), two new MAP loci (3p21.31 in MAP4 and 10q25.3 near ADRB1) and one locus associated with both of these traits (2q24.3 near FIGN) that has also recently been associated with SBP in east Asians. For three of the new PP loci, the estimated effect for SBP was opposite of that for DBP, in contrast to the majority of common SBP- and DBP-associated variants, which show concordant effects on both traits. These findings suggest new genetic pathways underlying blood pressure variation, some of which may differentially influence SBP and DBP.

Anti-HLA antibody repertoire after IVIg infusion in highly sensitized patients waiting for kidney transplantation.

Polyclonal intravenous immunoglobulin (IVIg) treatment reduces crossmatch positivity and increases rates of transplantation in highly sensitised patients (HS). We quantified the panel reactive antibody (PRA) by microlymphocytotoxicity (MLCC), and we analysed anti-HLA class I and class II IgG specific antibody repertoire by Luminex before and after IVIg infusion alone in HS patients awaiting kidney transplantation. Five patients received three monthly infusions of 1 g/kg of IVIg. Serum samples collected pre and post IVIg treatment were submitted for PRA analysis by MLCC. Anti-class I and anti-class II antibody specificities were then tested by Luminex. We focused on the anti-HLA class I and class II antibodies directed against HLA expressed by a previous graft. We also analysed the anti-HLA antibody repertoire in three patients who had not received IVIg infusion. The PRA level determined by MLCC decreased significantly in one of the five patients, dropping from 40% to 17%. The Luminex assay showed fluctuations of the anti-HLA antibody levels over time, but no significant longterm modifications of the anti-HLA antibody repertoire were observed, even in the patient with a strong and prolonged reduction of the PRA determined by MLCC. Our results show that IVIg at 1 g/kg is not sufficient to reduce PRA and does not modify the repertoire of specific anti-HLA antibody determined by Luminex.

Diurnal and exercise-related variability of haemoglobin and reticulocytes in athletes.

Haemoglobin (Hb) and Reticulocytes (Ret) are measured as indirect markers of doping in athletes. We studied the diurnal variation, the impact of exercise, fluid intake and ambient temperature in athletes on these parameters. Hourly venous blood samples were obtained from 36 male athletes of different disciplines (endurance (END) and non-endurance (NON-END)) over 12 h during a typical training day. Seven inactive subjects served as controls (CON). Hb and Ret were determined. A mixed model procedure was used to analyse the data. At baseline, Hb was similar for all groups, END showed lower Ret than NON-END and CON. Exercise showed a significant impact on Hb (+0.46 g/dl, p<0.001), the effect disappeared approximately 2 h after exercise. Hb decreased over the day by approximately 0.55 g/dl (p<0.01). There was no relevant effect on Ret. Fluid intake and ambient temperature had no significant effect. Hb shows significant diurnal- and exercise related variations. In an anti-doping context, most of these variations are in favour of the athlete. Blood samples taken after exercise might therefore provide reliable results and thus be used for the longitudinal monitoring of athletes if a timeframe for the re-equilibration of vascular volumes is respected.

Educational achievement in Swiss childhood cancer survivors compared with the general population.

BACKGROUND: The objective of this study was to describe educational achievements of childhood cancer survivors in Switzerland compared with the general population. In particular, the authors investigated educational problems during childhood, final educational achievement in adulthood, and its predictors. METHODS: Childhood cancer survivors who were aged <16 years at diagnosis from 1976 to 2003 who had survived for ≥5 years and were currently ages 20 to 40 years received a postal questionnaire during 2007 to 2009. Controls were respondents of the Swiss Health Survey ages 20 to 40 years. Educational achievement included compulsory schooling, vocational training, upper secondary schooling, and university degree. The analysis was weighted to optimize comparability of the populations. The authors analyzed the association between demographic and clinical predictors and educational achievement using multivariable logistic regression. Subgroup analyses focused on survivors aged ≥27 years. RESULTS: One-third of survivors encountered educational problems during schooling (30% repeated 1 year, and 35% received supportive tutoring). In the total sample, more survivors than controls achieved compulsory schooling only (8.7% vs 5.2%) and fewer acquired a university degree (7.3% vs 11%), but more survivors than controls achieved an upper secondary education (36.1 vs 24.1%). In those aged ≥27 years, differences in compulsory schooling and university education largely disappeared. In survivors and controls, sex, nationality, language region, and migration background were strong predictors of achievement. Survivors of central nervous system tumors or those who had a relapse had poorer outcomes (P < .05). CONCLUSIONS: Childhood cancer survivors encountered problems during schooling and completed professional education with some delay. However, with the exception of patients who had central nervous system tumors and those who experienced a relapse, the final educational achievement in survivors of child cancer was comparable to that of the general population.

Single nucleotide polymorphisms, apoptosis, and the development of severe late adverse effects after radiotherapy.

PURPOSE: Evidence has accumulated in recent years suggestive of a genetic basis for a susceptibility to the development of radiation injury after cancer radiotherapy. The purpose of this study was to assess whether patients with severe radiation-induced sequelae (RIS; i.e., National Cancer Institute/CTCv3.0 grade, > or =3) display both a low capacity of radiation-induced CD8 lymphocyte apoptosis (RILA) in vitro and possess certain single nucleotide polymorphisms (SNP) located in candidate genes associated with the response of cells to radiation. EXPERIMENTAL DESIGN: DNA was isolated from blood samples obtained from patients (n = 399) included in the Swiss prospective study evaluating the predictive effect of in vitro RILA and RIS. SNPs in the ATM, SOD2, XRCC1, XRCC3, TGFB1, and RAD21 genes were screened in patients who experienced severe RIS (group A, n = 16) and control subjects who did not manifest any evidence of RIS (group B, n = 18). RESULTS: Overall, 13 and 21 patients were found to possess a total of <4 and > or =4 SNPs in the candidate genes. The median (range) RILA in group A was 9.4% (5.3-16.5) and 94% (95% confidence interval, 70-100) of the patients (15 of 16) had > or =4 SNPs. In group B, median (range) RILA was 25.7% (20.2-43.2) and 33% (95% confidence interval, 13-59) of patients (6 of 18) had > or =4 SNPs (P < 0.001). CONCLUSIONS: The results of this study suggest that patients with severe RIS possess 4 or more SNPs in candidate genes and low radiation-induced CD8 lymphocyte apoptosis in vitro.

Comparison of aortic elasticity determined by cardiovascular magnetic resonance imaging in obese versus lean adults.

The vascular properties of large vessels in the obese have not been adequately studied. We used cardiovascular magnetic resonance imaging to quantify the cross-sectional area and elastic properties of the ascending thoracic and abdominal aorta in 21 clinically healthy obese young adult men and 25 men who were age-matched lean controls. Obese subjects had greater maximal cross-sectional area of the ascending thoracic aorta (984 +/- 252 vs 786 +/- 109 mm(2), p <0.01) and of the abdominal aorta (415 +/- 71 vs 374 +/- 51 mm(2), p <0.05). When indexed for height the differences persisted, but when indexed for body surface area, a significant difference between groups was found only for the maximal abdominal aortic cross-sectional area. The obese subjects also had decreased abdominal aortic elasticity, characterized by 24% lower compliance (0.0017 +/- 0.0004 vs 0.0021 +/- 0.0005 mm(2)/kPa/mm, p <0.01), 22% higher stiffness index beta (6.0 +/- 1.5 vs 4.9 +/- 0.7, p <0.005), and 41% greater pressure-strain elastic modulus (72 +/- 25 vs 51 +/- 9, p <0.005). At the ascending thoracic aorta, only the pressure-strain elastic modulus was different between obese and lean subjects (85 +/- 42 vs 65 +/- 26 kPa, respectively; p <0.05), corresponding to a 31% difference-but arterial compliance and stiffness index were not significantly different between groups. In clinically healthy young adult obese men, obesity is associated with increased cross-sectional aortic area and decreased aortic elasticity.

Type of alcoholic beverage and risk of head and neck cancer: a pooled analysis within the INHANCE Consortium

The authors pooled data from 15 case-control studies of head and neck cancer (9,107 cases, 14,219 controls) to investigate the independent associations with consumption of beer, wine, and liquor. In particular, they calculated associations with different measures of beverage consumption separately for subjects who drank beer only (858 cases, 986 controls), for liquor-only drinkers (499 cases, 527 controls), and for wine-only drinkers (1,021 cases, 2,460 controls), with alcohol never drinkers (1,124 cases, 3,487 controls) used as a common reference group. The authors observed similar associations with ethanol-standardized consumption frequency for beer-only drinkers (odds ratios (ORs) = 1.6, 1.9, 2.2, and 5.4 for < or =5, 6-15, 16-30, and >30 drinks per week, respectively; P(trend) < 0.0001) and liquor-only drinkers (ORs = 1.6, 1.5, 2.3, and 3.6; P < 0.0001). Among wine-only drinkers, the odds ratios for moderate levels of consumption frequency approached the null, whereas those for higher consumption levels were comparable to those of drinkers of other beverage types (ORs = 1.1, 1.2, 1.9, and 6.3; P < 0.0001). Study findings suggest that the relative risks of head and neck cancer for beer and liquor are comparable. The authors observed weaker associations with moderate wine consumption, although they cannot rule out confounding from diet and other lifestyle factors as an explanation for this finding. Given the presence of heterogeneity in study-specific results, their findings should be interpreted with caution.

Valproic acid monotherapy in pregnancy and major congenital malformations.

BACKGROUND: The use of valproic acid in the first trimester of pregnancy is associated with an increased risk of spina bifida, but data on the risks of other congenital malformations are limited. METHODS: We first combined data from eight published cohort studies (1565 pregnancies in which the women were exposed to valproic acid, among which 118 major malformations were observed) and identified 14 malformations that were significantly more common among the offspring of women who had received valproic acid during the first trimester. We then assessed the associations between use of valproic acid during the first trimester and these 14 malformations by performing a case-control study with the use of the European Surveillance of Congenital Anomalies (EUROCAT) antiepileptic-study database, which is derived from population-based congenital-anomaly registries. Registrations (i.e., pregnancy outcomes with malformations included in EUROCAT) with any of these 14 malformations were compared with two control groups, one consisting of infants with malformations not previously linked to valproic acid use (control group 1), and one consisting of infants with chromosomal abnormalities (control group 2). The data set included 98,075 live births, stillbirths, or terminations with malformations among 3.8 million births in 14 European countries from 1995 through 2005. RESULTS: Exposure to valproic acid monotherapy was recorded for a total of 180 registrations, with 122 registrations in the case group, 45 in control group 1, and 13 in control group 2. As compared with no use of an antiepileptic drug during the first trimester (control group 1), use of valproic acid monotherapy was associated with significantly increased risks for 6 of the 14 malformations under consideration; the adjusted odds ratios were as follows: spina bifida, 12.7 (95% confidence interval [CI], 7.7 to 20.7); atrial septal defect, 2.5 (95% CI, 1.4 to 4.4); cleft palate, 5.2 (95% CI, 2.8 to 9.9); hypospadias, 4.8 (95% CI, 2.9 to 8.1); polydactyly, 2.2 (95% CI, 1.0 to 4.5); and craniosynostosis, 6.8 (95% CI, 1.8 to 18.8). Results for exposure to valproic acid were similar to results for exposure to other antiepileptic drugs. CONCLUSIONS: The use of valproic acid monotherapy in the first trimester was associated with significantly increased risks of several congenital malformations, as compared with no use of antiepileptic drugs or with use of other antiepileptic drugs.

No longitudinal mitochondrial DNA sequence changes in HIV-infected individuals with and without lipoatrophy.

The potential for mitochondrial (mt) DNA mutation accumulation during antiretroviral therapy (ART), and preferential accumulation in patients with lipoatrophy compared with control participants, remains controversial. We sequenced the entire mitochondrial genome, both before ART and after ART exposure, in 29 human immunodeficiency virus (HIV)-infected Swiss HIV Cohort Study participants initiating a first-line thymidine analogue-containing ART regimen. No accumulation of mtDNA mutations or deletions was detected in 13 participants who developed lipoatrophy or in 16 control participants after significant and comparable ART exposure (median duration, 3.3 and 3.7 years, respectively). In HIV-infected persons, the development of lipoatrophy is unlikely to be associated with accumulation of mtDNA mutations detectable in peripheral blood.

Factors predisposing to dislocation after primary total hip arthroplasty: a multivariate analysis.

We conducted this study to determine the relative influence of various mechanical and patient-related factors on the incidence of dislocation after primary total hip asthroplasty (THA). Of 2,023 THAs, 21 patients who had at least 1 dislocation were compared with a control group of 21 patients without dislocation, matched for age, gender, pathology, and year of surgery. Implant positioning, seniority of the surgeon, American Society of Anesthesiologists (ASA) score, and diminished motor coordination were recorded. Data analysis included univariate and multivariate methods. The dislocation risk was 6.9 times higher if total anteversion was not between 40 degrees and 60 degrees and 10 times higher in patients with high ASA scores. Surgeons should pay attention to total anteversion (cup and stem) of THA. The ASA score should be part of the preoperative assessment of the dislocation risk.