Multiple hot-deck imputation for network inference from RNA sequencing data.

Motivation: Network inference provides a global view of the relations existing between gene expression in a given transcriptomic experiment (often only for a restricted list of chosen genes). However, it is still a challenging problem: even if the cost of sequencing techniques has decreased over the last years, the number of samples in a given experiment is still (very) small compared to the number of genes. Results: We propose a method to increase the reliability of the inference when RNA-seq expression data have been measured together with an auxiliary dataset that can provide external information on gene expression similarity between samples. Our statistical approach, hd-MI, is based on imputation for samples without available RNA-seq data that are considered as missing data but are observed on the secondary dataset. hd-MI can improve the reliability of the inference for missing rates up to 30% and provides more stable networks with a smaller number of false positive edges. On a biological point of view, hd-MI was also found relevant to infer networks from RNA-seq data acquired in adipose tissue during a nutritional intervention in obese individuals. In these networks, novel links between genes were highlighted, as well as an improved comparability between the two steps of the nutritional intervention. Availability and implementation: Software and sample data are available as an R package, RNAseqNet, that can be downloaded from the Comprehensive R Archive Network (CRAN). Contact: alyssa.imbert@inra.fr or nathalie.villa-vialaneix@inra.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Label-free quantitative urinary proteomics identifies the arginase pathway as a new player in congenital obstructive nephropathy.

Obstructive nephropathy is a frequently encountered situation in newborns. In previous studies, the urinary peptidome has been analyzed for the identification of clinically useful biomarkers of obstructive nephropathy. However, the urinary proteome has not been explored yet and should allow additional insight into the pathophysiology of the disease. We have analyzed the urinary proteome of newborns (n = 5/group) with obstructive nephropathy using label free quantitative nanoLC-MS/MS allowing the identification and quantification of 970 urinary proteins. We next focused on proteins exclusively regulated in severe obstructive nephropathy and identified Arginase 1 as a potential candidate molecule involved in the development of obstructive nephropathy, located at the crossroad of pro- and antifibrotic pathways. The reduced urinary abundance of Arginase 1 in obstructive nephropathy was verified in independent clinical samples using both Western blot and MRM analysis. These data were confirmed in situ in kidneys obtained from a mouse obstructive nephropathy model. In addition, we also observed increased expression of Arginase 2 and increased total arginase activity in obstructed mouse kidneys. mRNA expression analysis of the related arginase pathways indicated that the pro-fibrotic arginase-related pathway is activated during obstructive nephropathy. Taken together we have identified a new actor in the development of obstructive nephropathy in newborns using quantitative urinary proteomics and shown its involvement in an in vivo model of disease. The present study demonstrates the relevance of such a quantitative urinary proteomics approach with clinical samples for a better understanding of the pathophysiology and for the discovery of potential therapeutic targets.

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis.

During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

Functional and radiological outcome of subtalar arthroereisis for flexible pes planovalgus in children: A retrospective analysis.

INTRODUCTION: Pediatric idiopathic pes planovalgus can correct itself with growth. Otherwise, in the event of functional impairment and after failed conservative treatment, surgery can be considered. Based on a multicenter retrospective study, we report the functional and radiographic results obtained after subtalar arthroereisis. HYPOTHESIS: We hypothesized that this surgery improves functional and radiological parameters in childhood. MATERIAL AND METHOD: Forty-eight medical records of children (78 feet) operated on between 2010 and 2019 were studied. Functional (FAOS score) and radiographic (Djian angle, calcaneal slope, lateral talocalcaneal divergence and calcaneus/M5 alignment, talonavicular coverage measurement, AP talocalcaneal divergence) results were studied. The analysis of these different criteria was carried out between the preoperative period and the last follow-up. RESULTS: The functional outcome was satisfactory with an average FAOS questionnaire score of 95.5 out of 100 total points. All the radiographic parameters studied were significantly improved (p<0.001). The average age at the time of surgery was 11.3 years (range: 7 to 16) with a mean follow-up of 35 months (range: 18 months to 84). Spontaneous screw expulsion and subtalar pain were the main complications found. DISCUSSION: The results obtained are consistent with those in the literature. The age at the time of the surgery is an essential factor to consider with the goal of optimal correction. CONCLUSION: This technique is reliable and effective in the short term. It can be offered as first-line therapy in the management of symptomatic pes planovalgus in children. The follow-up is short, which necessitates longer term studies of this population. The ideal age for surgery remains to be determined. LEVEL OF EVIDENCE: IV.

The GMO90+ Project: Absence of Evidence for Biologically Meaningful Effects of Genetically Modified Maize-based Diets on Wistar Rats After 6-Months Feeding Comparative Trial.

The GMO90+ project was designed to identify biomarkers of exposure or health effects in Wistar Han RCC rats exposed in their diet to 2 genetically modified plants (GMP) and assess additional information with the use of metabolomic and transcriptomic techniques. Rats were fed for 6-months with 8 maize-based diets at 33% that comprised either MON810 (11% and 33%) or NK603 grains (11% and 33% with or without glyphosate treatment) or their corresponding near-isogenic controls. Extensive chemical and targeted analyses undertaken to assess each diet demonstrated that they could be used for the feeding trial. Rats were necropsied after 3 and 6 months. Based on the Organization for Economic Cooperation and Development test guideline 408, the parameters tested showed a limited number of significant differences in pairwise comparisons, very few concerning GMP versus non-GMP. In such cases, no biological relevance could be established owing to the absence of difference in biologically linked variables, dose-response effects, or clinical disorders. No alteration of the reproduction function and kidney physiology was found. Metabolomics analyses on fluids (blood, urine) were performed after 3, 4.5, and 6 months. Transcriptomics analyses on organs (liver, kidney) were performed after 3 and 6 months. Again, among the significant differences in pairwise comparisons, no GMP effect was observed in contrast to that of maize variety and culture site. Indeed, based on transcriptomic and metabolomic data, we could differentiate MON- to NK-based diets. In conclusion, using this experimental design, no biomarkers of adverse health effect could be attributed to the consumption of GMP diets in comparison with the consumption of their near-isogenic non-GMP controls.

Multidimensional reduction of multicentric cohort heterogeneity: An alternative method to increase statistical power and robustness.

Modern clinical research takes advantage of multicentric cohorts to increase sample size and gain in statistical power. However, combining individuals from different recruitment centers provides heterogeneity in the dataset that needs to be accounted for to obtain robust results. Sophisticated statistical multivariate models adjusting for center effect can be implemented, but they can become unstable and can be complex to interpret with the increasing number of covariates to consider. Here, we present a multidimensional reduction technique to identify heterogeneity in a French multicentric cohort of hematopoietic stem cell transplantations and characterize a homogeneous subgroup prior to performing simple statistical univariate analyses. The exclusion of outliers allowed the identification of two genetic factors associated with post-transplantation overall survival. We therefore provide proof-of-concept that a sample size reduction method can efficiently account for heterogeneity and center effect in multicentric cohorts while increasing statistical power and robustness for discovery of new association signals.

Proceedings: human leukocyte antigen haplo-homozygous induced pluripotent stem cell haplobank modeled after the california population: evaluating matching in a multiethnic and admixed population.

The development of a California-based induced pluripotent stem cell (iPSC) bank based on human leukocyte antigen (HLA) haplotype matching represents a significant challenge and a valuable opportunity for the advancement of regenerative medicine. However, previously published models of iPSC banks have neither addressed the admixed nature of populations like that of California nor evaluated the benefit to the population as a whole. We developed a new model for evaluating an iPSC haplobank based on demographic and immunogenetic characteristics reflecting California. The model evaluates haplolines or cell lines from donors homozygous for a single HLA-A, HLA-B, HLA-DRB1 haplotype. We generated estimates of the percentage of the population matched under various combinations of haplolines derived from six ancestries (black/African American, American Indian, Asian/Pacific Islander, Hispanic, and white/not Hispanic) and data available from the U.S. Census Bureau, the California Institute for Regenerative Medicine, and the National Marrow Donor Program. The model included both cis (haplotype-level) and trans (genotype-level) matching between a modeled iPSC haplobank and the recipient population following resampling simulations. We showed that serving a majority (>50%) of a simulated California population through cis matching would require the creation, redundant storage, and maintenance of almost 207 different haplolines representing the top 60 most frequent haplotypes from each ancestry group. Allowances for trans matching reduced the haplobank to fewer than 141 haplolines found among the top 40 most frequent haplotypes. Finally, we showed that a model optimized, custom haplobank was able to serve a majority of the California population with fewer than 80 haplolines.

Adipose tissue CIDEA is associated, independently of weight variation, to change in insulin resistance during a longitudinal weight control dietary program in obese individuals.

AIM: Weight loss reduces risk factors associated with obesity. However, long-term metabolic improvement remains a challenge. We investigated quantitative gene expression of subcutaneous adipose tissue in obese individuals and its relationship with low calorie diet and long term weight maintenance induced changes in insulin resistance. RESEARCH DESIGN: Three hundred eleven overweight and obese individuals followed a dietary protocol consisting of an 8-week low calorie diet followed by a 6-month ad libitum weight-maintenance diet. Individuals were clustered according to insulin resistance trajectories assessed using homeostasis model assessment of insulin resistance (HOMA-IR) index. Adipose tissue mRNA levels of 267 genes selected for regulation according to obesity, metabolic status and response to dieting was assessed using high throughput RT-qPCR. A combination of discriminant analyses was used to identify genes with regulation according to insulin resistance trajectories. Partial correlation was used to control for change in body mass index. RESULTS: Three different HOMA-IR profile groups were determined. HOMA-IR improved during low calorie diet in the 3 groups. At the end of the 6-month follow-up, groups A and B had reduced HOMA-IR by 50%. In group C, HOMA-IR had returned to baseline values. Genes were differentially expressed in the adipose tissue of individuals according to groups but a single gene, CIDEA, was common to all phases of the dietary intervention. Changes in adipose tissue CIDEA mRNA levels paralleled variations in insulin sensitivity independently of change in body mass index. Overall, CIDEA was up-regulated in adipose tissue of individuals with successful long term insulin resistance relapse and not in adipose tissue of unsuccessful individuals. CONCLUSION: The concomitant change in adipose tissue CIDEA mRNA levels and insulin sensitivity suggests a beneficial role of adipose tissue CIDEA in long term glucose homeostasis, independently of weight variation. TRIAL REGISTRATION: ClinicalTrials.gov NCT00390637.

Why statistics matter: limited inter-rater agreement prevents using the psoriasis area and severity index as a unique determinant of therapeutic decision in psoriasis.

Although the demand for evidence-based decisions is increasing in clinical practice, recent systematic reviews on the accuracy of existing psoriasis severity scales, including the Psoriasis Area and Severity Index (PASI), suggest that their validity is not fully characterized. We simulated the evaluation of PASI by two practitioners in 1,000 sets of 100 patients. PASI data from several practitioners who examined the same patients were used to generate PASI scores by two practitioners, in order to compare how well commonly used statistics assess the inter-rater agreement for the PASI. Because the PASI score has an asymmetric distribution, statistics such as Pearson's linear correlation coefficient "r" and Spearman's rank correlation coefficient overestimated the inter-rater agreement of PASI as compared with the intra-class correlation coefficient (ICC; r=0.8, ρ=0.7, ICC=0.5). When restricting the analysis to patients with a PASI <20, inter-rater agreement severely decreased (r=0.38, ρ=0.41, ICC=0.17), resulting in unacceptable therapeutic decision agreement (κ=0.38). Our study indicates that owing to the skewed distribution of the PASI its validity to influence therapeutic decisions is questionable. The ICC is preferable to the commonly used statistics (r and ρ) for assessing the inter-rater agreement reliability of asymmetrically distributed scores such as the PASI.